gistfile1.txt

mysim <- function(){
  des <- expand.grid(
    congruence=c(1,0),
    cuedresp=c("L", "R", "None"),
    trial=1:5,
    block=1:10,
    participant=1:20
  )
  samp <- des %>% rowwise() %>% mutate(y=rbinom(1,1,.03*congruence))
  modelfit <- tidy(lmer(y~congruence +(1|participant), data=samp))
  est <- modelfit %>% filter(term=="congruence")
  return(est$estimate/est$std.error)
}

sims <- data.frame(simindex=1:10) %>% rowwise() %>% mutate(e=mysim())
summary(sims)