2 xlabtrayylabseedling dry weight col21d1bpch cex2

Info iconThis preview shows page 1. Sign up to view the full content.

View Full Document Right Arrow Icon
This is the end of the preview. Sign up to access the rest of the document.

Unformatted text preview: 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 B B B B B B B B B B B B B B B B B B 6 6 6 6 7 7 7 7 7 7 8 8 8 8 8 8 8 8 4 5 6 7 1 2 3 4 5 6 1 2 3 4 5 6 7 8 8 13 9 11 12 16 17 15 15 15 9 6 8 8 13 9 9 10 2 4 plot(d[,2],d[,4]+rnorm(56,0,.2), xlab="Tray",ylab="Seedling Dry Weight", col=2*(1+(d[,1]=="B")),pch="-",cex=2) legend("topright",c("Genotype A","Genotype B"), fill=c(2,4),border=c(2,4)) 2 5 2 6 library(nlme) lme(SeedlingWeight~Genotype,random=~1|Tray, method="ML",data=d) Linear mixed-effects model fit by maximum likelihood Data: d Log-likelihood: -126.3709 Fixed: SeedlingWeight ~ Genotype (Intercept) GenotypeB 15.301832 -3.567017 Random effects: Formula: ~1 | Tray (Intercept) Residual StdDev: 2.932294 1.882470 Number of Observations: 56 Number of Groups: 8 2 7 library(lme4) lmer(SeedlingWeight~Genotype+(1|Tray),REML=F,data=d) Linear mixed model fit by maximum likelihood Formula: SeedlingWeight ~ Genotype + (1 | Tray) Data: d AIC BIC logLik deviance REMLdev 260.7 268.8 -126.4 252.7 247.4 Random effects: Groups Name Variance Std.Dev. Tray (Intercept) 8.5984 2.9323 Residual 3.5437 1.8825 Number of obs: 56, groups: Tray, 8 2 8 Fixed effects: Estimate Std. Error t value (Intercept) 15.302 1.510 10.14 GenotypeB -3.567 2.136 -1.67 Correlation of Fixed Effects: (Intr) GenotypeB -0.707 2 9 lme(SeedlingWeight~Genotype,random=~1|Tray,data=d) Linear mixed-effects model fit by REML Data: d Log-restricted-likelihood: -123.5705 Fixed: SeedlingWeight ~ Genotype (Intercept) GenotypeB 15.288838 -3.550201 Random effects: Formula: ~1 | Tray (Intercept) Residual StdDev: 3.414856 1.882230 Number of Observations: 56 Number of Groups: 8 3 0 lmer(SeedlingWeight~Genotype+(1|Tray),data=d) Linear mixed model fit by REML Formula: SeedlingWeight ~ Genotype + (1 | Tray) Data: d AIC BIC logLik deviance REMLdev 255.1 263.2 -123.6 253 247.1 Random effects: Groups Name Variance Std.Dev. Tray (Intercept) 11.6612 3.4149 Residual 3.5428 1.8822 Number of obs: 56, groups: Tray, 8 Fixed effects: Estimate Std. Error t value (Intercept) 15.289 1.745 8.762 GenotypeB -3.550 2.468 -1.438 Correlation of Fixed Effects: (Intr) GenotypeB -0.707 3 1...
View Full Document

This document was uploaded on 03/27/2014.

Ask a homework question - tutors are online