Modeling conditional reference regions: Application to glycemic markers