<html><head><style type='text/css'>p { margin: 0; }</style></head><body><div style='font-family: times new roman,new york,times,serif; font-size: 12pt; color: #000000'><p class="MsoNormal" style="margin: 0cm 0cm 0.0001pt; font-size: medium; font-family: Calibri, sans-serif;">Trata-se de um experimento deduzas temperaturas e seis períodos de avaliação, em DIC. Resolvi conforme mostrado a seguir.</p><p class="MsoNormal" style="margin: 0cm 0cm 0.0001pt; font-size: medium; font-family: Calibri, sans-serif;">Gostaria de saber se está correto. As analises ficaram boas.</p><p class="MsoNormal" style="margin: 0cm 0cm 0.0001pt; font-size: medium; font-family: Calibri, sans-serif;"><br></p><p class="MsoNormal" style="margin: 0cm 0cm 0.0001pt; font-size: medium; font-family: Calibri, sans-serif;"><span lang="EN-US"><br></span></p><p class="MsoNormal" style="margin: 0cm 0cm 0.0001pt; font-size: medium; font-family: Calibri, sans-serif;"><span lang="EN-US">pkg= c("lattice", "latticeExtra", "doBy", "multcomp", "reshape", "plyr", "nlme") <o:p></o:p></span></p><p class="MsoNormal" style="margin: 0cm 0cm 0.0001pt; font-size: medium; font-family: Calibri, sans-serif;"><span lang="EN-US"> </span></p><p class="MsoNormal" style="margin: 0cm 0cm 0.0001pt; font-size: medium; font-family: Calibri, sans-serif;"><span lang="EN-US">sapply(pkg, require, character.only=TRUE)<o:p></o:p></span></p><p class="MsoNormal" style="margin: 0cm 0cm 0.0001pt; font-size: medium; font-family: Calibri, sans-serif;"><span lang="EN-US"> </span></p><p class="MsoNormal" style="margin: 0cm 0cm 0.0001pt; font-size: medium; font-family: Calibri, sans-serif;"><span lang="EN-US">source("http://dl.dropboxusercontent.com/u/48140237/bandas.R")<o:p></o:p></span></p><p class="MsoNormal" style="margin: 0cm 0cm 0.0001pt; font-size: medium; font-family: Calibri, sans-serif;"><span lang="EN-US"> </span></p><p class="MsoNormal" style="margin: 0cm 0cm 0.0001pt; font-size: medium; font-family: Calibri, sans-serif;"><span lang="EN-US">sessionInfo()<o:p></o:p></span></p><p class="MsoNormal" style="margin: 0cm 0cm 0.0001pt; font-size: medium; font-family: Calibri, sans-serif;"><span lang="EN-US"> </span></p><p class="MsoNormal" style="margin: 0cm 0cm 0.0001pt; font-size: medium; font-family: Calibri, sans-serif;"><span lang="EN-US">trellis.device(color=FALSE)<o:p></o:p></span></p><p class="MsoNormal" style="margin: 0cm 0cm 0.0001pt; font-size: medium; font-family: Calibri, sans-serif;"><span lang="EN-US"> </span></p><p class="MsoNormal" style="margin: 0cm 0cm 0.0001pt; font-size: medium; font-family: Calibri, sans-serif;"><span lang="EN-US">ang=read.csv("ang.csv")<o:p></o:p></span></p><p class="MsoNormal" style="margin: 0cm 0cm 0.0001pt; font-size: medium; font-family: Calibri, sans-serif;"><span lang="EN-US">str(ang)<o:p></o:p></span></p><p class="MsoNormal" style="margin: 0cm 0cm 0.0001pt; font-size: medium; font-family: Calibri, sans-serif;"><span lang="EN-US">dh=melt(subset(ang, select=1:8), id.vars=1:2)<o:p></o:p></span></p><p class="MsoNormal" style="margin: 0cm 0cm 0.0001pt; font-size: medium; font-family: Calibri, sans-serif;"><span lang="EN-US">d=diff(as.Date(c("2017-01-01","2017-01-16", "2017-01-31","2017-02-15", "2017-03-02","2017-03-17"))); d<o:p></o:p></span></p><p class="MsoNormal" style="margin: 0cm 0cm 0.0001pt; font-size: medium; font-family: Calibri, sans-serif;">d =cumsum(c(15, d)); d<o:p></o:p></p><p class="MsoNormal" style="margin: 0cm 0cm 0.0001pt; font-size: medium; font-family: Calibri, sans-serif;">dh$dia=d[as.integer(dh$variable)]<o:p></o:p></p><p class="MsoNormal" style="margin: 0cm 0cm 0.0001pt; font-size: medium; font-family: Calibri, sans-serif;"><span lang="EN-US">levels(dh$variable)<o:p></o:p></span></p><p class="MsoNormal" style="margin: 0cm 0cm 0.0001pt; font-size: medium; font-family: Calibri, sans-serif;"><span lang="EN-US">xyplot(value~dia|Temperatura, groups=Temperatura, data=dh, type=c("p","a"), auto.key=TRUE)<o:p></o:p></span></p><p class="MsoNormal" style="margin: 0cm 0cm 0.0001pt; font-size: medium; font-family: Calibri, sans-serif;"><span lang="EN-US"> </span></p><p class="MsoNormal" style="margin: 0cm 0cm 0.0001pt; font-size: medium; font-family: Calibri, sans-serif;"><o:p> </o:p></p><p class="MsoNormal" style="margin: 0cm 0cm 0.0001pt; font-size: medium; font-family: Calibri, sans-serif;"><span lang="EN-US">dh1=transform(dh, D=factor(dia), parc=interaction(Rep, Temperatura), subp=interaction(Rep, Temperatura))<o:p></o:p></span></p><p class="MsoNormal" style="margin: 0cm 0cm 0.0001pt; font-size: medium; font-family: Calibri, sans-serif;"><span lang="EN-US">str(dh1)<o:p></o:p></span></p><p class="MsoNormal" style="margin: 0cm 0cm 0.0001pt; font-size: medium; font-family: Calibri, sans-serif;"><span lang="EN-US"> </span></p><p class="MsoNormal" style="margin: 0cm 0cm 0.0001pt; font-size: medium; font-family: Calibri, sans-serif;"><span lang="EN-US">m2=lme(value~Temperatura*D, random=~1|parc/subp, data=dh1, method="ML") <o:p></o:p></span></p><p class="MsoNormal" style="margin: 0cm 0cm 0.0001pt; font-size: medium; font-family: Calibri, sans-serif;"><span lang="EN-US">m3=lme(value~Temperatura*D, random=~1|parc/subp, data=dh1, method="REML") <o:p></o:p></span></p><br>-- <br><div><span name="x"></span>Andrea Vita Reis Mendonça<br>Profa. Associado I<br>UFRB/CCAAB<span name="x"></span><br></div></div>
<br><br></body></html>