Modelling the impact of protein-kinase R allelic variant on HIV biomarkers trajectories by means of latent class mixed models