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
|
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