The data were simulated like
y <- rpois(1, exp(fe + x1 + x2 + 2.5*x1.lead))
where fe, x1, and x2 are standard normal random variables.
fe varies only across id.
x1.lead is a one period lead of x1 which causes the violation
of the conditional mean assumption.