MortalitySmooth (version 2.1)

Mort2Dsmooth.estimate: Estimate 2D P-splines for two given lambdas

Description

This is an internal function of package MortalitySmooth which estimates coefficients and computes diagnostics for two-dimensional penalized B-splines for two given smoothing parameters within the function Mort2Dsmooth.

Usage

Mort2Dsmooth.estimate(x, y, Z, offset, psi2, wei,
                      Bx, By, nbx, nby, RTBx, RTBy,
                      lambdas, Px, Py, a.init, 
                      MON, TOL1, MAX.IT)

Arguments

x

vector for the abscissa of data.

y

vector for the ordinate of data.

Z

matrix of counts response.

offset

matrix with an a priori known component (optional).

psi2

overdispersion parameter.

wei

an optional matrix of weights to be used in the fitting process.

Bx

B-splines basis for the x-axis.

By

B-splines basis for the y-axis.

nbx

number of B-splines for the x-axis.

nby

number of B-splines for the y-axis.

RTBx

tensors product of B-splines basis for the x-axis.

RTBy

tensors product of B-splines basis for the y-axis.

lambdas

vector with the two smoothing parameters.

Px

penalty factor for the x-axis.

Py

penalty factor for the y-axis.

a.init

matrix with the initial coefficients.

MON

logical switch indicating if monitoring is required.

TOL1

the tolerance level in the IWLS algorithm.

MAX.IT

the maximum number of iterations.

Value

• A list with components:
• a
  fitted coefficients (in a matrix).
• h
  diagonal of the hat-matrix.
• df
  effective dimension of used degree of freedom.
• aic
  Akaike's Information Criterion.
• bic
  Bayesian Information Criterion.
• dev
  Poisson deviance.
• tol
  tolerance level.
• BWB
  inner product of basis and weights.
• P
  penalty matrix.

Details

Internal function used in Mort2Dsmooth for estimating coefficients and computing diagnostics.

See Also

Mort2Dsmooth.update, Mort2Dsmooth.