R 显示参数'data1[, i]'的种类(list)不对，求大神帮忙

16472

收藏 2017-03-15

drawbiomass <- function(...){
  i=2:5
  par(mfrow = c(2,2))
  data1[,1] <- factor( data1[,1], ordered = TRUE)
  fit <- aov(data1[, i] ~ data1[,1], data = data1)
  out2 <- duncan.test(fit,'data1[,1]',main = 'Aboveground biomass of Alfalfa')
  c <- out2$groups[order(out2$groups[,1]),]
  # draw pic of significance and error bar
  bar.group(c,ylim = c(0,4500), border='red',xaxt = 'n',yaxt = 'n')
  bargraph.CI(x.factor = data1[,1], response = data1[,i], data = data1,font = 2,bty = 2,add = TRUE)
}
drawbiomass(data1)

Error in model.frame.default(formula = data1[, i] ~ data1[, 1], data = data1,  :
  参数'data1[, i]'的种类(list)不对

扫码加我拉你入群

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

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

全部回复

bbslover

2017-3-15 20:23:57

data1是data.frame 还是矩阵？ data.frame不行的话，i直接换成2:5试试

扫码加我拉你入群

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

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

ａ学无止境

2017-3-16 09:01:44

bbslover 发表于 2017-3-15 20:23
data1是data.frame 还是矩阵？ data.frame不行的话，i直接换成2:5试试

data1 是数据框，i = 2:5，我试了下，也不行
TREAT 2013          2014          2015          2016
AN0P0 1651.88    2173.297 4049.292 783.6149
AN0P0 2222.004 2502.421 3768.736 1894.5498
AN0P0 2750.829 4392.023 4329.848 2006.5715
AN0P17 3247.651 4585.234 3902.159 1072.5105
AN0P17 3200.106 3574.604 3897.366 2604.3048
AN0P17 3005.354 3799.958 3906.952 1686.2499
AN0P34 3335.929 3932.378 4170.743 1549.9063
AN0P34 3258.475 4085.76    4341.844 2458.9928
AN0P34 3758.179 5483.095 4256.294 1168.6818
AN70P0 2252.605 2682.849 4092.044 2340.3207
AN70P0 2574.198 2877.207 3168.03 1720.8272
AN70P0 2799.737 3440.817 4507.321 1900.9414
AN70P17 2989.078 4028.316 4601.738 967.2613
AN70P17 3120.712 4303.842 4601.738 2963.3739
AN70P17 3943.287 4037.512 4601.738 1578.914
AN70P34 4271.586 5045.99 2645.101 1552.9302
AN70P34 2915.064 4920.284 3489.496 2446.5258
AN70P34 4278.096 4423.874 3832.361 2226.9163

扫码加我拉你入群

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

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

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

分享

扫码加好友，拉您进群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群