
*Marcelo,* * * *O Ivan, professor da UESC/BA, preparou o script anexo para as análises do experimento que informei anteriormente.* * * *E ele me autorizou a socializá-lo. Então, segue anexo, o script R que ele preparou. Se alguém utilizar, não se esqueçam de dar os créditos a ele.* * * *Aproveito para agradecer ao prof. Ivan, a autorização.* * * *Luiz Roberto* Luiz Roberto Martins Pinto Prof. Pleno/DCET/UESC luizroberto.uesc@gmail.com skype: lrmpinto http://lattes.cnpq.br/2732314327604831 Em 12 de março de 2013 11:32, Marcelo Laia <marcelolaia@gmail.com> escreveu:
Colega,
A dias que venho lendo a respeito de como realizar uma análise de medidas repetidas no tempo. Encontrei uma farta bibliografia e diversas sugestões de como fazer.
Eis um dado mínimo:
datafilename="http://dl.dropbox.com/u/34009642/Dados_Omissao.csv" data.min <- read.table(datafilename, header=T, sep="\t", dec=",") head(data.min) data.min <- within(data.min, { Clone <- factor(Clone) Tempo <- factor(Tempo) Trat <- factor(Trat) Planta <- factor(Planta) }) summary(data.min)
A partir desses dados e do que eu li, tentei:
dados_omissao_h.aov <- aov(h ~ Clone * Trat * Tempo + Error(Planta), data = data_min) # http://tinyurl.com/bjlas86
Esse http://pealco.net/2009/01/30/repeated-measures-anova-in-r.html dá outra abordagem.
Eu também rodei o código abaixo.
library(nlme) lme(h ~ Clone+Trat+Tempo,random=~1|Planta,data=data.min) #
Ainda li os textos seguintes:
http://egret.psychol.cam.ac.uk/statistics/R/anova.html#8.4
http://egret.psychol.cam.ac.uk/statistics/R/enteringdata.html#reshape
http://www.psych.yorku.ca/cribbie/6130/ex14_rscript.R
No entanto, sinceramente, não sei qual abordagem utilizar nos nossos dados. Desejamos verificar se há diferenças entre clones, entre tratamentos, a interação clone tratamento e a evolução ao longo do tempo das variáveis. Ainda, há plantas que morreram ao longo do experimento devido ao próprio tratamento. E ainda tem a dica do Wlames no ridiculas para dados binomial!
Alguma sugestão?
-- Marcelo Brasil (Brazil, for English Speakers) Linux user number 487797 _______________________________________________ R-br mailing list R-br@listas.c3sl.ufpr.br 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.