melt: Multiple Empirical Likelihood Tests

Description

Performs multiple empirical likelihood tests. It offers an easy-to-use interface and flexibility in specifying hypotheses and calibration methods, extending the framework to simultaneous inferences. The core computational routines are implemented using the 'Eigen' 'C++' library and 'RcppEigen' interface, with 'OpenMP' for parallel computation. Details of the testing procedures are provided in Kim, MacEachern, and Peruggia (2023) doi:10.1080/10485252.2023.2206919. A companion paper by Kim, MacEachern, and Peruggia (2024) doi:10.18637/jss.v108.i05 is available for further information. This work was supported by the U.S. National Science Foundation under Grants No. SES-1921523 and DMS-2015552.

Author(s)

Maintainer: Eunseop Kim markean@pm.me [copyright holder]

Other contributors:

• Steven MacEachern [contributor, thesis advisor]

• Mario Peruggia [contributor, thesis advisor]

• Pierre Chausse [reviewer]

• Alex Stringer [reviewer]

References

Kim E, MacEachern SN, Peruggia M (2023). “Empirical likelihood for the analysis of experimental designs.” Journal of Nonparametric Statistics, 35(4), 709–732. doi:10.1080/10485252.2023.2206919.

Kim E, MacEachern SN, Peruggia M (2024). “melt: Multiple Empirical Likelihood Tests in R.” Journal of Statistical Software, 108(5), 1–33. doi:10.18637/jss.v108.i05.

See Also

Useful links:

• https://docs.ropensci.org/melt/

• https://github.com/ropensci/melt

• Report bugs at https://github.com/ropensci/melt/issues

CEL class

Description

S4 class for constrained empirical likelihood. It inherits from EL class. Note that the optim slot has constrained optimization results with respect to the parameters, not the Lagrange multiplier.

Details

Let l(\theta) denote minus twice the empirical log-likelihood ratio function. We consider a linear hypothesis of the form

L\theta = r,

where the left-hand-side L is a q by p matrix and the right-hand-side r is a q-dimensional vector. Under some regularity conditions, l(\theta) converges in distribution to \chi^2_q under the constraint of hypothesis, i.e.,

\min_{\theta: L\theta = r} l(\theta) \to_d \chi^2_q .

Minimization of l(\theta) with respect to \theta is computationally expensive since it implicitly involves the evaluation step as described in EL. Further, depending on the form of g(X_i, \theta) and the constraint, the optimization problem can be nonconvex and have multiple local minima. For this reason, the package melt only considers linear hypotheses and performs local minimization of l(\theta) using projected gradient descent method. With the orthogonal projection matrix P and a step size \gamma, the algorithm updates \theta as

\theta^{(k + 1)} \leftarrow \theta^{(k)} - \gamma P \nabla l(\theta^{(k)}),

where \nabla l(\theta^{(k)}) denotes the gradient of l at \theta^{(k)}. The first order optimality condition is P \nabla l(\theta) = 0, which is used as the stopping criterion.

Slots

optim

A list of the following optimization results:

• par A numeric vector of the solution to the constrained optimization problem.

• lambda A numeric vector of the Lagrange multipliers of the dual problem corresponding to par.

• iterations A single integer for the number of iterations performed.

• convergence A single logical for the convergence status.

• cstr A single logical for whether constrained EL optimization is performed or not.

logp

A numeric vector of the log probabilities of the constrained empirical likelihood.

logl

A single numeric of the constrained empirical log-likelihood.

loglr

A single numeric of the constrained empirical log-likelihood ratio.

statistic

A single numeric of minus twice the constrained empirical log-likelihood ratio with an asymptotic chi-square distribution.

df

A single integer for the degrees of freedom of the statistic.

pval

A single numeric for the p-value of the statistic.

nobs

A single integer for the number of observations.

npar

A single integer for the number of parameters.

weights

A numeric vector of the re-scaled weights used for the model fitting.

coefficients

A numeric vector of the maximum empirical likelihood estimates of the parameters.

method

A single character for the method dispatch in internal functions.

data

A numeric matrix of the data for the model fitting.

control

An object of class ControlEL constructed by el_control().

References

Adimari G, Guolo A (2010). “A Note on the Asymptotic Behaviour of Empirical Likelihood Statistics.” Statistical Methods & Applications, 19(4), 463–476. doi:10.1007/s10260-010-0137-9.

Qin J, Lawless J (1995). “Estimating Equations, Empirical Likelihood and Constraints on Parameters.” Canadian Journal of Statistics, 23(2), 145–159. doi:10.2307/3315441.

Examples

showClass("CEL")

ConfregEL class

Description

S4 class for confidence region. It inherits from "matrix".

Slots

estimates

A numeric vector of length two for the parameter estimates.

level

A single numeric for the confidence level required.

cv

A single numeric for the critical value for calibration of empirical likelihood ratio statistic.

pnames

A character vector of length two for the name of parameters.

Examples

showClass("ConfregEL")

ControlEL class

Description

S4 class for computational details of empirical likelihood.

Slots

maxit

A single integer for the maximum number of iterations for the optimization with respect to \theta.

maxit_l

A single integer for the maximum number of iterations for the optimization with respect to \lambda.

tol

A single numeric for the convergence tolerance denoted by \epsilon. The iteration stops when

\|P \nabla l(\theta^{(k)})\| < \epsilon.

tol_l

A single numeric for the relative convergence tolerance denoted by \delta. The iteration stops when

\|\lambda^{(k)} - \lambda^{(k - 1)}\| < \delta\|\lambda^{(k - 1)}\| + \delta^2.

step

A single numeric for the step size \gamma for the projected gradient descent method.

th

A single numeric for the threshold for the negative empirical log-likelihood ratio.

verbose

A single logical for whether to print a message on the convergence status.

keep_data

A single logical for whether to keep the data used for fitting model objects.

nthreads

A single integer for the number of threads for parallel computation via OpenMP (if available).

seed

A single integer for the seed for random number generation.

an

A single numeric representing the scaling factor for adjusted empirical likelihood calibration.

b

A single integer for the number of bootstrap replicates.

m

A single integer for the number of Monte Carlo samples.

Examples

showClass("ControlEL")

EL class

Description

S4 class for empirical likelihood.

Details

Let X_i be independent and identically distributed p-dimensional random variable from an unknown distribution P for i = 1, \dots, n. We assume that P has a positive definite covariance matrix. For a parameter of interest \theta(F) \in {\rm{I\!R}}^p, consider a p-dimensional smooth estimating function g(X_i, \theta) with a moment condition

\textrm{E}[g(X_i, \theta)] = 0.

We assume that there exists an unique \theta_0 that solves the above equation. Given a value of \theta, the (profile) empirical likelihood ratio is defined by

R(\theta) = \max_{p_i}\left\{\prod_{i = 1}^n np_i : \sum_{i = 1}^n p_i g(X_i, \theta) = 0, p_i \geq 0, \sum_{i = 1}^n p_i = 1 \right\}.

The Lagrange multiplier \lambda \equiv \lambda(\theta) of the dual problem leads to

p_i = \frac{1}{n}\frac{1}{1 + \lambda^\top g(X_i, \theta)},

where \lambda solves

\frac{1}{n}\sum_{i = 1}^n \frac{g(X_i, \theta)} {1 + \lambda^\top g(X_i, \theta)} = 0.

Then the empirical log-likelihood ratio is given by

\log R(\theta) = -\sum_{i = 1}^n \log(1 + \lambda^\top g(X_i, \theta)).

This problem can be efficiently solved by the Newton-Raphson method when the zero vector is contained in the interior of the convex hull of \{g(X_i, \theta)\}_{i = 1}^n.

It is known that -2\log R(\theta_0) converges in distribution to \chi^2_p, where \chi^2_p has a chi-square distribution with p degrees of freedom. See the references below for more details.

Slots

optim

A list of the following optimization results:

• par A numeric vector of the specified parameters.

