> reg.len <- lm(chao1.ave ~ lg.std.len, b.div) # b.div is my data frame imported from a CSV file
> reg.len
Call:
lm(formula = chao1.ave ~ lg.std.len, data = b.div)
Coefficients:
(Intercept) lg.std.len
282.4 -115.7
> newx <- seq(0.6, 1.4, 0.01)
> prd.len <- predict(reg.len, newdata=data.frame(x=newx), interval="confidence", level=0.90, type="response")
Error in eval(expr, envir, enclos) : object 'lg.std.len' not found
我试过像这样做 lm: lm(b.div$chao1.ave ~ b.div$lg.std.len)
,但是随后predict()
警告说 newdata 和 variables 的长度不同。所以,我尝试了上面的方法,现在predict()
给出一个错误,说它无法识别对象。请问怎么修?