Modelos

# MODELOS

end=file.choose() # abrir o arquivo ic_medias
medias=read.table(end,sep=" ",header=T)


m0=lm(bm~1,medias)
m1=lm(bm~sp,medias)
anova(m0,m1)

m2=lm(bm~sp+light,medias)
anova(m1,m2)

m3=lm(bm~sp+light+bloco,medias)
anova(m2,m3)

# conclusão: o melhor modelo é o que inclui espécie, luz e bloco

# porém, não sei se o fato de os tratamentos estarem aninhado em bloco inviesa esse resultado