posterior_traj function

Estimate the subject-specific or marginal longitudinal trajectory