Re: [R-br] [R_STAT] Seleção de Quadrados Médios
Fernando, (agora mandei para a lista nova :P) O problema é que o comando aov() não aceita esse procedimento... Se eu faço uma Anova do fit ele da a seguinte mensagem de erro:
anova(fit) Erro em UseMethod("anova") : método não aplicável para 'anova' aplicado a um objeto de classe "c('aovlist', 'listof')"
2011/4/23 Fernando Henrique Toledo <fernandohtoledo@gmail.com>
Greice, bom dia.
Acesse:
anova(fit)[,3]
Isso resultará no vetor dos quadrados médios! a 1ª coluna é dos graus de liberdade, a segunda das somas de quadrados, 3ª dos quadrados médios ...
PS: Mudamos de lista, mude p/ a "r-br", esta em breve será desativada!
att, FH
2011/4/23 greice.laureano <greice.laureano@gmail.com>
Bom dia,
Gostaria de ajuda para extrair os quadrados médios de uma Anova. O banco de dados está abaixo assim como o código para fazer a Anova. Muito Obrigado;
Greice.
Dados:
sujeito aferidor tempo valor 1 1 1 -2 1 1 2 0 1 1 3 1 1 2 1 2 1 2 2 1 1 2 3 1 2 1 1 16 2 1 2 16 2 1 3 15 2 2 1 12 2 2 2 14 2 2 3 13 3 1 1 5 3 1 2 6 3 1 3 6 3 2 1 4 3 2 2 4 3 2 3 4 4 1 1 11 4 1 2 10 4 1 3 10 4 2 1 9 4 2 2 7 4 2 3 8 5 1 1 7 5 1 2 8 5 1 3 6 5 2 1 5 5 2 2 6 5 2 3 6 6 1 1 -7 6 1 2 -8 6 1 3 -8 6 2 1 -9 6 2 2 -10 6 2 3 -9 7 1 1 18 7 1 2 19 7 1 3 19 7 2 1 17 7 2 2 17 7 2 3 17 8 1 1 4 8 1 2 5 8 1 3 5 8 2 1 5 8 2 2 5 8 2 3 5 9 1 1 0 9 1 2 -3 9 1 3 -2 9 2 1 -7 9 2 2 -6 9 2 3 -5 10 1 1 0 10 1 2 0 10 1 3 -2 10 2 1 1 10 2 2 2 10 2 3 1 11 1 1 -3 11 1 2 -2 11 1 3 -2 11 2 1 -4 11 2 2 -3 11 2 3 -3 12 1 1 3 12 1 2 -1 12 1 3 1 12 2 1 -1 12 2 2 -2 12 2 3 1 13 1 1 7 13 1 2 9 13 1 3 9 13 2 1 4 13 2 2 4 13 2 3 2 14 1 1 -6 14 1 2 -7 14 1 3 -6 14 2 1 -8 14 2 2 -10 14 2 3 -9 15 1 1 1 15 1 2 1 15 1 3 0 15 2 1 -2 15 2 2 -2 15 2 3 -3 16 1 1 -13 16 1 2 -14 16 1 3 -14 16 2 1 -12 16 2 2 -12 16 2 3 -12 17 1 1 2 17 1 2 1 17 1 3 0 17 2 1 -1 17 2 2 0 17 2 3 0 18 1 1 4 18 1 2 4 18 1 3 3 18 2 1 7 18 2 2 6 18 2 3 4 19 1 1 -10 19 1 2 -9 19 1 3 -10 19 2 1 -10 19 2 2 -11 19 2 3 -10 20 1 1 8 20 1 2 9 20 1 3 8 20 2 1 2 20 2 2 8 20 2 3 8 21 1 1 7 21 1 2 6 21 1 3 7 21 2 1 8 21 2 2 7 21 2 3 7 22 1 1 -3 22 1 2 -2 22 1 3 -4 22 2 1 -5 22 2 2 -5 22 2 3 -5 23 1 1 -5 23 1 2 -5 23 1 3 -7 23 2 1 -6 23 2 2 -8 23 2 3 -7 24 1 1 5 24 1 2 5 24 1 3 5 24 2 1 3 24 2 2 4 24 2 3 4 25 1 1 0 25 1 2 -1 25 1 3 -1 25 2 1 -4 25 2 2 -3 25 2 3 -4 26 1 1 7 26 1 2 6 26 1 3 6 26 2 1 4 26 2 2 4 26 2 3 4 27 1 1 -8 27 1 2 -8 27 1 3 -8 27 2 1 -10 27 2 2 -11 27 2 3 -10 28 1 1 1 28 1 2 1 28 1 3 2 28 2 1 1 28 2 2 -1 28 2 3 0 29 1 1 -3 29 1 2 -3 29 1 3 -3 29 2 1 -5 29 2 2 -4 29 2 3 -5
sj=as.factor(dados$sujeito) af=as.factor(dados$aferidor) tempo=as.factor(dados$tempo)
fit <- aov(dados$valor ~ sj*af + sj:tempo:af + Error(sj/af)) summary(fit)
__._,_.___
<fernandohtoledo@gmail.com?subject=Res%3A%20Re%3A%20%5BR_STAT%5D%20Sele%E7%E3o%20de%20Quadrados%20M%E9dios>| através de email<R_STAT@yahoogrupos.com.br?subject=Res%3A%20Re%3A%20%5BR_STAT%5D%20Sele%E7%E3o%20de%20Quadrados%20M%E9dios>| Responder através da web<http://br.groups.yahoo.com/group/R_STAT/post;_ylc=X3oDMTJyNWprYWVzBF9TAzk3NDkwNDM3BGdycElkAzExOTIzMjc1BGdycHNwSWQDMjEzNzExMTYwNQRtc2dJZAMxNjA0MARzZWMDZnRyBHNsawNycGx5BHN0aW1lAzEzMDM1Njc5NDA-?act=reply&messageNum=16040>| Adicionar um novo tópico<http://br.groups.yahoo.com/group/R_STAT/post;_ylc=X3oDMTJmNDRncmJ1BF9TAzk3NDkwNDM3BGdycElkAzExOTIzMjc1BGdycHNwSWQDMjEzNzExMTYwNQRzZWMDZnRyBHNsawNudHBjBHN0aW1lAzEzMDM1Njc5NDA-> Mensagens neste tópico<http://br.groups.yahoo.com/group/R_STAT/message/16039;_ylc=X3oDMTM3N2Q4Y2ViBF9TAzk3NDkwNDM3BGdycElkAzExOTIzMjc1BGdycHNwSWQDMjEzNzExMTYwNQRtc2dJZAMxNjA0MARzZWMDZnRyBHNsawN2dHBjBHN0aW1lAzEzMDM1Njc5NDAEdHBjSWQDMTYwMzk->( 2) Atividade nos últimos dias:
Visite seu Grupo<http://br.groups.yahoo.com/group/R_STAT;_ylc=X3oDMTJmbDdmMWM3BF9TAzk3NDkwNDM3BGdycElkAzExOTIzMjc1BGdycHNwSWQDMjEzNzExMTYwNQRzZWMDdnRsBHNsawN2Z2hwBHN0aW1lAzEzMDM1Njc5NDA-> Não deixe de ir atrás da pessoa dos seus sonhos!<http://global.ard.yahoo.com/SIG=15me74fii/M=758712.14532720.14422460.12960164/D=brclubs/S=2137111605:MKP1/Y=BR/EXP=1303575141/L=b473c190-6db3-11e0-860f-2b0ed2bf8b7f/B=f6ybi0wNO6Q-/J=1303567941524925/K=bkkHq4.Jix.KKZcJHGPNcw/A=6351151/R=0/id=mkp1/SIG=135g0plhh/*http://tracking.parperfeito.com.br/ppbanner/bannerTracker?originId=1&identifierId=461809&actionId=1> [image: Yahoo! Grupos]<http://br.groups.yahoo.com/;_ylc=X3oDMTJlc2RmczJyBF9TAzk3NDkwNDM1BGdycElkAzExOTIzMjc1BGdycHNwSWQDMjEzNzExMTYwNQRzZWMDZnRyBHNsawNnZnAEc3RpbWUDMTMwMzU2Nzk0MA--> Trocar para: Só Texto<R_STAT-traditional@yahoogrupos.com.br?subject=Mudar+Formato+de+Envio:+Tradicional>, Resenha Diária<R_STAT-digest@yahoogrupos.com.br?subject=Envio+de+email:+Resenha>• Sair do grupo <R_STAT-unsubscribe@yahoogrupos.com.br?subject=Sair+do+grupo> • Termos de uso <http://br.yahoo.com/info/utos.html> .
__,_._,___
Agora que vi... Você está usando o argumento Error(), o que altera a saída para uma lista! se você efetivamente está interessada em resgatar os QM pode ajustar o mesmo modelo sem esse Error() e extrair como te expliquei, afinal com o Error() especificado o que muda são os testes F e os p-valor, os QM permanecem os mesmos! Desculpe pelo confusão! Vale lembrar que ainda existe um meio de se extrair nesse caso, com o aovlist! Vou ver aqui se acho e te envio em tempo. att, FH 2011/4/23 Greice Laureano <greice.laureano@gmail.com>
Fernando, (agora mandei para a lista nova :P)
O problema é que o comando aov() não aceita esse procedimento... Se eu faço uma Anova do fit ele da a seguinte mensagem de erro:
anova(fit) Erro em UseMethod("anova") : método não aplicável para 'anova' aplicado a um objeto de classe "c('aovlist', 'listof')"
2011/4/23 Fernando Henrique Toledo <fernandohtoledo@gmail.com>
Greice, bom dia.
Acesse:
anova(fit)[,3]
Isso resultará no vetor dos quadrados médios! a 1ª coluna é dos graus de liberdade, a segunda das somas de quadrados, 3ª dos quadrados médios ...
PS: Mudamos de lista, mude p/ a "r-br", esta em breve será desativada!
att, FH
2011/4/23 greice.laureano <greice.laureano@gmail.com>
Bom dia,
Gostaria de ajuda para extrair os quadrados médios de uma Anova. O banco de dados está abaixo assim como o código para fazer a Anova. Muito Obrigado;
Greice.
Dados:
sujeito aferidor tempo valor 1 1 1 -2 1 1 2 0 1 1 3 1 1 2 1 2 1 2 2 1 1 2 3 1 2 1 1 16 2 1 2 16 2 1 3 15 2 2 1 12 2 2 2 14 2 2 3 13 3 1 1 5 3 1 2 6 3 1 3 6 3 2 1 4 3 2 2 4 3 2 3 4 4 1 1 11 4 1 2 10 4 1 3 10 4 2 1 9 4 2 2 7 4 2 3 8 5 1 1 7 5 1 2 8 5 1 3 6 5 2 1 5 5 2 2 6 5 2 3 6 6 1 1 -7 6 1 2 -8 6 1 3 -8 6 2 1 -9 6 2 2 -10 6 2 3 -9 7 1 1 18 7 1 2 19 7 1 3 19 7 2 1 17 7 2 2 17 7 2 3 17 8 1 1 4 8 1 2 5 8 1 3 5 8 2 1 5 8 2 2 5 8 2 3 5 9 1 1 0 9 1 2 -3 9 1 3 -2 9 2 1 -7 9 2 2 -6 9 2 3 -5 10 1 1 0 10 1 2 0 10 1 3 -2 10 2 1 1 10 2 2 2 10 2 3 1 11 1 1 -3 11 1 2 -2 11 1 3 -2 11 2 1 -4 11 2 2 -3 11 2 3 -3 12 1 1 3 12 1 2 -1 12 1 3 1 12 2 1 -1 12 2 2 -2 12 2 3 1 13 1 1 7 13 1 2 9 13 1 3 9 13 2 1 4 13 2 2 4 13 2 3 2 14 1 1 -6 14 1 2 -7 14 1 3 -6 14 2 1 -8 14 2 2 -10 14 2 3 -9 15 1 1 1 15 1 2 1 15 1 3 0 15 2 1 -2 15 2 2 -2 15 2 3 -3 16 1 1 -13 16 1 2 -14 16 1 3 -14 16 2 1 -12 16 2 2 -12 16 2 3 -12 17 1 1 2 17 1 2 1 17 1 3 0 17 2 1 -1 17 2 2 0 17 2 3 0 18 1 1 4 18 1 2 4 18 1 3 3 18 2 1 7 18 2 2 6 18 2 3 4 19 1 1 -10 19 1 2 -9 19 1 3 -10 19 2 1 -10 19 2 2 -11 19 2 3 -10 20 1 1 8 20 1 2 9 20 1 3 8 20 2 1 2 20 2 2 8 20 2 3 8 21 1 1 7 21 1 2 6 21 1 3 7 21 2 1 8 21 2 2 7 21 2 3 7 22 1 1 -3 22 1 2 -2 22 1 3 -4 22 2 1 -5 22 2 2 -5 22 2 3 -5 23 1 1 -5 23 1 2 -5 23 1 3 -7 23 2 1 -6 23 2 2 -8 23 2 3 -7 24 1 1 5 24 1 2 5 24 1 3 5 24 2 1 3 24 2 2 4 24 2 3 4 25 1 1 0 25 1 2 -1 25 1 3 -1 25 2 1 -4 25 2 2 -3 25 2 3 -4 26 1 1 7 26 1 2 6 26 1 3 6 26 2 1 4 26 2 2 4 26 2 3 4 27 1 1 -8 27 1 2 -8 27 1 3 -8 27 2 1 -10 27 2 2 -11 27 2 3 -10 28 1 1 1 28 1 2 1 28 1 3 2 28 2 1 1 28 2 2 -1 28 2 3 0 29 1 1 -3 29 1 2 -3 29 1 3 -3 29 2 1 -5 29 2 2 -4 29 2 3 -5
sj=as.factor(dados$sujeito) af=as.factor(dados$aferidor) tempo=as.factor(dados$tempo)
fit <- aov(dados$valor ~ sj*af + sj:tempo:af + Error(sj/af)) summary(fit)
__._,_.___
<fernandohtoledo@gmail.com?subject=Res%3A%20Re%3A%20%5BR_STAT%5D%20Sele%E7%E3o%20de%20Quadrados%20M%E9dios>| através de email<R_STAT@yahoogrupos.com.br?subject=Res%3A%20Re%3A%20%5BR_STAT%5D%20Sele%E7%E3o%20de%20Quadrados%20M%E9dios>| Responder através da web<http://br.groups.yahoo.com/group/R_STAT/post;_ylc=X3oDMTJyNWprYWVzBF9TAzk3NDkwNDM3BGdycElkAzExOTIzMjc1BGdycHNwSWQDMjEzNzExMTYwNQRtc2dJZAMxNjA0MARzZWMDZnRyBHNsawNycGx5BHN0aW1lAzEzMDM1Njc5NDA-?act=reply&messageNum=16040>| Adicionar um novo tópico<http://br.groups.yahoo.com/group/R_STAT/post;_ylc=X3oDMTJmNDRncmJ1BF9TAzk3NDkwNDM3BGdycElkAzExOTIzMjc1BGdycHNwSWQDMjEzNzExMTYwNQRzZWMDZnRyBHNsawNudHBjBHN0aW1lAzEzMDM1Njc5NDA-> Mensagens neste tópico<http://br.groups.yahoo.com/group/R_STAT/message/16039;_ylc=X3oDMTM3N2Q4Y2ViBF9TAzk3NDkwNDM3BGdycElkAzExOTIzMjc1BGdycHNwSWQDMjEzNzExMTYwNQRtc2dJZAMxNjA0MARzZWMDZnRyBHNsawN2dHBjBHN0aW1lAzEzMDM1Njc5NDAEdHBjSWQDMTYwMzk->( 2) Atividade nos últimos dias:
Visite seu Grupo<http://br.groups.yahoo.com/group/R_STAT;_ylc=X3oDMTJmbDdmMWM3BF9TAzk3NDkwNDM3BGdycElkAzExOTIzMjc1BGdycHNwSWQDMjEzNzExMTYwNQRzZWMDdnRsBHNsawN2Z2hwBHN0aW1lAzEzMDM1Njc5NDA-> Não deixe de ir atrás da pessoa dos seus sonhos!<http://global.ard.yahoo.com/SIG=15me74fii/M=758712.14532720.14422460.12960164/D=brclubs/S=2137111605:MKP1/Y=BR/EXP=1303575141/L=b473c190-6db3-11e0-860f-2b0ed2bf8b7f/B=f6ybi0wNO6Q-/J=1303567941524925/K=bkkHq4.Jix.KKZcJHGPNcw/A=6351151/R=0/id=mkp1/SIG=135g0plhh/*http://tracking.parperfeito.com.br/ppbanner/bannerTracker?originId=1&identifierId=461809&actionId=1> [image: Yahoo! Grupos]<http://br.groups.yahoo.com/;_ylc=X3oDMTJlc2RmczJyBF9TAzk3NDkwNDM1BGdycElkAzExOTIzMjc1BGdycHNwSWQDMjEzNzExMTYwNQRzZWMDZnRyBHNsawNnZnAEc3RpbWUDMTMwMzU2Nzk0MA--> Trocar para: Só Texto<R_STAT-traditional@yahoogrupos.com.br?subject=Mudar+Formato+de+Envio:+Tradicional>, Resenha Diária<R_STAT-digest@yahoogrupos.com.br?subject=Envio+de+email:+Resenha>• Sair do grupo <R_STAT-unsubscribe@yahoogrupos.com.br?subject=Sair+do+grupo> • Termos de uso <http://br.yahoo.com/info/utos.html> .
__,_._,___
_______________________________________________ R-br mailing list R-br@listas.c3sl.ufpr.br https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
Greice, achei... Tente esse! Acho que o que possa dar erro é o que deixei colorido, não sei se tem que ser sj/af ou sj:af! summary(anova)$"Error: sj/af"[[1]]$"Mean Sq" Qualquer problema fico a disposição. att, FH 2011/4/23 Fernando Henrique Toledo <fernandohtoledo@gmail.com>
Agora que vi... Você está usando o argumento Error(), o que altera a saída para uma lista! se você efetivamente está interessada em resgatar os QM pode ajustar o mesmo modelo sem esse Error() e extrair como te expliquei, afinal com o Error() especificado o que muda são os testes F e os p-valor, os QM permanecem os mesmos!
Desculpe pelo confusão!
Vale lembrar que ainda existe um meio de se extrair nesse caso, com o aovlist! Vou ver aqui se acho e te envio em tempo.
att, FH
2011/4/23 Greice Laureano <greice.laureano@gmail.com>
Fernando, (agora mandei para a lista nova :P)
O problema é que o comando aov() não aceita esse procedimento... Se eu faço uma Anova do fit ele da a seguinte mensagem de erro:
anova(fit) Erro em UseMethod("anova") : método não aplicável para 'anova' aplicado a um objeto de classe "c('aovlist', 'listof')"
2011/4/23 Fernando Henrique Toledo <fernandohtoledo@gmail.com>
Greice, bom dia.
Acesse:
anova(fit)[,3]
Isso resultará no vetor dos quadrados médios! a 1ª coluna é dos graus de liberdade, a segunda das somas de quadrados, 3ª dos quadrados médios ...
PS: Mudamos de lista, mude p/ a "r-br", esta em breve será desativada!
att, FH
2011/4/23 greice.laureano <greice.laureano@gmail.com>
Bom dia,
Gostaria de ajuda para extrair os quadrados médios de uma Anova. O banco de dados está abaixo assim como o código para fazer a Anova. Muito Obrigado;
Greice.
Dados:
sujeito aferidor tempo valor 1 1 1 -2 1 1 2 0 1 1 3 1 1 2 1 2 1 2 2 1 1 2 3 1 2 1 1 16 2 1 2 16 2 1 3 15 2 2 1 12 2 2 2 14 2 2 3 13 3 1 1 5 3 1 2 6 3 1 3 6 3 2 1 4 3 2 2 4 3 2 3 4 4 1 1 11 4 1 2 10 4 1 3 10 4 2 1 9 4 2 2 7 4 2 3 8 5 1 1 7 5 1 2 8 5 1 3 6 5 2 1 5 5 2 2 6 5 2 3 6 6 1 1 -7 6 1 2 -8 6 1 3 -8 6 2 1 -9 6 2 2 -10 6 2 3 -9 7 1 1 18 7 1 2 19 7 1 3 19 7 2 1 17 7 2 2 17 7 2 3 17 8 1 1 4 8 1 2 5 8 1 3 5 8 2 1 5 8 2 2 5 8 2 3 5 9 1 1 0 9 1 2 -3 9 1 3 -2 9 2 1 -7 9 2 2 -6 9 2 3 -5 10 1 1 0 10 1 2 0 10 1 3 -2 10 2 1 1 10 2 2 2 10 2 3 1 11 1 1 -3 11 1 2 -2 11 1 3 -2 11 2 1 -4 11 2 2 -3 11 2 3 -3 12 1 1 3 12 1 2 -1 12 1 3 1 12 2 1 -1 12 2 2 -2 12 2 3 1 13 1 1 7 13 1 2 9 13 1 3 9 13 2 1 4 13 2 2 4 13 2 3 2 14 1 1 -6 14 1 2 -7 14 1 3 -6 14 2 1 -8 14 2 2 -10 14 2 3 -9 15 1 1 1 15 1 2 1 15 1 3 0 15 2 1 -2 15 2 2 -2 15 2 3 -3 16 1 1 -13 16 1 2 -14 16 1 3 -14 16 2 1 -12 16 2 2 -12 16 2 3 -12 17 1 1 2 17 1 2 1 17 1 3 0 17 2 1 -1 17 2 2 0 17 2 3 0 18 1 1 4 18 1 2 4 18 1 3 3 18 2 1 7 18 2 2 6 18 2 3 4 19 1 1 -10 19 1 2 -9 19 1 3 -10 19 2 1 -10 19 2 2 -11 19 2 3 -10 20 1 1 8 20 1 2 9 20 1 3 8 20 2 1 2 20 2 2 8 20 2 3 8 21 1 1 7 21 1 2 6 21 1 3 7 21 2 1 8 21 2 2 7 21 2 3 7 22 1 1 -3 22 1 2 -2 22 1 3 -4 22 2 1 -5 22 2 2 -5 22 2 3 -5 23 1 1 -5 23 1 2 -5 23 1 3 -7 23 2 1 -6 23 2 2 -8 23 2 3 -7 24 1 1 5 24 1 2 5 24 1 3 5 24 2 1 3 24 2 2 4 24 2 3 4 25 1 1 0 25 1 2 -1 25 1 3 -1 25 2 1 -4 25 2 2 -3 25 2 3 -4 26 1 1 7 26 1 2 6 26 1 3 6 26 2 1 4 26 2 2 4 26 2 3 4 27 1 1 -8 27 1 2 -8 27 1 3 -8 27 2 1 -10 27 2 2 -11 27 2 3 -10 28 1 1 1 28 1 2 1 28 1 3 2 28 2 1 1 28 2 2 -1 28 2 3 0 29 1 1 -3 29 1 2 -3 29 1 3 -3 29 2 1 -5 29 2 2 -4 29 2 3 -5
sj=as.factor(dados$sujeito) af=as.factor(dados$aferidor) tempo=as.factor(dados$tempo)
fit <- aov(dados$valor ~ sj*af + sj:tempo:af + Error(sj/af)) summary(fit)
__._,_.___
<fernandohtoledo@gmail.com?subject=Res%3A%20Re%3A%20%5BR_STAT%5D%20Sele%E7%E3o%20de%20Quadrados%20M%E9dios>| através de email<R_STAT@yahoogrupos.com.br?subject=Res%3A%20Re%3A%20%5BR_STAT%5D%20Sele%E7%E3o%20de%20Quadrados%20M%E9dios>| Responder através da web<http://br.groups.yahoo.com/group/R_STAT/post;_ylc=X3oDMTJyNWprYWVzBF9TAzk3NDkwNDM3BGdycElkAzExOTIzMjc1BGdycHNwSWQDMjEzNzExMTYwNQRtc2dJZAMxNjA0MARzZWMDZnRyBHNsawNycGx5BHN0aW1lAzEzMDM1Njc5NDA-?act=reply&messageNum=16040>| Adicionar um novo tópico<http://br.groups.yahoo.com/group/R_STAT/post;_ylc=X3oDMTJmNDRncmJ1BF9TAzk3NDkwNDM3BGdycElkAzExOTIzMjc1BGdycHNwSWQDMjEzNzExMTYwNQRzZWMDZnRyBHNsawNudHBjBHN0aW1lAzEzMDM1Njc5NDA-> Mensagens neste tópico<http://br.groups.yahoo.com/group/R_STAT/message/16039;_ylc=X3oDMTM3N2Q4Y2ViBF9TAzk3NDkwNDM3BGdycElkAzExOTIzMjc1BGdycHNwSWQDMjEzNzExMTYwNQRtc2dJZAMxNjA0MARzZWMDZnRyBHNsawN2dHBjBHN0aW1lAzEzMDM1Njc5NDAEdHBjSWQDMTYwMzk->( 2) Atividade nos últimos dias:
Visite seu Grupo<http://br.groups.yahoo.com/group/R_STAT;_ylc=X3oDMTJmbDdmMWM3BF9TAzk3NDkwNDM3BGdycElkAzExOTIzMjc1BGdycHNwSWQDMjEzNzExMTYwNQRzZWMDdnRsBHNsawN2Z2hwBHN0aW1lAzEzMDM1Njc5NDA-> Não deixe de ir atrás da pessoa dos seus sonhos!<http://global.ard.yahoo.com/SIG=15me74fii/M=758712.14532720.14422460.12960164/D=brclubs/S=2137111605:MKP1/Y=BR/EXP=1303575141/L=b473c190-6db3-11e0-860f-2b0ed2bf8b7f/B=f6ybi0wNO6Q-/J=1303567941524925/K=bkkHq4.Jix.KKZcJHGPNcw/A=6351151/R=0/id=mkp1/SIG=135g0plhh/*http://tracking.parperfeito.com.br/ppbanner/bannerTracker?originId=1&identifierId=461809&actionId=1> [image: Yahoo! Grupos]<http://br.groups.yahoo.com/;_ylc=X3oDMTJlc2RmczJyBF9TAzk3NDkwNDM1BGdycElkAzExOTIzMjc1BGdycHNwSWQDMjEzNzExMTYwNQRzZWMDZnRyBHNsawNnZnAEc3RpbWUDMTMwMzU2Nzk0MA--> Trocar para: Só Texto<R_STAT-traditional@yahoogrupos.com.br?subject=Mudar+Formato+de+Envio:+Tradicional>, Resenha Diária<R_STAT-digest@yahoogrupos.com.br?subject=Envio+de+email:+Resenha>• Sair do grupo <R_STAT-unsubscribe@yahoogrupos.com.br?subject=Sair+do+grupo> • Termos de uso <http://br.yahoo.com/info/utos.html> .
__,_._,___
_______________________________________________ R-br mailing list R-br@listas.c3sl.ufpr.br https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
Oi Fernando, Muito Obrigado pela ajuda! Tentei rodar o teu comando mas deu erro aqui... summary(anova)$"Error: sj/af"[[1]]$"Mean Sq" Erro em object[[i]] : objeto de tipo 'closure' não possível dividir em subconjuntos Em ultimo caso posso tentar realizar o teste sem o Erro, no entanto gostaria de conseguir retirar essas estimativas do modelo geral, pois posteriormente precisarei dessas informações de diversos modelos... Att, Greice. 2011/4/23 Fernando Henrique Toledo <fernandohtoledo@gmail.com>
Greice, achei...
Tente esse!
Acho que o que possa dar erro é o que deixei colorido, não sei se tem que ser sj/af ou sj:af!
summary(anova)$"Error: sj/af"[[1]]$"Mean Sq"
Qualquer problema fico a disposição.
att, FH
2011/4/23 Fernando Henrique Toledo <fernandohtoledo@gmail.com>
Agora que vi... Você está usando o argumento Error(), o que altera a saída para uma lista! se você efetivamente está interessada em resgatar os QM pode ajustar o mesmo modelo sem esse Error() e extrair como te expliquei, afinal com o Error() especificado o que muda são os testes F e os p-valor, os QM permanecem os mesmos!
Desculpe pelo confusão!
Vale lembrar que ainda existe um meio de se extrair nesse caso, com o aovlist! Vou ver aqui se acho e te envio em tempo.
att, FH
2011/4/23 Greice Laureano <greice.laureano@gmail.com>
Fernando, (agora mandei para a lista nova :P)
O problema é que o comando aov() não aceita esse procedimento... Se eu faço uma Anova do fit ele da a seguinte mensagem de erro:
anova(fit) Erro em UseMethod("anova") : método não aplicável para 'anova' aplicado a um objeto de classe "c('aovlist', 'listof')"
2011/4/23 Fernando Henrique Toledo <fernandohtoledo@gmail.com>
Greice, bom dia.
Acesse:
anova(fit)[,3]
Isso resultará no vetor dos quadrados médios! a 1ª coluna é dos graus de liberdade, a segunda das somas de quadrados, 3ª dos quadrados médios ...
PS: Mudamos de lista, mude p/ a "r-br", esta em breve será desativada!
att, FH
2011/4/23 greice.laureano <greice.laureano@gmail.com>
Bom dia,
Gostaria de ajuda para extrair os quadrados médios de uma Anova. O banco de dados está abaixo assim como o código para fazer a Anova. Muito Obrigado;
Greice.
Dados:
sujeito aferidor tempo valor 1 1 1 -2 1 1 2 0 1 1 3 1 1 2 1 2 1 2 2 1 1 2 3 1 2 1 1 16 2 1 2 16 2 1 3 15 2 2 1 12 2 2 2 14 2 2 3 13 3 1 1 5 3 1 2 6 3 1 3 6 3 2 1 4 3 2 2 4 3 2 3 4 4 1 1 11 4 1 2 10 4 1 3 10 4 2 1 9 4 2 2 7 4 2 3 8 5 1 1 7 5 1 2 8 5 1 3 6 5 2 1 5 5 2 2 6 5 2 3 6 6 1 1 -7 6 1 2 -8 6 1 3 -8 6 2 1 -9 6 2 2 -10 6 2 3 -9 7 1 1 18 7 1 2 19 7 1 3 19 7 2 1 17 7 2 2 17 7 2 3 17 8 1 1 4 8 1 2 5 8 1 3 5 8 2 1 5 8 2 2 5 8 2 3 5 9 1 1 0 9 1 2 -3 9 1 3 -2 9 2 1 -7 9 2 2 -6 9 2 3 -5 10 1 1 0 10 1 2 0 10 1 3 -2 10 2 1 1 10 2 2 2 10 2 3 1 11 1 1 -3 11 1 2 -2 11 1 3 -2 11 2 1 -4 11 2 2 -3 11 2 3 -3 12 1 1 3 12 1 2 -1 12 1 3 1 12 2 1 -1 12 2 2 -2 12 2 3 1 13 1 1 7 13 1 2 9 13 1 3 9 13 2 1 4 13 2 2 4 13 2 3 2 14 1 1 -6 14 1 2 -7 14 1 3 -6 14 2 1 -8 14 2 2 -10 14 2 3 -9 15 1 1 1 15 1 2 1 15 1 3 0 15 2 1 -2 15 2 2 -2 15 2 3 -3 16 1 1 -13 16 1 2 -14 16 1 3 -14 16 2 1 -12 16 2 2 -12 16 2 3 -12 17 1 1 2 17 1 2 1 17 1 3 0 17 2 1 -1 17 2 2 0 17 2 3 0 18 1 1 4 18 1 2 4 18 1 3 3 18 2 1 7 18 2 2 6 18 2 3 4 19 1 1 -10 19 1 2 -9 19 1 3 -10 19 2 1 -10 19 2 2 -11 19 2 3 -10 20 1 1 8 20 1 2 9 20 1 3 8 20 2 1 2 20 2 2 8 20 2 3 8 21 1 1 7 21 1 2 6 21 1 3 7 21 2 1 8 21 2 2 7 21 2 3 7 22 1 1 -3 22 1 2 -2 22 1 3 -4 22 2 1 -5 22 2 2 -5 22 2 3 -5 23 1 1 -5 23 1 2 -5 23 1 3 -7 23 2 1 -6 23 2 2 -8 23 2 3 -7 24 1 1 5 24 1 2 5 24 1 3 5 24 2 1 3 24 2 2 4 24 2 3 4 25 1 1 0 25 1 2 -1 25 1 3 -1 25 2 1 -4 25 2 2 -3 25 2 3 -4 26 1 1 7 26 1 2 6 26 1 3 6 26 2 1 4 26 2 2 4 26 2 3 4 27 1 1 -8 27 1 2 -8 27 1 3 -8 27 2 1 -10 27 2 2 -11 27 2 3 -10 28 1 1 1 28 1 2 1 28 1 3 2 28 2 1 1 28 2 2 -1 28 2 3 0 29 1 1 -3 29 1 2 -3 29 1 3 -3 29 2 1 -5 29 2 2 -4 29 2 3 -5
sj=as.factor(dados$sujeito) af=as.factor(dados$aferidor) tempo=as.factor(dados$tempo)
fit <- aov(dados$valor ~ sj*af + sj:tempo:af + Error(sj/af)) summary(fit)
__._,_.___
<fernandohtoledo@gmail.com?subject=Res%3A%20Re%3A%20%5BR_STAT%5D%20Sele%E7%E3o%20de%20Quadrados%20M%E9dios>| através de email<R_STAT@yahoogrupos.com.br?subject=Res%3A%20Re%3A%20%5BR_STAT%5D%20Sele%E7%E3o%20de%20Quadrados%20M%E9dios>| Responder através da web<http://br.groups.yahoo.com/group/R_STAT/post;_ylc=X3oDMTJyNWprYWVzBF9TAzk3NDkwNDM3BGdycElkAzExOTIzMjc1BGdycHNwSWQDMjEzNzExMTYwNQRtc2dJZAMxNjA0MARzZWMDZnRyBHNsawNycGx5BHN0aW1lAzEzMDM1Njc5NDA-?act=reply&messageNum=16040>| Adicionar um novo tópico<http://br.groups.yahoo.com/group/R_STAT/post;_ylc=X3oDMTJmNDRncmJ1BF9TAzk3NDkwNDM3BGdycElkAzExOTIzMjc1BGdycHNwSWQDMjEzNzExMTYwNQRzZWMDZnRyBHNsawNudHBjBHN0aW1lAzEzMDM1Njc5NDA-> Mensagens neste tópico<http://br.groups.yahoo.com/group/R_STAT/message/16039;_ylc=X3oDMTM3N2Q4Y2ViBF9TAzk3NDkwNDM3BGdycElkAzExOTIzMjc1BGdycHNwSWQDMjEzNzExMTYwNQRtc2dJZAMxNjA0MARzZWMDZnRyBHNsawN2dHBjBHN0aW1lAzEzMDM1Njc5NDAEdHBjSWQDMTYwMzk->( 2) Atividade nos últimos dias:
Visite seu Grupo<http://br.groups.yahoo.com/group/R_STAT;_ylc=X3oDMTJmbDdmMWM3BF9TAzk3NDkwNDM3BGdycElkAzExOTIzMjc1BGdycHNwSWQDMjEzNzExMTYwNQRzZWMDdnRsBHNsawN2Z2hwBHN0aW1lAzEzMDM1Njc5NDA-> Não deixe de ir atrás da pessoa dos seus sonhos!<http://global.ard.yahoo.com/SIG=15me74fii/M=758712.14532720.14422460.12960164/D=brclubs/S=2137111605:MKP1/Y=BR/EXP=1303575141/L=b473c190-6db3-11e0-860f-2b0ed2bf8b7f/B=f6ybi0wNO6Q-/J=1303567941524925/K=bkkHq4.Jix.KKZcJHGPNcw/A=6351151/R=0/id=mkp1/SIG=135g0plhh/*http://tracking.parperfeito.com.br/ppbanner/bannerTracker?originId=1&identifierId=461809&actionId=1> [image: Yahoo! Grupos]<http://br.groups.yahoo.com/;_ylc=X3oDMTJlc2RmczJyBF9TAzk3NDkwNDM1BGdycElkAzExOTIzMjc1BGdycHNwSWQDMjEzNzExMTYwNQRzZWMDZnRyBHNsawNnZnAEc3RpbWUDMTMwMzU2Nzk0MA--> Trocar para: Só Texto<R_STAT-traditional@yahoogrupos.com.br?subject=Mudar+Formato+de+Envio:+Tradicional>, Resenha Diária<R_STAT-digest@yahoogrupos.com.br?subject=Envio+de+email:+Resenha>• Sair do grupo <R_STAT-unsubscribe@yahoogrupos.com.br?subject=Sair+do+grupo>• Termos de uso <http://br.yahoo.com/info/utos.html> .
__,_._,___
_______________________________________________ R-br mailing list R-br@listas.c3sl.ufpr.br https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
_______________________________________________ R-br mailing list R-br@listas.c3sl.ufpr.br https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
Greice, Segue o exemplo reproduzível de como acessar de diversas formas os valores do quadro de análise de variância para um modelo com Error() na fórmula. da <- expand.grid(bloco=gl(4,1,la="b"), parce=gl(3,1,la="p"), subpa=gl(5,1,la="s")) da$y <- rnorm(nrow(da)) m0 <- aov(y~bloco+parce*subpa+Error(bloco:parce), data=da) summary(m0) str(m0) # apresenta a estrutura do objeto, você aprende como acessá-lo summ <- summary(m0) str(summ) # o que é? o que tem? como acesso? names(summ) summ[1] # porção parcela da anova summ[2] # porção subparcela da anova summ$"Error: bloco:parce" # acesso pelo nome do slot summ$"Error: Within" # acesso pelo nome do slot summ[[1]][[1]] # acesso pelo número do slot summ[[2]][[1]] # acesso pelo número do slot summ$"Error: bloco:parce"[[1]][,"Df"] # acesso por nome summ$"Error: bloco:parce"[[1]]["bloco","Df"] # acesso por nome summ[[1]][[1]][,"Sum Sq"] # acesso por nome summ[[2]][[1]]["Residuals",] # acesso por nome summ[[2]][[1]]["Residuals","Mean Sq"] # acesso por nome summ$"Error: bloco:parce"[[1]][,1] # acesso por nome summ$"Error: bloco:parce"[[1]][1,1] # acesso por nome À disposição. Walmes. ========================================================================== Walmes Marques Zeviani LEG (Laboratório de Estatística e Geoinformação, 25.450418 S, 49.231759 W) Departamento de Estatística - Universidade Federal do Paraná fone: (+55) 41 3361 3573 VoIP: (3361 3600) 1053 1173 e-mail: walmes@ufpr.br twitter: @walmeszeviani homepage: http://www.leg.ufpr.br/~walmes linux user number: 531218 ==========================================================================
Walmes, Muito obrigado! Funcionou perfeitamente =D Att, Greice. 2011/4/23 Walmes Zeviani <walmeszeviani@gmail.com>
Greice,
Segue o exemplo reproduzível de como acessar de diversas formas os valores do quadro de análise de variância para um modelo com Error() na fórmula.
da <- expand.grid(bloco=gl(4,1,la="b"), parce=gl(3,1,la="p"), subpa=gl(5,1,la="s")) da$y <- rnorm(nrow(da))
m0 <- aov(y~bloco+parce*subpa+Error(bloco:parce), data=da) summary(m0)
str(m0) # apresenta a estrutura do objeto, você aprende como acessá-lo summ <- summary(m0) str(summ) # o que é? o que tem? como acesso? names(summ)
summ[1] # porção parcela da anova summ[2] # porção subparcela da anova
summ$"Error: bloco:parce" # acesso pelo nome do slot summ$"Error: Within" # acesso pelo nome do slot
summ[[1]][[1]] # acesso pelo número do slot summ[[2]][[1]] # acesso pelo número do slot
summ$"Error: bloco:parce"[[1]][,"Df"] # acesso por nome summ$"Error: bloco:parce"[[1]]["bloco","Df"] # acesso por nome
summ[[1]][[1]][,"Sum Sq"] # acesso por nome summ[[2]][[1]]["Residuals",] # acesso por nome summ[[2]][[1]]["Residuals","Mean Sq"] # acesso por nome
summ$"Error: bloco:parce"[[1]][,1] # acesso por nome summ$"Error: bloco:parce"[[1]][1,1] # acesso por nome
À disposição. Walmes.
========================================================================== Walmes Marques Zeviani LEG (Laboratório de Estatística e Geoinformação, 25.450418 S, 49.231759 W) Departamento de Estatística - Universidade Federal do Paraná fone: (+55) 41 3361 3573 VoIP: (3361 3600) 1053 1173 e-mail: walmes@ufpr.br twitter: @walmeszeviani homepage: http://www.leg.ufpr.br/~walmes linux user number: 531218 ==========================================================================
_______________________________________________ R-br mailing list R-br@listas.c3sl.ufpr.br https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
participantes (3)
-
Fernando Henrique Toledo -
Greice Laureano -
Walmes Zeviani