I am using the wbm() from the panelr package by Prof Long (https://panelr.jacob-long.com/articles/wbm.html) to fit a within-between (mixed-effects) model.
wbm()