Displaying 1 result from an estimated 1 matches for "cat_sum1".
2010 Feb 17
1
Replicating output from a function
....21)
mod1<-factor(c(1,2,2, rep(c(1,2,3),3)))
mod2<-c(1,2,15,rep(3,9))
datas<-data.frame(id,n,r,mod1,mod2)
# intermediate level fuction (courtesy of Hadley Wickham):
pick_one <- function(x) {
if (length(x) == 1) return(x)
sample(x, 1)
}
# Function that I want replicated 100 times:
cat_sum1 <- function(meta, mod) {
m <- meta
m$mod <- mod
meta <- ddply(m, .(id), summarize, r = mean(r), n=mean(n), mod =
pick_one(mod))
meta$z <- 0.5*log((1 + meta$r)/(1-meta$r))
meta$var.z <- 1/(meta$n-3)
meta$wi <- 1/meta$var.z
return(meta)
}
# output from 1 run:
cat_...