[R-br] Informação sobre efeito marginal para Modelos Logit Ordenado Generalizado

[Abel Brasil Ramos da Silva]

Olá a todos,

Como solicitado fiz uma simulação dos dados e reproduzi o meu problema:

####### R GRUPO ######
require("VGAM")
set.seed(279)
##### Simulação dos dados ####
y=sample(1:5,100,rep=T) #Exemplo: 1:Péssimo 2:Ruim 3:Razoável 4:Bom
5:Excelente
x3=rnorm(100,50,22)
x4=rpois(100, lambda=24)
#-----------------------------------------------------------------------------------
#### Estimando o modelo utilizando a função vglm() ###
y=ordered(y,levels=c(1,2,3,4,5))
mod=vglm(y~x3+x4,family = cumulative(parallel=T,reverse=F))
summary(mod)  ## aqui são as estimativas dos parametros em logit(P(Y<j)) e
não em probabilidade
#---------------------------------------------------------
##### Extraindo os efeitos marginais (probabilidade) através da função
margeff() ##
efeito.marg=margeff(mod)
#--------------------------------------------------------------
##### Aplicando a média para a categoria por variável ####
efeitomarg=matrix(c(rowMeans(efeito.marg['(Intercept)',,]),
                    rowMeans(efeito.marg['x3',,]),
                    rowMeans(efeito.marg['x4',,])
),ncol=5)
rownames(efeitomarg) = c("Intercepto","Var1","Var2")
colnames(efeitomarg)=c("Péssimo","Ruim","Razoável","Bom","Excelente")
#-----------------------------------------------------------------------
### efeito marginal da probabilidade###
round(efeitomarg,7)  ## porém não tenho informação sobre o erro-padrão e etc
Então, o procedimento que realizo é esse... Como já salientei, tenho apenas
a informação média do efeito marginal da probabilidade... Gostaria de saber
se existe um modo ou pacote auxiliar que possa extrair as informações mais
completa sobre isso.

Agredeço pela ajuda de todos,

P.S: Caro, Leonardo, rsrsrsrs, valeu pela nota... olha eu sempre observo as
msg neste grupo e nunca vi a sigla CMR. Eu fui atrás de saber sobre a
sigla, rsrsrsrs, no google, como o que estou executando é algo da análise
de dados categorizados e o Walmes mencionou sobre efeito aleatório... ai
pensei que fosse isso Categorical Model Random... (paguei mico em rede
nacional, mas tudo bem!) Enfim, na proxima vez eu ajudo os futuros colegas
com esse problema do CMR... rsrsrsrs..

Abraços
Em 15 de outubro de 2012 20:44, Leonard de Assis
<assis.leonard em gmail.com>escreveu:

>  CMR = Código mínimo Reproduzivel
>
> Mas leva um 8 pela criatividade, hehehe
>
> []s
> Leonard de Assis
> assis <dot> leonard <at> gmail <dot> com
>
> Em 15/10/2012 14:50, Abel Brasil Ramos da Silva escreveu:
>
> Prezado Walmes,
>
> Inicialmente, não pretendo utilizar o modelo com efeito aleatório... CMR
> seria categorical model random? Na verdade estou fazendo uma estimação do
> modelo, a função margeff() extrai o efeito marginal para cada individuo em
> particular (o que não quero), então eu utilizo o recurso do rowMeans() para
> extrair a média do efeito individual para cada variável. Então, gostaria de
> obter uma função ou pacote auxiliar que possa gerar as estimativas do
> efeito marginal para cada variável bem como seus erro padrões... Estive
> lendo que o metodo delta é usada para calcular erros padrões para modelos
> ordenados (não tenho certeza sobre isso!), também não sei onde faço isso no
> R.
>
> Agradeço pelo ajuda...
>
> Att.
>
>
>
> Em 15 de outubro de 2012 10:31, Walmes Zeviani <walmeszeviani em gmail.com>escreveu:
>
>> É um modelo de efeito aleatório? É uma estimação do efeito de indivíduo
>> ou uma predição? Que tal um CMR?
>>
>> À 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 em ufpr.br
>> skype: walmeszeviani
>> twitter: @walmeszeviani
>> homepage: http://www.leg.ufpr.br/~walmes
>> linux user number: 531218
>> ==========================================================================
>>
>> _______________________________________________
>> 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 listR-br em listas.c3sl.ufpr.brhttps://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/2ce2cef7/attachment.html>