• lambda A numeric vector of the Lagrange multipliers of the dual problem corresponding to par.

• iterations A single integer for the number of iterations performed.

• convergence A single logical for the convergence status.

• cstr A single logical for whether constrained EL optimization is performed or not.

logp

A numeric vector of the log probabilities of the empirical likelihood.

logl

A single numeric of the empirical log-likelihood.

loglr

A single numeric of the empirical log-likelihood ratio.

statistic

A single numeric of minus twice the empirical log-likelihood ratio with an asymptotic chi-square distribution.

df

A single integer for the degrees of freedom of the statistic.

pval

A single numeric for the p-value of the statistic.

nobs

A single integer for the number of observations.

npar

A single integer for the number of parameters.

weights

A numeric vector of the re-scaled weights used for the model fitting.

coefficients

A numeric vector of the maximum empirical likelihood estimates of the parameters.

method

A single character for the method dispatch in internal functions.

data

A numeric matrix of the data for the model fitting.

control

An object of class ControlEL constructed by el_control().

References

Owen A (2001). Empirical Likelihood. Chapman & Hall/CRC. doi:10.1201/9781420036152.

Qin J, Lawless J (1994). “Empirical Likelihood and General Estimating Equations.” The Annals of Statistics, 22(1), 300–325. doi:10.1214/aos/1176325370.

Examples

showClass("EL")

ELD class

Description

S4 class for empirical likelihood displacement. It inherits from "numeric".

Examples

showClass("ELD")

ELMT class

Description

S4 class for empirical likelihood multiple tests.

Slots

estimates

A list of numeric vectors of the estimates of the linear hypotheses.

statistic

A numeric vector of minus twice the (constrained) empirical log-likelihood ratios with asymptotic chi-square distributions.

df

An integer vector of the marginal degrees of freedom of the statistic.

pval

A numeric vector for the multiplicity adjusted p-values.

cv

A single numeric for the multiplicity adjusted critical value.

rhs

A numeric vector for the right-hand sides of the hypotheses.

lhs

A numeric matrix for the left-hand side of the hypotheses.

alpha

A single numeric for the overall significance level.

calibrate

A single character for the calibration method used.

weights

A numeric vector of the re-scaled weights used for the model fitting.

coefficients

A numeric vector of the maximum empirical likelihood estimates of the parameters.

method

A single character for the method dispatch in internal functions.

data

A numeric matrix of the data for the model fitting.

control

An object of class ControlEL constructed by el_control().

Examples

showClass("ELMT")

ELT class

Description

S4 class for empirical likelihood test.

Slots

optim

A list of the following optimization results:

• par A numeric vector of the solution to the (constrained) optimization problem.

• lambda A numeric vector of the Lagrange multipliers of the dual problem corresponding to par.

• iterations A single integer for the number of iterations performed.

• convergence A single logical for the convergence status.

• cstr A single logical for whether constrained EL optimization is performed or not.

logp

A numeric vector of the log probabilities of the (constrained) empirical likelihood.

logl

A single numeric of the (constrained) empirical log-likelihood.

loglr

A single numeric of the (constrained) empirical log-likelihood ratio.

statistic

A single numeric of minus twice the (constrained) empirical log-likelihood ratio with an asymptotic chi-square distribution.

df

A single integer for the chi-square degrees of freedom of the statistic.

pval

A single numeric for the (calibrated) p-value of the statistic.

cv

A single numeric for the critical value.

rhs

A numeric vector for the right-hand side of the hypothesis.

lhs

A numeric matrix for the left-hand side of the hypothesis.

alpha

A single numeric for the significance level.

calibrate

A single character for the calibration method used.

control

An object of class ControlEL constructed by el_control().

Examples

showClass("ELT")

GLM class

Description

S4 class for generalized linear models. It inherits from LM class.

Details

The overall test involves a constrained optimization problem. All the parameters except for the intercept are constrained to zero. The optim slot contains the results. When there is no intercept, all parameters are set to zero, and the results need to be understood in terms of EL class since no constrained optimization is involved. Once the solution is found, the log probabilities (logp) and the (constrained) empirical likelihood values (logl, loglr, statistic) readily follow, along with the degrees of freedom (df) and the p-value (pval). The significance tests for each parameter also involve constrained optimization problems where only one parameter is constrained to zero. The sigTests slot contains the results.

Slots

family

A family object used.

dispersion

A single numeric for the estimated dispersion parameter.

sigTests

A list of the following results of significance tests:

• statistic A numeric vector of minus twice the (constrained) empirical log-likelihood ratios with asymptotic chi-square distributions.

• iterations An integer vector for the number of iterations performed for each parameter.

• convergence A logical vector for the convergence status of each parameter.

• cstr A single logical for whether constrained EL optimization is performed or not.

call

A matched call.

terms

A terms object used.

misc

A list of various outputs obtained from the model fitting process. They are used in other generics and methods.

optim

A list of the following optimization results:

• par A numeric vector of the solution to the (constrained) optimization problem.

• lambda A numeric vector of the Lagrange multipliers of the dual problem corresponding to par.

• iterations A single integer for the number of iterations performed.

• convergence A single logical for the convergence status.

logp

A numeric vector of the log probabilities of the (constrained) empirical likelihood.

logl

A single numeric of the (constrained) empirical log-likelihood.

loglr

A single numeric of the (constrained) empirical log-likelihood ratio.

statistic

A single numeric of minus twice the (constrained) empirical log-likelihood ratio with an asymptotic chi-square distribution.

df

A single integer for the degrees of freedom of the statistic.

pval

A single numeric for the p-value of the statistic.

nobs

A single integer for the number of observations.

npar

A single integer for the number of parameters.

weights

A numeric vector of the re-scaled weights used for the model fitting.

coefficients

A numeric vector of the maximum empirical likelihood estimates of the parameters.

method

A single character for the method dispatch in internal functions.

data

A numeric matrix of the data for the model fitting.

control

An object of class ControlEL constructed by el_control().

Examples

showClass("GLM")

LM class

Description

S4 class for linear models with empirical likelihood. It inherits from CEL class.

Details

The overall test involves a constrained optimization problem. All the parameters except for the intercept are constrained to zero. The optim slot contains the results. When there is no intercept, all parameters are set to zero, and the results need to be understood in terms of EL class since no constrained optimization is involved. Once the solution is found, the log probabilities (logp) and the (constrained) empirical likelihood values (logl, loglr, statistic) readily follow, along with the degrees of freedom (df) and the p-value (pval). The significance tests for each parameter also involve constrained optimization problems where only one parameter is constrained to zero. The sigTests slot contains the results.

Methods (by generic)

• formula(LM): Extracts the symbolic model formula used in el_lm() or el_glm().

Slots

sigTests

A list of the following results of significance tests:

• statistic A numeric vector of minus twice the (constrained) empirical log-likelihood ratios with asymptotic chi-square distributions.

• iterations An integer vector for the number of iterations performed for each parameter.

• convergence A logical vector for the convergence status of each parameter.

call

A matched call.

terms

A terms object used.

misc

A list of various outputs obtained from the model fitting process. They are used in other generics and methods.

optim

A list of the following optimization results:

• par A numeric vector of the solution to the (constrained) optimization problem.

• lambda A numeric vector of the Lagrange multipliers of the dual problem corresponding to par.

• iterations A single integer for the number of iterations performed.

• convergence A single logical for the convergence status.

logp

A numeric vector of the log probabilities of the (constrained) empirical likelihood.

logl

A single numeric of the (constrained) empirical log-likelihood.

loglr

A single numeric of the (constrained) empirical log-likelihood ratio.

statistic

A single numeric of minus twice the (constrained) empirical log-likelihood ratio with an asymptotic chi-square distribution.

df

A single integer for the degrees of freedom of the statistic.

pval

A single numeric for the p-value of the statistic.

nobs

A single integer for the number of observations.

npar

