[R] using 'lrm' for logistic regression

From: Utkarsh Singhal <utkarshs_at_ambaresearch.com>
Date: Mon, 03 Mar 2008 18:53:18 +0530


Hi R,  

I am getting this error while trying to use 'lrm' function with nine independent variables:  

> res =

lrm(y1994~WC08301+WC08376+WC08316+WC08311+WC01001+WC08221+WC08106+WC0810 1+WC08231,data=y)  

singular information matrix in lrm.fit (rank= 8 ). Offending variable(s):

WC08101 WC08221

Error in j:(j + params[i] - 1) : NA/NaN argument  

Now, if I take choose only four independent variables then there is no error.  

> res = lrm(y1994~ WC08221+WC08106+WC08101+WC08231,data=y)
   

But strangely, if I use 'glm', with the family as binomial(logit) and with the same dataset, it is working perfectly fine.  

> res =

glm(y1994~WC08301+WC08376+WC08316+WC08311+WC01001+WC08221+WC08106+WC0810 1+WC08231,data=y,family=binomial(logit))  

Any ideas..?  

Regards

Utkarsh  

This e-mail may contain confidential and/or privileged i...{{dropped:13}}



R-help_at_r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code. Received on Mon 03 Mar 2008 - 13:31:35 GMT

Archive maintained by Robert King, hosted by the discipline of statistics at the University of Newcastle, Australia.
Archive generated by hypermail 2.2.0, at Mon 03 Mar 2008 - 14:30:18 GMT.

Mailing list information is available at https://stat.ethz.ch/mailman/listinfo/r-help. Please read the posting guide before posting to the list.

list of date sections of archive