
Olá prezados, boa tarde! Estou com uma dificuldade para efetuar os contrastes de tratamentos de um ensaio em DBC 3*3 + 2. Minha situação assemelha-se a do código: https://ridiculas.wordpress.com/tag/fatorial/ Entretanto quando chego na parte de apresentar os contrastes utilizando da função cld() aparece o seguinte erro: Trecho do código:
#----------------------------------------------------------------------------- # fazendo as comparações com a glht, método para correção de p-valor # fdr, o single-step demora muito c0 <- summary(glht(m6, linfct=Xc), test=adjusted(type="fdr")) c0
Simultaneous Tests for General Linear Hypotheses Fit: lm(formula = y ~ -1 + mm, data = fu) Linear Hypotheses: Estimate Std. Error t value Pr(>|t|) benlate:puro-benlate:antes == 0 0.008761 0.070134 0.125 0.901153 benlate:puro-benlate:misturado == 0 -0.100584 0.070134 -1.434 0.249216 benlate:puro-captam:puro == 0 -0.364723 0.070134 -5.200 1.71e-05 *** benlate:puro-captam:antes == 0 -0.055481 0.070134 -0.791 0.554020 benlate:puro-captam:misturado == 0 -0.302970 0.070134 -4.320 0.000201 *** benlate:puro-derosal:puro == 0 0.027171 0.070134 0.387 0.755244 benlate:puro-derosal:antes == 0 0.040517 0.070134 0.578 0.649000 benlate:puro-derosal:misturado == 0 0.072273 0.070134 1.031 0.436194 benlate:puro-polimero:polimero == 0 -0.375107 0.070134 -5.348 1.19e-05 *** benlate:puro-testemunha:testemunha == 0 -0.417647 0.070134 -5.955 2.98e-06 *** benlate:antes-benlate:misturado == 0 -0.109345 0.070134 -1.559 0.204047 benlate:antes-captam:puro == 0 -0.373485 0.070134 -5.325 1.20e-05 *** benlate:antes-captam:antes == 0 -0.064242 0.070134 -0.916 0.496458 benlate:antes-captam:misturado == 0 -0.311731 0.070134 -4.445 0.000147 *** benlate:antes-derosal:puro == 0 0.018409 0.070134 0.262 0.839986 benlate:antes-derosal:antes == 0 0.031756 0.070134 0.453 0.718215 benlate:antes-derosal:misturado == 0 0.063512 0.070134 0.906 0.496458 benlate:antes-polimero:polimero == 0 -0.383869 0.070134 -5.473 9.13e-06 *** benlate:antes-testemunha:testemunha == 0 -0.426409 0.070134 -6.080 2.36e-06 *** benlate:misturado-captam:puro == 0 -0.264140 0.070134 -3.766 0.001033 ** benlate:misturado-captam:antes == 0 0.045103 0.070134 0.643 0.625923 benlate:misturado-captam:misturado == 0 -0.202386 0.070134 -2.886 0.011843 * benlate:misturado-derosal:puro == 0 0.127755 0.070134 1.822 0.129457 benlate:misturado-derosal:antes == 0 0.141101 0.070134 2.012 0.092362 . benlate:misturado-derosal:misturado == 0 0.172857 0.070134 2.465 0.033547 * benlate:misturado-polimero:polimero == 0 -0.274524 0.070134 -3.914 0.000685 *** benlate:misturado-testemunha:testemunha == 0 -0.317064 0.070134 -4.521 0.000121 *** captam:puro-captam:antes == 0 0.309242 0.070134 4.409 0.000158 *** captam:puro-captam:misturado == 0 0.061754 0.070134 0.881 0.502026 captam:puro-derosal:puro == 0 0.391894 0.070134 5.588 6.78e-06 *** captam:puro-derosal:antes == 0 0.405241 0.070134 5.778 4.37e-06 *** captam:puro-derosal:misturado == 0 0.436997 0.070134 6.231 1.70e-06 *** captam:puro-polimero:polimero == 0 -0.010384 0.070134 -0.148 0.899323 captam:puro-testemunha:testemunha == 0 -0.052924 0.070134 -0.755 0.568124 captam:antes-captam:misturado == 0 -0.247489 0.070134 -3.529 0.002019 ** captam:antes-derosal:puro == 0 0.082652 0.070134 1.178 0.364093 captam:antes-derosal:antes == 0 0.095998 0.070134 1.369 0.271962 captam:antes-derosal:misturado == 0 0.127755 0.070134 1.822 0.129457 captam:antes-polimero:polimero == 0 -0.319626 0.070134 -4.557 0.000113 *** captam:antes-testemunha:testemunha == 0 -0.362166 0.070134 -5.164 1.81e-05 *** captam:misturado-derosal:puro == 0 0.330141 0.070134 4.707 7.29e-05 *** captam:misturado-derosal:antes == 0 0.343487 0.070134 4.898 4.13e-05 *** captam:misturado-derosal:misturado == 0 0.375243 0.070134 5.350 1.19e-05 *** captam:misturado-polimero:polimero == 0 -0.072138 0.070134 -1.029 0.436194 captam:misturado-testemunha:testemunha == 0 -0.114678 0.070134 -1.635 0.181928 derosal:puro-derosal:antes == 0 0.013347 0.070134 0.190 0.882020 derosal:puro-derosal:misturado == 0 0.045103 0.070134 0.643 0.625923 derosal:puro-polimero:polimero == 0 -0.402278 0.070134 -5.736 4.53e-06 *** derosal:puro-testemunha:testemunha == 0 -0.444818 0.070134 -6.342 1.45e-06 *** derosal:antes-derosal:misturado == 0 0.031756 0.070134 0.453 0.718215 derosal:antes-polimero:polimero == 0 -0.415625 0.070134 -5.926 2.98e-06 *** derosal:antes-testemunha:testemunha == 0 -0.458165 0.070134 -6.533 1.45e-06 *** derosal:misturado-polimero:polimero == 0 -0.447381 0.070134 -6.379 1.45e-06 *** derosal:misturado-testemunha:testemunha == 0 -0.489921 0.070134 -6.986 6.60e-07 *** polimero:polimero-testemunha:testemunha == 0 -0.042540 0.070134 -0.607 0.640414 --- Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1 (Adjusted p values reported -- fdr method)
c0$focus <- "comparacoes" # para poder usar a cld() cld(c0) Error in contrMat(table(x), type = "Tukey") : less than two groups fdr(c0) Error: could not find function "fdr"
Achei que poderia ser algo relacionado com atualização da função cld(), porém, buscando em outros CMR não encontrei uma resposta pertinente. Na tentantiva de uma solução, conto com a ajuda de vcs. Att., Marcos Paulo fatorial | Ridículas<https://ridiculas.wordpress.com/tag/fatorial/> ridiculas.wordpress.com Posts sobre fatorial escritos por Walmes Zeviani and fernandohtoledo Técnico Agrícola em Zootecnia Bacharel em Agronomia Mestrando em Produção Vegetal - Universidade Federal de Goiás Estagiário de Pós-Graduação da Embrapa Arroz e Feijão/CNPAF-Santo Antônio de Goiás-GO Contato: (62) 85075783 http://lattes.cnpq.br/4322347592884852

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.
participantes (2)
-
marcos paulo
-
Walmes Zeviani