Displaying 1 result from an estimated 1 matches for "mlmobject".
Did you mean:
lmobject
2011 Jan 19
1
buglet in weighted.residuals(mlmObject)
When weighted.residuals() is given a fitted model
object with several responses (class mlm) and some
zero weights it returns a vector instead of a matrix.
It looks like it is doing
resids[ weights != 0 ]
instead of
resids[ weights != 0, , drop=FALSE]
in the multi-response case.
E.g.,
> d4 <- data.frame(y1=1:4, y2=2^(0:3), wt=log(1:4),
fac=LETTERS[c(1,1,2,2)])
> fit <-