library(frontier)
fit<-sfa(log(N)~log(D),data=mydata)
fit
summary(fit)
拟合结果如下:
fit
Call:
sfa(formula = log(N) ~ log(D), data = mydata)
Maximum likelihood estimates
(Intercept) log(D) sigmaSq gamma
8.730e+00 -4.214e-01 2.079e-03 2.015e-05
> summary(fit)
Error Components Frontier (see Battese & Coelli 1992)
Inefficiency decreases the endogenous variable (as in a production function)
The dependent variable is logged
Iterative ML estimation terminated after 47 iterations:
log likelihood values and parameters of two successive iterations
are within the tolerance limit