attach(regrowth)
names(regrowth)
plot(Root,Fruit,pch=16,col=c("blue","red")[as.numeric(Grazing)])
levels(Grazing)
abline(lm(Fruit[Grazing=="Grazed"]~Root[Grazing=="Grazed"]),col="blue")
abline(lm(Fruit[Grazing=="Ungrazed"]~Root[Grazing=="Ungrazed"]),col="red")
tapply(Fruit,Grazing, mean)
t.test(Fruit~Grazing)
sum(Root);sum(Rootˆ2)
sum(Fruit);sum(Fruitˆ2)
sum(Root*Fruit)
sum(Root[Grazing=="Grazed"]);sum(Root[Grazing=="Grazed"]ˆ2)
sum(Root[Grazing=="Ungrazed"]);sum(Root[Grazing=="Ungrazed"]ˆ2)
sum(Fruit[Grazing=="Grazed"]);sum(Fruit[Grazing=="Grazed"]ˆ2)
sum(Fruit[Grazing=="Ungrazed"]);sum(Fruit[Grazing=="Ungrazed"]ˆ2)
sum(Root[Grazing=="Grazed"]*Fruit[Grazing=="Grazed"])
sum(Root[Grazing=="Ungrazed"]*Fruit[Grazing=="Ungrazed"])
ancova <- lm(Fruit~Grazing*Root)
summary(ancova)
anova(ancova)
ancova2 <- update(ancova, ~ . - Grazing:Root)
anova(ancova,ancova2)
ancova3 <- update(ancova2, ~ . - Grazing)
anova(ancova2,ancova3)
summary(ancova2)
anova(ancova2)

本案例使用R语言进行了一次详细的数据分析过程,包括从数据读取到可视化展示、假设检验等多个环节。通过分析不同放牧条件下植物根部与果实生长的关系,展示了如何利用线性模型和方差分析等统计方法来探索变量间的相互作用。
1614

被折叠的 条评论
为什么被折叠?



