[R] Manly's randomization analysis of multiple regression

Frank E Harrell Jr fharrell at virginia.edu
Thu May 15 12:38:10 CEST 2003

On Wed, 14 May 2003 22:52:33 -0400
Nurnberg-LaZerte <mail at fwr.on.ca> wrote:

> My wife has been using a diagnostic from Manley (1991; "Randomization and MonteCarlo Methods in Biology") that compares a normal multiple regression's performance with that using random predicted variables. 
> Is there something like this already available in R? 
> If not, the "boot" package looks like a good place to start looking for methods, no?
> Thanks in advance 
> Bruce L.

This is related to that: Look at the validate function in the Design package (http://hesweb1.med.virginia.edu/biostat/s/Design.html) which has a "randomize" option to get the apparent and validated performance of models fitted by randomly permuting the vector of response variable values.  This is not a high-precision way to estimate bias in your original R^2 etc. but rather is a teaching tool.  For efficient estimation of future model performance use the validation bootstrap option (the default for validate( )).
Frank E Harrell Jr              Prof. of Biostatistics & Statistics
Div. of Biostatistics & Epidem. Dept. of Health Evaluation Sciences
U. Virginia School of Medicine  http://hesweb1.med.virginia.edu/biostat

More information about the R-help mailing list