<div>Olá a todos,</div>
<div> </div>
<div>Como solicitado fiz uma simulação dos dados e reproduzi o meu problema:</div>
<div> </div>
<div>####### R GRUPO ######<br>require("VGAM")<br>set.seed(279)</div>
<div>##### Simulação dos dados #### <br>y=sample(1:5,100,rep=T) #Exemplo: 1:Péssimo 2:Ruim 3:Razoável 4:Bom 5:Excelente<br>x3=rnorm(100,50,22)<br>x4=rpois(100, lambda=24)<br>#-----------------------------------------------------------------------------------</div>
<div>#### Estimando o modelo utilizando a função vglm() ###<br>y=ordered(y,levels=c(1,2,3,4,5))<br>mod=vglm(y~x3+x4,family = cumulative(parallel=T,reverse=F))<br>summary(mod) ## aqui são as estimativas dos parametros em logit(P(Y<j)) e não em probabilidade<br>
#---------------------------------------------------------</div>
<div>##### Extraindo os efeitos marginais (probabilidade) através da função margeff() ##<br>efeito.marg=margeff(mod)<br>#--------------------------------------------------------------</div>
<div>##### Aplicando a média para a categoria por variável ####<br>efeitomarg=matrix(c(rowMeans(efeito.marg['(Intercept)',,]),<br> rowMeans(efeito.marg['x3',,]),<br> rowMeans(efeito.marg['x4',,])<br>
),ncol=5)</div>
<div>rownames(efeitomarg) = c("Intercepto","Var1","Var2")<br>colnames(efeitomarg)=c("Péssimo","Ruim","Razoável","Bom","Excelente")<br>#-----------------------------------------------------------------------</div>
<div>### efeito marginal da probabilidade### <br>round(efeitomarg,7) ## porém não tenho informação sobre o erro-padrão e etc<br></div>
<div>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. </div>
<div> </div>
<div>Agredeço pela ajuda de todos,</div>
<div> </div>
<div>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..</div>
<div> </div>
<div>Abraços<br></div>
<div class="gmail_quote">Em 15 de outubro de 2012 20:44, Leonard de Assis <span dir="ltr"><<a href="mailto:assis.leonard@gmail.com" target="_blank">assis.leonard@gmail.com</a>></span> escreveu:<br>
<blockquote style="BORDER-LEFT:#ccc 1px solid;MARGIN:0px 0px 0px 0.8ex;PADDING-LEFT:1ex" class="gmail_quote">
<div bgcolor="#FFFFFF" text="#000000">
<div>CMR = Código mínimo Reproduzivel<br><br>Mas leva um 8 pela criatividade, hehehe<br><pre cols="72">[]s
Leonard de Assis
assis <dot> leonard <at> gmail <dot> com
</pre>Em 15/10/2012 14:50, Abel Brasil Ramos da Silva escreveu:<br></div>
<div>
<div class="h5">
<blockquote type="cite">
<div>Prezado Walmes,</div>
<div> </div>
<div>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. </div>
<div> </div>
<div>Agradeço pelo ajuda...</div>
<div> </div>
<div>Att.<br><br><br><br></div>
<div class="gmail_quote">Em 15 de outubro de 2012 10:31, Walmes Zeviani <span dir="ltr"><<a href="mailto:walmeszeviani@gmail.com" target="_blank">walmeszeviani@gmail.com</a>></span> escreveu:<br>
<blockquote style="BORDER-LEFT:#ccc 1px solid;MARGIN:0px 0px 0px 0.8ex;PADDING-LEFT:1ex" class="gmail_quote"><font face="trebuchet ms,sans-serif">É um modelo de efeito aleatório? É uma estimação do efeito de indivíduo ou uma predição? Que tal um CMR?<br>
<br>À disposição.<br>Walmes.<br><br clear="all"></font><span style="FONT-FAMILY:trebuchet ms,sans-serif">==========================================================================</span><br style="FONT-FAMILY:trebuchet ms,sans-serif">
<span style="FONT-FAMILY:trebuchet ms,sans-serif">Walmes Marques Zeviani</span><br style="FONT-FAMILY:trebuchet ms,sans-serif"><span style="FONT-FAMILY:trebuchet ms,sans-serif">LEG (Laboratório de Estatística e Geoinformação, 25.450418 S, 49.231759 W)</span><br style="FONT-FAMILY:trebuchet ms,sans-serif">
<span style="FONT-FAMILY:trebuchet ms,sans-serif">Departamento de Estatística - Universidade Federal do Paraná</span><br style="FONT-FAMILY:trebuchet ms,sans-serif"><span style="FONT-FAMILY:trebuchet ms,sans-serif">fone: <a href="tel:%28%2B55%29%2041%203361%203573" target="_blank" value="+554133613573">(+55) 41 3361 3573</a></span><br style="FONT-FAMILY:trebuchet ms,sans-serif">
<span style="FONT-FAMILY:trebuchet ms,sans-serif">VoIP: (3361 3600) 1053 1173</span><br style="FONT-FAMILY:trebuchet ms,sans-serif"><span style="FONT-FAMILY:trebuchet ms,sans-serif">e-mail: <a href="mailto:walmes@ufpr.br" target="_blank">walmes@ufpr.br</a><br>
skype: walmeszeviani<br style="FONT-FAMILY:trebuchet ms,sans-serif"></span><span style="FONT-FAMILY:trebuchet ms,sans-serif">twitter: @walmeszeviani</span><br style="FONT-FAMILY:trebuchet ms,sans-serif"><span style="FONT-FAMILY:trebuchet ms,sans-serif">homepage: <a href="http://www.leg.ufpr.br/~walmes" target="_blank">http://www.leg.ufpr.br/~walmes</a></span><br style="FONT-FAMILY:trebuchet ms,sans-serif">
<span style="FONT-FAMILY:trebuchet ms,sans-serif">linux user number: 531218</span><br style="FONT-FAMILY:trebuchet ms,sans-serif"><span style="FONT-FAMILY:trebuchet ms,sans-serif">==========================================================================</span><br>
<br>_______________________________________________<br>R-br mailing list<br><a href="mailto:R-br@listas.c3sl.ufpr.br" target="_blank">R-br@listas.c3sl.ufpr.br</a><br><a href="https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br" target="_blank">https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br</a><br>
Leia o guia de postagem (<a href="http://www.leg.ufpr.br/r-br-guia" target="_blank">http://www.leg.ufpr.br/r-br-guia</a>) e forneça código mínimo reproduzível.<br></blockquote></div><br><br>
<fieldset></fieldset> <br><pre>_______________________________________________
R-br mailing list
<a href="mailto:R-br@listas.c3sl.ufpr.br" target="_blank">R-br@listas.c3sl.ufpr.br</a>
<a href="https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br" target="_blank">https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br</a>
Leia o guia de postagem (<a href="http://www.leg.ufpr.br/r-br-guia" target="_blank">http://www.leg.ufpr.br/r-br-guia</a>) e forneça código mínimo reproduzível.</pre></blockquote><br></div></div></div><br>_______________________________________________<br>
R-br mailing list<br><a href="mailto:R-br@listas.c3sl.ufpr.br">R-br@listas.c3sl.ufpr.br</a><br><a href="https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br" target="_blank">https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br</a><br>
Leia o guia de postagem (<a href="http://www.leg.ufpr.br/r-br-guia" target="_blank">http://www.leg.ufpr.br/r-br-guia</a>) e forneça código mínimo reproduzível.<br></blockquote></div><br>