https://github.com/cran/aster
Tip revision: c45193827d47cd090900035e43f3188cd3fbc795 authored by Charles J. Geyer on 05 May 2013, 00:00:00 UTC
version 0.8-23
version 0.8-23
Tip revision: c451938
reaster1.Rout.save
R version 2.15.1 (2012-06-22) -- "Roasted Marshmallows"
Copyright (C) 2012 The R Foundation for Statistical Computing
ISBN 3-900051-07-0
Platform: x86_64-unknown-linux-gnu (64-bit)
R is free software and comes with ABSOLUTELY NO WARRANTY.
You are welcome to redistribute it under certain conditions.
Type 'license()' or 'licence()' for distribution details.
R is a collaborative project with many contributors.
Type 'contributors()' for more information and
'citation()' on how to cite R or R packages in publications.
Type 'demo()' for some demos, 'help()' for on-line help, or
'help.start()' for an HTML browser interface to help.
Type 'q()' to quit R.
>
> library(aster)
Loading required package: trust
>
> data(radish)
>
> options(digits=4) # avoid rounding differences
>
> pred <- c(0,1,2)
> fam <- c(1,3,2)
>
> rout <- reaster(resp ~ varb + fit : (Site * Region), ~ 0 + fit : Block,
+ pred, fam, varb, id, root, data = radish)
> summary(rout)
Call:
reaster.formula(fixed = resp ~ varb + fit:(Site * Region), random = ~0 +
fit:Block, pred = pred, fam = fam, varvar = varb, idvar = id,
root = root, data = radish)
Fixed Effects:
Estimate Std. Error z value Pr(>|z|)
(Intercept) -471.6346 1.7869 -263.94 < 2e-16 ***
varbFlowers 478.5295 1.7892 267.45 < 2e-16 ***
varbFruits 470.5709 1.7938 262.33 < 2e-16 ***
fit:SitePoint Reyes -0.1139 0.1852 -0.62 0.54
fit:RegionS -0.0856 0.0121 -7.06 1.7e-12 ***
fit:SiteRiverside:RegionS 0.4334 0.0157 27.67 < 2e-16 ***
---
Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
Square Roots of Variance Components (P-values are one-tailed):
Estimate Std. Error z value Pr(>|z|)/2
sigma1 0.2922 0.0655 4.46 4e-06 ***
---
Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
> summary(rout, stand = FALSE)
Call:
reaster.formula(fixed = resp ~ varb + fit:(Site * Region), random = ~0 +
fit:Block, pred = pred, fam = fam, varvar = varb, idvar = id,
root = root, data = radish)
Fixed Effects:
Estimate Std. Error z value Pr(>|z|)
(Intercept) -471.6346 1.7869 -263.94 < 2e-16 ***
varbFlowers 478.5295 1.7892 267.45 < 2e-16 ***
varbFruits 470.5709 1.7938 262.33 < 2e-16 ***
fit:SitePoint Reyes -0.1139 0.1852 -0.62 0.54
fit:RegionS -0.0856 0.0121 -7.06 1.7e-12 ***
fit:SiteRiverside:RegionS 0.4334 0.0157 27.67 < 2e-16 ***
---
Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
Variance Components (P-values are one-tailed):
Estimate Std. Error z value Pr(>|z|)/2
nu1 0.0854 0.0383 2.23 0.013 *
---
Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
>
>
> proc.time()
user system elapsed
1.155 0.040 1.324