A single integer for the number of parameters.

weights

A numeric vector of the re-scaled weights used for the model fitting.

coefficients

A numeric vector of the maximum empirical likelihood estimates of the parameters.

method

A single character for the method dispatch in internal functions.

data

A numeric matrix of the data for the model fitting.

control

An object of class ControlEL constructed by el_control().

Examples

showClass("LM")

QGLM class

Description

S4 class for generalized linear models with quasi-likelihood methods. It inherits from GLM class.

Details

The overall test involves a constrained optimization problem. All the parameters except for the intercept are constrained to zero. The optim slot contains the results. When there is no intercept, all parameters are set to zero, and the results need to be understood in terms of EL class since no constrained optimization is involved. Once the solution is found, the log probabilities (logp) and the (constrained) empirical likelihood values (logl, loglr, statistic) readily follow, along with the degrees of freedom (df) and the p-value (pval). The significance tests for each parameter also involve constrained optimization problems where only one parameter is constrained to zero. The sigTests slot contains the results.

Slots

family

A family object used.

dispersion

A single numeric for the estimated dispersion parameter.

sigTests

A list of the following results of significance tests:

• statistic A numeric vector of minus twice the (constrained) empirical log-likelihood ratios with asymptotic chi-square distributions.

• iterations An integer vector for the number of iterations performed for each parameter.

• convergence A logical vector for the convergence status of each parameter.

• cstr A single logical for whether constrained EL optimization is performed or not.

call

A matched call.

terms

A terms object used.

misc

A list of various outputs obtained from the model fitting process. They are used in other generics and methods.

optim

A list of the following optimization results:

• par A numeric vector of the solution to the (constrained) optimization problem.

• lambda A numeric vector of the Lagrange multipliers of the dual problem corresponding to par.

• iterations A single integer for the number of iterations performed.

• convergence A single logical for the convergence status.

logp

A numeric vector of the log probabilities of the (constrained) empirical likelihood.

logl

A single numeric of the (constrained) empirical log-likelihood.

loglr

A single numeric of the (constrained) empirical log-likelihood ratio.

statistic

A single numeric of minus twice the (constrained) empirical log-likelihood ratio with an asymptotic chi-square distribution.

df

A single integer for the degrees of freedom of the statistic.

pval

A single numeric for the p-value of the statistic.

nobs

A single integer for the number of observations.

npar

A single integer for the number of parameters.

weights

A numeric vector of the re-scaled weights used for the model fitting.

coefficients

A numeric vector of the maximum empirical likelihood estimates of the parameters.

method

A single character for the method dispatch in internal functions.

data

A numeric matrix of the data for the model fitting.

control

An object of class ControlEL constructed by el_control().

Examples

showClass("QGLM")

SD class

Description

S4 class for standard deviation. It inherits from EL class.

Slots

optim

A list of the following optimization results:

• par A numeric vector of the specified parameters.

• lambda A numeric vector of the Lagrange multipliers of the dual problem corresponding to par.

• iterations A single integer for the number of iterations performed.

• convergence A single logical for the convergence status.

• cstr A single logical for whether constrained EL optimization is performed or not.

logp

A numeric vector of the log probabilities of the empirical likelihood.

logl

A single numeric of the empirical log-likelihood.

loglr

A single numeric of the empirical log-likelihood ratio.

statistic

A single numeric of minus twice the empirical log-likelihood ratio with an asymptotic chi-square distribution.

df

A single integer for the degrees of freedom of the statistic.

pval

A single numeric for the p-value of the statistic.

nobs

A single integer for the number of observations.

npar

A single integer for the number of parameters.

weights

A numeric vector of the re-scaled weights used for the model fitting.

coefficients

A numeric vector of the maximum empirical likelihood estimates of the parameters.

method

A single character for the method dispatch in internal functions.

data

A numeric matrix of the data for the model fitting.

control

An object of class ControlEL constructed by el_control().

Examples

showClass("SD")

SummaryEL class

Description

S4 class for a summary of EL objects.

Slots

optim

A list of the following optimization results:

• par A numeric vector of the specified parameters.

• lambda A numeric vector of the Lagrange multipliers of the dual problem corresponding to par.

• iterations A single integer for the number of iterations performed.

• convergence A single logical for the convergence status.

• cstr A single logical for whether constrained EL optimization is performed or not.

logl

A single numeric of the empirical log-likelihood.

loglr

A single numeric of the empirical log-likelihood ratio.

statistic

A single numeric of minus twice the empirical log-likelihood ratio with an asymptotic chi-square distribution.

df

A single integer for the degrees of freedom of the statistic.

pval

A single numeric for the p-value of the statistic.

nobs

A single integer for the number of observations.

npar

A single integer for the number of parameters.

weighted

A single logical for whether the data are weighted or not.

coefficients

A numeric vector of the maximum empirical likelihood estimates of the parameters.

method

A single character for the method dispatch in internal functions.

control

An object of class ControlEL constructed by el_control().

Examples

showClass("SummaryEL")

SummaryELMT class

Description

S4 class for a summary of ELMT objects.

Slots

aliased

A named logical vector showing if the original coefficients are aliased.

Examples

showClass("SummaryELMT")

SummaryELT class

Description

S4 class for a summary of ELT objects.

Slots

optim

A list of the following optimization results:

• par A numeric vector of the solution to the (constrained) optimization problem.

• lambda A numeric vector of the Lagrange multipliers of the dual problem corresponding to par.

• iterations A single integer for the number of iterations performed.

• convergence A single logical for the convergence status.

• cstr A single logical for whether constrained EL optimization is performed or not.

logl

A single numeric of the (constrained) empirical log-likelihood.

loglr

A single numeric of the (constrained) empirical log-likelihood ratio.

statistic

A single numeric of minus twice the (constrained) empirical log-likelihood ratio with an asymptotic chi-square distribution.

df

A single integer for the chi-square degrees of freedom of the statistic.

pval

A single numeric for the (calibrated) p-value of the statistic.

cv

A single numeric for the critical value.

rhs

A numeric vector for the right-hand side of the hypothesis.

lhs

A numeric matrix for the left-hand side of the hypothesis.

alpha

A single numeric for the significance level.

calibrate

A single character for the calibration method used.

control

An object of class ControlEL constructed by el_control().

Examples

showClass("SummaryELT")

SummaryGLM class

Description

S4 class for a summary of GLM objects. It inherits from SummaryLM class.

Slots

family

A family object used.

dispersion

A single numeric for the estimated dispersion parameter.

coefficients

A numeric matrix of the results of significance tests.

intercept

A single logical for whether the given model has an intercept term or not.

na.action

Information returned by model.frame on the special handling of NAs.

call

A matched call.

terms

A terms object used.

aliased

A named logical vector showing if the original coefficients are aliased.

optim

A list of the following optimization results:

• par A numeric vector of the solution to the (constrained) optimization problem.

• lambda A numeric vector of the Lagrange multipliers of the dual problem corresponding to par.

• iterations A single integer for the number of iterations performed.

• convergence A single logical for the convergence status.

• cstr A single logical for whether constrained EL optimization is performed or not.

logl

A single numeric of the empirical log-likelihood.

loglr

A single numeric of the empirical log-likelihood ratio.

statistic

A single numeric of minus twice the (constrained) empirical log-likelihood ratio for the overall test.

df

A single integer for the degrees of freedom of the statistic.

pval

A single numeric for the p-value of the statistic.

nobs

A single integer for the number of observations.

npar

A single integer for the number of parameters.

weighted

A single logical for whether the data are weighted or not.

method

A single character for the method dispatch in internal functions.

control

An object of class ControlEL constructed by el_control().

Examples

showClass("SummaryGLM")

SummaryLM class

Description

S4 class for a summary of LM objects.

Slots

coefficients

A numeric matrix of the results of significance tests.

intercept

A single logical for whether the given model has an intercept term or not.

na.action

