[R-br] sugestão análise experimento no tempo

[Wagner Bonat]

O que vc pretende fazer com a estrutura multivariada?
Vc esta assumindo que os tempos e as medidas repetidas são independentes?

Talvez realmente sejam, mas seria interessante testar. Uma estrutura
simples tipo compound symmetry, talvez já seja suficiente pra verificar se
existe alguma indicação de dependencia.


Em 25 de outubro de 2015 23:21, Maurício Lordêlo <mslordelo em gmail.com>
escreveu:

> Caros,
> Recebi os dados de um experimento sem a orientação prévia de um
> profissional
> da área de Estatística.
>
> Foram seis grupos de plantas, sendo que cada grupo tinha inicialmente 15
> plantas.
>
> Um desses grupos foi o controle e em cada um dos outros cinco grupos foi
>  aplicado um regulador de crescimento (hormônio).
>
> O esquema de aplicação/avaliação  de cada tratamento foi o seguinte:
>
> 1. Tempo zero: aplicação de cada um dos reguladores nos grupos;
>
> 2. Após 7 dias, sorteio de  5 plantas por grupo e realizadas as medições
>
> 3. 7 dias depois, aplicação dos hormônios nas 10 plantas restantes
>
> 4. 7 dias depois, sorteio de 5 plantas por grupo e realizadas as medições
>
> 5. 7 dias depois, aplicação dos hormônios nas 5 plantas restantes
>
> 6. 7 dias depois, medições nas 5 plantas restantes.
>
> Em cada avaliação foram analisadas as variáveis resposta: altura,
> comprimento da raiz, peso de folhas, ...
> O objetivo é avaliar qual dos hormônios produziu melhores resultados ao
> longo do tempo.
> Fiz alguns gráficos e ajustei um modelo considerando um fatorial duplo.
> Devido a forma como foi realizado o delineamento, não sei se o fatorial é
> adequado.
> Alguém tem sugestões de como analisar estes dados?
>
> https://www.dropbox.com/s/lkhzb87gh1vtgas/dados.txt?dl=0
>
> dados = read.table("dados.txt",sep="",header=TRUE)
>
> attach(dados)
> summary(dados)
>
> plot(altura~hormonio)
> plot(altura~tempo)
>
> require(latticeExtra)
> xyplot(altura~hormonio|tempo)
> xyplot(altura~tempo|hormonio)
>
> modelo = aov(altura ~ hormonio*tempo)
> summary(modelo)
> plot(modelo)
> residuos = modelo$residuals
> shapiro.test(residuos)
>
>
>
> _______________________________________________
> R-br mailing list
> R-br em 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.
>



-- 
Wagner Hugo Bonat
----------------------------------------------------------------------------------------------
Department of Mathematics and Computer Science (IMADA)
University of Southern Denmark (SDU) and
Laboratório de Estatística e Geoinformação (LEG)
Universidade Federal do Paraná (UFPR)
-------------- Próxima Parte ----------
Um anexo em HTML foi limpo...
URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20151026/378f9e7e/attachment.html>