marginal | R Documentation |
This function returns marginal effects of the inefficiency drivers
from classic or latent class stochastic frontier models estimated with
sfacross
or lcmcross.
## S3 method for class 'sfacross' marginal(object, ...) ## S3 method for class 'lcmcross' marginal(object, ...)
object |
A classic or latent class stochastic frontier model returned
by |
... |
Currently ignored. |
marginal
operates in the presence of exogenous variables that
explain inefficiency, namely the inefficiency drivers (uhet = ~ Z_u or muhet = ~
Z_{mu}).
Two components are computed for each variable: the marginal effects on the expected inefficiency (\frac{\partial E[u]}{\partial Z_{(m)u}}) and the marginal effects on the variance of inefficiency (\frac{\partial V[u]}{\partial Z_{(m)u}}).
The model also allows the Wang (2002) parametrization of μ and σ_u^2 by the same vector of exogenous variables. This double parameterization accounts for non-monotonic relationships between the inefficiency and its drivers.
marginal
returns a data frame containing the marginal effects of the Z_u variables
on the expected inefficiency (each variable has the prefix "Eu_"
)
and on the variance of the inefficiency (each variable has the prefix
"Vu_"
) is returned.
In the case of the latent class model (LCM), each variable terminates with
"_c#"
where "#"
is the class number.
K Hervé Dakpo, Yann Desjeux and Laure Latruffe
Wang, H.J. 2002. Heteroscedasticity and non-monotonic efficiency effects of a stochastic frontier model. Journal of Productivity Analysis, 18:241–253.
sfacross
, for the stochastic frontier analysis model fitting
function.
lcmcross
, for the latent class stochastic frontier analysis
model fitting function.
## Using data on fossil fuel fired steam electric power generation plants in the U.S. # Translog SFA (cost function) truncated normal with scaling property tl_u_ts <- sfacross(formula = log(tc/wf) ~ log(y) + I(1/2 * (log(y))^2) + log(wl/wf) + log(wk/wf) + I(1/2 * (log(wl/wf))^2) + I(1/2 * (log(wk/wf))^2) + I(log(wl/wf) * log(wk/wf)) + I(log(y) * log(wl/wf)) + I(log(y) * log(wk/wf)), udist = "tnormal", muhet = ~ regu + wl, uhet = ~ regu + wl, data = utility, S = -1, scaling = TRUE, method = "mla") marg.tl_u_ts <- marginal(tl_u_ts) summary(marg.tl_u_ts) ## Using data on eighty-two countries production (DGP) # LCM Cobb Douglas (production function) half normal distribution cb_2c_h <- lcmcross(formula = ly ~ lk + ll + yr, udist = "hnormal", data = worldprod, uhet = ~ initStat + h, S = 1, method = "mla") marg.cb_2c_h <- marginal(cb_2c_h) summary(marg.cb_2c_h)
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.