Information returned by model.frame on the special handling of NAs.

call

A matched call.

terms

A terms object used.

aliased

A named logical vector showing if the original coefficients are aliased.

optim

A list of the following optimization results:

• par A numeric vector of the solution to the (constrained) optimization problem.

• lambda A numeric vector of the Lagrange multipliers of the dual problem corresponding to par.

• iterations A single integer for the number of iterations performed.

• convergence A single logical for the convergence status.

• cstr A single logical for whether constrained EL optimization is performed or not.

logl

A single numeric of the empirical log-likelihood.

loglr

A single numeric of the empirical log-likelihood ratio.

statistic

A single numeric of minus twice the (constrained) empirical log-likelihood ratio for the overall test.

df

A single integer for the degrees of freedom of the statistic.

pval

A single numeric for the p-value of the statistic.

nobs

A single integer for the number of observations.

npar

A single integer for the number of parameters.

weighted

A single logical for whether the data are weighted or not.

method

A single character for the method dispatch in internal functions.

control

An object of class ControlEL constructed by el_control().

Examples

showClass("SummaryLM")

SummaryQGLM class

Description

S4 class for a summary of QGLM objects. It inherits from SummaryGLM class.

Slots

family

A family object used.

dispersion

A single numeric for the estimated dispersion parameter.

coefficients

A numeric matrix of the results of significance tests.

intercept

A single logical for whether the given model has an intercept term or not.

na.action

Information returned by model.frame on the special handling of NAs.

call

A matched call.

terms

A terms object used.

aliased

A named logical vector showing if the original coefficients are aliased.

optim

A list of the following optimization results:

• par A numeric vector of the solution to the (constrained) optimization problem.

• lambda A numeric vector of the Lagrange multipliers of the dual problem corresponding to par.

• iterations A single integer for the number of iterations performed.

• convergence A single logical for the convergence status.

• cstr A single logical for whether constrained EL optimization is performed or not.

logl

A single numeric of the empirical log-likelihood.

loglr

A single numeric of the empirical log-likelihood ratio.

statistic

A single numeric of minus twice the (constrained) empirical log-likelihood ratio for the overall test.

df

A single integer for the degrees of freedom of the statistic.

pval

A single numeric for the p-value of the statistic.

nobs

A single integer for the number of observations.

npar

A single integer for the number of parameters.

weighted

A single logical for whether the data are weighted or not.

method

A single character for the method dispatch in internal functions.

control

An object of class ControlEL constructed by el_control().

Examples

showClass("SummaryQGLM")

Chi-square statistic

Description

Extracts the chi-square statistic from a model.

Usage

## S4 method for signature 'EL'
chisq(object, ...)

## S4 method for signature 'ELMT'
chisq(object, ...)

## S4 method for signature 'ELT'
chisq(object, ...)

## S4 method for signature 'SummaryEL'
chisq(object, ...)

## S4 method for signature 'SummaryELMT'
chisq(object, ...)

## S4 method for signature 'SummaryELT'
chisq(object, ...)

## S4 method for signature 'SummaryLM'
chisq(object, ...)

Arguments

Value

The form of the value returned by chisq() depends on the class of its argument.

Methods (by class)

• chisq(EL): Extracts the chi-square statistic.

• chisq(ELMT): Extracts the vector of chi-square statistics.

• chisq(ELT): Extracts the chi-square statistic.

• chisq(SummaryEL): Extracts the chi-square statistic.

• chisq(SummaryELMT): Extracts the vector of chi-square statistics.

• chisq(SummaryELT): Extracts the chi-square statistic.

• chisq(SummaryLM): Extracts the chi-square statistic for the overall test of the model.

See Also

EL, ELMT, ELT, pVal()

Examples

data("precip")
fit <- el_mean(precip, par = 40)
chisq(fit)

Clothianidin concentration in maize plants

Description

A dataset summarizing field experiments result of seed treatments on clothianidin concentration.

Usage

data("clothianidin")

Format

A data frame with 102 observations and 3 variables:

blk

New blocks constructed from original data. The format is 'days post planting_original block_year'.

trt

Seed treatment.

clo

Log transformed clothianidin concentration (µg).

Details

The original data is provided by Alford and Krupke (2017). Only some of the shoot region observations are taken from the original data and processed for illustration.

Source

Alford A, Krupke CH (2017). “Translocation of the Neonicotinoid Seed Treatment Clothianidin in Maize.” PLOS ONE, 12(3), 1–19. doi:10.1371/journal.pone.0173836.

Examples

data("clothianidin")
clothianidin

Model coefficients

Description

Extracts the maximum empirical likelihood estimates from a model.

Usage

## S4 method for signature 'EL'
coef(object, ...)

## S4 method for signature 'ELMT'
coef(object, ...)

## S4 method for signature 'SummaryEL'
coef(object, ...)

## S4 method for signature 'SummaryLM'
coef(object, ...)

Arguments

Value

The form of the value returned by coef() depends on the class of its argument.

Methods (by class)

• coef(EL): Extracts the numeric vector of the maximum empirical likelihood estimates.

• coef(ELMT): Extracts the list of numeric vectors of the maximum empirical likelihood estimates. Each element of the list corresponds to a distinct hypothesis.

• coef(SummaryEL): Extracts the numeric vector of the maximum empirical likelihood estimates.

• coef(SummaryLM): Extracts a matrix with the results of significance tests.

See Also

EL, ELMT

Examples

data("mtcars")
fit <- el_lm(mpg ~ wt, data = mtcars)
coef(fit)

Confidence interval for model parameters

Description

Computes confidence intervals for one or more parameters in a model.

Usage

## S4 method for signature 'EL'
confint(object, parm, level = 0.95, cv = NULL, control = NULL)

## S4 method for signature 'ELMT'
confint(object, cv = NULL, control = NULL)

Arguments

Value

A matrix with columns giving lower and upper confidence limits for each parameter. In contrast to other methods that rely on studentization, the lower and upper limits obtained from empirical likelihood do not correspond to the (1 - level) / 2 and 1 - (1 - level) / 2 in %, respectively.

References

Owen A (1990). “Empirical Likelihood Ratio Confidence Regions.” The Annals of Statistics, 18(1), 90–120. doi:10.1214/aos/1176347494.

See Also

EL, ELMT, confreg(), elt(), el_control()

Examples

data("mtcars")
fit <- el_lm(mpg ~ ., data = mtcars)
confint(fit, parm = c(2, 3))

Confidence region for model parameters

Description

Computes boundary points of a two-dimensional confidence region for model parameters.

Usage

## S4 method for signature 'EL'
confreg(object, parm, level = 0.95, cv = NULL, npoints = 50L, control = NULL)

Arguments

Value

An object of class ConfregEL.

References

Owen A (1990). “Empirical Likelihood Ratio Confidence Regions.” The Annals of Statistics, 18(1), 90–120. doi:10.1214/aos/1176347494.

See Also

EL, confint(), elt(), plot(), el_control()

Examples

data("mtcars")
fit <- el_lm(mpg ~ wt + qsec, data = mtcars)
cr <- confreg(fit, parm = c(2, 3), cv = qchisq(0.90, 2))
plot(cr)

Convergence check

Description

Extracts the convergence status from a model.

Usage

## S4 method for signature 'CEL'
conv(object, ...)

## S4 method for signature 'EL'
conv(object, ...)

## S4 method for signature 'ELT'
conv(object, ...)

## S4 method for signature 'SummaryEL'
conv(object, ...)

## S4 method for signature 'SummaryELT'
conv(object, ...)

## S4 method for signature 'SummaryLM'
conv(object, ...)

Arguments

Value

A single logical.

Methods (by class)

• conv(CEL): Extracts the convergence status of the model with respect to the parameter.

• conv(EL): Extracts the convergence status of the model with respect to the Lagrange multiplier.

• conv(ELT): Extracts the convergence status of the test with respect to the parameter (or the Lagrange multiplier if the argument lhs is NULL).

