Abstract Background This study aimed to jointly model HIV disease progression patterns based on viral load (VL) among adult ART patients adjusting for the time-varying “incremental transients states” variable, and the CD4 cell counts orthogonal variable in a single 5-stage time-homogenous multistate Markov model. We further jointly mapped the relative risks of HIV disease progression outcomes (detectable VL (VL ≥ 50copies/uL) and immune deterioration (CD4