[R-br] Informação sobre efeito marginal para Modelos Logit Ordenado Generalizado
Fabio Mathias Corrêa
fabio.ufla em yahoo.com.br
Terça Outubro 16 11:58:06 BRT 2012
Abel,
Veja se é isto o que vc quer.
Se vc tem a média do efeito marginal, você pode ter também o desvio-padrão e outras estimativas.
# Média do efeito marginal
apply(efeito.marg['(Intercept)',,],1,mean) # isto é igual ao rowMeans que vc usa.
1 2 3 4 5
0.10965817 0.29555932 0.18275531 0.02128801 -0.60926081
apply(efeito.marg['(Intercept)',,],1,sd) # Desvio-padrão da estimativa.
1 2 3 4 5
0.02903666 0.02472514 0.05574690 0.13221473 0.14065723
Valeu!
Fábio Mathias Corrêa
Universidade Estadual de Santa Cruz
Departamento de Ciências Exatas e da Terra - DCET
Campus Soane Nazaré de Andrade, km 16 Rodovia Ilhéus-Itabuna
CEP 45662-900. Ilhéus-Bahia
Tel.: 73-3680-5076
________________________________
De: Abel Brasil Ramos da Silva <abelbrasil88 em gmail.com>
Para: r-br em listas.c3sl.ufpr.br; Ivan Bezerra Allaman <ivanalaman em yahoo.com.br>
Enviadas: Terça-feira, 16 de Outubro de 2012 11:37
Assunto: Re: [R-br] Informação sobre efeito marginal para Modelos Logit Ordenado Generalizado
Olá Ivan,
Agradeço a ajuda...
Mas, não é isso o que quero... O fitted() da a probabilidade do indivíduo pertencer a categoria, somente... isso já é uma predição do modelo...
Note, vou tentar exemplicar com os comandos de forma mais simples...
O comando summary(mod) mostra as estimativas dos parametros do modelo, bem como as estimativas dos threshold... Porém, a interpretação das estimativas é sobre logit(P(Y<j)) = a + b*X3 + c*X4 + erro, ou seja, a interpretação da estimativas dos parametros é em cima deste modelo (a interpretação no odds). O efeito marginal de b é na mudança no logit(P(Y<j)).
Então, o efeito marginal na probabilidade, cuja modelo vem da distribuição acumulada da logistica, não é extraido no summary(). Para isso, existe esta função margeff() que faz isso. Porém, não tenho informação do erro-padrão e etc, como vem no summary().
Att.
Em 16 de outubro de 2012 11:10, Ivan Bezerra Allaman <ivanalaman em yahoo.com.br> escreveu:
fitted(mod) e terás o que quer!
>
>
>(s,f,p)
>Allaman
>
>
>
>
>\begin{signature}
><<>>=
>Prof. Dr. Ivan Bezerra Allaman
>Universidade Estadual de Santa Cruz
>Departamento de Ciências Exatas e Tecnológicas
>Ilhéus/BA - Brasil
>Fone: +55 73 3680-5596
>E-mail: ivanalaman em yahoo.com.br/ivanalaman em gmail.com
>@
>\end{signature}
>_______________________________________________
>R-br mailing list
>R-br em listas.c3sl.ufpr.br
>https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
>Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça código mínimo reproduzível.
>
_______________________________________________
R-br mailing list
R-br em listas.c3sl.ufpr.br
https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça código mínimo reproduzível.
-------------- Próxima Parte ----------
Um anexo em HTML foi limpo...
URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20121016/b1cfb046/attachment.html>
Mais detalhes sobre a lista de discussão R-br