changed order of parameters for avg_mu so it could be used with
DeltaMethod function in mrds
created compute_Nhat for se(Nhat) via DeltaMethod; pulled this code from
plotfit
Addition of code from paper in demo; added std err calculation for N