diferença minina significativa
Pessoal mais uma vez procurando ajuda de vc's, tenho um DIC e estou precissando estimar DMS no R, alguem tem algum material que possa me ajudar. Obrigado -- Eng. Agr.Gilson Sánchez Chia, MSc. Laboratório de Fisiologia Vegetal Embrapa Amazônia Ocidental Fone: (92) 3303-7841 --- Antes de imprimir, pense na sua responsabilidade com o Meio Ambiente. Please consider the environment before printing this email. Renew, Reduce and Recycle. The Planet Knows.
Se for balanceado, use a função agricolae::HSD.test() que solta essa estimativa na saída da função, veja
x <- gl(5,5) y <- rnorm(x) agricolae::HSD.test(y, x, DFerror=20, MSError=1)
Study: HSD Test for y Mean Square Error: 1 x, means y std.err replication 1 -0.3230810 0.1327155 5 2 0.3595503 0.5492514 5 3 -0.1703578 0.3377758 5 4 -0.1483504 0.5129000 5 5 -0.1101464 0.2727455 5 alpha: 0.05 ; Df Error: 20 Critical Value of Studentized Range: 4.231857 Honestly Significant Difference: 1.892544 <------------------------------------------ DMS Means with the same letter are not significantly different. Groups, Treatments and means a 2 0.3595503 a 5 -0.1101464 a 4 -0.1483504 a 3 -0.1703578 a 1 -0.323081
========================================================================== Walmes Marques Zeviani LEG (Laboratório de Estatística e Geoinformação, 25.450418 S, 49.231759 W) Departamento de Estatística - Universidade Federal do Paraná fone: (+55) 41 3361 3573 VoIP: (3361 3600) 1053 1173 e-mail: walmes@ufpr.br twitter: @walmeszeviani homepage: http://www.leg.ufpr.br/~walmes linux user number: 531218 ==========================================================================
Estava olhando a discussão, e não consegui rodar esta função agricolae! De que pacote ela é? Se for balanceado, use a função agricolae::HSD.test() que solta essa estimativa na saída da função, veja
x <- gl(5,5) y <- rnorm(x) agricolae::HSD.test(y, x, DFerror=20, MSError=1)
Study: HSD Test for y Mean Square Error: 1 x, means y std.err replication 1 -0.3230810 0.1327155 5 2 0.3595503 0.5492514 5 3 -0.1703578 0.3377758 5 4 -0.1483504 0.5129000 5 5 -0.1101464 0.2727455 5 alpha: 0.05 ; Df Error: 20 Critical Value of Studentized Range: 4.231857 Honestly Significant Difference: 1.892544 <------------------------------------------ DMS Means with the same letter are not significantly different. Groups, Treatments and means a 2 0.3595503 a 5 -0.1101464 a 4 -0.1483504 a 3 -0.1703578 a 1 -0.323081
========================================================================== Walmes Marques Zeviani LEG (Laboratório de Estatística e Geoinformação, 25.450418 S, 49.231759 W) Departamento de Estatística - Universidade Federal do Paraná fone: (+55) 41 3361 3573 VoIP: (3361 3600) 1053 1173 e-mail: walmes@ufpr.br twitter: @walmeszeviani homepage: http://www.leg.ufpr.br/~walmes linux user number: 531218
Mauro, agricolae é o pacote. A função HSD.test() poderá ser encontrada nela. att Peternelli 2011/9/20 Mauro Sznelwar <sznelwar@uol.com.br>
** Estava olhando a discussão, e não consegui rodar esta função agricolae! De que pacote ela é?
Se for balanceado, use a função agricolae::HSD.test() que solta essa estimativa na saída da função, veja
x <- gl(5,5) y <- rnorm(x) agricolae::HSD.test(y, x, DFerror=20, MSError=1)
Study:
HSD Test for y
Mean Square Error: 1
x, means
y std.err replication 1 -0.3230810 0.1327155 5 2 0.3595503 0.5492514 5 3 -0.1703578 0.3377758 5 4 -0.1483504 0.5129000 5 5 -0.1101464 0.2727455 5
alpha: 0.05 ; Df Error: 20 Critical Value of Studentized Range: 4.231857
Honestly Significant Difference: 1.892544 <------------------------------------------ DMS
Means with the same letter are not significantly different.
Groups, Treatments and means a 2 0.3595503 a 5 -0.1101464 a 4 -0.1483504 a 3 -0.1703578 a 1 -0.323081
========================================================================== Walmes Marques Zeviani LEG (Laboratório de Estatística e Geoinformação, 25.450418 S, 49.231759 W) Departamento de Estatística - Universidade Federal do Paraná fone: (+55) 41 3361 3573 VoIP: (3361 3600) 1053 1173 e-mail: walmes@ufpr.br twitter: @walmeszeviani homepage: http://www.leg.ufpr.br/~walmes linux user number: 531218
_______________________________________________ R-br mailing list R-br@listas.c3sl.ufpr.br 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.
-- ***************************************************************** Prof. Luiz A. Peternelli, Professor Associado II Departamento de Estatística Universidade Federal de Vicosa (UFV) Vicosa, MG, Brasil, 36570-000 55 (31) 3899-1787 - http://www.det.ufv.br/~peternelli *****************************************************************
participantes (4)
-
Gilson Sanchez -
Luiz Alexandre Peternelli -
Mauro Sznelwar -
Walmes Zeviani