An embedded and charset-unspecified text was scrubbed... Name: not available Url: https://stat.ethz.ch/pipermail/r-help/attachments/20070511/5402cdb5/attachment.pl
PRESS criterion in leaps
2 messages · Andrew Smith, Thomas Lumley
2 days later
On Fri, 11 May 2007, Andrew Smith wrote:
I thought it would be simplest to build on already existing functions like regsubsets in package leaps. It's easy enough to calculate the PRESS criterion for a fitted lm object, but I'm having trouble deciphering the structure of the regsubsets objects that leaps works with. Is there any way to calculate press from a regsubsets? Or, to put it another way, can I get the residual vector and the diagonal entries of the hat matrix from a regsubsets object? In fact, if the hat matrix is never calculated explicitly, the columns of Q from the QR factorization would suffice.
Not only is the hat matrix never calculated explicitly, the Q matrix isn't calculated either. The code forms R and Q^TY directly (the same code is used in the biglm package to provide bounded-space linear regression). -thomas