View source: R/snqProfitImposeConvexity.R
snqProfitImposeConvexity | R Documentation |
Imposing Convexity on a Symmetric Normalized Quadratic (SNQ) Profit function.
snqProfitImposeConvexity( estResult, rankReduction = 0, start = 10, optimMethod = "BFGS", control = list( maxit=5000 ), stErMethod = "none", nRep = 1000, verbose = 0 )
estResult |
object returned by |
rankReduction |
an integer specifying the reduction of the rank of the β matrix. |
start |
starting values of the triangular Cholesky matrix. |
optimMethod |
method to be used by |
control |
list of control parameters passed to |
stErMethod |
method to compute standard errors, either 'none', 'resample', 'jackknife' or 'coefSim' (see details). |
nRep |
number of replications to compute the standard errors
if |
verbose |
an integer idicating the verbose level. |
The procedure proposed by Koebel, Falk and Laisney (2000, 2003) is
applied to impose convexity in prices on an estimated symmetric
normalized quadratic (SNQ) profit function.
The standard errors of the restricted coefficients can be either
calculated by
bootstrap resampling ('resampling'),
jackknife ('jacknife') or
by simulating the distribution of the unrestricted coefficients
using its variance covariance matrix ('coefSim').
a list of class snqProfitImposeConvexity
containing the same
objects as an object of class snqProfitEst
and additionally the objects:
mindist |
object returned by |
sim |
results of the simulation to obtain the standard errors of the estimated coefficients. |
Arne Henningsen
Koebel, B., M. Falk and F. Laisney (2000), Imposing and Testing Curvature Conditions on a Box-Cox Cost Function. Discussion Paper No. 00-70, ZEW, Mannheim, https://madoc.bib.uni-mannheim.de/515/1/dp0070.pdf.
Koebel, B., M. Falk and F. Laisney (2003), Imposing and Testing Curvature Conditions on a Box-Cox Cost Function. Journal of Business and Economic Statistics, 21, p. 319-335.
snqProfitEst
.
if( requireNamespace( 'micEcon', quietly = TRUE ) ) { data( germanFarms, package = "micEcon" ) germanFarms$qOutput <- germanFarms$vOutput / germanFarms$pOutput germanFarms$qVarInput <- -germanFarms$vVarInput / germanFarms$pVarInput germanFarms$qLabor <- -germanFarms$qLabor priceNames <- c( "pOutput", "pVarInput", "pLabor" ) quantNames <- c( "qOutput", "qVarInput", "qLabor" ) estResult <- snqProfitEst( priceNames, quantNames, "land", data = germanFarms ) estResult # Note: it is NOT convex in netput prices estResultConvex <- snqProfitImposeConvexity( estResult ) estResultConvex # now it is convex }
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.