C:\sta6167>C:\Progra~1\R\R-2.3.1\bin\Rterm --vanilla R : Copyright 2006, The R Foundation for Statistical Computing Version 2.3.1 (2006-06-01) ISBN 3-900051-07-0 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. Natural language support but running in an English locale 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. > pdf("wintwheat.pdf") > > ww0 <- read.fwf("C:\\data\\wintwheat.dat", width=c(8,8,8,8), + col.names=c("idnum", "variety", "yield", "moist")) > > ww1 <- data.frame(yield=ww0$yield, moist=ww0$moist, variety=factor(ww0$variety)) > > attach(ww1) > > library(nlme) > > yield.rcr1 <- lme(yield ~ moist , random =~1+moist | variety) > > summary(yield.rcr1) Linear mixed-effects model fit by REML Data: NULL AIC BIC logLik 198.1463 210.5089 -93.07314 Random effects: Formula: ~1 + moist | variety Structure: General positive-definite, Log-Cholesky parametrization StdDev Corr (Intercept) 4.34679573 (Intr) moist 0.04893021 -0.342 Residual 0.59334279 Fixed effects: yield ~ moist Value Std.Error DF t-value p-value (Intercept) 33.43388 1.3984776 49 23.90734 0 moist 0.66166 0.0167837 49 39.42240 0 Correlation: (Intr) moist -0.376 Standardized Within-Group Residuals: Min Q1 Med Q3 Max -1.67798841 -0.56317596 -0.05373476 0.58221036 1.97151334 Number of Observations: 60 Number of Groups: 10 > > ranef(yield.rcr1) (Intercept) moist 1 0.9577951 -0.04921124 2 -2.2842787 -0.06669719 3 -0.4081180 0.06722274 4 0.6960210 -0.02330617 5 1.1159073 -0.01990371 6 4.6391475 0.02388878 7 -10.7300453 0.05642355 8 2.4011668 0.02243373 9 -0.1762116 0.02335677 10 3.7886159 -0.03420725 > > dev.off() null device 1 >