<div dir="ltr">Obrigado Rodrigo.<div>A função bild apresenta um problema ao inserir efeito aleatório no modelo,. Segundo o help, faz-se necessário modificar limites de integração para verossimilhança e gradiente em cada banco de dados. </div>
<div>Para um número de simulações pequeno, ele solta a seguinte mensagem</div><div><div>Mensagens de aviso:</div><div>In sqrt(prob * (1 - prob)) : NaNs produzidos</div></div><div>mas emite os resultados.</div><div>Quando uso um número de simulações é muito grande(que é o que preciso), sai seguinte mensagem</div>
<div><div>Erro em solve.default(Info) : </div><div> sistema é computacionalmente singular: condição recíproca número = 3.18165e-19</div></div><div>e os resultado não são emitidos.</div><div>Isto porque eu teria que modificar os limites para cada uma das respostas geradas. E isto é inviável.</div>
<div>A ideia do "replicate" dentro do "for" foi o que surgiu como alternativa para que a função geral() (que extrai os coeficientes) seja rodada uma vez e cada resultado seja armazenado em um objeto para que, por fim, seja obtida uma média.</div>
<div>Entendo que esta opção não é a forma prática, mas foi o que eu pensei como alternativa para fugir deste problema. Qualquer outra sugestão será muito bem vinda.<br></div><div>Grato,</div><div>Maurício</div><div class="gmail_extra">
<br><br><div class="gmail_quote">Em 26 de outubro de 2013 10:08, Rodrigo Coster <span dir="ltr"><<a href="mailto:rcoster@gmail.com" target="_blank">rcoster@gmail.com</a>></span> escreveu:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
<div dir="ltr">replicate() dentro do for não faz sentido, ou tu usa ou outro (nesse caso). Note que tu ta mandando repetir uma vez (o 1o parametro do replicate) para cada iteração do for, e sobrescrevendo o resultado...<div>
<br></div><div><br></div><div>Tente (não testei, mas deve funcionar):</div><div><span style="font-family:arial,sans-serif;font-size:13px">repeticoes <- replicate(ns, geral())</span><br></div><div><span style="font-family:arial,sans-serif;font-size:13px">colMeans(repeticoes)</span></div>
</div><div class="HOEnZb"><div class="h5"><div class="gmail_extra"><br><br><div class="gmail_quote">2013/10/26 Maurício Lordêlo <span dir="ltr"><<a href="mailto:mslordelo@gmail.com" target="_blank">mslordelo@gmail.com</a>></span><br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
<div dir="ltr">Obrigado pela atenção Manoel. Quanto a instalação do pacote, o mesmo foi desenvolvido para versões >= 2.10.0.<div>Fiz o que recomendou mas </div><div><div>Erro em colMeans(repeticoes, na.rm = TRUE) : </div>
<div> 'x' deve ser um array de pelo menos duas dimensões</div></div><div><br></div><div>A saída de "repeticoes" é mostrada assim:</div><div><div>[[1]]</div><div>[1] 0.11926720 -0.06925387 -0.49590355</div>
<div><br></div><div>[[1]]</div><div>[1] -0.1505574 -0.2059867 0.6151096</div><div><br></div><div>[[1]]</div><div>[1] -0.2860369 -0.1084762 -0.1997825</div><div><br></div><div>[[1]]</div><div>[1] 0.4159422 -0.9583708 -0.4319101</div>
<div><br></div><div>[[1]]</div><div>[1] 0.3784151 -0.4214649 -0.1431581</div></div><div><br></div><div><br></div><div>Fiz </div><div>as.array(repeticoes)<br></div><div>e ele me retorna apenas a última linha do objeto "repeticoes".</div>
<div><div>[[1]]</div><div>[1] 0.3784151 -0.4214649 -0.1431581</div></div><div>Tentei também</div><div><div>coefs=do.call(rbind, repeticoes)</div><div>media_coefs=apply(coefs, 2, mean)</div><div>media_coefs</div></div><div>
[1] 0.3784151 -0.4214649 -0.1431581<br></div><div>Que também retorna a última linha.</div><div>Maurício</div><div><br></div></div><div class="gmail_extra"><br><br><div class="gmail_quote">Em 26 de outubro de 2013 08:12, Manoel Galdino <span dir="ltr"><<a href="mailto:mcz.fea@gmail.com" target="_blank">mcz.fea@gmail.com</a>></span> escreveu:<div>
<div><br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr"><div><div>Eu não consegui instalar a biblioteca build no R 3.0<br><br></div>De todo modo, pelo que entendi, colMeans não resolve seu problema?<br>
<br></div>M<br></div><div class="gmail_extra"><br><br><div class="gmail_quote">
2013/10/26 Maurício Lordêlo <span dir="ltr"><<a href="mailto:mslordelo@gmail.com" target="_blank">mslordelo@gmail.com</a>></span><br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
<div><div>
<div dir="ltr">Caros, <div>Devido a um problema com a função "bild" (pertencente ao pacote de mesmo nome), tive que modificar um script de simulação e não estou conseguindo finalizá-lo.</div><div>No CRM abaixo, a função "geral" ajusta o modelo e extrai os coeficientes de interesse. Em seguida, coloco um "for" com o número de simulações que desejo (coloquei 5 mas pretendo 2000). Ao usar o "replicate" obtenho os resultados. Preciso obter uma média por coluna destes resultados. Pelo que percebi, o objeto que denominei de "repeticoes" armazena um objeto em cada uma de suas linhas. Alguém saberia uma forma de solucionar isso, ou seja, encontrar uma média por coluna do objeto "repeticoes"? Caso não seja esta a opção adequada, aceito sugestões.</div>
<div>Agradeço desde já.</div><div><br></div><div><br></div><div><div>library(bild)</div><div>nt <- 8; ntr <- 3 ; nr <- 5</div><div>n <- ntr*nr*nt; nid <- ntr * nr </div><div>id <- rep(1:nid, each = nt) </div>
<div>tempo <- rep(1:nt, times = nid) </div><div>trat <- rep(factor(LETTERS[0:ntr]), each = nr*nt) </div><div><br></div><div>dados=data.frame(id,tempo,trat)</div><div><br></div><div>geral <- function() </div><div>
{</div><div> z <- numeric(n)</div><div> for(i in 1:n) </div><div> z[i] <- rbinom(1,1,0.5)</div><div> dat <- data.frame(dados, y=z)</div><div> (bild(y ~ trat+tempo, data = dat, time="tempo",start=NULL,dependence="MC1R"))@coefficients[1:3]</div>
<div>}</div><div><br></div><div>ns=5 ### número de simulações</div><div>for(i in 1:ns) </div><div>{</div><div>repeticoes <- replicate(1, geral(), simplify = FALSE)</div><div>print(repeticoes)</div><div>}</div></div><div>
<br></div></div>
<br></div></div>_______________________________________________<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.<span><font color="#888888"><br></font></span></blockquote>
</div><span><font color="#888888"><br><br clear="all"><br>-- <br>Manoel Galdino<br>
<a href="https://sites.google.com/site/galdinomcz/" target="_blank">https://sites.google.com/site/galdinomcz/</a><br>
</font></span></div>
<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></div></div><br></div>
<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></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></div></div>