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)不对


雷达卡





京公网安备 11010802022788号







