PSQR: A Stable and Efficient Penalized Spline Algorithm

Simple item page
dc.contributor.authorCalderon, Christopher P.
dc.contributor.authorMartinez, Josue G.
dc.contributor.authorCarroll, Raymond J.
dc.contributor.authorSorensen, Danny C.
dc.date.accessioned2018-06-19T17:45:07Z
dc.date.available2018-06-19T17:45:07Z
dc.date.issued2009-05
dc.date.noteMay 2009
dc.description.abstractWe introduce an algorithm for reliably computing quantities associated with several types of semiparametric mixed models in situations where the condition number on the random effects matrix is large. The algorithm is numerically stable and efficient. It was designed to process penalized spline (P-spline) models without making unnecessary numerical approximations. The algorithm, PSQR (P-splines via QR), is formulated in terms of QR decompositions. PSQR can treat both exactly rank deficient and ill-conditioned matrices. The latter situation often arises in large scale mixed models and/or when a P-spline is estimated using a basis with poor numerical properties, e.g. a truncated power function (TPF) basis. We provide concrete examples where unnecessary numerical approximations introduce both subtle and dramatic errors that would likely go undetected, thus demonstrating the importance of using this reliable numerical algorithm. Simulation results studying a univariate function and a longitudinal data set are used to demonstrate the algorithm. Extensions and the utility of the method in more general semiparametric regression applications are briefly discussed. MATLAB scripts demonstrating implementation are provided in the Supplemental Materials.
dc.format.extent24 pp
dc.identifier.citationCalderon, Christopher P., Martinez, Josue G., Carroll, Raymond J., et al.. "PSQR: A Stable and Efficient Penalized Spline Algorithm." (2009) <a href="https://hdl.handle.net/1911/102120">https://hdl.handle.net/1911/102120</a>.
dc.identifier.digitalTR09-15
dc.identifier.urihttps://hdl.handle.net/1911/102120
dc.language.isoeng
dc.titlePSQR: A Stable and Efficient Penalized Spline Algorithm
dc.typeTechnical report
dc.type.dcmiText
Files
Original bundle
Now showing 1 - 1 of 1
Thumbnail Image
Name:
TR09-15.pdf
Size:
505.42 KB
Format:
Adobe Portable Document Format
Download
Collections
CAAM Technical Reports