Não foi possível enviar o arquivo. Será algum problema com as permissões?
Diferenças
Aqui você vê as diferenças entre duas revisões dessa página.
| — |
pessoais:codigocpi [2010/12/02 14:55] (atual) mmerces criada |
||
|---|---|---|---|
| Linha 1: | Linha 1: | ||
| + | ########################################################################################## | ||
| + | ########################################################################################## | ||
| + | ##### ANÁLISE DE EXPERIMENTO FATORIAL FRACIONADO 2^(5-1) ##### | ||
| + | ########################################################################################## | ||
| + | ########################################################################################## | ||
| + | #Leitura dos dados: | ||
| + | A<-c(,-1,1,-1,1,-1,1,-1,1,-1,1,-1,1,-1,1,-1,1,-1,1,-1,1,-1,1,-1,1,-1,1,-1,1,-1,1,-1,1) | ||
| + | B<-c(-1,-1,1,1,-1,-1,1,1,-1,-1,1,1,-1,-1,1,1,-1,-1,1,1,-1,-1,1,1,-1,-1,1,1,-1,-1,1,1) | ||
| + | C<-c(-1,-1,-1,-1,1,1,1,1,-1,-1,-1,-1,1,1,1,1,-1,-1,-1,-1,1,1,1,1,-1,-1,-1,-1,1,1,1,1) | ||
| + | D<-c(-1,-1,-1,-1,-1,-1,-1,-1,1,1,1,1,1,1,1,1,-1,-1,-1,-1,-1,-1,-1,-1,1,1,1,1,1,1,1,1) | ||
| + | E<-c(-1,-1,-1,-1,-1,-1,-1,-1,-1,-1,-1,-1,-1,-1,-1,-1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1) | ||
| + | VR<-c(10,5,10,5,10,18,11,18,16,7,15,7,16,35,16,18,30,9,30,10,11,28,11,28,25,10,25,10,25,35,25,25) | ||
| + | ABCDE<-c(-1,1,1,-1,1,-1,-1,1,1,-1,-1,1,-1,1,1,-1,1,-1,-1,1,-1,1,1,-1,-1,1,1,-1,1,-1,-1,1) | ||
| + | dados<-data.frame(A,B,C,D,E,ABCDE,VR) | ||
| + | dados<-as.data.frame(dados) | ||
| + | |||
| + | #Ajuste do modelo saturado; | ||
| + | ajust1<-lm(VR~A*B*C*D*E) | ||
| + | summary(ajust1) | ||
| + | anova(ajust1) | ||
| + | |||
| + | #Ajuste do modelo com interações duplas; | ||
| + | ajust2<-lm(VR~A+B+C+D+E+A*B+A*C+A*D+A*E+B*C+B*D+B*E+C*D+C*E+D*E,data=dados) | ||
| + | summary(ajust2) | ||
| + | anova(ajust2) | ||
| + | |||
| + | #Ajuste do modelo com interações duplas e com efeito de bloco; | ||
| + | ajust3<-lm(VR~A+B+C+D+E+ABCDE + A*B+A*C+A*D+A*E+B*C+B*D+B*E+C*D+C*E+D*E,data=dados) | ||
| + | summary(ajust3) | ||
| + | anova(ajust3) | ||
| + | |||
| + | |||
| + | #Ajuste do modelo selecionado | ||
| + | |||
| + | #Interações do modelo selecionado | ||
| + | interaction.plot() | ||
