daily pastebin goal
66%
SHARE
TWEET

Untitled

a guest Dec 9th, 2018 65 Never
Not a member of Pastebin yet? Sign Up, it unlocks many cool features!
  1. g1rfinal = lmer(log(trmLatency) ~ WMC + order + distance + DomSide + (1 + BlockType |
  2. pno) + (1 + firstBlType | pno) + BlockType + (1 + distance | pno) + WMC:order + WMC:distance + WMC:DomSide + order:DomSide +
  3. distance:DomSide + WMC:BlockType + distance:BlockType + DomSide:BlockType +
  4. order:BlockType, data = g1Study)
  5.    
  6. g1Study = cbind(g1Study, predict = exp(predict(g1rfinal)))
  7. g1test = g1Study %>% group_by(BlockType, order) %>% summarise(avgrt = mean(predict))
  8.    
  9. ggplot() + geom_point(data = g1test, aes(x = order, y = avgrt, color = BlockType)) + geom_line(data = as.data.frame(effect("order : BlockType", g1rfinal)), aes(x = order, y =exp(fit), color = BlockType))
RAW Paste Data
We use cookies for various purposes including analytics. By continuing to use Pastebin, you agree to our use of cookies as described in the Cookies Policy. OK, I Understand
 
Top