# NOT RUN {
# just a stupid simple example
snqProfitFixEla( matrix(1:6/6,3,2 ), matrix(4:1/4,2 ), c(1,1,1), c(1,1),
c(0.4,0.3,0.3) )
# now with real data
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
germanFarms$time <- c( 0:19 )
priceNames <- c( "pOutput", "pVarInput", "pLabor" )
quantNames <- c( "qOutput", "qVarInput", "qLabor" )
fixNames <- c( "land", "time" )
estResult <- snqProfitEst( priceNames, quantNames, fixNames, data=germanFarms )
estResult$fixEla # price elasticities at mean quantities of netputs
# and fixed factors
# fixed factor elasticities at the last observation (1994/95)
snqProfitFixEla( estResult$coef$delta, estResult$coef$gamma,
estResult$data[ 20, quantNames ], estResult$data[ 20, fixNames ],
estResult$weights, estResult$scalingFactors )
}
# }
Run the code above in your browser using DataLab