MortSmooth_BWB: Calculates the inner product of a matrix and a sparse weight...
In MortalitySmooth: Smoothing and Forecasting Poisson Counts with P-Splines

Description Usage Arguments Details Value Author(s) References See Also

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 within the function Mort2Dsmooth. It construct the LHD of the IWLS algorithm following the idea proposed in the Generalized Linear Array Models (see references).

1	MortSmooth_BWB(RTBx, RTBy, nbx, nby, W)

`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.
`W`	matrix of weights.

The function employs an arithmetic of arrays which allows to define the LHS of the Poisson system of equations as a sequence of nested matrix operations. Such way of operating with arrays leads to low storage and high speed computation when data are structure as array.

A matrix of inner product of a matrix and a sparse weight matrix.

Carlo G Camarda

Eilers, P. H. C., I. D. Currie, and M. Durban (2006). Fast and compact smoothing on large multidimensional grids. Computational Statistics & Data Analysis. 50, 61-76.

Currie, I. D., M. Durban, and P. H. C. Eilers (2006). Generalized linear array models with applications to multidimentional smoothing. Journal of the Royal Statistical Society. Series B. 68, 259-280.

Mort2Dsmooth, Mort2Dsmooth_estimate, Mort2Dsmooth_update.

MortalitySmooth documentation built on May 2, 2019, 6:07 a.m.