[R-br] dic parcela subdividida no tempo

[ANDREA VITA REIS MENDONCA]

Trata-se de um experimento deduzas temperaturas e seis períodos de avaliação, em DIC. Resolvi conforme mostrado a seguir. 
Gostaria de saber se está correto. As analises ficaram boas. 




pkg= c("lattice", "latticeExtra", "doBy", "multcomp", "reshape", "plyr", "nlme") 

sapply(pkg, require, character.only=TRUE) 

source("http://dl.dropboxusercontent.com/u/48140237/bandas.R") 

sessionInfo() 

trellis.device(color=FALSE) 

ang=read.csv("ang.csv") 
str(ang) 
dh=melt(subset(ang, select=1:8), id.vars=1:2) 
d=diff(as.Date(c("2017-01-01","2017-01-16", "2017-01-31","2017-02-15", "2017-03-02","2017-03-17"))); d 
d =cumsum(c(15, d)); d 
dh$dia=d[as.integer(dh$variable)] 
levels(dh$variable) 
xyplot(value~dia|Temperatura, groups=Temperatura, data=dh, type=c("p","a"), auto.key=TRUE) 


dh1=transform(dh, D=factor(dia), parc=interaction(Rep, Temperatura), subp=interaction(Rep, Temperatura)) 
str(dh1) 

m2=lme(value~Temperatura*D, random=~1|parc/subp, data=dh1, method="ML") 
m3=lme(value~Temperatura*D, random=~1|parc/subp, data=dh1, method="REML") 
-- 

Andrea Vita Reis Mendonça 
Profa. Associado I 
UFRB/CCAAB 

-------------- Próxima Parte ----------
Um anexo em HTML foi limpo...
URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20190920/0200beef/attachment.html>