我为CES生产功能编写了以下代码:
df <- data.frame(K = CESdata$Kap, L=CESdata$Lab, Q=CESdata$Out)
library(minpack.lm)
ceS<-nlsLM(df$Q~gamma * (delta*K^(-rho) + (1-delta)*L^(-rho))^(-nu/rho), data=df,
start=list(gamma=0.07, delta=0.5, rho=0.25, nu=01.5),
trace = TRUE, control = list( maxiter=1000, reltol = 1e-5, minFactor = 1/1024,printEval = FALSE, warnOnly = FALSE),
algorithm="LM", lower=c(0,0,-1,0), upper=c(Inf,1,Inf,Inf))
summary(ceS)
我尝试使用不同的初始参数,但是报告的结果并没有太大不同:Reason stopped: Number of calls to 'fcn' has reached or exceeded 'maxfev' == 500
。
附件是任何协助的数据。
Data