Description Usage Arguments Details Value Note Author(s) References See Also Examples

View source: R/tenonradialbc.R

Routine `tenonradialbc`

performs bias correction of the nonradial Russell input- or output-based measure of technical efficiency, computes bias and constructs confidence intervals via bootstrapping techniques.

1 2 3 4 5 6 |

`formula` |
an object of class “formula” (or one that can be coerced to that class): a symbolic description of the model. The details of model specification are given under ‘Details’. |

`data` |
an optional data frame containing the variables in the model. If not found in data, the variables are taken from environment ( |

`subset` |
an optional vector specifying a subset of observations for which technical efficiency is to be computed. |

`rts` |
character or numeric. string: first letter of the word “c” for constant, “n” for non-increasing, or “v” for variable returns to scale assumption. numeric: 3 for constant, 2 for non-increasing, or 1 for variable returns to scale assumption. |

`base` |
character or numeric. string: first letter of the word “o” for computing output-based or “i” for computing input-based technical efficiency measure. string: 2 for computing output-based or 1 for computing input-based technical efficiency measure |

`ref` |
an object of class “formula” (or one that can be coerced to that class): a symbolic description of inputs and outputs that are used to define the technology reference set. The details of technology reference set specification are given under ‘Details’. If reference is not provided, the technical efficiency measures for data points are computed relative to technology based on data points themselves. |

`data.ref` |
an optional data frame containing the variables in the technology reference set. If not found in |

`subset.ref` |
an optional vector specifying a subset of observations to define the technology reference set. |

`smoothed` |
logical. If TRUE, the reference set is bootstrapped with smoothing; if FALSE, the reference set is bootstrapped with subsampling. |

`homogeneous` |
logical. Relevant if |

`kappa` |
relevant if |

`reps` |
specifies the number of bootstrap replications to be performed. The default is 999. The minimum is 100. Adequate estimates of confidence intervals using bias-corrected methods typically require 1,000 or more replications. |

`level` |
sets confidence level for confidence intervals; default is |

`show.progress` |
logical. Relevant if |

`print.level` |
numeric. 0 - nothing is printed; 1 - print summary of the model and data. 2 - print summary of technical efficiency measures. 3 - print estimation results observation by observation. Default is 1. |

`seed` |
numeric. The seed (for replication purposes). |

Routine `tenonradialbc`

performs bias correction of the nonradial Russell input- or output-based measure of technical efficiency, computes bias and constructs confidence intervals via bootstrapping techniques (see Badunenko and Mozharovskyi (2020), doi: 10.1080/01605682.2019.1599778).

Models for `tenonradialbc`

are specified symbolically. A typical model has the form `outputs ~ inputs`

, where `outputs`

(`inputs`

) is a series of (numeric) terms which specifies outputs (inputs). The same goes for reference set. Refer to the examples.

Results can be summarized using `summary.npsf`

.

`tenonradialbc`

returns a list of class `npsf`

containing the following elements:

`K` |
numeric: number of data points. |

`M` |
numeric: number of outputs. |

`N` |
numeric: number of inputs. |

`rts` |
string: RTS assumption. |

`base` |
string: base for efficiency measurement. |

`reps` |
numeric: number of bootstrap replications. |

`level` |
numeric: confidence level for confidence intervals. |

`te` |
numeric: radial measure (Russell) of technical efficiency. |

`tebc` |
numeric: bias-corrected radial measures of technical efficiency. |

`biasboot` |
numeric: bootstrap bias estimate for the original radial measures of technical efficiency. |

`varboot` |
numeric: bootstrap variance estimate for the radial measures of technical efficiency. |

`biassqvar` |
numeric: one-third of the ratio of bias squared to variance for radial measures of technical efficiency. |

`realreps` |
numeric: actual number of replications used for statistical inference. |

`telow` |
numeric: lower bound estimate for radial measures of technical efficiency. |

`teupp` |
numeric: upper bound estimate for radial measures of technical efficiency. |

`teboot` |
numeric: |

`esample` |
logical: returns TRUE if the observation in user supplied data is in the estimation subsample and FALSE otherwise. |

Before specifying option `homogeneous`

it is advised to preform the test of independence (see `nptestind`

). Routine `nptestrts`

may help deciding regarding option `rts`

.

Results can be summarized using `summary.npsf`

.

Oleg Badunenko <oleg.badunenko@brunel.ac.uk>, Pavlo Mozharovskyi <pavlo.mozharovskyi@telecom-paris.fr>

Badunenko, O. and Mozharovskyi, P. (2020), Statistical inference for the Russell measure of technical efficiency, *Journal of the Operational Research Society*, **71**3, 517–527, doi: 10.1080/01605682.2019.1599778

Färe, R., Grosskopf, S. and Lovell, C. A. K. (1994), *Production Frontiers*, Cambridge U.K.: Cambridge University Press, doi: 10.1017/CBO9780511551710

`teradial`

, `tenonradial`

, `teradialbc`

, `nptestrts`

, `nptestind`

, `sf`

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 | ```
## Not run:
data( ccr81 )
head( ccr81 )
# Subsampling
t9 <- tenonradialbc(y1 + y2 + y3 ~ x1 + x2 + x3 + x4 + x5, data = ccr81,
ref = NULL, data.ref = NULL, subset.ref = NULL,
rts = "v", base = "i",
homogeneous = FALSE, smoothed = TRUE, kappa = .6,
reps = 999, level = 95,
print.level = 1, show.progress = TRUE, seed = NULL)
# display the results
cbind(te = t9$te, telow = t9$telow, tebc = t9$tebc, teupp = t9$teupp,
biasboot = t9$biasboot, varboot = t9$varboot, biassqvar = t9$biassqvar)
# Smoothing
t10 <- tenonradialbc(y1 + y2 + y3 ~ x1 + x2 + x3 + x4 + x5, data = ccr81,
ref = NULL, data.ref = NULL, subset.ref = NULL,
rts = "v", base = "i",
homogeneous = TRUE, smoothed = TRUE, kappa = .6,
reps = 999, level = 95,
print.level = 1, show.progress = TRUE, seed = NULL)
# display the results
cbind(te = t10$te, telow = t10$telow, tebc = t10$tebc, teupp = t10$teupp,
biasboot = t10$biasboot, varboot = t10$varboot, biassqvar = t10$biassqvar)
## End(Not run)
``` |

Embedding an R snippet on your website

Add the following code to your website.

For more information on customizing the embed code, read Embedding Snippets.