Informação sobre efeito marginal para Modelos Logit Ordenado Generalizado

Olá a todos, Estou atualmente trabalhando com modelos logit ordenados generalizados, estou trabalhando no R atráves do pacote VGAM (vector generalized linear and additive model) na função vglm(). Bom, eu consigo estimar o modelo com ou sem hipótese de regressão paralela (proportional odds model). Meu objetivo é estimar o efeito marginal da probabilidade, para isso, utilizo a função margeff() contida no pacote. Porém, a função executa o efeito para cada indivíduo, sem erro-padrão e estatística de teste. Então, gostaria de saber se existe outro modo (ou pacote) de estimar tais efeitos marginais de forma completa. P.S: já utilizei as funções clm (pacote ordinal) e polr() (pacote MASS) e não obtive sucesso. Agradeço a ajuda. Abel Brasil Estatístico - UFC

É 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@ufpr.br skype: walmeszeviani twitter: @walmeszeviani homepage: http://www.leg.ufpr.br/~walmes linux user number: 531218 ==========================================================================

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@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@ufpr.br skype: walmeszeviani twitter: @walmeszeviani homepage: http://www.leg.ufpr.br/~walmes linux user number: 531218 ==========================================================================
_______________________________________________ R-br mailing list R-br@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.

Caro Abel, CMR significa Código Mínimo Reproduzível. Veja o seguinte link: http://r-br.2285057.n4.nabble.com/R-br-Produzindo-um-Codigo-Minimo-Reproduzi... Att, Eliardo. Em 15 de outubro de 2012 14:50, Abel Brasil Ramos da Silva < abelbrasil88@gmail.com> 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@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@ufpr.br skype: walmeszeviani twitter: @walmeszeviani homepage: http://www.leg.ufpr.br/~walmes linux user number: 531218 ==========================================================================
_______________________________________________ R-br mailing list R-br@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@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.

Caro Eliardo, Vou preparar um documento com os comandos na qual estou utilizando e enviarei brevemente. Agredeço pela ajuda... Att. Em 15 de outubro de 2012 15:06, Eliardo Costa <eliardocosta@gmail.com>escreveu:
Caro Abel,
CMR significa Código Mínimo Reproduzível. Veja o seguinte link:
http://r-br.2285057.n4.nabble.com/R-br-Produzindo-um-Codigo-Minimo-Reproduzi...
Att, Eliardo.
Em 15 de outubro de 2012 14:50, Abel Brasil Ramos da Silva < abelbrasil88@gmail.com> 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@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@ufpr.br skype: walmeszeviani twitter: @walmeszeviani homepage: http://www.leg.ufpr.br/~walmes linux user number: 531218
==========================================================================
_______________________________________________ R-br mailing list R-br@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@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@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.

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@gmail.com <mailto:walmeszeviani@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 <tel:%28%2B55%29%2041%203361%203573> VoIP: (3361 3600) 1053 1173 e-mail: walmes@ufpr.br <mailto:walmes@ufpr.br> skype: walmeszeviani twitter: @walmeszeviani homepage: http://www.leg.ufpr.br/~walmes <http://www.leg.ufpr.br/%7Ewalmes> linux user number: 531218 ==========================================================================
_______________________________________________ R-br mailing list R-br@listas.c3sl.ufpr.br <mailto:R-br@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@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.

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@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@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@ufpr.br skype: walmeszeviani twitter: @walmeszeviani homepage: http://www.leg.ufpr.br/~walmes linux user number: 531218 ==========================================================================
_______________________________________________ R-br mailing list R-br@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@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@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.
participantes (4)
-
Abel Brasil Ramos da Silva
-
Eliardo Costa
-
Leonard de Assis
-
Walmes Zeviani