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)
### 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..
[]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.
É 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 ==========================================================================