taojiea1014的博客

regrowth <- read.table("c:\\temp\\ipomopsis.txt",header=T)attach(regrowth)names(regrowth)plot(Root,Fruit,pch=16,col=c("blue","red")[as.numeric(Grazing)])levels(Grazing)abline(lm(Fruit[Grazing=="Gra...

Ancovacontrasts <- read.table("c:\\temp\\Ancovacontrasts.txt",header=T)attach(Ancovacontrasts)names(Ancovacontrasts)lm(weight[sex=="male"]~age[sex=="male"])lm(weight~age,subset=(sex=="female"))lm(w...

Gain <- read.table("c:\\temp\\Gain.txt",header=T)attach(Gain)names(Gain)m1 <- lm(Weight~Sex*Age*Genotype)summary(m1)m2 <- step(m1)summary(m2)newGenotype <- Genotypelevels(newGenotype)level...

summary.aov(lm(weight~sex*age))summary.aov(lm(weight~age*sex))summary.aov(lm(Fruit~Grazing*Root))summary.aov(lm(Fruit~Root*Grazing))tapply(Root,Grazing, mean)summary(lm(Fruit~Root*Grazing))summary(lm(...