- rd(t)~N(t) plots

-theta model per macrounit (does not work yet)