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 <[email protected]> > 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 <[email protected]> > >> 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 <[email protected]> >> >>> 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 <[email protected]> >>> >>>> >>>> >>>> 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 <[email protected]> >>>> >>>> >>>>> >>>>> 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) >>>>> >>>>> >>>> __._,_.___ >>>> >>>> <[email protected]?subject=Res%3A%20Re%3A%20%5BR_STAT%5D%20Sele%E7%E3o%20de%20Quadrados%20M%E9dios>| >>>> através >>>> de >>>> email<[email protected]?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: >>>> >>>> - Novos >>>> usuários<http://br.groups.yahoo.com/group/R_STAT/members;_ylc=X3oDMTJnNmk2bWJoBF9TAzk3NDkwNDM3BGdycElkAzExOTIzMjc1BGdycHNwSWQDMjEzNzExMTYwNQRzZWMDdnRsBHNsawN2bWJycwRzdGltZQMxMzAzNTY3OTQw?o=6> >>>> 2 >>>> >>>> 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<[email protected]?subject=Mudar+Formato+de+Envio:+Tradicional>, >>>> Resenha >>>> Diária<[email protected]?subject=Envio+de+email:+Resenha>• >>>> Sair >>>> do grupo <[email protected]?subject=Sair+do+grupo>• >>>> Termos >>>> de uso <http://br.yahoo.com/info/utos.html> >>>> . >>>> >>>> __,_._,___ >>>> >>> >>> >>> _______________________________________________ >>> R-br mailing list >>> [email protected] >>> https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br >>> >>> >> > > _______________________________________________ > R-br mailing list > [email protected] > https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br > >
_______________________________________________ R-br mailing list [email protected] https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
