Afer having created an lmer model:
my_m = lmer(y ~ x1*x2 + (1|id) + (1|tp), data=my_data)
I now want to extract the conditional modes of the