View source: R/hypothesis_tests.R
| poolability_test | R Documentation |
Tests the null hypothesis that all groups share a common frontier (i.e., the metafrontier coincides with all group frontiers) against the alternative that group-specific frontiers differ. For SFA-based metafrontiers a likelihood ratio test is used; for DEA-based metafrontiers a permutation test is used.
poolability_test(object, B = 199, seed = NULL, ...)
object |
a fitted |
B |
integer. Number of permutation replicates for the DEA permutation test (default 199). Ignored for SFA objects. |
seed |
integer or |
... |
additional arguments (currently unused). |
Likelihood ratio test (SFA). The LR statistic is:
LR = -2 [LL_{pooled} - \sum_j LL_j]
where LL_{pooled} is the log-likelihood of the pooled (single
frontier) model and LL_j are the group-specific
log-likelihoods. Under H0, the statistic follows a chi-squared
distribution with degrees of freedom equal to
df = k_{groups} - k_{pooled}, where k_{groups} is the
total number of parameters across all group-specific models and
k_{pooled} is the number of parameters in the pooled model.
For J groups each with p frontier parameters plus
distributional parameters, this equals
(J - 1) \times p_{total} where p_{total} includes
frontier coefficients, \sigma_v, and \sigma_u
(and \mu for truncated-normal). This test requires a
likelihood and is therefore only available for SFA-based
metafrontiers.
Permutation test (DEA). DEA has no likelihood, so the
poolability hypothesis is assessed by a permutation test. Under the
null of a single pooled technology, group labels are exchangeable:
reassigning observations to groups at random should not
systematically change the distance between the group frontiers and
the metafrontier. The observed statistic is the mean technology gap,
S_{obs} = \mathrm{mean}(1 - TGR_i), and its null distribution
is approximated by refitting the metafrontier on B random
permutations of the group labels. The p-value is
(1 + \#\{S_b \ge S_{obs}\}) / (B + 1), following the
aggregate-efficiency inference logic of Simar and Zelenyuk (2007).
The smoothed subsampling approach of Kneip, Simar, and Wilson (2016)
is the asymptotically rigorous alternative for testing hypotheses in
nonparametric production models; the permutation test offered here
is a computationally simple approximation. The default B = 199
is a pragmatic choice; p-values have resolution 1/(B + 1), so
increase B for finer resolution.
A list of class "htest" with components:
the test statistic (LR statistic for SFA; the
mean technology gap, \bar{S} = \mathrm{mean}(1 - TGR),
for DEA)
degrees of freedom (SFA) or the effective number of permutation replicates (DEA)
p-value of the test
description of the test
Simar, L. and Zelenyuk, V. (2007). Statistical inference for aggregates of Farrell-type efficiencies. Journal of Applied Econometrics, 22(7), 1367–1394. \Sexpr[results=rd]{tools:::Rd_expr_doi("10.1002/jae.991")}
Kneip, A., Simar, L. and Wilson, P.W. (2016). Testing hypotheses in nonparametric models of production. Journal of Business & Economic Statistics, 34(3), 435–456. \Sexpr[results=rd]{tools:::Rd_expr_doi("10.1080/07350015.2015.1049747")}
set.seed(42)
sim <- simulate_metafrontier(n_groups = 2, n_per_group = 200,
tech_gap = c(0, 0.5))
fit <- metafrontier(log_y ~ log_x1 + log_x2,
data = sim$data, group = "group")
poolability_test(fit)
# DEA permutation test
fit_dea <- metafrontier(log_y ~ log_x1 + log_x2,
data = sim$data, group = "group",
method = "dea")
poolability_test(fit_dea, B = 99, seed = 1)
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.