Abstract
The medication possession ratio (MPR) method is commonly used for the determination of antiretroviral medication adherence. However, different ways of calculating MPR and methodological issues hinder the interpretation of the results and the reproducibility of the method. Thus, this study used three different models of MPR calculation and aimed to identify the one that best represents the situation of patient adherence. The results show that there was a statistically significant difference between the adherence rates determined by the three models, which indicates the need to specify the parameters used for calculation in the MPR method. However, the models individually were found to be related to viral suppression, but none of them had a greater effect than the other in this regard. The model that used residual medication (RM) and a fixed period of analysis allowed for a more precise identification of the number of doses that the patient used when compared to the others. Health services should avoid the application of the model using a variable analysis period. This study found that RM and the period of analysis considered are the main influencing factors in the accuracy of adherence results when the MPR method is used.
Introduction
The main goals of antiretroviral therapy (ART) are to reduce the morbidity and mortality of people living with HIV(PLWH), to increase life expectancy and quality of life, and to decrease the chance of virus transmission. In this context, adherence to treatment is one of the greatest challenges in caring for these patients.1–3 Adherence to ART is a multifactorial and complex process with no consensus regarding the optimal method for its determination. The assessment of ART adherence includes self-reports, pharmacy records, electronic monitoring devices, therapeutic serum level measurements and daily drug registration.3–12
Pharmacy records are widely used because the data can generally be accessed easily by the health service, incurs a low cost for execution and provides more objective information.7,13,14 To ensure best operation of this technique, certain aspects should be considered. The locality from which the patient obtains medication is one of them, because the information is more reliable when a single establishment or system records the withdrawal dates of antiretrovirals (ARVs).13,15 However, when the residual medication (RM) of previous withdrawals is not considered, the patient may eventually be incorrectly classified as non-adherent. 13 There are different ways of calculating adherence from pharmacy records. One of these is the medication possession ratio (MPR), which measures the amount of time the individual is in possession of ARVs. 7
Although pharmacy records are correlated to virological and immunological measures, some studies show that patients with high adherence rates do not achieve viral suppression, while others with low adherence maintain an undetectable viral load (VL). These results may suggest, among other reasons, that there are methodological problems in determining adherence using this method. 7 In this sense, this research aimed to recommend, after assessing three models of adherence calculation by the MPR method, the model that satisfactorily represents the real situation of the patient regarding the use of ARVs.
Methods
This study was carried out at the Testing and Counselling Centre (TCC) located in Curitiba, Paraná State, Brazil. The study was conducted from November 2016 to July 2017. Volunteers, both genders, over 18 years of age, not pregnant and with at least one year and one month of treatment with ARVs, were randomly selected at the TCC during the monthly ARV withdrawal. The results are discussed from the perspective of the statistical difference between the percentages of adherence determined by the three models and their associations with viral suppression, as well as the classification of patients as adherent or not, at distinct cut-off points.
Data collection and adherence measurement
The adherence was determined by using the MPR method. The ARV withdrawal dates of each participant, retrospective over one year, as well as the most recent VL within that period, were obtained from the Logistic Control System of Medicines – a national unified system.
MPR is the ratio between the number of days the patient has medication (numerator) and the number of days of a predetermined period for the study (denominator). The result, when multiplied by 100, indicates the percentage of days that patient would have taken the drug within that time interval. This study used three different ways to calculate MPR 14 that differed by the use (or not) of RM and by the analysis period applied in the denominator. RM can be defined as the doses that the patient had in their possession of a drug withdrawal prior to the first withdrawal within the study period.
The three models exemplified in Table 1 consider a study period from 1 January to 31 December 2017 (365 days). In order to evaluate the presence or absence of RM, the last ARV withdrawal in 2016 is also presented. At each withdrawal, the patient receives medication for 30 days of treatment. Figure 1 shows the similarities between each model.

