Learn R Programming

WALS (version 0.2.5)

computeGammaUnSVD: Internal function: Computes unrestricted one-step ML estimator for transformed regressors in walsNB

Description

Computes one-step ML estimator for the unrestricted model in walsNB (coefs of transformed regressors \(\bar{Z}\)) by using SVD on entire transformed design matrix \(\bar{Z}\). The matrix \(\bar{Z}\) should have full column rank.

Usage

computeGammaUnSVD(
  U,
  V,
  singularVals,
  ellStart,
  gStart,
  epsilonStart,
  qStart,
  y0Start,
  tStart,
  psiStart
)

Arguments

U

Left singular vectors of \(\bar{Z}\) or \(\bar{Z}_{1}\) from svd.

V

Right singular vectors of \(\bar{Z}\) or \(\bar{Z}_{1}\) from svd.

singularVals

Singular values of \(\bar{Z}\) or \(\bar{Z}_{1}\) from svd.

ellStart

Vector \(\bar{\ell}\) see details.

gStart

Derivative of dispersion parameter \(\rho\) of NB2 with respect to \(\alpha = \log(\rho)\) evaluated at starting values of one-step ML. gStart is a scalar. See section "ML estimation" of huynhwalsnb;textualWALS.

epsilonStart

Scalar \(\bar{\epsilon}\), see section "One-step ML estimator" of huynhwalsnb;textualWALS for definition.

qStart

Vector \(\bar{q}\), see section "One-step ML estimator" of huynhwalsnb;textualWALS for definition.

y0Start

Vector \(\bar{y}_0\), see section "One-step ML estimator" of huynhwalsnb;textualWALS for definition.

tStart

Scalar \(\bar{t}\), see section "One-step ML estimator" of huynhwalsnb;textualWALS for definition.

psiStart

Diagonal matrix \(\bar{\Psi}\), see section "One-step ML estimator" of huynhwalsnb;textualWALS for definition.

Details

See section "Simplification for computing \(\tilde{\gamma}_{u}\)" in the appendix of huynhwals;textualWALS for details of the implementation and for the definitions of argument ellStart.

All parameters that contain "start" feature the starting values for the one-step ML estimation of submodels. See section "One-step ML estimator" of huynhwalsnb;textualWALS for details.

Uses svdLSplus under-the-hood.

References