[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