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?
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)