[R-br] dic parcela subdividida no tempo

classic Classic list List threaded Threaded
2 messages Options
Reply | Threaded
Open this post in threaded view
|

[R-br] dic parcela subdividida no tempo

R-br mailing list

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
[hidden email]
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.
Reply | Threaded
Open this post in threaded view
|

Re: [R-br] dic parcela subdividida no tempo

R-br mailing list
Seu exemplo não é reproduzível.
Sua pergunta não está clara.

À disposição.
Walmes.


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