Eu não consegui instalar a biblioteca build no R 3.0 De todo modo, pelo que entendi, colMeans não resolve seu problema?
M 2013/10/26 Maurício Lordêlo <[email protected]> > Caros, > 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. > 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. > Agradeço desde já. > > > library(bild) > nt <- 8; ntr <- 3 ; nr <- 5 > n <- ntr*nr*nt; nid <- ntr * nr > id <- rep(1:nid, each = nt) > tempo <- rep(1:nt, times = nid) > trat <- rep(factor(LETTERS[0:ntr]), each = nr*nt) > > dados=data.frame(id,tempo,trat) > > geral <- function() > { > z <- numeric(n) > for(i in 1:n) > z[i] <- rbinom(1,1,0.5) > dat <- data.frame(dados, y=z) > (bild(y ~ trat+tempo, data = dat, > time="tempo",start=NULL,dependence="MC1R"))@coefficients[1:3] > } > > ns=5 ### número de simulações > for(i in 1:ns) > { > repeticoes <- replicate(1, geral(), simplify = FALSE) > print(repeticoes) > } > > > _______________________________________________ > R-br mailing list > [email protected] > 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. > -- Manoel Galdino https://sites.google.com/site/galdinomcz/
_______________________________________________ R-br mailing list [email protected] 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.
