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.
