I'm interested in using the user-written traj command (link below) to identify latent trajectories of change in patients' BMIs. The command has useful features like joint trajectory modeling and accounting for non-random attrition.

However, patients in my dataset are nested within physicians; I have unique physician identifiers for each physician.

Questions:

1. Is there a method or workaround that would allow traj to account for hierarchically nested data?
2. If not, to what degree would traj be robust to violation of the assumption that patients are independent of each other?


https://www.andrew.cmu.edu/user/bjones/