g1rfinal = lmer(log(trmLatency) ~ WMC + order + distance + DomSide + (1 + BlockType |
pno) + (1 + firstBlType | pno) + BlockType + (1 + distance | pno) + WMC:order + WMC:distance + WMC:DomSide + order:DomSide +
distance:DomSide + WMC:BlockType + distance:BlockType + DomSide:BlockType +
order:BlockType, data = g1Study)

g1Study = cbind(g1Study, predict = exp(predict(g1rfinal)))
g1test = g1Study %>% group_by(BlockType, order) %>% summarise(avgrt = mean(predict))

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))