Mort2Dsmooth.se {MortalitySmooth}R Documentation

Compute a 2D standard errors

Description

This is an internal function of package MortalitySmooth which calculates the inner product of a matrix (from a Kronecker product) and a sparse weight matrix in order to obtain standard errors. It uses the same idea employed in MortSmooth.BWB and the elements after the IWLS converged, including the penalty term.

Usage

Mort2Dsmooth.se(RTBx, RTBy, nbx, nby, BWB.P1)

Arguments

RTBx Tensors product of B-splines basis for the x-axis.
RTBy Tensors product of 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.
BWB.P1 Inverse of the LHS of the Poisson system of equations.

Details

This function is only used within predict.Mort2Dsmooth when standard errors are required. The arguments BWB.P1 is the LHS after convergence is reached and smoothing parameter selected. The standard errors as given in the function are computed for the linear predictor term and simple computation is needed to obtain standard errors for the Poisson counts. Anyway predict.Mort2Dsmooth takes care of such differences.

The Generalized Linear Array Models setting is explained in the reference in MortSmooth.BWB and Mort2Dsmooth.

Value

A matrix of standard errors for the linear predictor term.

Author(s)

Carlo G Camarda

See Also

Mort2Dsmooth, MortSmooth.BWB, predict.Mort2Dsmooth.


[Package MortalitySmooth version 1.0 Index]