• conv(SummaryEL): Extracts the convergence status of the model with respect to the Lagrange multiplier.

• conv(SummaryELT): Extracts the convergence status of the test with respect to the parameter (or the Lagrange multiplier if the argument lhs is NULL).

• conv(SummaryLM): Extracts the convergence status of the model. See the documentation of EL and CEL.

See Also

CEL, EL, ELT, getOptim()

Examples

## Convergence check for the overall model test
data("mtcars")
fit <- el_lm(mpg ~ ., data = mtcars)
conv(fit)

Critical value

Description

Extracts the critical value from a model.

Usage

## S4 method for signature 'ELMT'
critVal(object, ...)

## S4 method for signature 'ELT'
critVal(object, ...)

## S4 method for signature 'SummaryELMT'
critVal(object, ...)

## S4 method for signature 'SummaryELT'
critVal(object, ...)

Arguments

Value

A single numeric.

See Also

ELMT, ELT

Examples

## F-calibrated critical value
data("precip")
fit <- el_mean(precip, 30)
elt <- elt(fit, rhs = 34, calibrate = "f")
critVal(elt)

Control parameters for computation

Description

Specifies computational details of (constrained) empirical likelihood.

Usage

el_control(
  maxit = 200L,
  maxit_l = 25L,
  tol = 1e-06,
  tol_l = 1e-06,
  step = NULL,
  th = NULL,
  verbose = FALSE,
  keep_data = TRUE,
  nthreads,
  seed = NULL,
  an = NULL,
  b = 10000L,
  m = 1000000L
)

Arguments

Value

An object of class of ControlEL.

See Also

el_eval(), elt()

Examples

optcfg <- el_control(maxit = 300, step = 0.01, th = 200, nthreads = 1)

Empirical likelihood for general estimating functions

Description

Computes empirical likelihood with general estimating functions.

Usage

el_eval(g, weights = NULL, control = el_control())

Arguments

Details

Let X_i be independent and identically distributed p-dimensional random variable from an unknown distribution P for i = 1, \dots, n. We assume that P has a positive definite covariance matrix. For a parameter of interest \theta(F) \in {\rm{I\!R}}^p, consider a p-dimensional smooth estimating function g(X_i, \theta) with a moment condition

\textrm{E}[g(X_i, \theta)] = 0.

We assume that there exists an unique \theta_0 that solves the above equation. Given a value of \theta, the (profile) empirical likelihood ratio is defined by

R(\theta) = \max_{p_i}\left\{\prod_{i = 1}^n np_i : \sum_{i = 1}^n p_i g(X_i, \theta) = 0, p_i \geq 0, \sum_{i = 1}^n p_i = 1 \right\}.

el_mean() computes the empirical log-likelihood ratio statistic -2\log R(\theta) with the n by p numeric matrix g, whose ith row is g(X_i, \theta). Since the estimating function can be arbitrary, el_eval() does not return an object of class EL, and the associated generics and methods are not applicable.

Value

A list of the following optimization results:

• optim A list with the following optimization results:

  • lambda A numeric vector of the Lagrange multipliers of the dual problem.

  • iterations A single integer for the number of iterations performed.

  • convergence A single logical for the convergence status.

• logp A numeric vector of the log probabilities of the empirical likelihood.

• logl A single numeric of the empirical log-likelihood.

• loglr A single numeric of the empirical log-likelihood ratio.

• statistic A single numeric of minus twice the empirical log-likelihood ratio with an asymptotic chi-square distribution.

• df A single integer for the degrees of freedom of the statistic.

• pval A single numeric for the p-value of the statistic.

• nobs A single integer for the number of observations.

• npar A single integer for the number of parameters.

• weights A numeric vector of the re-scaled weights used for the model fitting.

References

Qin J, Lawless J (1994). “Empirical Likelihood and General Estimating Equations.” The Annals of Statistics, 22(1), 300–325. doi:10.1214/aos/1176325370.

See Also

EL, el_control()

Examples

set.seed(123526)
mu <- 0
sigma <- 1
x <- rnorm(100)
g <- matrix(c(x - mu, (x - mu)^2 - sigma^2), ncol = 2)
el_eval(g, weights = rep(c(1, 2), each = 50))

Empirical likelihood for generalized linear models

Description

Fits a generalized linear model with empirical likelihood.

Usage

el_glm(
  formula,
  family = gaussian,
  data,
  weights = NULL,
  na.action,
  start = NULL,
  etastart = NULL,
  mustart = NULL,
  offset,
  control = el_control(),
  ...
)

Arguments

Details

Suppose that we observe n independent random variables {Z_i} \equiv {(X_i, Y_i)} from a common distribution, where X_i is the p-dimensional covariate (including the intercept if any) and Y_i is the response. A generalized linear model specifies that {\textrm{E}(Y_i | X_i)} = {\mu_i}, {G(\mu_i)} = {X_i^\top \theta}, and {\textrm{Var}(Y_i | X_i)} = {\phi V(\mu_i)}, where \theta = (\theta_0, \dots, \theta_{p-1}) is an unknown p-dimensional parameter, \phi is an optional dispersion parameter, G is a known smooth link function, and V is a known variance function.

With H denoting the inverse link function, define the quasi-score

{g_1(Z_i, \theta)} = \left\{ H^\prime(X_i^\top \theta) \left(Y_i - H(X_i^\top \theta)\right) / \left(\phi V\left(H(X_i^\top \theta)\right)\right) \right\} X_i.

Then we have the estimating equations \sum_{i = 1}^n g_1(Z_i, \theta) = 0. When \phi is known, the (profile) empirical likelihood ratio for a given \theta is defined by

R_1(\theta) = \max_{p_i}\left\{\prod_{i = 1}^n np_i : \sum_{i = 1}^n p_i g_1(Z_i, \theta) = 0,\ p_i \geq 0,\ \sum_{i = 1}^n p_i = 1 \right\}.

With unknown \phi, we introduce another estimating function based on the squared residuals. Let {\eta} = {(\theta, \phi)} and

{g_2(Z_i, \eta)} = \left(Y_i - H(X_i^\top \theta)\right)^2 / \left(\phi^2 V\left(H(X_i^\top \theta)\right)\right) - 1 / \phi.

Now the empirical likelihood ratio is defined by

R_2(\eta) = \max_{p_i}\left\{\prod_{i = 1}^n np_i : \sum_{i = 1}^n p_i g_1(Z_i, \eta) = 0,\ \sum_{i = 1}^n p_i g_2(Z_i, \eta) = 0,\ p_i \geq 0,\ \sum_{i = 1}^n p_i = 1 \right\}.

el_glm() first computes the parameter estimates by calling glm.fit() (with ... if any) with the model.frame and model.matrix obtained from the formula. Note that the maximum empirical likelihood estimator is the same as the the quasi-maximum likelihood estimator in our model. Next, it tests hypotheses based on asymptotic chi-square distributions of the empirical likelihood ratio statistics. Included in the tests are overall test with

H_0: \theta_1 = \theta_2 = \cdots = \theta_{p-1} = 0,

and significance tests for each parameter with

H_{0j}: \theta_j = 0,\ j = 0, \dots, p-1.

The available families and link functions are as follows:

• gaussian: "identity", "log", and "inverse".

• binomial: "logit", "probit", and "log".

• poisson: "log", "identity", and "sqrt".

• quasipoisson: "log", "identity", and "sqrt".

Value

An object of class of GLM.

References

Chen SX, Cui H (2003). “An Extended Empirical Likelihood for Generalized Linear Models.” Statistica Sinica, 13(1), 69–81.

Kolaczyk ED (1994). “Empirical Likelihood for Generalized Linear Models.” Statistica Sinica, 4(1), 199–218.

See Also

EL, GLM, el_lm(), elt(), el_control()

Examples

