https://github.com/cran/aster
Raw File
Tip revision: c45193827d47cd090900035e43f3188cd3fbc795 authored by Charles J. Geyer on 05 May 2013, 00:00:00 UTC
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 
back to top