[R-br] Particionar soma de quadrado

[Andre Oliveira]

Boa noite,

 dados que houve a internação entre fatores, preciso particionar a soma de quadrado. É apenas um exercício que estou fazendo co dados que está na web. Evidentemente que posso construir os contraste e testar o que me interessa e também poderia usar a ExpDes para isto. Alguém poderia ajudar com um script para tal desdobramento?
Grato
Pequei este exemplo! 

 da <- read.table("http://www.leg.ufpr.br/~walmes/cursoR/anovareg.txt", header=TRUE, sep="\t")
 str(da)

 attach(da)
# Este exemplo testa regressão! 

fit=aov(indice ~ cultivar*ordered(dose))
summary(fit)
summary.lm(fit)


# Gostaria de ajuda para desdobrar e testar doses para este modelo abaixo! 

 fit=aov(indice ~ cultivar*factor(dose)+Error(bloco:(cultivar+dose))

 summary(fit)
Grato


-------------- Próxima Parte ----------
Um anexo em HTML foi limpo...
URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20190108/1b0d0316/attachment.html>