Re: [R] singular information matrix in lrm.fit

2008-10-12 Thread Gad Abraham
Prof Brian Ripley wrote: I believe lrm has a criterion appropriate to single-precision calculations (as S-PLUS used to use). Try reducing 'tol' from its default of 1e-7. But your design matrix *is* near singular kappa(cbind(1,x)) [1] 557390.5 so try centring/scaling your variables. Than

Re: [R] singular information matrix in lrm.fit

2008-10-11 Thread Prof Brian Ripley
I believe lrm has a criterion appropriate to single-precision calculations (as S-PLUS used to use). Try reducing 'tol' from its default of 1e-7. But your design matrix *is* near singular kappa(cbind(1,x)) [1] 557390.5 so try centring/scaling your variables. On Sun, 12 Oct 2008, Gad Abraham

[R] singular information matrix in lrm.fit

2008-10-11 Thread Gad Abraham
Hi, I'm trying to do binary logistic regression on 10 covariables, comparing glm to lrm from Harrell's Design package. They don't seem to agree on whether the data is collinear: > library(Design) > load(url("http://www.csse.unimelb.edu.au/~gabraham/data.Rdata";)) > lrm(y ~ X1 + X2 + X3 + X4 +

Re: [R] Singular information matrix in lrm.fit

2008-10-09 Thread Charles C. Berry
On Thu, 9 Oct 2008, useR wrote: Hi R helpers, I'm fitting large number of single factor logistic regression models as a way to immediatly discard factor which are insignificant. Everything works fine expect that for some factors I get error message "Singular information matrix in lrm.fit" whi

[R] Singular information matrix in lrm.fit

2008-10-09 Thread useR
Hi R helpers, I'm fitting large number of single factor logistic regression models as a way to immediatly discard factor which are insignificant. Everything works fine expect that for some factors I get error message "Singular information matrix in lrm.fit" which breaks whole execution loop... h