data("warpbreaks")
fit <- el_glm(wool ~ .,
  family = binomial, data = warpbreaks, weights = NULL, na.action = na.omit,
  start = NULL, etastart = NULL, mustart = NULL, offset = NULL
)
summary(fit)

Empirical likelihood for linear models

Description

Fits a linear model with empirical likelihood.

Usage

el_lm(
  formula,
  data,
  weights = NULL,
  na.action,
  offset,
  control = el_control(),
  ...
)

Arguments

Details

Suppose that we observe n independent random variables {Z_i} \equiv {(X_i, Y_i)} from a common distribution, where X_i is the p-dimensional covariate (including the intercept if any) and Y_i is the response. We consider the following linear model:

Y_i = X_i^\top \theta + \epsilon_i,

where \theta = (\theta_0, \dots, \theta_{p-1}) is an unknown p-dimensional parameter and the errors \epsilon_i are independent random variables that satisfy \textrm{E}(\epsilon_i | X_i) = 0. We assume that the errors have finite conditional variances. Then the least square estimator of \theta solves the following estimating equations:

\sum_{i = 1}^n(Y_i - X_i^\top \theta)X_i = 0.

Given a value of \theta, let {g(Z_i, \theta)} = {(Y_i - X_i^\top \theta)X_i} and the (profile) empirical likelihood ratio is defined by

R(\theta) = \max_{p_i}\left\{\prod_{i = 1}^n np_i : \sum_{i = 1}^n p_i g(Z_i, \theta) = \theta,\ p_i \geq 0,\ \sum_{i = 1}^n p_i = 1 \right\}.

el_lm() first computes the parameter estimates by calling lm.fit() (with ... if any) with the model.frame and model.matrix obtained from the formula. Note that the maximum empirical likelihood estimator is the same as the the quasi-maximum likelihood estimator in our model. Next, it tests hypotheses based on asymptotic chi-square distributions of the empirical likelihood ratio statistics. Included in the tests are overall test with

H_0: \theta_1 = \theta_2 = \cdots = \theta_{p-1} = 0,

and significance tests for each parameter with

H_{0j}: \theta_j = 0,\ j = 0, \dots, p-1.

Value

An object of class of LM.

References

Owen A (1991). “Empirical Likelihood for Linear Models.” The Annals of Statistics, 19(4), 1725–1747. doi:10.1214/aos/1176348368.

See Also

EL, LM, el_glm(), elt(), el_control()

Examples

## Linear model
data("thiamethoxam")
fit <- el_lm(fruit ~ trt, data = thiamethoxam)
summary(fit)

## Weighted data
wfit <- el_lm(fruit ~ trt, data = thiamethoxam, weights = visit)
summary(wfit)

## Missing data
fit2 <- el_lm(fruit ~ trt + scb, data = thiamethoxam,
  na.action = na.omit, offset = NULL
)
summary(fit2)

Empirical likelihood for the mean

Description

Computes empirical likelihood for the mean.

Usage

el_mean(x, par, weights = NULL, control = el_control())

Arguments

Details

Let X_i be independent and identically distributed p-dimensional random variable from an unknown distribution P for i = 1, \dots, n. We assume that {\textrm{E}[X_i]} = {\theta_0} \in {\rm{I\!R}}^p and that P has a positive definite covariance matrix. Given a value of \theta, the (profile) empirical likelihood ratio is defined by

R(\theta) = \max_{p_i}\left\{\prod_{i = 1}^n np_i : \sum_{i = 1}^n p_i X_i = \theta,\ p_i \geq 0,\ \sum_{i = 1}^n p_i = 1 \right\}.

el_mean() computes the empirical log-likelihood ratio statistic -2\log R(\theta), along with other values in EL.

Value

An object of class EL.

References

Owen A (1990). “Empirical Likelihood Ratio Confidence Regions.” The Annals of Statistics, 18(1), 90–120. doi:10.1214/aos/1176347494.

See Also

EL, elt(), el_eval(), el_control()

Examples

## Scalar mean
data("precip")
fit <- el_mean(precip, 30)
fit
summary(fit)

## Vector mean
data("faithful")
fit2 <- el_mean(faithful, par = c(3.5, 70))
summary(fit2)

## Weighted data
w <- rep(c(1, 2), each = nrow(faithful) / 2)
fit3 <- el_mean(faithful, par = c(3.5, 70), weights = w)
summary(fit3)

Empirical likelihood for the standard deviation

Description

Computes empirical likelihood for the standard deviation.

Usage

el_sd(x, mean, sd, weights = NULL, control = el_control())

Arguments

Details

Let X_i be independent and identically random variable from an unknown distribution P for i = 1, \dots, n. We assume that {\textrm{E}[X_i]} = {\mu_0} is known and that P has a variance \sigma_0^2. Given a value of \sigma, the (profile) empirical likelihood ratio is defined by

R(\sigma) = \max_{p_i}\left\{\prod_{i = 1}^n np_i : \sum_{i = 1}^n p_i (X_i - \mu_0)^2 = \sigma^2,\ p_i \geq 0,\ \sum_{i = 1}^n p_i = 1 \right\}.

el_sd() computes the empirical log-likelihood ratio statistic -2\log R(\sigma), along with other values in SD.

Value

An object of class SD.

See Also

EL, SD, el_mean(), elt(), el_control()

Examples

data("women")
x <- women$height
w <- women$weight
fit <- el_sd(x, mean = 65, sd = 5, weights = w)
fit
summary(fit)

Empirical likelihood displacement

Description

Computes empirical likelihood displacement for model diagnostics and outlier detection.

Usage

## S4 method for signature 'EL'
eld(object, control = NULL)

## S4 method for signature 'GLM'
eld(object, control = NULL)

Arguments

Details

Let L(\theta) be the empirical log-likelihood function based on the full sample with n observations. The maximum empirical likelihood estimate is denoted by \hat{\theta}. Consider a reduced sample with the ith observation deleted and the corresponding estimate \hat{\theta}_{(i)}. The empirical likelihood displacement is defined by

\textrm{ELD}_i = 2\{L(\hat{\theta}) - L(\hat{\theta}_{(i)})\}.

If \textrm{ELD}_i is large, then the ith observation is an influential point and can be inspected as a possible outlier. eld computes \textrm{ELD}_i for i = 1, \dots, n .

Value

An object of class ELD.

References

Lazar NA (2005). “Assessing the Effect of Individual Data Points on Inference From Empirical Likelihood.” Journal of Computational and Graphical Statistics, 14(3), 626–642. doi:10.1198/106186005X59568.

Zhu H, Ibrahim JG, Tang N, Zhang H (2008). “Diagnostic Measures for Empirical Likelihood of General Estimating Equations.” Biometrika, 95(2), 489–507. doi:10.1093/biomet/asm094.

See Also

EL, ELD, el_control(), plot()

Examples

data("precip")
fit <- el_mean(precip, par = 30)
eld <- eld(fit)
plot(eld)

Empirical likelihood multiple tests

Description

Tests multiple linear hypotheses simultaneously.

Usage

## S4 method for signature 'EL'
elmt(object, rhs = NULL, lhs = NULL, alpha = 0.05, control = NULL)

Arguments

Details

elmt() tests multiple hypotheses simultaneously. Each hypothesis corresponds to the constrained empirical likelihood ratio described in CEL. rhs and lhs cannot be both NULL. The right-hand side and left-hand side of each hypothesis must be specified as described in elt().

For specifying linear contrasts more conveniently, rhs and lhs also take a numeric vector and a numeric matrix, respectively. Each element of rhs and each row of lhs correspond to a contrast (hypothesis).

The vector of empirical likelihood ratio statistics asymptotically follows a multivariate chi-square distribution under the complete null hypothesis. The multiple testing procedure asymptotically controls the family-wise error rate at the level alpha. Based on the distribution of the maximum of the test statistics, the adjusted p-values are estimated by Monte Carlo simulation.

Value

An object of class of ELMT.

References

