Displaying 1 result from an estimated 1 matches for "uasgn".
Did you mean:
asgn
2006 Mar 13
1
anova.mlm (single-model case) does not handle factors? (PR#8679)
...e problem lies in the computation of the 'ss' terms (line 237 in the
file mlm.R in the source code). Changing this (and the following line) by
something similar to what is done in summary.manova seems to resolve the
problem:
comp <- as.matrix(fit$effects)[seq(along=asgn), ,drop=FALSE]
uasgn <- unique(asgn)
nterms <- length(uasgn)
ss <- list(nterms)
df <- numeric(nterms)
for(i in seq(nterms)) {
ai <- (asgn == uasgn[i])
ss[[i]] <- crossprod(comp[ai, ,drop=FALSE])
df[i] <- sum(ai)
}