[R-br] [R_STAT] Seleção de Quadrados Médios

Sábado Abril 23 17:49:59 BRT 2011

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 em 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 em 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 em 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 em 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 em 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 em gmail.com?subject=Res%3A%20Re%3A%20%5BR_STAT%5D%20Sele%E7%E3o%20de%20Quadrados%20M%E9dios>| através
>>>> de email<R_STAT em 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:
>>>>
>>>>    - 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<R_STAT-traditional em yahoogrupos.com.br?subject=Mudar+Formato+de+Envio:+Tradicional>,
>>>> Resenha DiÃ¡ria<R_STAT-digest em yahoogrupos.com.br?subject=Envio+de+email:+Resenha>• Sair
>>>> do grupo <R_STAT-unsubscribe em yahoogrupos.com.br?subject=Sair+do+grupo>• Termos
>>>> de uso <http://br.yahoo.com/info/utos.html>
>>>>    .
>>>>
>>>> __,_._,___
>>>>
>>>
>>>
>>> _______________________________________________
>>> R-br mailing list
>>> R-br em listas.c3sl.ufpr.br
>>> https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
>>>
>>>
>>
>
> _______________________________________________
> R-br mailing list
> R-br em listas.c3sl.ufpr.br
> https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
>
>
-------------- Próxima Parte ----------
Um anexo em HTML foi limpo...
URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20110423/e5d69555/attachment.html>