Kim E, MacEachern SN, Peruggia M (2023). “Empirical likelihood for the analysis of experimental designs.” Journal of Nonparametric Statistics, 35(4), 709–732. doi:10.1080/10485252.2023.2206919.

Kim E, MacEachern SN, Peruggia M (2024). “melt: Multiple Empirical Likelihood Tests in R.” Journal of Statistical Software, 108(5), 1–33. doi:10.18637/jss.v108.i05.

See Also

EL, ELMT, elt(), el_control()

Examples

## Bivariate mean (list `rhs` & no `lhs`)
set.seed(143)
data("women")
fit <- el_mean(women, par = c(65, 135))
rhs <- list(c(64, 133), c(66, 140))
elmt(fit, rhs = rhs)

## Pairwise comparison (no `rhs` & list `lhs`)
data("clothianidin")
fit2 <- el_lm(clo ~ -1 + trt, clothianidin)
lhs2 <- list(
  "trtNaked - trtFungicide",
  "trtFungicide - trtLow",
  "trtLow - trtHigh"
)
elmt(fit2, lhs = lhs2)

## Arbitrary hypotheses (list `rhs` & list `lhs`)
data("mtcars")
fit3 <- el_lm(mpg ~ wt + qsec, data = mtcars)
lhs3 <- list(c(1, 4, 0), rbind(c(0, 1, 0), c(0, 0, 1)))
rhs3 <- list(0, c(-6, 1))
elmt(fit3, rhs = rhs3, lhs = lhs3)

Empirical likelihood test

Description

Tests a linear hypothesis with various calibration options.

Usage

## S4 method for signature 'EL'
elt(
  object,
  rhs = NULL,
  lhs = NULL,
  alpha = 0.05,
  calibrate = "chisq",
  control = NULL
)

Arguments

Details

elt() performs the constrained minimization of l(\theta) described in CEL. rhs and lhs cannot be both NULL. For non-NULL lhs, it is required that lhs have full row rank q \leq p and p be equal to the number of parameters in the object.

Depending on the specification of rhs and lhs, we have the following three cases:

1. If both rhs and lhs are non-NULL, the constrained minimization is performed with the right-hand side r and the left-hand side L as

   \inf_{\theta: L\theta = r} l(\theta).

2. If rhs is NULL, r is set to the zero vector as \inf_{\theta: L\theta = 0} l(\theta).

3. If lhs is NULL, L is set to the identity matrix and the problem reduces to evaluating at r as l(r).

calibrate specifies the calibration method used. Four methods are available: "ael" (adjusted empirical likelihood calibration), "boot" (bootstrap calibration), "chisq" (chi-square calibration), and "f" (F calibration). When lhs is not NULL, only "chisq" is available. "f" is applicable only to the mean parameter. The an slot in control applies specifically to "ael", while the nthreads, seed, and B slots apply to "boot".

Value

An object of class of ELT. If lhs is non-NULL, the optim slot corresponds to that of CEL. Otherwise, it corresponds to that of EL.

References

Adimari G, Guolo A (2010). “A Note on the Asymptotic Behaviour of Empirical Likelihood Statistics.” Statistical Methods & Applications, 19(4), 463–476. doi:10.1007/s10260-010-0137-9.

Chen J, Variyath AM, Abraham B (2008). “Adjusted Empirical Likelihood and Its Properties.” Journal of Computational and Graphical Statistics, 17(2), 426–443.

Kim E, MacEachern SN, Peruggia M (2024). “melt: Multiple Empirical Likelihood Tests in R.” Journal of Statistical Software, 108(5), 1–33. doi:10.18637/jss.v108.i05.

Qin J, Lawless J (1995). “Estimating Equations, Empirical Likelihood and Constraints on Parameters.” Canadian Journal of Statistics, 23(2), 145–159. doi:10.2307/3315441.

See Also

EL, ELT, elmt(), el_control()

Examples

## Adjusted empirical likelihood calibration
data("precip")
fit <- el_mean(precip, 32)
elt(fit, rhs = 100, calibrate = "ael")

## Bootstrap calibration
elt(fit, rhs = 32, calibrate = "boot")

## F calibration
elt(fit, rhs = 32, calibrate = "f")

## Test of no treatment effect
data("clothianidin")
contrast <- matrix(c(
  1, -1, 0, 0,
  0, 1, -1, 0,
  0, 0, 1, -1
), byrow = TRUE, nrow = 3)
fit2 <- el_lm(clo ~ -1 + trt, clothianidin)
elt(fit2, lhs = contrast)

## A symbolic description of the same hypothesis
elt(fit2, lhs = c(
  "trtNaked - trtFungicide",
  "trtFungicide - trtLow",
  "trtLow - trtHigh"
))

Degrees of freedom

Description

Extracts the degrees of freedom from a model.

Usage

## S4 method for signature 'EL'
getDF(object)

## S4 method for signature 'ELMT'
getDF(object)

## S4 method for signature 'ELT'
getDF(object)

## S4 method for signature 'SummaryEL'
getDF(object)

## S4 method for signature 'SummaryELMT'
getDF(object)

## S4 method for signature 'SummaryLM'
getDF(object)

Arguments

Value

An integer vector.

Methods (by class)

• getDF(EL): Extracts the degrees of freedom.

• getDF(ELMT): Extracts the vector of marginal degrees of freedom.

• getDF(ELT): Extracts the (chi-square) degrees of freedom.

• getDF(SummaryEL): Extracts the degrees of freedom.

• getDF(SummaryELMT): Extracts the vector of marginal degrees of freedom.

• getDF(SummaryLM): Extracts the degrees of freedom.

See Also

EL, ELMT, ELT

Examples

data("faithful")
fit <- el_mean(faithful, par = c(3.5, 70))
getDF(fit)

Optimization results

Description

Extracts the optimization results from a model.

Usage

## S4 method for signature 'EL'
getOptim(object, ...)

## S4 method for signature 'ELT'
getOptim(object, ...)

## S4 method for signature 'SummaryEL'
getOptim(object, ...)

## S4 method for signature 'SummaryELT'
getOptim(object, ...)

## S4 method for signature 'SummaryLM'
getOptim(object, ...)

Arguments

Value

A list with the following optimization results:

• par A numeric vector of the parameter value. See the documentation of EL and CEL.

• lambda A numeric vector of the Lagrange multipliers.

• iterations A single integer for the number of iterations performed.

• convergence A single logical for the convergence status.

See Also

EL, ELT, sigTests()

Examples

data("precip")
fit <- el_mean(precip, par = 40)
getOptim(fit)

Empirical log-likelihood

Description

Extracts the empirical log-likelihood from a model.

Usage

## S4 method for signature 'EL'
logL(object, ...)

## S4 method for signature 'ELT'
logL(object, ...)

## S4 method for signature 'SummaryEL'
logL(object, ...)

## S4 method for signature 'SummaryELT'
logL(object, ...)

## S4 method for signature 'SummaryLM'
logL(object, ...)

Arguments

Value

A single numeric.

References

Baggerly KA (1998). “Empirical Likelihood as a Goodness-of-Fit Measure.” Biometrika, 85(3), 535–547. doi:10.1093/biomet/85.3.535.

See Also

EL, ELT

Examples

data("precip")
fit <- el_mean(precip, par = 40)
logL(fit)

Empirical log-likelihood ratio

Description

Extracts the empirical log-likelihood ratio from a model.

Usage

## S4 method for signature 'EL'
logLR(object, ...)

## S4 method for signature 'ELT'
logLR(object, ...)

## S4 method for signature 'SummaryEL'
logLR(object, ...)

## S4 method for signature 'SummaryELT'
logLR(object, ...)

## S4 method for signature 'SummaryLM'
logLR(object, ...)

Arguments

Value

A single numeric.

References

Baggerly KA (1998). “Empirical Likelihood as a Goodness-of-Fit Measure.” Biometrika, 85(3), 535–547. doi:10.1093/biomet/85.3.535.

