这是我当前的代码:
dat = data.frame(matrix(c(rnorm(20), rep(c(1:5),10), rep(seq(1,4),5)), ncol=3))
colnames(dat)<- c('var1','var2','var3')
ggplot(dat, aes(y=var1,x=var2,colour=var3)) +
geom_point() + geom_smooth(method='aov',se=F)
我想为每种颜色显示 aov 线。我怎样才能做到这一点 ?
此外,我想提取summary(aov())
文本文件中显示的行。并且还提取这些假设的测试aov()
......换句话说,我想提取fligner.test(dat$var1~var2)
andfligner.test(dat$var1~var3)
和shapiro.test(residuals(aov(var1~var2*var3)))
.
有没有办法通过利用 ggplot 包来完成所有这些工作,还是我必须通过 for 循环来操作所有数据?