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?
Thank you,
Annie
Please read earlier posts. Avoid it by not having an overdescribed
problem which will result in overfitting.
Frank
On Fri, Sep 4, 2009 at 2:44 PM, Frank E Harrell Jr
<f.harr...@vanderbilt.edu <mailto:f.harr...@vanderbilt.edu>> wrote:
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 CovA CovB CovC CovD CovE
2 1 0 1 1 0 -0.005575280 2
2 1 0 1 0 1 -0.001959580 2
3 0 0 0 1 0 -0.004725880 2
0 0 0 1 0 0 -0.005504850 2
2 1 1 0 0 0 -0.003880170 1
2 1 0 0 1 0 -0.006074230 2
2 1 0 0 1 1 -0.003963920 2
2 1 0 0 1 0 0.000658230 2
2 1 0 1 0 1 -0.002033610 2
0 0 0 0 0 1 -0.007761770 2
0 0 0 0 1 1 -0.000951784 2
Note: In the complete data, there are NAs. But the error
messages showed in both the complete and above sets.
attach(mydata)
lrm(RC~RS+Sex+CovA+CovB+CovC+CovD+CovE,data=mydata,na.action=na.pass)
singular information matrix in lrm.fit (rank= 8 ). Offending
variable(s):
CovE Error in j:(j + params[i] - 1) : NA/NaN argument
detach(mydata)
lrm(mydata$RC~ mydata$RS+ mydata$Sex+ mydata$CovA+
mydata$CovB+ mydata$CovC+ mydata$CovD+
mydata$CovE,data=mydata,na.action=na.pass)
Please read the documentation. Use data=. Don't use $ in variables
in the model. I don't know what na.pass is. lrm does not allow NAs
to be used in the attempted model fit.
Frank
Error in if (!length(fname) || !any(fname == zname)) { :
missing value where TRUE/FALSE needed
I have debugged for a while, but no success. Does anyone know
what happens?
Thanks,
Mike
______________________________________________
R-help@r-project.org <mailto: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
<http://www.r-project.org/posting-guide.html>
and provide commented, minimal, self-contained, reproducible code.
______________________________________________
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.
