testing dynamic regimes code

working for no weights, 1 beta, Y~1, 1 time point (gives same estimate
for any number of equal regimens - old code didn't)