Similarities and differences between the three models. NRMF – not used RM in a 365-day fixed analysis period; NRMV – not used RM in a variable analysis period; RMF – used RM in a 365-day fixed analysis period.
Description of the three models of MPR calculation.
MPR: medication possession ratio; MW: medication withdrawal; RM: residual medication.
a11 MW × 30 days.
Sample size
The sample was determined based on the number of patients withdrawing ARVs in the TCC in 2016 (4322). Considering a 95% confidence interval, a margin of error of 5%, a non-adherence rate to ART of 30% 6 and a margin of 10% to compensate for losses, the minimum sample required was 331 patients.
Statistical analysis
To compare the three models, a mixed linear regression model was fitted in which the errors followed a normal distribution, including a random effect to accommodate repeated measures in the same patient. 16 The likelihood ratio test (LRT) was used to evaluate the effect of the three proposed models on the withdrawal record, with a significance level of 5%. The results were reported in terms of the adjusted means, standard error, 95% overall confidence intervals and Wald test to test the contrasts of the variable in the selected model, with the Bonferroni correction. 17
To evaluate the association between the three models and viral suppression, two procedures were performed through multiple logistic regression. In the first, the logistic regression was adjusted for each model separately to assess if they presented statistical significance, demonstrating the results in terms of Akaike Information Criterion (AIC), odds ratios and p-values. Then, it was determined if all models had the same correlation with viral suppression using the Vuong’s LRT 18 for strictly non-nested models, which is based on AIC.
R software 19 was used to perform the analysis, using the Ime4, 20 ImerTest 21 and Nonnest2 22 packages (May 2019, Vienna, Austria).
Ethical approval
The study was approved by the Ethics Committee of the Health Sciences Sector of the Federal University of Paraná (protocol 1739379) and the Municipal Health Department of Curitiba (protocol 1620414), Brazil.
Results
A total of 349 volunteers were included in the study. In the sampled population, there was a predominance of men (70.8%) of white ethnicity (69.9%), heterosexual (49%), with no fixed partner (51%) or unmarried (52.2%), without children (57%), and with a mean age of 44.9 years. The majority had more than 11 years of education (77.9%), were Catholic (45.6%) and had more than five years of ARV use (60.2%). Most of the participants used the Unified Health System, provided by the government, for their medical care (59.6%) and 6.3% had a detectable VL in the interview period.
The comparison of the mean adherence rates determined by the three models is presented in Table 2. Although the nominal values between RMF and NRMV were similar, there was a statistically significant difference in the results, as well as between RMF versus NRMF and NRMV versus NRMF, as shown in Table 3.
Comparison of means of adherence between the three models.
CI: confidence interval; NRMF: not used residual medication in a 365-day fixed analysis period; NRMV: not used residual medication in a variable analysis period; RMF: used residual medication in a 365-day fixed analysis period.
Difference between the adjusted means of the three models.
CI: confidence interval; NRMF: not used residual medication in a 365-day fixed analysis period; NRMV: not used residual medication in a variable analysis period; RMF: used residual medication in a 365-day fixed analysis period.
aWald test.
It is important that patient adherence rates, regardless of the methodology used for calculation, correlate with virological results in order to be relevant. This is demonstrated in Table 4, based on logistic regression considering each model individually. All models were found to be related to viral suppression, but none of them had a greater effect than the other in this regard (Table 5), demonstrated by the p values and because of the small difference in AIC between the two models (Table 4).
Correlation of adherence rates with viral suppression.
AIC: Akaike information criterion; CI: confidence interval; NRMF: not used residual medication in a 365-day fixed analysis period; NRMV: not used residual medication in a variable analysis period; RMF: used residual medication in a 365-day fixed analysis period.
aWald test.
Difference in correlation with viral suppression.
AIC: Akaike information criterion; CI: confidence interval; NRMF: not used residual medication in a 365-day fixed analysis period; NRMV: not used residual medication in a variable analysis period; RMF: used residual medication in a 365-day fixed analysis period.
aVueng test.
Analysing the differences in adherence obtained by the models, the adherence rates were categorised as ≥95%, ≥90%, ≥85% and ≥80%. The models were compared to each other by evaluating the occurrence of divergence in patient classification with respect to adherence and cut-off points. Table 6 presents these results, but does not include the comparison between the NRMF model and the other two because no patients with disagreement in the classification of adherence were identified.
Comparison between models in relation to the different adherence cut-off points.
A/NA: number of patients considered adherent by the first model and not adherent by the second one, considering each cut-off point; NRMF: not used residual medication in a 365-day fixed analysis period; NRMV: not used residual medication in a variable analysis period; RMF: used residual medication in a 365-day fixed analysis period; T: total number of participants.
It was observed that the divergence between the models was always higher with larger cut-offs, regardless of the comparison being made. We highlight the differences in the classification of adherence between the 90 and 95% cut-off points, especially for the comparison between RMF versus NRMF and NRMV versus NRMF, in which case approximately half of the individuals at the cut-off point of 95% had an erroneous categorisation of adherence.
Discussion
Among the most widely used methods to determine adherence to ART, the use of pharmacy records can be considered the most adequate in relation to the reliability, availability and facility of data collection, since techniques involving self-reports, for example, may overestimate adherence. 23 This study demonstrates that a detailed description of data collection considering withdrawal dates and periods of analysis considered is important to present the results safely. Thus, in the case of analyses that involve the use of MPR, one should describe how the numerator and denominator of the equation were defined.
This study revealed that the three models give statistically different adherence rates (Table 3). Patient adherence when assessed by the RMF and NRMV models resulted in values, on average, 4.7 and 5.6% higher than with the NRMF model, respectively. This means, for example, that for a one-year period of treatment, the adherence rate calculated by NRMV would indicate approximately 20 more days of therapy than by NRMF.
It was observed that NRMV gave a higher average adherence value than the others (Table 2), which can be explained by the fact that the analysis period was less than 365 days in most cases. This is a limiting factor in the use of this model since it considers the period of analysis starting at the date of the first ARV withdrawal within the study range. For example, if the study period is from 1 January to 31 December and drug withdrawals within this range occurred only on 1 August, September, October, November and December (150 tablets in 153 days), the MPR determined by NRMV will be 98% ([150/153] × 100), indicating adequate adherence to ART. However, if the study period of one year (RMF or NRMF models) is considered, adherence will be 41%, which is a better approximation of reality, since the patient remained without medication for seven months. Therefore, the use of variable denominator models must be carefully evaluated on a case-by-case basis so that patients with inadequate adherence are not misclassified as adherent. These outcomes corroborate the findings of Kozma et al., 24 who suggest the use of a fixed denominator for the assessment of adherence to long-term therapies such as ART.
Observing the adherence determined by RMF and NRMF, which differed only in the use of RM, it was expected that the values would be lower with NRMF. Also, if the patient does not have RM, the adherence rate calculated by the two models will be the same. However, the use of NRMF can misclassify individuals as non-adherent when in fact they are adherent. For a cut-off point of 95%, for example, it was observed that approximately half of the patients (47.5%) were mistakenly classified (Table 6), which occurred similarly between NRMV and NRMF; this is a relevant finding when the researcher needs to define which model will be used in a study. Thus, the RMF model indicates adherence to ART more reliably, which corroborates the findings of Boer et al. 13
It should be emphasised that, in the health services, standardised, cost-effective and routine adherence monitoring is necessary, since adherence to chronic treatment is dynamic and oscillates over time. The identification of patients having difficulties with ART allows for intervention actions by the health team, with a view to promoting an adequate therapeutic efficiency of ARVs.7,8,23 In this sense, it would be possible to use NRMF in the scope of professional practice, as it is simpler with regard to data collection; it is not necessary to seek withdrawals prior to the observational period. In this situation, given the 4.7% increase in the adherence rate obtained by the NRMF model, the result is close to what would be found by RMF, as recommended in this study (Table 3). Therefore, the use of NRMF would only be feasible with the use of this factor, which was achieved in a TCC with approximately 4000 registered patients and could be extrapolated to other services with similar characteristics.
Regarding the period of analysis considered in the studies, it is noted that the use of short periods for adherence calculation using pharmacy records could give imprecise results regarding the association between adherence and viral suppression. 8 Studies that contemplate longer time intervals, such as the 12-month period considered in this research, tend to provide results with greater accuracy. It is important to point out that these considerations are relevant for patients who have been on ART for some time before joining a study. For those starting treatment, as long as the study period coincides with the date of the user’s first ARV withdrawal, model selection will not make a significant difference because there will be no RM.
Another consideration that must be made concerns the cut-off point of adherence. Although patient retention in HIV care services is associated with viral suppression, adherence to ARVs is the factor that best relates to this condition. 25 Inadequate treatment follow-up is a strong indicator of therapeutic failure, development of drug resistance, disease progression and death. 15 The 95% cut-off value of adherence is the gold standard for older therapies, while studies point out that use of at least 80% of the prescribed doses is necessary to obtain an effective therapeutic response, especially for newer ARVs.4,15,26 It has been noted that rates of 80% or less do not necessarily mean that the individual has a detectable VL. However, these cases should represent an alert to the health service, since these individuals may be at the limit of virus detection or there may be possible difficulties in the follow-up. For this reason, the cut-off point defined to classify patients as non-adherent should guide the choice of the calculation model to be used. It is enhanced here in the non-recommendation of the NRMF model, mainly for cut-offs above 90%, when the discrepancies are more remarkable (Table 6). It is evident that, with smaller cut-offs, the differences in this classification are less impactful, but no less important to the scientific aspects of the publication of results.
In practice, what is the consequence of misclassifying non-adherents? Especially for the management of the health service, this situation may result in patient rescue procedures that require the time and resources of different health professionals, both of which are generally tight in public services. As a result, there is a waste of resources with patients who are under adequate therapeutic follow-up, leaving those who may actually need assistance without intervention. In addition, individuals with adequate adherence to ART achieve such low VL levels that there is little chance of transmitting the virus. Thus, the correct identification of patients with adherence problems is of social relevance, since this group must be monitored by health professionals and encouraged to return to treatment. It is also emphasised that adequate therapeutic follow-up and clinical monitoring of VL are strategic components of combined HIV prevention, which is designed to make joint use of biomedical, behavioural and structural interventions applied to the individual, their relationships and social groups to which they belong.
It has been shown in this study that the adherence rates obtained by the three models are statistically different, to a greater or lesser degree. The NRMV and NRMF models are simpler in relation to data collection, but the researcher should know that the adherences rates calculated by NRMF will be smaller when compared to the others and may not represent the reality of the studied individuals. Likewise, NRMV may mistakenly indicate individuals as adherent when in fact they are not. In this context, the importance of the use of a scientific methodology with clearly defined validity, safety and precision is pointed out, in order to ensure that the results can be used in the practice of health services with the main objective of adequately identifying patients with adherence problems.
Conclusion
In this study, we identified that the three proposed models showed statistically significant differences in adherence rates, but all of them correlated with viral suppression. The RMF model allowed for a more precise identification of the number of doses that the patient used because it considered RM. In addition, this model reduces the possibility of errors in the classification of adherence. Furthermore, in the RMF model, the time interval previously defined in the study is in fact considered, in contrast to the NRMV model. Thus, application of the NRMV model must be carefully evaluated, given the possibility of mistakenly classifying adherent individuals. Health services should be therefore advised that, when using the MPR method to define adherence rates, the RMF model should be applied, avoiding the use of NRMV in particular.
In studies that use the MPR calculation, it is necessary to describe whether RM was used and the period of analysis considered, since these elements are the main influencing factors regarding the accuracy and safety of the results.
Footnotes
Declaration of conflicting interests
The authors declared no potential conflicts of interest with respect to the research, authorship, and/or publication of this article.
Funding
The authors received no financial support for the research, authorship, and/or publication of this article.
