Dear All,
I'm using S-PLUS 6.0 on WindowsXp.
I met a problem on lm.fit.chol.
Since this function has "singular.ok" argument,
I've tried;
> data <- data.frame(x=1:10, y=11:20, z=101:110)
> mat <- model. matrix(z ~ x+ y, data=data)
> lm.fit.chol(y=data$z, x=mat, singular.ok=T)
Problem in dimnames(.A1) <- .A0: Cannot have dimnames for nonarray
Use traceback() to see the call stack
So, I edit some fraction of lm.fit.chol, but I have not
received correct result yet.
Does anybody have an experience to modify this function?
Thank you,
Mika
Mika Nakazono
Mathematical Systems Inc. Japan
|