Nothing
summary.frontier <- function( object, extraPar = FALSE, effic = FALSE,
logDepVar = TRUE, effMinusU = farrell, farrell = TRUE, ... ) {
# check argument 'extraPar'
if( length( extraPar ) != 1 || !is.logical( extraPar[1] ) ) {
stop( "argument 'extraPar' must be a single logical value" )
}
# save variable 'logDepVar'
object$logDepVar <- logDepVar
# calculate efficiency estimates
object$effic <- efficiencies( object, logDepVar = logDepVar,
minusU = effMinusU )
# matrix of OLS estimates, their standard errors, t-values and P-values
olsParam <- matrix( NA, length( object$olsParam ) , 4 )
rownames( olsParam ) <- names( object$olsParam )
colnames( olsParam ) <- c( "Estimate", "Std. Error", "t value",
"Pr(>|t|)" )
olsParam[ , 1 ] <- object$olsParam
olsParam[ , 2 ] <- c( object$olsStdEr, NA )
olsParam[ , 3 ] <- olsParam[ , 1 ] / olsParam[ , 2 ]
df <- object$nob - length( object$olsParam )
olsParam[ , 4 ] <- 2 * pt( abs( olsParam[ , 3 ] ), df, lower.tail = FALSE )
object$olsParam <- olsParam
# matrix of ML estimates, their standard errors, t-values and P-values
mlePar <- coef( object, extraPar = extraPar )
mleParam <- matrix( NA, nrow = length( mlePar ), ncol = 4 )
rownames( mleParam ) <- names( mlePar )
colnames( mleParam ) <- c( "Estimate", "Std. Error", "z value",
"Pr(>|z|)" )
mleParam[ , 1 ] <- mlePar
mleParam[ , 2 ] <- diag( vcov( object, extraPar = extraPar ) )^0.5
mleParam[ , 3 ] <- mleParam[ , 1 ] / mleParam[ , 2 ]
df <- object$nob - length( object$mleParam )
mleParam[ , 4 ] <- 2 * pnorm( abs( mleParam[ , 3 ] ), lower.tail = FALSE )
object$mleParam <- mleParam
object$printEffic <- effic
if( ncol( object$effic ) > 1 ) {
object$efficYearMeans <- rep( NA, object$nt )
resid <- residuals( object )
for( i in 1:object$nt ) {
object$efficYearMeans[ i ] <-
mean( object$effic[ !is.na( resid[ , i ] ), i ] )
}
names( object$efficYearMeans ) <- colnames( object$effic )
}
if( object$modelType == 1 && !object$timeEffect ) {
object$efficMean <- mean( object$effic )
} else {
object$efficMean <- mean( object$effic[ !is.na( residuals( object ) ) ] )
}
class( object ) <- "summary.frontier"
return( object )
}
Any scripts or data that you put into this service are public.
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.