R语言中 如何实现不同因素在公式中循环，见代码

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

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

立即领取

感谢您参与论坛问题回答

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

+2 论坛币

1. # 数据
   
2. head(npk)
   
3. 
   
4. # 问题，我想写个循环，分别对N，P，K三个因素单独进行方差分析
   
5. #类似
   
6. aov(yield ~ N, data=npk)
   
7. aov(yield ~ P, data=npk)
   
8. aov(yield ~ K, data=npk)
   
9. 
   
10. # 实施：
    
11. c <- c("N","P","K")
    
12. for (i in 1:3) print(c[i])
    
13. # [1] "N"
    
14. # [1] "P"
    
15. # [1] "K"
    
16. 
    
17. # 放到循环中
    
18. for (i in 1:3) aov(yield ~ c[i] ,data=npk)
    
19. # Error in model.frame.default(formula = yield ~ c[i], data = npk, drop.unused.levels = TRUE) :
    
20. #   变数的长度不一样('c[i]')

复制代码

看到报错后，一脸懵逼相，求大师帮助！

扫码加我 拉你入群

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

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

分享0 收藏0 回帖

关键词：如何实现 R语言 yield Data 方差分析

for (i in 1:3){
  aov(paste0('yield ~ ', c[i]), data=npk)
}

Eric1028 发表于 2017-4-12 20:38
for (i in 1:3){
  aov(paste0('yield ~ ', c), data=npk)
}

谢谢你的回答，不过没有解决，报错如下：

1. > for (i in 1:3){
   
2. +     aov(paste0('yield ~ ', c[i]), data=npk)
   
3. +   }
   
4. Error: $ operator is invalid for atomic vectors

复制代码

data(npk)
c <- c("N","P","K")
mod <- list(length = 3)
for (i in 1:3) {
    mod[[i]] <- aov(formula(paste0('yield ~ ', c[i])), data = npk)
}
mod