AR(1) in the error matrix rather than in the linear model

nlme uses a slightly different paramerization which leads to sigma being no longer sigma but something like sigma / sqrt(1 - ar^2) or something like that (hence higher sigma being reported for higher ar). I am at my phone so please don’t take my word for the formula.

brms parametizes correctly so that sigma is in fact what we expect it to be in this case.