Nothing
library( frontier )
suppressWarnings( RNGversion( "3.5.3" ) )
## *****************************
## Testing front41WriteInput
data( front41Data )
front41Data$logOutput <- log( front41Data$output )
front41Data$logCapital <- log( front41Data$capital )
front41Data$logLabour <- log( front41Data$labour )
front41Ins <- front41WriteInput( front41Data, "firm", yName = "logOutput",
xNames = c( "logCapital", "logLabour" ), path = tempdir() )
print( front41Ins[ names( front41Ins ) != "path" ] )
print( readLines( file.path( tempdir(), "front41.ins" ) ) )
print( readLines( file.path( tempdir(), "front41.dta" ) ) )
print( readLines( file.path( tempdir(), "front41.000" ) ) )
# irregular firm (cross section) identifier
set.seed( 20061705 )
front41Data$firm <- sample( c( 1:( nrow( front41Data ) + 20 ) ) )[ 1:nrow( front41Data ) ]
front41Ins <- front41WriteInput( front41Data, "firm", yName = "logOutput",
xNames = c( "logCapital", "logLabour" ), path = tempdir() )
print( front41Ins[ names( front41Ins ) != "path" ] )
print( readLines( file.path( tempdir(), "front41.ins" ) ) )
print( readLines( file.path( tempdir(), "front41.dta" ) ) )
print( readLines( file.path( tempdir(), "front41.000" ) ) )
## *****************************
## Testing front41ReadOutput
outFile <- system.file( "front41/EG1.OUT", package = "frontier" )
sfa <- front41ReadOutput( outFile )
print( coef( sfa, which = "OLS" ) )
print( coef( sfa, which = "GRID" ) )
print( coef( sfa ) )
print( summary( sfa ) )
print( summary( sfa ), efficiencies = TRUE )
print( coef( summary( sfa ), which = "OLS" ) )
print( coef( summary( sfa ), which = "GRID" ) )
print( coef( summary( sfa ) ) )
print( vcov( sfa ) )
print( sfa )
print( sfa, efficiencies = TRUE )
class( sfa ) <- NULL
print( sfa )
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.