adea_hierarchical: Selection of an optimal subset of variables for DEA analysis

View source: R/adea_hierarchical.R

adea_hierarchicalR Documentation

Selection of an optimal subset of variables for DEA analysis

Description

The function returns a list of DEA models dropping one variable in each step following ADEA methodology.

Usage

adea_hierarchical(
  input,
  output,
  orientation = c("input", "output"),
  load.orientation = c("inoutput", "input", "output"),
  name = "",
  direction = c("backward", "backward/input", "backward/output"),
  verbose = 0
)

Arguments

input

A matrix or a data frame with the inputs of units to be evaluated, one row for each DMU and one column for each input.

output

A matrix or a data frame with the outputs of units to be evaluated, one row for each DMU and one column for each output.

orientation

Use "input" for input orientation or use "output" for output orientation in DEA model.

load.orientation

It allows the selection of variables to be included in load analysis. Its default value is "inoutput" which means that all input and all output variables will be included. Use "input" or "output" to include only input or output variables in load analysis.

name

An optional descriptive name for the model. It will be shown in print and summary results.

direction

The direction in which the variables go in and out of the model. Until now, only backward option is implemented.

verbose

Use 0 for minimal output, only final model. 1 or more to get detailed information of each step. This option affects only to printed output but not the result.

Details

This procedure provides a list of all DEA models for all nested sets of variables.

In each model, the variable with lower value of load is dropped. In this case, the load of a new model can be lower than a previous one.

See examples for more details.

Value

A list with all computed models. So, you can retrieve each model individually.

See Also

adea_parametric

Examples

# Read data
data('cardealers4')
input <- cardealers4[, 1:2]
output <- cardealers4[, 3:4]

# Compute all dea models in hierarchical way
adea_hierarchical(input, output)
#        Load #Efficient Factors #Inputs #Outputs                  Inputs              Outputs
# 4 0.6666667          2       4       2        2 Employees, Depreciation CarsSold, WorkOrders
# 3 0.9575672          2       3       1        2            Depreciation CarsSold, WorkOrders
# 2 1.0000000          1       2       1        1            Depreciation CarsSold

# Compute again and store result in sol.ah
sol.ah <- adea_hierarchical(input, output)

# Summary the model with 3 variables
summary(sol.ah$models[[3]])
# Model name:
# Orientation is input
# Inputs: Depreciation
# Outputs: CarsSold WorkOrders
# Load: 0.9575672
# Input loads: 1
# Output loads: 1.042433 0.9575672
# Summary of efficiencies:
# Mean        sd      Min.   1st Qu.    Median   3rd Qu.      Max.
# 0.9002232 0.1351949 0.6515044 0.8722527 0.9422250 0.9978982 1.0000000
# Efficiencies:
# Dealer A  Dealer B  Dealer C  Dealer D  Dealer E  Dealer F
# 0.9915929 1.0000000 0.8928571 0.8653846 1.0000000 0.6515044

# Get efficiencies for the model with 3 variables
sol.ah$models[[3]]$eff
# Dealer A  Dealer B  Dealer C  Dealer D  Dealer E  Dealer F
# 0.9915929 1.0000000 0.8928571 0.8653846 1.0000000 0.6515044 

adea documentation built on March 18, 2022, 7:24 p.m.