[R] lrm in Design package--missing value where TRUE/FALSE needed

Hi, A error message arose while I was trying to fit a ordinal model with lrm() I am using R 2.8 with Design package. Here is a small set of mydata: RC RS Sex CovACovBCovCCovDCovE 2 1 0 1 1 0 -0.0055752802 2 1 0

Re: [R] lrm in Design package--missing value where TRUE/FALSE needed

The idea of estimating 7 parameters from 11 cases using mostly binary variables just seems begging for a singular matrix. On Sep 4, 2009, at 3:59 PM, Ping-Hsun Hsieh wrote: Hi, A error message arose while I was trying to fit a ordinal model with lrm() I am using R 2.8 with Design package.

Re: [R] lrm in Design package--missing value where TRUE/FALSE needed

Ping-Hsun Hsieh wrote: Hi, A error message arose while I was trying to fit a ordinal model with lrm() I am using R 2.8 with Design package. Here is a small set of mydata: RC RS Sex CovACovBCovCCovDCovE 2 1 0 1 1 0 -0.005575280

Re: [R] lrm in Design package--missing value where TRUE/FALSE needed

annie Zhang wrote: Hi, Frank, I met the same problem. My data does not have NA, when I run fit - lrm(Y_t~.,data=X) The error message is: singular information matrix in lrm.fit (rank= 35 ). Offending variable(s): X35 Error in j:(j + params[i] - 1) : NA/NaN argument How can I avoid this?

Re: [R] lrm in Design package--missing value where TRUE/FALSE needed

Hi, Frank, I met the same problem. My data does not have NA, when I run fit - lrm(Y_t~.,data=X) The error message is: singular information matrix in lrm.fit (rank= 35 ). Offending variable(s): X35 Error in j:(j + params[i] - 1) : NA/NaN argument How can I avoid this? Thank you, Annie On Fri,

Re: [R] lrm in Design

Thank you for your reply, it helped. That thread is about glm, I take it that lrm behaves in the same way?!? See: https://stat.ethz.ch/pipermail/r-help/2009-July/204192.html and also other posts in that thread. -- View this message in context:

Re: [R] lrm in Design

Presumably the same idea would apply to lrm if appropriately adapted. On Sun, Aug 30, 2009 at 5:19 AM, loch1l...@gmx.li wrote: Thank you for your reply, it helped. That thread is about glm, I take it that lrm behaves in the same way?!? See:

[R] lrm in Design

Hello everybody, I am trying to do a logistic regression model with lrm() from the design package. I am comparing to groups with different medical outcome which can either be good or bad. In the help file it says that lrm codes al responses to 0,1,2,3, etc. internally and does so in alphabetical

Re: [R] lrm in Design

See: https://stat.ethz.ch/pipermail/r-help/2009-July/204192.html and also other posts in that thread. On Sat, Aug 29, 2009 at 6:02 PM, loch1l...@gmx.li wrote: Hello everybody, I am trying to do a logistic regression model with lrm() from the design package. I am comparing to groups with

[R] lrm in Design package, missing value problem

__ R-help@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.

[R] lrm in Design package, missing value problem

Hi, This concerns lrm in the Design package, which I need for the clustered errors option with panel data. I have created a variable (using sign()) that has values -1, 0, and 1 (and some NAs). call this variable x3. lrm does not like t = lrm(y ~ x1 + x2 + I(x3==-1) + I(x3==1),

Re: [R] lrm in Design package, missing value problem

James Fearon wrote: Hi, This concerns lrm in the Design package, which I need for the clustered errors option with panel data. I have created a variable (using sign()) that has values -1, 0, and 1 (and some NAs). call this variable x3. lrm does not like t = lrm(y ~ x1 + x2 + I(x3==-1) +