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