Alexandre,
No seu CMR tem-se um fator apenas. Você não precisa de "a" e "a.int", eles
são iguais. Será que não é a versão da apc que você tem? Ela tem sido
aprimorada, está no meu pacote wzRfun, disponível do
www.github.com/walmes/wzRfun. A mais recente tá assim:
## wzRfun::apc
apc <- function (lfm, lev = NULL)
{
nlev <- nrow(lfm)
rn <- rownames(lfm)
a <- attr(lfm, "grid")
if (is.null(lev)) {
if (!is.null(a)) {
lev <- apply(a, 1, paste, collapse = ":")
}
else if (!is.null(rn)) {
lev <- rn
}
else {
lev <- as.character(1:nlev)
}
}
cbn <- combn(seq_along(lev), 2)
M <- lfm[cbn[1, ], ] - lfm[cbn[2, ], ]
if (is.vector(M))
dim(M) <- c(1, length(M))
rownames(M) <- paste(lev[cbn[1, ]], lev[cbn[2, ]], sep = "-")
return(M)
}
À disposição.
Walmes.
_______________________________________________
R-br mailing list
[email protected]
https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça código
mínimo reproduzível.
