Obtaining (orthogonal) basis for null space and range of the penalty matrix
INTERNAL function to obtain (orthogonal) basis for the null space and range space of the penalty, and obtain actual null space dimension components are roughly rescaled to avoid any dominating.
totalPenaltySpace(S, H, off, p)
S |
a list of penalty matrices, in packed form. |
H |
the coefficient matrix of an user supplied fixed quadratic penalty on the parameters of the GAM. |
off |
a vector where the i-th element is the offset for the i-th matrix. |
p |
total number of parameters. |
A list of matrix square roots such that S[[i]]=B[[i]]%*%t(B[[i]]).
Simon N. Wood <simon.wood@r-project.org>.
Please choose more modern alternatives, such as Google Chrome or Mozilla Firefox.