emplik: Self-concordant empirical likelihood for a vector mean

In mev: Modelling of Extreme Values

Self-concordant empirical likelihood for a vector mean

Description

Self-concordant empirical likelihood for a vector mean

Usage

emplik(
  dat,
  mu = rep(0, ncol(dat)),
  lam = rep(0, ncol(dat)),
  eps = 1/nrow(dat),
  M = 1e+30,
  thresh = 1e-30,
  itermax = 100
)

Arguments

dat	n by d matrix of d-variate observations
mu	d vector of hypothesized mean of dat
lam	starting values for Lagrange multiplier vector, default to zero vector
eps	lower cutoff for -\log, with default 1/nrow(dat)
M	upper cutoff for -\log.
thresh	convergence threshold for log likelihood (default of 1e-30 is aggressive)
itermax	upper bound on number of Newton steps.

Value

a list with components

• logelr log empirical likelihood ratio.

• lam Lagrange multiplier (vector of length d).

• wts n vector of observation weights (probabilities).

• conv boolean indicating convergence.

• niter number of iteration until convergence.

• ndec Newton decrement.

• gradnorm norm of gradient of log empirical likelihood.

Author(s)

Art Owen, C++ port by Leo Belzile

References

Owen, A.B. (2013). Self-concordance for empirical likelihood, Canadian Journal of Statistics, 41(3), 387–397.

mev documentation built on Sept. 11, 2024, 8:14 p.m.