

Youk AO, Mundy LM, St. Laurent S, Lann MF and Bowlin SJ. Assessment of risk in longitudinal observational studies with time-varying drug exposures for treatment of a chronic disease. Epidemiol Rep. 2015; 3:2. http://dx.doi.org/10.7243/2054-9911-3-2
Ada O. Youk1†*, Linda M. Mundy2†, Samantha St. Laurent2†, Michael F. Lann1† and Steven J. Bowlin3†
*Correspondence: Ada O. Youk ayouk@pitt.edu
†These authors contributed equally to this work.
1. Department of Biostatistics, University of Pittsburgh, Graduate School of Public Health Pittsburgh, PA 15261, USA.
2. GlaxoSmithKline, Collegeville, PA, 19426, USA.
3. Shire, Bedminster, NJ, 07921, USA.
This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/3.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.
Background: Drug exposures can be continuous or intermittent with several competing exposures which
can influence the analysis of exposure-response relationships in longitudinal observational data.
Purpose: We report the in-depth application of research methods employed in a nested case-control study to assess risk for fracture among persons with human immunodeficiency virus (HIV) infection treated with antiretroviral (ARV) therapy. These pharmacoepidemiological research methods, extrapolated from occupational epidemiology, permitted detailed exposure-response analysis in longitudinal observational data sources.
Methods: Multidrug administrative prescription claims data were reformatted to create time-varying, nonoverlapping drug exposure histories that preserved the information in the original prescriptions. We applied research methods from occupational health to estimate drug-class (N=5) and drug-specific (N=29) exposure summary measures associated with fracture events.
Results: Using pharmacy claims from the Clinformatics™ DataMart Multiplan (IMPACT) database, a product of OptumInsight Life Sciences, Inc., 1,929,067 prescription histories from 30,405 ARV-treated HIV subjects were restructured retaining gaps and removing overlaps and encompassments in the prescription histories. Most subjects were dispensed at least three different ARV drugs; over 75% had prescription histories with at least one gap in prescriptions and 23.5% had over 65 prescription claims. Drug exposure summary metrics were computed from the restructured claims to determine time-varying exposureresponse relationships for individual ARV drug exposures and fractures.
Conclusions: We report a research method to assess drug exposure using administrative prescription claims in observational data that goes beyond dichotomous or cumulative exposures measurements of HIV drug classes and specific ARV drugs. The time-varying measures derived here are easily applicable to other observational data sources such as drug exposures for other chronic diseases where prescription histories are complex and evolve over time.
Keywords: Pharmacovigilance, pharmacoepidemiology, antiretroviral (ARV), prescription claims data, benefit-risk analysis
Advancing age, comorbid conditions, and cumulative drug exposures contribute to long-term health outcomes among persons with chronic diseases [1-3]. Historically, observational studies using administrative pharmacy claims characterize drug exposure as either dichotomized variables or as cumulative exposure for drug class, combination regimens, or specific drugs [4]. Measurement bias in prescribed drug exposures measured in administrative claims data may under- or overestimate risk identified in exposure-response relationships of health outcomes and in pharmacovigilance, and is of concern in this era of comparative effectiveness research [5,6]. While randomized trials provide enhanced understanding of treatment efficacy, often well-designed studies with adequate person-years of observation in epidemiology studies are needed for longer-term assessments of health outcomes, drug safety, and benefit-risk balance [7,8].
Recently, we used time-varying exposure metrics using data from an administrative claims database and reported the differential low-impact fracture risk of 22 antiretroviral (ARV) drug exposures among persons treated for human immunodeficiency virus (HIV) infection [9]. The research methodology for this nested case-control study was used because of the time-varying and complex ARV exposures in our observational data source [9]. The analytical process included the reformat of the administrative drug prescription claims data to drugexposure histories for computation of time-varying summary measures of exposure utilizing the Occupational Mortality Analysis Program (OCMAP-Plus®-http://www.ibridgenetwork.org/university-of-pittsburgh/ocmap-software). This structured approach permitted the extrapolation of exposures and risks, as defined in occupational epidemiology, to our assessment of drug exposures and event risk in pharmacoepidemiology [10-14]. The detailed prescription histories for longitudinal drug exposures are often as complex as occupational exposure histories [10-14]. We computed typical occupational exposure metrics (duration of exposure, cumulative exposure, and average intensity of exposure) in our nested case-control study, along with the assessment of drug lag periods and latency, defined as time since first drug exposure [10-14]. By reformatting prescription claims data and using these metrics, we explored drug-exposure relationships in HIV treatment and risk of fracture. The methods, however, are applicable to any disease where partitioning the disease and complex treatment exposures are needed to further inform drug exposure response relationships. In this report, we illustrate in detail the application of these occupational exposure methods for a nested case-control study of ARV treatment exposure and response defined as risk for fracture. This framework offers application in other pharmacoepidemiology settings, comparative effective research methods, and in-depth assessments of treatment benefit-risk balance.
Study population
The case-control study population was selected from a healthinsured
cohort of 59,584 adults with HIV infection enrolled in
the Clinformatics™ DataMart Multiplan (IMPACT) database, a
product of OptumInsight Life Sciences, Inc, between January 1,
1997 and March 31, 2008. The database was Health Insurance Portability and Accountability Act (HIPAA) compliant; HIV
infection was identified by ICD-9-CM codes 042 or v08 and
other enrollment criteria, definitions of cases and controls,
and other study methodology were previously reported [9].
Initial assessment of exposure-response relationships: ARV
drug class prescription histories utilized as dichotomous
exposures
We performed an initial feasibility assessment of the source
data using years 1997-2006 to explore ARV prescription histories
and their relationship to risk of fracture. By class, there
were 29 nucleoside reverse transcriptase inhibitors (NRTI),
nine non-nucleoside reverse transcriptase inhibitors (NNRTI),
27 protease inhibitors (PI), and one fusion inhibitor (FI). Given
the broad distribution of ARV exposures, the exposures were
grouped by a two-step categorization for each ARV drug class.
First, subject-level exposure to each of the four drug classes
was categorized as a dichotomous exposure (exposed or not
exposed across the whole prescription history). We further
partitioned the NRTI drugs as either nucleoside reverse
transcriptase inhibitors (NsRTI) or nucleotide reverse transcriptase
inhibitors, (NtRTI). Then, the subjects' were grouped
by cumulative combinations of ARV drug exposures based
on national treatment guidelines for the entire observation
period (Table 1). Fracture incidence rate ratios were computed
but did not include time-varying exposure measures and
were strictly static measures of exposure computed across
the total prescription history.
Table 1 : Initial assessment of antiretroviral drug prescription claims among 23,207 subjects with HIV infection during calendar years 1997-2006: Groupings of exposures based on National treatment guidelines.
Assessment of exposure-response relationships using a
time-varying strategy
Step 1: Convert claims into prescription history records
The final drug dataset contained 1,929,067 single drug ARV
prescription history records after the 435,319 original combination
ARV drug claims were converted to 958,509 single ARV
drug claims and combined with the 970,558 original single
drug claims for years 1997 quarter 1 through 2008 quarter 2.
The standard components of each claim included sales date,
drug name, strength, days supplied, and metric quantity
dispensed. Initially, there were 40,105 ARV claims (2.1%) with
one or more missing standard components of the prescription.
In addition, 16,379 prescriptions (0.8%) had zero for days
supplied, and 107,348 prescriptions (5.6%) had a cost-driven
vendor-derived imputation for metric quantity. For the 56,484
(2.9%) prescription claims with missing, or zero, for days supplied,
we imputed the mode of the frequency distribution for
days supplied for each specific ARV drug-strength.
Step 2: Deriving the daily-dose estimates of drug exposure
To quantify the time-varying drug exposures for the 68 different
drug and drug strength combinations, estimates of
daily dose were needed. We calculated the drug-specific
daily doses as drug strength and a multiple of pills per day
(PPD) which is a function of metric quantity. For those drugs that were missing metric quantity, we assessed the vendor
supplied cost-derived imputations from the prescription
claims. The estimated values of metric quantity were not
equivalent with either clinically-derived values or datadriven
estimates for metric quantity. Therefore, we used an
alternative method to assign daily-dose estimates for those
drug-strengths with missing metric quantity that combined
the prescription claims data and standard dosing practices.
When the calculated daily dose for >80% of the prescription
claims for an ARV drug-strength was equivalent to the
standard adult daily dose defined from the drug label, that
daily dose was assigned to all of the prescription claims for
that ARV drug-strength. The ARV drug-strength claims that
met this >80% threshold for daily-dose assignment primarily
included drugs that were either once-daily regimens or
single-strength formulations. For drug-strengths that did
not meet the >80% equivalency threshold for assigning a
daily dose, we assessed the drug-strength prescriptions by formulation (elixir, power, tablets, etc.) and created groupings
of daily-dose estimates for imputation. Daily-dose estimates
were assigned to 22 drug-strength combinations that did not
meet the >80% threshold (Appendix 1).
Step 3: Restructure the prescription claims data
To utilize the OCMAP-Plus® program for computing the ARV
drug exposure metrics, it was necessary to restructure the
pharmacy history records to remove encompassed and
overlapping prescriptions with the same drug-strength.
While dispensed drugs most likely overlapped as real time
prescription sales, the assumption was that drug ingestion
did not overlap. An overlapped record occurred when
the prescription start date of the current prescription was
before the previous prescription end date and the current
prescription end date was after the previous prescription's
end date. An encompassed record occurred when the current
prescription start and stop dates were before the previous
prescription end date. An ARV drug gap in the prescription
record was one or more days' interval with no drug claim for
the same ARV drug-strength exposure. A continuous exposure
record had the same ARV drug-strength prescriptions with
no overlaps or encompassments. When a prescription for
an ARV drug was overlapping with another prescription of
the same drug-strength, the overlap was removed by taking
the days of supply that constituted the overlap and adding
the days to the end of the first prescription, pushing out the
second prescription the same number of days as the overlap;
this removed the overlap but preserved the days of supply.
When a prescription for an ARV drug was totally encompassed
by another prescription of the same drug-strength, the last
date of the first prescription was extended forward in time
the same amount of days as the encompassed prescription
to ensure the correct number of days supplied was retained.
Gaps between prescription claims were retained when they
occurred in the prescription history.
A systematic, quality-control check of the restructured data was performed by stratifying the subset of the cohort exposed to ARV drugs into approximate quartiles of number of prescriptions per subject (<15, 15-34, 35-64, >65 prescriptions). Within each quartile a random sample of 25 subjects was selected for a manual review of the restructure process. Each of these 100 records were manually restructured and then compared to the computerized restructured record to check for accuracy. No errors were found during this quality-control assessment. We used PERL language in a MySQL database to perform the computerized restructure of the pharmacy claims.
Step 4a: Construct risk sets for exposure-response modeling
Traditional analyses of exposure-response typically include
the use of Poisson regression of the internal cohort rates. Relative
risk regression (RRR) is an alternative to the traditional
analyses and was used in our nested case-control study of
ARV drug use and bone fracture risk [9]. For each case, a risk set was defined as the case plus all cohort members alive
and at risk for fracture but not yet having had a fracture (noncases)
who were the same age of the case. We used age as
the time dimension to construct these full risk sets (using the
RISKSET module in OCMAP-Plus®) [10] from the HIV cohort.
Risk sets were further matched on gender and year of birth
(within five years) to control for birth cohort effects because
individuals of the same age were not necessarily born in the
same year. We randomly selected four non-case subjects
(controls) from each risk set. For each non-case in each risk
set, each exposure variable was evaluated as of the age of
the corresponding case with fracture.
Relative risk regression (RRR) was used to investigate the dependence of internal cohort fracture rates (modeled as time to diagnosis) on the duration of time exposed to a drug and separately to the cumulative exposure of the drug with adjustment for the potential confounding factors. Multiplicative RRR models of the form λ(t)=λ0(t) exp{x(t)β} were fit to the internal cohort rates using stratified conditional logistic regression [15-17]. In this model, λ0(t) is the hazard of an event at time t for an individual with baseline levels of all covariates, x(t) is a vector of covariates (exposures and potential confounding variables evaluated as of each event time t), and β is the corresponding parameter vector estimated by partial likelihood [15-17]. The RRR modeling was performed using the conditional logistic regression module in Stata [18]. The risk sets were explicitly constructed for the objectives of the nested case-control study of ARV exposures and fracture events, and the conditional logistic regression algorithm to compute the estimates of relative risk was equivalent to fitting a Cox model. When there were multiple time-varying covariates (i.e., ARV drug exposures), the most feasible way to fit these models was to explicitly construct risk sets and to estimate parameters using a conditional logistic regression [15-17].
Step 4b: Summary measures of ARV drug exposure
The restructured prescription data were used for the computation
of time-varying summary measures of drug exposure
using the RISKSET program module in OCMAP-Plus® [10].
An ARV drug exposure period was the entire time period of
exposure to the specific ARV drug in the person's prescription
history. The cumulative measure of exposure for a drug
(Drug_Cumj) computed over Nj prescriptions during drug
exposure period j was given by:
(1)
where Timei and Exp_dosei represented the duration of the prescription (in days) and daily dose of the drug (in milligrams), respectively, of the ith prescription during drug exposure period j.
The duration of exposure for a drug (Drug_Durj) computed over Nj prescriptions during drug exposure period j was given by:
(2)
where Timei represented the duration of prescription (in days) and Exp_YNi was equal to 1 if exposed to that drug and 0 otherwise. These quantities were specific for the ith prescription during drug exposure period j.
The average intensity of exposure for a drug (Drug_AIE) computed over Nj prescriptions during drug exposure period j was:
(3)
Step 4c (if warranted): Assess drug exposure lags and latency
The internal event rates relative to cumulative drug exposures
in the regression models assumed that an event
during an observation period was related to the cumulative
drug exposure from study entry up to the event. The actual
exposure-response relationship may change over time and
may be negligible during parts of the drug history. To explore
for this possible changing exposure-response relationship,
alternative characterizations of the above summary measures
were computed by incorporating a lag period [10,19]. Lag is
the period of time before the event during which drug exposures
are not counted. This is based on the assumption that a
latent period exists between when an individual is exposed
and when the exposure may possibly impact the disease [19].
In occupational epidemiology, the choice of the lag period
depends on a known or hypothesized exposure-response
and latency (time since first exposure). For the assessment of
ARV drug latency and fracture, latency (defined as the time
from first exposure to development of fracture) remained
exploratory so we used lag periods of 3, 6 and 12 months.
Initial assessment of exposure-response relationships: ARV
drug class prescription histories utilized as dichotomous
exposures
The initial analysis of the cohort was based on 23,207 subjects
with HIV infection from calendar years 1997-2006. In
this feasibility assessment, the ARV prescription claims were
grouped into four categorical combination regimens (Table 1).
The majority of the claims were for NsRTI/PI, NsRTI/NNRTI, and
NsRTI/NNRTI/PI respectively. There were 1,096 subjects with
fracture, for a fracture incidence rate of 2.68 per 100 person
years (95% CI: 2.46-2.91) for subjects without ARV exposure
(reference group). When comparing the subjects with ARV
exposure to those no ARV exposure, we estimated incidence
rate ratios that ranged from 0.46-1.19 for subjects with various
combinations of drug class exposures (data not shown).
When comparing those with no ARV exposure to regimens
with and without tenofovir, we estimated an incidence rate
ratio of 2.57 (95% CI: 2.32-2.83) among subjects on combination
regimens without tenofovir, 1.16 (95% CI: 0.94-1.38) among subjects on combination regimens with tenofovir,
and 3.19 (2.34-4.04) among subjects on other ARV regimens
(data not shown). The structure of the prescription histories
did not permit ARV-class or ARV-drug specific assessments
of risk for fracture.
Assessment of exposure-response relationship
An example of the time-varying drug exposure restructuring
and computation of summary measures of ARV drugs
The prescription data for one study participant showing the
original overlapped and encompassed ARV prescriptions, as
well as gaps, are displayed in Figure 1a. The dispensed drugs
were a combination regimen of stavudine, lamivudine, and
ritonavir. The 1st ritonavir prescription was filled on October
27, with the 1st stavudine and lamivudine prescriptions filled
on November 9. The 2nd stavudine and lamivudine prescriptions
on December 4 overlapped with the 3rd stavudine and
lamivudine prescriptions filled on December 29, while the 2nd
and 3rd ritonavir prescriptions occurred one day apart with
the 3rd totally encompassed within the 2nd, after a gap of 20
days from the 1st prescription of ritonavir. The restructure
of the original pharmacy claims, inclusive of the prescribed
dose, resulted in daily, multidrug exposure intervals while
preserving the original number of days supplied and drug
strength for each prescription (Figure 1b).
Figure 1 : Antiretroviral administrative prescription claims data from original pharmacy file by date prescribed and number of days supplied (a) and as a restructured pharmacy history (b) for assessment of risk in the Occupational Mortality Analysis Program (OCMAP-Plus®).
The details for the computation of the summary exposure measures for this subject are displayed as the original prescription claims data (Table 2a) restructured to create non-overlapping daily, drug-specific exposures (Table 2b). Using the drug's metric quantity and days supplied from the original pharmacy claim (Table 2a) and the derived daily-dose estimates (Appendix 1), Table 2b shows the daily-dose was 1200 mg (prescriptions 1, 2, 6-10 and 15-17) or 800 mg for ritonavir, 80 mg for stavudine, and 300 mg for lamivudine (prescription 12). As noted in Table 2b, the subject received 100 days of ritonavir (90 days at 1200 mg and 10 days at 800 mg), 90 days of stavudine (80 mg), and 90 days of lamivudine (300 mg). The ARV drug exposure summary measures of duration of exposure, cumulative exposure, and average intensity of exposure were computed across the total prescription history (Table 2c).
Table 2 : Original pharmacy prescription claims data from the Clinformatics™ DataMart Multiplan (IMPACT) database for a single study subject (Table 2a), restructured as non-overlapping prescriptions of daily ARV drugspecific exposure (Table 2b), and as summary measures of drug exposures computed across the subject's total prescription history (Table 2c).
ARV-exposed HIV study cohort and their prescription claims
There were 30,405 subjects with HIV infection and one or
more ARV drug prescription claims. The median age was
41.9 years, 15.5% of the subjects were women, 41.4% were
from the northeastern United States (US) census region, the
median duration of observation was 2.0 years, and the majority
(91.2%) was dispensed three or more different ARV drugs
(Table 3). Most subjects (98.9%) were dispensed one or more
NRTI, 58.1% were dispensed one or more NNRTI, 55.8% were
dispensed one or more PI, and few subjects were dispensed
an FI or EI. Over 75% of the subjects had a prescription history
with at least one gap, 25.3% of subjects had less than 15 prescription claims, and 23.5% had 65 or more prescription
claims. Gaps in the prescription history were evident after the
restructure of the original claims. In the restructure of 1,929,067
ARV prescription history records, 144,258 (13%) prescription
history records had one or more days in prescription gaps.
Table 3 : Demographic and prescription-level data for 30,405 subjects with human immunodeficiency virus (HIV) infection prescribed antiretroviral (ARV) therapy between January 1, 1997 and March 31, 2008.
Within each drug class, most subjects had prescriptions with no overlap or encompassed records (Table 3). For subjects with NRTIs, 32.3% had at least one prescription that overlapped, less than 1% had at least one prescription encompassed within another, and 3% had NRTIs that both overlapped and encompassed other prescriptions. For subjects with NNRTIs, 3.5% had records that overlapped or encompassed another NNRTI. For subjects with PIs, 28.1% had overlap, less than 1% had encompassed prescriptions, and 4.4% had prescriptions which overlapped and encompassed other PIs. All subjects with FI and EI exposures had no prescriptions with overlap or encompassment.
Summarizing ARV drug prescriptions
Each of the 1,929,067 ARV prescription history records were for
a single ARV drug and associated strength as well as a start and
end date. By class, 62.1% of the records were for NRTIs, 14.3%
were for NNRTIs, 23.2% were for PIs, and less than 1% were for
FIs or EIs. Subjects' single drug prescription history records
were combined into intervals containing multiple drugs when prescription start and end dates were the same for separate
single-agent prescription history records. This compilation of
prescription records resulted in 914,821 different, time-varying
prescription history records, with 542,418 (59.3%) overlapped
and 9,379 (1.0%) encompassed prescriptions (Table 4).
Among these time-varying multidrug prescription history
records, 72.9% contained the NRTIs, 30.1% NNRTIs, 34.5% PIs,
and less than 1% FIs or EIs. Within drug classes, prescription
gaps for one or more days occurred 17-31% of the time. From
the restructured prescription history records, we calculated the summary ARV exposure metrics for the entire exposed
cohort. Distributions of the ARV drug and class-specific summary
measures are summarized in Appendix 2.
Table 4 : Description of original antiretroviral (ARV) single drug prescriptions (N=1,929,067) and time-varying multiple ARV drug prescriptions (N=914,821) among 30,405 persons with human immunodeficiency virus (HIV) infection.
Using the summary measures for ARV drug exposure
The use of the cumulative exposure metric (in milligrams
months) for fracture risk is illustrated for ritonavir, stavudine,
and lamivudine (Table 5). There were no increased risks for
fracture with ritonavir or stavudine with increasing levels of
cumulative exposure in unadjusted and adjusted models. In
contrast, lamivudine showed statistically significant reduced
risk with increasing levels of cumulative exposure in both
unadjusted and adjusted models.
Table 5 : Cumulative exposure-response findings of ritonavir, stavudine, and lamivudine for fracture in a nested case-control study of 11,621 persons with human immunodeficiency virus infection.
Assessment of risk by drug class using lagged exposures
The use of lagging with ARV drug class duration of exposure
is illustrated in Table 6. Duration of exposure (in months) for each drug class was categorized into approximate quartiles
using the unlagged exposures (zero lag). Risk estimates were
minimally different across lag periods and NRTI and NNRTI
showed reductions of risk with increased duration of exposure.
Exposure to the PI class showed a null effect that became
slightly protective with the longest duration of exposure
quartile across each lag period. No effect was noted for exposure
to the FI in a small number of cases with fracture and
there were no exposures to EIs among cases with fracture.
Table 6 : Lagged antiretroviral (ARV) drug class duration of exposures (in months) and risk for fracture among 11,621 persons with HIV infection by different lag periods.
The drug exposure summary methods in this report introduce an alternative approach for pharmacovigilance and drug safety studies when using pharmacy claims from an administrative database. Traditional analyses of pharmacoepidemiology studies typically do not use the complete, and often complex, drug histories found within administrative pharmacy claims data for patients. This may lead to incomplete descriptions of drug benefits and risks when studied using observational claims data. Historically, prescription claims in observational assessments of benefits and untoward events utilize either dichotomous or cumulative exposure to a single or few drugs, drug classes, or multidrug regimens. The restructure of the prescription claims into non-overlapping prescription histories in a defined study cohort permitted analysis of the data in ways that are not commonly executed in pharmacoepidemiology. The approaches we describe led to the quantifiable assessment of time-varying, ARV drug-specific exposureresponse relationships. These exposure-response relationships employed duration of exposure, cumulative exposure, and average intensity of exposure, as is often done in occupational epidemiological research [10-14]. Additionally, pharmacoepidemiological studies can also readily incorporate exposure lags and latency into complex exposure-response assessments and well-designed comparative effectiveness research.
These methods can be readily applied to time-varying drug exposures for chronic diseases where prescription histories are complex and evolve over time. The restructure of prescriptions into multidrug, daily-dose exposures allows capture of time-varying, multidrug prescription histories similar to occupational exposures that accumulate during an individuals' job history [9]. The restructure enabled the temporal re-alignment of same drug prescription overlaps and encompassments while giving proper attention to the gaps in the prescription history. Future investigative work can readily be expanded to include analysis of many concurrent exposures similar to analyses in occupational epidemiology [10]. Generalizing these methods for use with other pharmacy claims data-based studies of health outcomes and benefitrisk assessments is straightforward. The major objectives of the occupational studies that have employed such methods were to identify the health consequences of workplace exposure, provide data for setting standards for protection of workers to potentially toxic substances, make projections of risk to less-exposed populations, and elucidate exposureresponse relationships [11-14,20-23]. The exposures may be either continuous or intermittent and there may be several competing exposures and covariables that change over time. Phamacoepidemiology studies, particularly those involving chronic diseases that study aging, adverse drug events, cumulative comorbid conditions, and time-varying multidrug exposures have similar objectives for hypothesis generation and hypothesis testing of exposure-response relationships. To date, traditional analyses of pharmacoepidemiology studies often have not made complete use of drug exposure information contained within the complex drug histories found in pharmacy claims data. We have shown a method of restructuring the data that incorporates these complex drug histories more easily into exposure response analyses [10].
Models that include the accumulation of comorbid conditions and additional drug exposures over time will more reliably determine exposure-response estimates for untoward outcomes in chronic diseases such as HIV infection, heart disease, chronic obstructive pulmonary disease, and diabetes mellitus. We do, however, acknowledge limitations and biases associated with retrospective studies that use administrative pharmacy and medical claims to assess drugexposure response relationships. These include information bias, measurement error, and misclassification of exposures and outcome. Known limitations of measuring drug exposure that may not be documented in claims data include imperfect medication adherence by patients, stretching of drug supply, pill splitting, free samples, and double supplies to cover absences. While the methods we propose do not eliminate these limitations and biases, our methods do allow a more complete summarization of complex drug exposures over time.
We report detailed methodology for the creation of nonoverlapping, time-varying prescription histories from an administrative claims database. The adoption of occupational epidemiology methods to pharmacoepidemiology provided a quantified assessment of varied cumulative, combination ARV exposures for HIV treatment to the risk of fracture in a nested case-control study. These methods offer an alternative approach to the assessment of long-term untoward health outcomes among persons with complex drug exposures for chronic diseases.
Additional files
Appendix 1
Appendix 2
The authors declare that they have no competing interests.
Authors' contributions | AOY | LMM | SSL | MFL | SJB |
Research concept and design | √ | √ | -- | -- | √ |
Collection and/or assembly of data | √ | -- | √ | √ | -- |
Data analysis and interpretation | √ | -- | -- | √ | -- |
Writing the article | √ | √ | -- | -- | √ |
Critical revision of the article | √ | √ | √ | -- | √ |
Final approval of article | √ | √ | √ | √ | √ |
Statistical analysis | √ | -- | -- | -- | -- |
This study was supported by a sponsored research contract between the University of Pittsburgh and GlaxoSmithKline (GSK). LMM, SSL and SJB worked on this manuscript while employed with GSK. SJB is currently employed by Shire and the opinions, assertions and/or findings are those of the authors and are not to be construed as reflecting the views of Shire or any of its affiliates. GSK nor Shire had no influence in the analysis, writing or conclusions of this manuscript.
Editor: Kwok-ming Ho, Royal Perth Hospital, Australia.
Received: 27-Mar-2015 Final Revised: 14-May-2015
Accepted: 25-May-2015 Published: 03-Jun-2015
Marsh GM, Buchanich JM and Youk AO. Fiber glass exposure and human respiratory system cancer risk: lack of evidence persists since 2001 IARC re-evaluation. Regul Toxicol Pharmacol. 2011; 60:84-92. | Article | PubMed
Youk AO, Mundy LM, St. Laurent S, Lann MF and Bowlin SJ. Assessment of risk in longitudinal observational studies with time-varying drug exposures for treatment of a chronic disease. Epidemiol Rep. 2015; 3:2. http://dx.doi.org/10.7243/2054-9911-3-2
Copyright © 2015 Herbert Publications Limited. All rights reserved.