[R-br] Atualização da função cld()
Walmes Zeviani
walmeszeviani em gmail.com
Quinta Dezembro 29 14:11:39 BRST 2016
Marcos,
Eu não sei te dizer o que mudou na mulcomp::cld() mas na época dessa
matéria ela devia funcionar. Para situações semelhantes a essa do
fatorial, eu acabei fazendo uma cópia da cld() e modifiquei algumas coisas
para que pudesse usá-la. Ela está disponível no fonte de um pacote onde
mantenho uma miscelânea de funções: github.com/walmes/wzRfun. O código
abaixo usa a cld2() que é a cópia que mencionei.
# source("
https://raw.githubusercontent.com/walmes/wzRfun/master/R/pairwise.R")
# cld2(c0)
# Função disponível no pacote wzRfun: github.com/walmes/wzRfun.
cld2 <- function(object, level = 0.05) {
lvl_order <- unique(unlist(
strsplit(rownames(object$linfct), "-")))
signif <- (object$test$pvalues < level)
ret <- list()
ret$signif <- signif
ret$comps <- do.call(rbind,
strsplit(rownames(object$linfct), "-"))
# Modificação feita aqui.
ret$mcletters <-
multcomp:::insert_absorb(x = signif,
decreasing = FALSE,
comps = ret$comps,
lvl_order = lvl_order)
class(ret) <- "cld"
return(ret)
}
cld2(c0)
result$cld <- cld2(c0)$mcletters$Letters
result[with(result, order(Estimate)), ]
ordered_cld <- function(x) {
s <- strsplit(x, "")
ul <- unique(unlist(s))
UL <- toupper(sort(ul, decreasing = TRUE))
l <- sapply(s, FUN = function(i) {
tolower(paste(sort(UL[match(i, table = ul)]), collapse = ""))
})
return(l)
}
# Ordena o vetor de letras pelas médias dos tratamentos.
v <- result[with(result, order(Estimate)), ]$cld
result$cld <- ordered_cld(v)
À disposição.
Walmes.
-------------- Próxima Parte ----------
Um anexo em HTML foi limpo...
URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20161229/a7fff636/attachment.html>
Mais detalhes sobre a lista de discussão R-br