我正在学习r为我的项目编码,并且由于这一错误,有几天我没有进步:
model.frame.default中的错误(公式=〜out + K + L,数据= df, warnings = FALSE):可变长度不同(在“(警告)”中找到)
control = list( maxit = 1000, reltol = 1e-5 )
,但没有用。nls.lm
在参数上设置边界可以有所帮助,但无济于事。在我面临最后期限时,有人可以帮助我吗?谢谢您。
列表项:
#Generate data frame with regressor and AR(1) error. The error term is
set.seed(122)
#\eps_t = 0.5 * \eps_{t-1} + v_t
df <- data.frame(K = rnorm(100,0,1), L=rnorm(100,0,1), err=filter(rnorm(100,0,1)/5,filter=0.3))
gamm=1; delt=0.4; rh=0.5; nu1=0; error=df$err
df$out=(gamm*(delt*K^(-rh)+(1-delt)*L^(-rh))^(-nu1/rh)+error)
library(minpack.lm)
ceS<-nls(out~gamma * (delta * K^(-rho) + (1-delta) * L^(-rho))^(-nu/rho), data=df, start = c(gamma=0.5, delta=0.6, rho=0.25, nu=0), trace = TRUE, control= list( maxit = 1000, reltol = 1e-5 ),warnings = FALSE )
summary(ceS)