Hi all,
I am looking to fit a logistic regression using the lrm function from the Design library. I am interested in this function because I would like to obtain "pseudo-R2" values (see http://tolstoy.newcastle.edu.au/R/help/02b/1011.html).
Can anyone help me with the syntax?
If I fit the model using the stats library, the code looks like this:
model <- glm(x$trait ~ x$PC1 + I((x$PC1)^2) + I((x$PC1)^3), family = binomial)
What would be the equivalent syntax for the lrm function?
Thanks very much in advance,
-----------------------------------
Josh Banta, Ph.D
Center for Genomics and Systems Biology
New York University
100 Washington Square East
New York, NY 10003
Tel: (212) 998-8465
http://plantevolutionaryecology.org
[[alternative HTML version deleted]]