rkpk: Numerical Engine for ssanova and gssanova

rkpkR Documentation

Numerical Engine for ssanova and gssanova

Description

Perform numerical calculations for the ssanova and gssanova suites.

Usage

sspreg1(s, r, q, y, wt, method, alpha, varht, random)
mspreg1(s, r, id.basis, y, wt, method, alpha, varht, random, skip.iter)
ngreg1(family, s, r, id.basis, y, wt, offset, method, varht, alpha, nu, random, skip.iter)

sspreg91(s, r, q, y, cov, method, alpha, varht)
mspreg91(s, r, id.basis, y, cov, method, alpha, varht, skip.iter)

sspngreg(family, s, r, q, y, wt, offset, alpha, nu, random)
mspngreg(family, s, r, id.basis, y, wt, offset, alpha, nu, random, skip.iter)
ngreg(dc, family, sr, q, y, wt, offset, nu, alpha)

regaux(s, r, q, nlambda, fit)

ngreg.proj(dc, family, sr, q, y0, wt, offset, nu)

Arguments

family

Description of the error distribution. Supported are exponential families "binomial", "poisson", "Gamma", and "nbinomial". Also supported are accelerated life model families "weibull", "lognorm", and "loglogis".

s

Unpenalized terms evaluated at data points.

r

Basis of penalized terms evaluated at data points.

q

Penalty matrix.

id.basis

Index of observations to be used as "knots."

y

Response vector.

wt

Model weights.

cov

Input for covariance function for correlated data.

offset

Model offset.

method

"v" for GCV, "m" for GML, or "u" for Mallows' CL.

alpha

Parameter modifying GCV or Mallows' CL scores for smoothing parameter selection.

nu

Optional argument for future support of nbinomial, weibull, lognorm, and loglogis families.

varht

External variance estimate needed for method="u".

random

Input for parametric random effects in nonparametric mixed-effect models.

skip.iter

Flag indicating whether to use initial values of theta and skip theta iteration.

nlambda

Smoothing parameter in effect.

fit

Fitted model.

dc

Coefficients of fits.

sr

cbind(s,r).

y0

Components of the fit to be projected.

Details

sspreg1 is used by ssanova to compute regression estimates with a single smoothing parameter. mspreg1 is used by ssanova to compute regression estimates with multiple smoothing parameters.

ssngpreg is used by gssanova to compute non-Gaussian regression estimates with a single smoothing parameter. mspngreg is used by gssanova to compute non-Gaussian regression estimates with multiple smoothing parameters. ngreg is used by ssngpreg and mspngreg to perform Newton iteration with fixed smoothing parameters and to calculate cross-validation scores on return.

regaux is used by sspreg1, mspreg1, ssngpreg, and mspngreg to obtain auxiliary information needed by predict.ssanova for standard error calculation.

ngreg.proj is used by project.gssanova to calculate the Kullback-Leibler projection for non-Gaussian regression.

References

Gu, C. (1992), Cross validating non Gaussian data. Journal of Computational and Graphical Statistics, 1, 169–179.

Kim, Y.-J. and Gu, C. (2004), Smoothing spline Gaussian regression: more scalable computation via efficient approximation. Journal of the Royal Statistical Society, Ser. B, 66, 337–356.


gss documentation built on Aug. 16, 2023, 9:07 a.m.

Related to rkpk in gss...