[R-br] média em um "replicate"

Manoel Galdino mcz.fea em gmail.com
Sábado Outubro 26 09:12:13 BRST 2013


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 <mslordelo em gmail.com>

> 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
> R-br em listas.c3sl.ufpr.br
> 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/
-------------- Próxima Parte ----------
Um anexo em HTML foi limpo...
URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20131026/eacc9d6a/attachment.html>


Mais detalhes sobre a lista de discussão R-br