DIC e parcelas subdivididas

Ver o tópico anterior Ver o tópico seguinte Ir em baixo

DIC e parcelas subdivididas

Mensagem  Ítalo L F em Ter Jul 03, 2012 11:19 pm

Olá,

Tenho um delineamento inteiramente casualizado (DIC) num experimento em parcelas subdivididas. Entretanto meu quadro ANOVA está errado. Bem provavelmente por inserção dos comandos errados.

Meu comandos são:

### Inserção dos dados ###
Embalagem<-factor(rep(1:2,each=20),labels=c("pote","saco"))
Tempo<-factor(rep(1:5,each=4,2),labels=c("0d","20d","40d","60d","80d"))
Repeticao<-factor(rep(1:4,10),labels=c("r1","r2","r3","r4"))
Resposta<-c(4.93, 4.77, 4.84, 4.83, 5.08, 5.03, 5.05, 5.05, 4.99, 4.99, 5.01, 4.98, 4.85, 4.81, 4.85, 4.8, 4.92, 4.95, 4.9, 4.96, 4.72, 4.84, 4.65, 4.9, 4.98, 5.04, 5.09, 5.12, 4.95, 4.97, 4.95, 5.01, 4.72, 4.93, 5, 4.84, 4.85, 4.85, 4.9, 4.9)
pH<-data.frame(Embalagem,Tempo,Repeticao,Resposta)
attach(pH)

### ANOVA ###
av<-aov(Resposta~Embalagem+Tempo+Embalagem*Tempo+Error(Repeticao:Embalagem),pH)
summary(av)
plot(av$fit,av$res, xlab="Valores Ajustados",ylab="Resíduos")
qqnorm(av$res,xlab="Quantil da Normal",ylab="Resíduos")
qqline(av$res)
bartlett.test(av$res,Embalagem)
bartlett.test(av$res,Tempo)
shapiro.test(av$res)

1 - A ANOVA apresenta um error de modelo singular, destacados na figura;
2 - O quadro anova apresenta SQ Erros errados, apesar de ter os graus de liberdade ok. Destaque amarelo na figura; A tabela do SISVAR (outro programa estatístico) está correta;
3 - Preciso executar os teste de normalidade dos resíduos e homogeneidade das variâncias, tbm não estou conseguindo fazer isso.

Ver imagem em: http://imageshack.us/photo/my-images/839/anova.jpg/


Ítalo L F

Mensagens : 1
Data de inscrição : 03/07/2012

Ver perfil do usuário

Voltar ao Topo Ir em baixo

Re: DIC e parcelas subdivididas

Mensagem  gustavoreis em Qua Jul 04, 2012 11:48 pm

Boa noite Ítalo!

Por enquanto só pude analisar até o summary(av), ainda não descobri o motivo da diferença, apenas verifiquei alguns comandos que você pode eliminar:

pH<-data.frame(Embalagem,Tempo,Repeticao,Resposta) # não vejo necessidade de criar esse data.frame
attach(pH) # Mesmo que você crie o data.frame não há necessidade de fazer esse attach, pois as variáveis já existem no R (é o que diz na msg que aparece quando utiliza-se esse comando)

Eliminando os dois comandos acima, o modelo ficaria assim:

av<-aov(Resposta~Embalagem*Tempo+Error(Repeticao:Embalagem))

Veja que também eliminei Embalagem + Tempo, pois:

Embalagem*Tempo = Embalagem + Tempo + Embalagem:Tempo


Todas as modificações que fiz não irão alterar nada no resultado final, é somente uma forma mais "limpa" de se trabalhar.

A dica que posso dar é que verifique item por item na variável Repeticao, é possível que tenha algum dado diferente do que você utilizou no outro programa.

gustavoreis
Admin

Mensagens : 70
Data de inscrição : 01/03/2011
Idade : 31

Ver perfil do usuário http://estatisticanor.forumeiros.com

Voltar ao Topo Ir em baixo

Ver o tópico anterior Ver o tópico seguinte Voltar ao Topo

- Tópicos similares

 
Permissão deste fórum:
Você não pode responder aos tópicos neste fórum