<html>
<head>
<meta http-equiv="Content-Type" content="text/html; charset=iso-8859-1">
<style type="text/css" style="display:none;"><!-- P {margin-top:0;margin-bottom:0;} --></style>
</head>
<body dir="ltr">
<div id="divtagdefaultwrapper" style="font-size:12pt;color:#000000;font-family:'Times New Roman', Times, serif;" dir="ltr">
<p>Bom dia pessoal, </p>
<p>Mais uma vez venho recorrer a vocês em relação a uma dúvida...</p>
<p><br>
</p>
<p>Estou trabalhando com um experimento em DBC em faixas, porém tenho encontrado poucos exemplos em relação a comparação de médias no caso de interação. </p>
<p><br>
</p>
<p>Encontrei um exemplo em uma aula na internet e um exemplo com os comandos no pacote agricolae (<a href="http://127.0.0.1:25234/library/agricolae/html/strip.plot.html" class="OWAAutoLink" id="LPlnk17600" previewremoved="true">http://127.0.0.1:25234/library/agricolae/html/strip.plot.html</a>),
porém mesmo usando os mesmos valores a comparação de médias não está batendo entre o exemplo da aula e quando rodo usando os comandos do exemplo do agricolae. </p>
<p><br>
</p>
<p>Alguém ai poderia conferir se esses comandos que envio logo a baixo estão corretos? Principalmente na parte referente a comparação de médias.</p>
<p><br>
</p>
<p>(Seria um experimento em blocos completos no esquema em faixas.)</p>
<p><br>
</p>
<p>Desde já muito obrigado!</p>
<p><br>
</p>
<p>Abraços</p>
<p><br>
</p>
<p><span style="font-size: 12pt;"><br>
</span></p>
<p><span style="font-size: 12pt;"><br>
</span></p>
<p><span style="font-size: 12pt;"># Yield</span><br>
</p>
<p></p>
<pre>library(agricolae)
data(huasahuasi)
YIELD<-huasahuasi$YIELD
market <- YIELD$y1da + YIELD$y2da
non_market <- YIELD$y3da
yield <- market + non_market
model<-with(YIELD,strip.plot(block, clon, trt, yield))
comparison<-with(YIELD,LSD.test(yield,clon,model$gl.a,model$Ea))
comparison<-with(YIELD,LSD.test(yield,trt,model$gl.b,model$Eb))
# simple effects
A<-model$data
a<-nlevels(A$clon)
b<-nlevels(A$trt)
r<-nlevels(A$block)
Ea<-model$Ea; Eb<-model$Ec; Ec<-model$Ec;
gla<-model$gl.a; glb<-model$gl.b; glc<-model$gl.c;
B <-tapply.stat(A[,4],A[,2:3],mean)
std<-tapply.stat(A[,4],A[,2:3],function(x) sd(x)/sqrt(length(x)))
B<-data.frame(B[,1:2],yield=B[,3],std=std[,3])
cmab<-(b-1)*Ec + Ea
cmba<-(a-1)*Ec + Eb
# order.group
# Tukey
ta <- qtukey(0.95,a,gla)
tb <- qtukey(0.95,b,glb)
tc <- qtukey(0.95,a*b,glc)
tab<- ((b-1)*Ec*tc + Ea*ta)/cmab
tba<- ((a-1)*Ec*tc + Eb*tb)/cmba
# Comparison of clones by treatment
groups<-by(B,B[,2], function(x) order.group(x$clon,x$yield,N=b*r, cmab,tab,
std.err=x$std,parameter=0.5,console=FALSE))
groups
# Comparison of treatments by clon
groups<-by(B,B[,1], function(x) order.group(x$trt,x$yield,N=a*r,cmba,tba,
std.err=x$std,parameter=0.5,console=FALSE))
groups</pre>
<br>
<p></p>
<p><br>
</p>
<div id="Signature">
<div id="divtagdefaultwrapper" style="font-size:12pt; color:#000000; background-color:#FFFFFF; font-family:'Times New Roman',Times,serif">
<font face="Times New Roman" size="3" style="font-size:12pt"><span class="ecxyiv1220415869ecxApple-style-span" style="line-height:20px; background-color:rgb(255,255,255)">
<div style="line-height:17px; text-align:center"><font style="line-height:normal"><b style="line-height:17px">..................</b></font></div>
<div style="line-height:17px; text-align:center"><font style="line-height:normal"><b style="line-height:17px">Me. Eng. Agr. Maurício Sangiogo</b></font></div>
<div style="line-height:17px; text-align:center"><font class="ecxyiv1220415869ecxApple-style-span" style="line-height:normal">TEL:(55) 96822030</font></div>
<div style="line-height:17px; text-align:center"><span class="ecxyiv1220415869ecxApple-style-span" style="font-size:12pt"> </span><span class="ecxyiv1220415869ecxApple-style-span" style="font-size:12pt"></span></div>
</span></font></div>
</div>
</div>
</body>
</html>