如何批量进行单因素方差分析

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

送您一个全额奖学金名额~ !

立即领取

感谢您参与论坛问题回答

经管之家送您两个论坛币！

+2 论坛币

各位老师好，我想用for循环量对多个变量进行单因素的方差分析，但循环体是字符型的变量，而aov(x~A,data=mydata)中的x肯定不能是字符型的，请问应该怎么转化呢？谢谢。我拿iris作为案例，试了下，果然不成功，求助各位老师，谢谢

1. head(iris)
   
2. fit.out=list()
   
3. for(i in c('Sepal.Length',"Sepal.Width","Petal.Length","Petal.Width")){
   
4.   fit<-aov(i~Species,data=iris)
   
5.   pair.test<-TukeyHSD(fit)
   
6.   fit.out<-c(fit.out,list(pair.test[[1]]))
   
7. }
   
8. fit.out

复制代码

扫码加我 拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

分享0 收藏0 回帖

我知道可以用向量化操作：
lapply(iris[1:4],function(x) TukeyHSD(aov(x~Species,data=iris))[[1]])
但是此类问题我碰到过好多次了，有时感觉lapply也不太好用，我就是想知道如果用for循环的话，怎么弄呢？？？有的强迫症哈，希望各位老师帮忙。

千车麦穗 发表于 2020-4-18 11:36
各位老师好，我想用for循环量对多个变量进行单因素的方差分析，但循环体是字符型的变量，而aov(x~A,data=my ...

第四句里的i改成get(i)就可以了

drunkfish69 发表于 2020-4-18 13:39
第四句里的i改成get(i)就可以了

嗯嗯，果然好使，谢谢…。ps：我之前以为assign和get得配套使用呢，看来思维束缚了

千车麦穗 发表于 2020-4-18 13:43
嗯嗯，果然好使，谢谢…。ps：我之前以为assign和get得配套使用呢，看来思维束缚了

请问如果使用LSD的比较方法，怎么返回直接标有显著性关系的那组数而不是LSD数值。谢谢

drunkfish69 发表于 2020-4-18 13:39
第四句里的i改成get(i)就可以了

请问如果使用LSD的比较方法，怎么返回直接标有显著性关系的那组数而不是LSD数值。谢谢

赵斯特 发表于 2022-1-5 22:59
请问如果使用LSD的比较方法，怎么返回直接标有显著性关系的那组数而不是LSD数值。谢谢

不好意思，我现在用的很少了（毕业了），我建议你新开一个帖子，问下群里的各位老师。确实不好意思了