Displaying 1 result from an estimated 1 matches for "leaflmm".
2003 May 21
1
Code Help
...s.txt attached. Comma seperations
Code:
options(digits=3)
clover <- read.table("Lines.txt",header=T,sep=",")
vnames <- names(clover);nv <- length(vnames)
flevels <- levels(clover$Line)
par(oma=c(0,6,4,6),mfrow=c(2,1),mar=c(4,4,1,0))
clover$Y <- clover[,"LeafLmm"]
mn.tab <- tapply(clover$Y,list(clover$Line),mean)
sd.tab <- sqrt(tapply(clover$Y,list(clover$Line),var))
om <- order(mn.tab)
plot(mn.tab,sd.tab,las=1,xlab="mean",ylab="s.d.")
lines(lowess(mn.tab,sd.tab))
model.i <- glm(Y ~ C(clover$Line,Cline),data=clover,...