https://github.com/javierbarbero/DataEnvelopmentAnalysis.jl
Tip revision: a9ad709c106fcaa6c21d03edad7c83fa02a3fbf0 authored by Javier Barbero on 14 September 2022, 17:40:23 UTC
Version 0.8.0
Version 0.8.0
Tip revision: a9ad709
radialboot.md
```@meta
CurrentModule = DataEnvelopmentAnalysis
```
# Bootstrap Radial Model
The bootstrap radial DEA model (Simar and Wilson, 1998) can be calculated with the `deaboot` function, indicating the number of bootstrap replications in the `nreps` parameter. The other parameters work the same as in the radial DEA model.
A random number generator can be specified in the `rng` parameter for reproducibility.
```@example radialboot
using DataEnvelopmentAnalysis
using StableRNGs
X = [2, 4, 3, 5, 6]
Y = [1, 2, 3, 4, 5]
ioboot = deaboot(X, Y, orient = :Input, rts = :VRS, nreps = 200, rng = StableRNG(1234567))
```
!!! warning "Number of bootstrap replications"
The example above uses 200 bootstrap replications for illustrative purposes. In practice, at least 1000 replications are recommended.
Bias-corrected efficiency scores are returned with the `efficiency` function:
```@example radialboot
efficiency(ioboot)
```
The bias, calculated as the difference between the reference efficiency score and the bias-corrected efficiency score, is returned with the `bias` function:
```@example radialboot
bias(ioboot)
```
Confidence intervals at the $95\%$, or any other desired level, are calculated with the `confint` function:
```@example radialboot
confint(ioboot, level = 0.95)
```
The optimal bandwidth computed for the model is returned with the `bandwidth` function:
```@example radialboot
bandwidth(ioboot)
```
### deaboot Function Documentation
```@docs
deaboot
bias(::BootstrapRadialDEAModel)
confint(::BootstrapRadialDEAModel)
bandwidth(::BootstrapRadialDEAModel)
```