我正在拟合随机生产边界以计算公司的效率。我正在使用边境包裹。但是,我有一个错误-
optim中的错误(pars,fn = maxlik,y = y,X = X,sc = sc,粗麻布= TRUE,: 非有限差分值
我的代码是
library(readxl)
d<- read_excel("Downloads/data_spf.xls",sheet = "2000 total")
k<-d$k
l<-d$l
y<-d$y
k2<-k*k
l2<-l*l
lk<-l*k
df<-data.frame(y,k,l,k2,l2,lk)
spf<-sfa(y~k+l+k2+l2+lk,data = df,intercept = TRUE,fun="hnormal",form = "production",method = "Nelder-Mead")