<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 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>_______________________________________________<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><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>
</div>