Displaying 2 results from an estimated 2 matches for "la_chol2inv".
2008 Apr 03
1
Lapack error in Design:::ols
...action 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(fit$qr$qr) : 'size' cannot exceed nrow(x) = 20
> traceback()
6: .Call("La_chol2inv", x, size, PACKAGE = "base")
5: chol2inv(fit$qr$qr)
4: ols(y ~ x1 * x2, data = x)
3: eval.with.vis(expr, envir, enclos)
2: eval.with.vis(ei, envir)
1: source("t.R", echo = TRUE)
Any ideas?
Thanks,
Gad
> sessionInfo()
R version 2.6.2 (2008-02-08)
x86_64-pc-linux-gnu...
2012 Dec 14
1
Found explanation for R-2.15.2 slowdown in one case; caution for any users of La_chol
2 days ago, I posted my long message about the observed slowdown in a
package between R-2.15.0 and R-2.15.2.
Uwe Ligges urged me to make a self-contained R example. That was the
encouragement I needed. I tracked the problem down to a failing use of
a LAPACK routine.
R's LAPACK C interface changed one variable in one function. But it
turned out to be an important change. In case others have