See Also

EL, ELT

Examples

data("precip")
fit <- el_mean(precip, par = 40)
logLR(fit)

Log probabilities

Description

Extracts log probabilities of empirical likelihood from a model.

Usage

## S4 method for signature 'EL'
logProb(object, ...)

## S4 method for signature 'ELT'
logProb(object, ...)

Arguments

Value

A numeric vector.

See Also

EL, ELT

Examples

data("precip")
fit <- el_mean(precip, par = 40)
logProb(fit)

Number of observations in a model

Description

Extracts the number of observations from a model.

Usage

## S4 method for signature 'EL'
nobs(object, ...)

## S4 method for signature 'SummaryEL'
nobs(object, ...)

## S4 method for signature 'SummaryLM'
nobs(object, ...)

Arguments

Value

A single integer.

See Also

EL

Examples

data("precip")
fit <- el_mean(precip, par = 40)
nobs(fit)

p-value

Description

Extracts the p-value from a model.

Usage

## S4 method for signature 'EL'
pVal(object, ...)

## S4 method for signature 'ELMT'
pVal(object, ...)

## S4 method for signature 'ELT'
pVal(object, ...)

## S4 method for signature 'SummaryEL'
pVal(object, ...)

## S4 method for signature 'SummaryELT'
pVal(object, ...)

## S4 method for signature 'SummaryELMT'
pVal(object, ...)

## S4 method for signature 'SummaryLM'
pVal(object, ...)

Arguments

Value

The form of the value returned by pVal() depends on the class of its argument.

Methods (by class)

• pVal(EL): Extracts the p-value.

• pVal(ELMT): Extracts the multiplicity adjusted p-values.

• pVal(ELT): Extracts the p-value.

• pVal(SummaryEL): Extracts the p-value.

• pVal(SummaryELT): Extracts the p-value.

• pVal(SummaryELMT): Extracts the multiplicity adjusted p-values.

• pVal(SummaryLM): Extracts the p-value.

See Also

EL, ELMT, ELT, chisq()

Examples

data("precip")
fit <- el_mean(precip, par = 40)
pVal(fit)

Plot methods

Description

Provides plot methods for objects.

Usage

## S4 method for signature 'ConfregEL'
plot(x, y, ...)

## S4 method for signature 'EL'
plot(x, y, ...)

## S4 method for signature 'ELD'
plot(x, y, ...)

Arguments

Value

No return value, called for side effects.

Methods (by class)

• plot(ConfregEL): Plots a two-dimensional confidence region for model parameters.

• plot(EL): Plots empirical likelihood displacement values versus observation index. eld() is called implicitly.

• plot(ELD): Plots empirical likelihood displacement values versus observation index.

See Also

ConfregEL, EL, ELD, confreg(), eld()

Examples

## Model
data("mtcars")
fit <- el_lm(hp ~ wt, data = mtcars)

## Confidence region
out1 <- confreg(fit, npoints = 500)
plot(out1)

## Empirical likelihood displacement
out2 <- eld(fit)
plot(out2)

## A shortcut to `ELD`
plot(fit)

Print methods

Description

Provides print methods for objects.

Usage

## S4 method for signature 'EL'
print(x, digits = max(3L, getOption("digits") - 3L), ...)

## S4 method for signature 'ELMT'
print(x, digits = max(3L, getOption("digits") - 3L), ...)

## S4 method for signature 'ELT'
print(x, digits = max(3L, getOption("digits") - 3L), ...)

## S4 method for signature 'LM'
print(x, digits = max(3L, getOption("digits") - 3L), ...)

## S4 method for signature 'SummaryEL'
print(x, digits = max(3L, getOption("digits") - 3L), ...)

## S4 method for signature 'SummaryELMT'
print(
  x,
  digits = max(3L, getOption("digits") - 3L),
  signif.stars = getOption("show.signif.stars"),
  ...
)

## S4 method for signature 'SummaryELT'
print(x, digits = max(3L, getOption("digits") - 3L), ...)

## S4 method for signature 'SummaryGLM'
print(
  x,
  digits = max(3L, getOption("digits") - 3L),
  signif.stars = getOption("show.signif.stars"),
  ...
)

## S4 method for signature 'SummaryLM'
print(
  x,
  digits = max(3L, getOption("digits") - 3L),
  signif.stars = getOption("show.signif.stars"),
  ...
)

Arguments

Value

The argument x (invisibly).

See Also

EL, ELMT, ELT, LM

Examples

data("precip")
fit <- el_mean(precip, par = 40)
print(fit)

Significance tests

Description

Extracts the results of significance tests from a model.

Usage

## S4 method for signature 'LM'
sigTests(object, ...)

Arguments

Value

The form of the value returned by sigTests() depends on the class of its argument.

Methods (by class)

• sigTests(LM): Extracts a list with the optimization results of significance tests.

See Also

LM, getOptim()

Examples

data("mtcars")
fit <- el_lm(mpg ~ ., data = mtcars)
sigTests(fit)

Summary methods

Description

Provides summary methods for objects.

Usage

## S4 method for signature 'EL'
summary(object, ...)

## S4 method for signature 'ELMT'
summary(object, ...)

## S4 method for signature 'ELT'
summary(object, ...)

## S4 method for signature 'GLM'
summary(object, ...)

## S4 method for signature 'LM'
summary(object, ...)

## S4 method for signature 'QGLM'
summary(object, ...)

Arguments

Value

The form of the value returned by summary() depends on the class of its argument.

Methods (by class)

• summary(EL): Summarizes the test results of the specified parameters.

• summary(ELMT): Summarizes the multiple testing results.

• summary(ELT): Summarizes the hypothesis test results.

• summary(GLM): Summarizes the results of the overall model test and the significance tests for coefficients. The dispersion parameter is extracted for display.

• summary(LM): Summarizes the results of the overall model test and the significance tests for coefficients.

• summary(QGLM): Summarizes the results of the overall model test and the significance tests for coefficients. The estimated dispersion parameter is extracted for display.

See Also

EL, ELMT, ELT, GLM, LM, QGLM,

Examples

data("faithful")
fit <- el_mean(faithful, par = c(3.5, 70))
summary(fit)

data("mtcars")
fit2 <- el_lm(mpg ~ wt, data = mtcars)
summary(fit2)

Thiamethoxam applications in squash crops

Description

A dataset on the effect of the thiamethoxam application method and plant variety on bees.

Usage

data("thiamethoxam")

Format

A data frame with 165 observations and 11 variables:

trt

Treatment.

var

Variety.

rep

Replicate.

fruit

Average fruit number per plant.

avg_mass

Individual Fruit mass average (g).

mass

Fruit mass per plant (g).

yield

Yield (4 plants).

visit

Bee visits per plot.

foliage

Proportion of foliage consumed by striped cucumber beetle.

scb

Striped cucumber beetle per plant.

defoliation

Defoliation percentage.

Source

Obregon D, Pederson G, Taylor A, Poveda K (2022). “The Pest Control and Pollinator Protection Dilemma: The Case of Thiamethoxam Prophylactic Applications in Squash Crops.” PLOS ONE, 17(5), 1–18. doi:10.1371/journal.pone.0267984.

Examples

data("thiamethoxam")
thiamethoxam

Model weights

Description

Extracts weights from model objects. The weights are re-scaled to up to the total number of observations in the fitting procedure.

Usage

## S4 method for signature 'EL'
weights(object, ...)

Arguments

Value

A numeric vector of the re-scaled weights.

References

Glenn N, Zhao Y (2007). “Weighted Empirical Likelihood Estimates and Their Robustness Properties.” Computational Statistics & Data Analysis, 51(10), 5130–5141. doi:10.1016/j.csda.2006.07.032.

See Also

EL

Examples

data("airquality")
x <- airquality$Wind
w <- airquality$Day
fit <- el_mean(x, par = 10, weights = w)
weights(fit)