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)