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
_______________________________________________
R-br mailing list
[email protected]
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.
