<font face="trebuchet ms,sans-serif">Se for balanceado, use a função agricolae::HSD.test() que solta essa estimativa na saída da função, veja<br><br><span style="font-family: courier new,monospace;">> x <- gl(5,5)</span><br style="font-family: courier new,monospace;">
<span style="font-family: courier new,monospace;">> y <- rnorm(x)</span><br style="font-family: courier new,monospace;"><span style="font-family: courier new,monospace;">> agricolae::HSD.test(y, x, DFerror=20, MSError=1)</span><br style="font-family: courier new,monospace;">
<br style="font-family: courier new,monospace;"><span style="font-family: courier new,monospace;">Study:</span><br style="font-family: courier new,monospace;"><br style="font-family: courier new,monospace;"><span style="font-family: courier new,monospace;">HSD Test for y </span><br style="font-family: courier new,monospace;">
<br style="font-family: courier new,monospace;"><span style="font-family: courier new,monospace;">Mean Square Error: 1 </span><br style="font-family: courier new,monospace;"><br style="font-family: courier new,monospace;">
<span style="font-family: courier new,monospace;">x, means</span><br style="font-family: courier new,monospace;"><br style="font-family: courier new,monospace;"><span style="font-family: courier new,monospace;"> y std.err replication</span><br style="font-family: courier new,monospace;">
<span style="font-family: courier new,monospace;">1 -0.3230810 0.1327155 5</span><br style="font-family: courier new,monospace;"><span style="font-family: courier new,monospace;">2 0.3595503 0.5492514 5</span><br style="font-family: courier new,monospace;">
<span style="font-family: courier new,monospace;">3 -0.1703578 0.3377758 5</span><br style="font-family: courier new,monospace;"><span style="font-family: courier new,monospace;">4 -0.1483504 0.5129000 5</span><br style="font-family: courier new,monospace;">
<span style="font-family: courier new,monospace;">5 -0.1101464 0.2727455 5</span><br style="font-family: courier new,monospace;"><br style="font-family: courier new,monospace;"><span style="font-family: courier new,monospace;">alpha: 0.05 ; Df Error: 20 </span><br style="font-family: courier new,monospace;">
<span style="font-family: courier new,monospace;">Critical Value of Studentized Range: 4.231857 </span><br style="font-family: courier new,monospace;"><br style="font-family: courier new,monospace;"><span style="font-family: courier new,monospace;">Honestly Significant Difference: 1.892544 <------------------------------------------ DMS</span><br style="font-family: courier new,monospace;">
<br style="font-family: courier new,monospace;"><span style="font-family: courier new,monospace;">Means with the same letter are not significantly different.</span><br style="font-family: courier new,monospace;"><br style="font-family: courier new,monospace;">
<span style="font-family: courier new,monospace;">Groups, Treatments and means</span><br style="font-family: courier new,monospace;"><span style="font-family: courier new,monospace;">a 2 0.3595503 </span><br style="font-family: courier new,monospace;">
<span style="font-family: courier new,monospace;">a 5 -0.1101464 </span><br style="font-family: courier new,monospace;"><span style="font-family: courier new,monospace;">a 4 -0.1483504 </span><br style="font-family: courier new,monospace;">
<span style="font-family: courier new,monospace;">a 3 -0.1703578 </span><br style="font-family: courier new,monospace;"><span style="font-family: courier new,monospace;">a 1 -0.323081 </span><br style="font-family: courier new,monospace;">
<span style="font-family: courier new,monospace;">> </span><br style="font-family: courier new,monospace;"><br clear="all"></font><span style="font-family: trebuchet ms,sans-serif;">==========================================================================</span><br style="font-family: trebuchet ms,sans-serif;">
<span style="font-family: trebuchet ms,sans-serif;">Walmes Marques Zeviani</span><br style="font-family: trebuchet ms,sans-serif;"><span style="font-family: trebuchet ms,sans-serif;">LEG (Laboratório de Estatística e Geoinformação, 25.450418 S, 49.231759 W)</span><br style="font-family: trebuchet ms,sans-serif;">
<span style="font-family: trebuchet ms,sans-serif;">Departamento de Estatística - Universidade Federal do Paraná</span><br style="font-family: trebuchet ms,sans-serif;"><span style="font-family: trebuchet ms,sans-serif;">fone: (+55) 41 3361 3573</span><br style="font-family: trebuchet ms,sans-serif;">
<span style="font-family: trebuchet ms,sans-serif;">VoIP: (3361 3600) 1053 1173</span><br style="font-family: trebuchet ms,sans-serif;"><span style="font-family: trebuchet ms,sans-serif;">e-mail: <a href="mailto:walmes@ufpr.br" target="_blank">walmes@ufpr.br</a></span><br style="font-family: trebuchet ms,sans-serif;">
<span style="font-family: trebuchet ms,sans-serif;">twitter: @walmeszeviani</span><br style="font-family: trebuchet ms,sans-serif;"><span style="font-family: trebuchet ms,sans-serif;">homepage: <a href="http://www.leg.ufpr.br/%7Ewalmes" target="_blank">http://www.leg.ufpr.br/~walmes</a></span><br style="font-family: trebuchet ms,sans-serif;">
<span style="font-family: trebuchet ms,sans-serif;">linux user number: 531218</span><br style="font-family: trebuchet ms,sans-serif;"><span style="font-family: trebuchet ms,sans-serif;">==========================================================================</span><br>