Re: [R] Lapack error in Design:::ols

2008-04-05 Thread Gad Abraham
Prof Brian Ripley wrote: > This is not an LAPACK error, but an R function chol2inv complaining that > it is being misused. > > The problem is that you have a vastly overparametrized model (using > recover()) > > Browse[1]> dim(X) > [1] 20 143 > > and ols() does not detect that, whereas lm() w

Re: [R] Lapack error in Design:::ols

2008-04-03 Thread Prof Brian Ripley
This is not an LAPACK error, but an R function chol2inv complaining that it is being misused. The problem is that you have a vastly overparametrized model (using recover()) Browse[1]> dim(X) [1] 20 143 and ols() does not detect that, whereas lm() would. I'd suggest ols() is to blame here.

[R] Lapack error in Design:::ols

2008-04-03 Thread Gad Abraham
Hi, I'm trying to use Frank Harrell's Design:::ols function to do regression of y (numeric) on the interaction of two factors (x1 and x2), but Lapack throws an error: > library(Design) ... > load(url("http://www.csse.unimelb.edu.au/~gabraham/x";)) > ols(y ~ x1 * x2, data=x) Error in chol2inv