dat <- Multiplyr (G = rep(c("A", "B"), each=50),
m = rep(c(5, 10), each=50),
alloc=1)
dat %>% group_by (G) %>% mutate (x=rnorm(length(m), mean=m))
dat %>% within_group ({
mdl <- lm (x ~ 1)
})
dat %>% summarise (x.mean = coef(mdl)[[1]])
dat %>% shutdown()
Run the code above in your browser using DataLab