Tree-based subgroup discovery using electronic health record data: heterogeneity of treatment effects for DTG-containing therapies.

The rich longitudinal individual level data available from electronic health records (EHRs) can be used to examine treatment effect heterogeneity. However, estimating treatment effects using EHR data poses several challenges, including time-varying confounding, repeated and temporally non-aligned measurements of covariates, treatment assignments and outcomes, and loss-to-follow-up due to dropout. Here, we develop the subgroup discovery for longitudinal data algorithm, a tree-based algorithm for discovering subgroups with heterogeneous treatment effects using longitudinal data by combining the generalized interaction tree algorithm, a general data-driven method for subgroup discovery, with longitudinal targeted maximum likelihood estimation. We apply the algorithm to EHR data to discover subgroups of people living with human immunodeficiency virus who are at higher risk of weight gain when receiving dolutegravir (DTG)-containing antiretroviral therapies (ARTs) versus when receiving non-DTG-containing ARTs.

Validity of Self-Report for Ascertaining HIV Status Among Circular Migrants and Permanent Residents in South Africa: A Cross-Sectional, Population-Based Analysis.

While expanded HIV testing is needed in South Africa, increasing accurate self-report of HIV status is an essential parallel goal in this highly mobile population. If self-report can ascertain true HIV-positive status, persons with HIV (PWH) could be linked to life-saving care without the existing delays required by producing medical records or undergoing confirmatory testing, which are especially burdensome for the country's high prevalence of circular migrants. We used Wave 1 data from The Migration and Health Follow-Up Study, a representative adult cohort, including circular migrants and permanent residents, randomly sampled from the Agincourt Health and Demographic Surveillance System in a rural area of Mpumalanga Province. Within the analytic sample (n = 1,918), sensitivity, specificity, positive predictive value (PPV), and negative predictive value (NPV) of self-report were calculated with dried blood spot (DBS) HIV test results as the standard. Among in-person participants (n = 2,468), 88.8% consented to DBS-HIV testing. HIV prevalence was 25.3%. Sensitivity of self-report was 43.9% (95% CI: 39.5-48.5), PPV was 93.4% (95% CI: 89.5-96.0); specificity was 99.0% (95% CI: 98.3-99.4) and NPV was 83.9% (95% CI: 82.8-84.9). Self-report of an HIV-positive status was predictive of true status for both migrants and permanent residents in this high-prevalence setting. Persons who self-reported as living with HIV were almost always truly positive, supporting a change to clinical protocol to immediately connect persons who say they are HIV-positive to ART and counselling. However, 56% of PWH did not report as HIV-positive, highlighting the imperative to address barriers to disclosure.

A flexible approach for causal inference with multiple treatments and clustered survival outcomes.

When drawing causal inferences about the effects of multiple treatments on clustered survival outcomes using observational data, we need to address implications of the multilevel data structure, multiple treatments, censoring, and unmeasured confounding for causal analyses. Few off-the-shelf causal inference tools are available to simultaneously tackle these issues. We develop a flexible random-intercept accelerated failure time model, in which we use Bayesian additive regression trees to capture arbitrarily complex relationships between censored survival times and pre-treatment covariates and use the random intercepts to capture cluster-specific main effects. We develop an efficient Markov chain Monte Carlo algorithm to draw posterior inferences about the population survival effects of multiple treatments and examine the variability in cluster-level effects. We further propose an interpretable sensitivity analysis approach to evaluate the sensitivity of drawn causal inferences about treatment effect to the potential magnitude of departure from the causal assumption of no unmeasured confounding. Expansive simulations empirically validate and demonstrate good practical operating characteristics of our proposed methods. Applying the proposed methods to a dataset on older high-risk localized prostate cancer patients drawn from the National Cancer Database, we evaluate the comparative effects of three treatment approaches on patient survival, and assess the ramifications of potential unmeasured confounding. The methods developed in this work are readily available in the R $$ \mathsf{R}\kern.15em $$ package riAFTBART $$ \mathsf{riAFTBART} $$ .

Seroprevalence of SARS-CoV-2 Antibodies in Rhode Island From a Statewide Random Sample.

Objectives. To characterize statewide seroprevalence and point prevalence of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) in Rhode Island.Methods. We conducted a cross-sectional survey of randomly selected households across Rhode Island in May 2020. Antibody-based and polymerase chain reaction (PCR)-based tests for SARS-CoV-2 were offered. Hispanics/Latinos and African Americans/Blacks were oversampled to ensure adequate representation. Seroprevalence estimations accounted for test sensitivity and specificity and were compared according to age, race/ethnicity, gender, housing environment, and transportation mode.Results. Overall, 1043 individuals from 554 households were tested (1032 antibody tests, 988 PCR tests). The estimated seroprevalence of SARS-CoV-2 antibodies was 2.1% (95% credible interval [CI] = 0.6, 4.1). Seroprevalence was 7.5% (95% CI = 1.3, 17.5) among Hispanics/Latinos, 3.8% (95% CI = 0.0, 15.0) among African Americans/Blacks, and 0.8% (95% CI = 0.0, 2.4) among non-Hispanic Whites. Overall PCR-based prevalence was 1.5% (95% CI = 0.5, 3.1).Conclusions. Rhode Island had low seroprevalence relative to other settings, but seroprevalence was substantially higher among African Americans/Blacks and Hispanics/Latinos. Rhode Island sits along the highly populated northeast corridor, making our findings broadly relevant to this region of the country. Continued monitoring via population-based sampling is needed to quantify these impacts going forward.

Modeling the Impact of Social Determinants of Health on HIV.

There is growing evidence for the key role of social determinants of health (SDOH) in understanding morbidity and mortality outcomes globally. Factors such as stigma, racism, poverty or access to health and social services represent complex constructs that affect population health via intricate relationships to individual characteristics, behaviors and disease prevention and treatment outcomes. Modeling the role of SDOH is both critically important and inherently complex. Here we describe different modeling approaches and their use in assessing the impact of SDOH on HIV/AIDS. The discussion is thematically divided into mechanistic models and statistical models, while recognizing the overlap between them. To illustrate mechanistic approaches, we use examples of compartmental models and agent-based models; to illustrate statistical approaches, we use regression and statistical causal models. We describe model structure, data sources required, and the scope of possible inferences, highlighting similarities and differences in formulation, implementation, and interpretation of different modeling approaches. We also indicate further needed research on representing and quantifying the effect of SDOH in the context of models for HIV and other health outcomes in recognition of the critical role of SDOH in achieving the goal of ending the HIV epidemic and improving overall population health.

Low Seroprevalence of SARS-CoV-2 in Rhode Island blood donors during may 2020 as determined using multiple serological assay formats.

BACKGROUND: Epidemic projections and public health policies addressing Coronavirus disease (COVID)-19 have been implemented without data reporting on the seroconversion of the population since scalable antibody testing has only recently become available. METHODS: We measured the percentage of severe acute respiratory syndrome- Coronavirus-2 (SARS-CoV-2) seropositive individuals from 2008 blood donors drawn in the state of Rhode Island (RI). We utilized multiple antibody testing platforms, including lateral flow immunoassays (LFAs), enzyme-linked immunosorbent assays (ELISAs) and high throughput serological assays (HTSAs). To estimate seroprevalence, we utilized the Bayesian statistical method to adjust for sensitivity and specificity of the commercial tests used. RESULTS: We report than an estimated seropositive rate of RI blood donors of approximately 0.6% existed in April-May of 2020. Daily new case rates peaked in RI in late April 2020. We found HTSAs and LFAs were positively correlated with ELISA assays to detect antibodies specific to SARS-CoV-2 in blood donors. CONCLUSIONS: These data imply that seroconversion, and thus infection, is likely not widespread within this population. We conclude that IgG LFAs and HTSAs are suitable to conduct seroprevalence assays in random populations. More studies will be needed using validated serological tests to improve the precision and report the kinetic progression of seroprevalence estimates.

Fixed choice design and augmented fixed choice design for network data with missing observations.

The statistical analysis of social networks is increasingly used to understand social processes and patterns. The association between social relationships and individual behaviors is of particular interest to sociologists, psychologists, and public health researchers. Several recent network studies make use of the fixed choice design (FCD), which induces missing edges in the network data. Because of the complex dependence structure inherent in networks, missing data can pose very difficult problems for valid statistical inference. In this article, we introduce novel methods for accounting for the FCD censoring and introduce a new survey design, which we call the augmented fixed choice design (AFCD). The AFCD adds considerable information to analyses without unduly burdening the survey respondent, resulting in improvements over the FCD, and other existing estimators. We demonstrate this new method through simulation studies and an analysis of alcohol use in a network of undergraduate students living in a residence hall.

Causal comparative effectiveness analysis of dynamic continuous-time treatment initiation rules with sparsely measured outcomes and death.

Evidence supporting the current World Health Organization recommendations of early antiretroviral therapy (ART) initiation for adolescents is inconclusive. We leverage a large observational data and compare, in terms of mortality and CD4 cell count, the dynamic treatment initiation rules for human immunodeficiency virus-infected adolescents. Our approaches extend the marginal structural model for estimating outcome distributions under dynamic treatment regimes, developed in Robins et al. (2008), to allow the causal comparisons of both specific regimes and regimes along a continuum. Furthermore, we propose strategies to address three challenges posed by the complex data set: continuous-time measurement of the treatment initiation process; sparse measurement of longitudinal outcomes of interest, leading to incomplete data; and censoring due to dropout and death. We derive a weighting strategy for continuous-time treatment initiation, use imputation to deal with missingness caused by sparse measurements and dropout, and define a composite outcome that incorporates both death and CD4 count as a basis for comparing treatment regimes. Our analysis suggests that immediate ART initiation leads to lower mortality and higher median values of the composite outcome, relative to other initiation rules.

Classification using ensemble learning under weighted misclassification loss.

Binary classification rules based on covariates typically depend on simple loss functions such as zero-one misclassification. Some cases may require more complex loss functions. For example, individual-level monitoring of HIV-infected individuals on antiretroviral therapy requires periodic assessment of treatment failure, defined as having a viral load (VL) value above a certain threshold. In some resource limited settings, VL tests may be limited by cost or technology, and diagnoses are based on other clinical markers. Depending on scenario, higher premium may be placed on avoiding false-positives, which brings greater cost and reduced treatment options. Here, the optimal rule is determined by minimizing a weighted misclassification loss/risk. We propose a method for finding and cross-validating optimal binary classification rules under weighted misclassification loss. We focus on rules comprising a prediction score and an associated threshold, where the score is derived using an ensemble learner. Simulations and examples show that our method, which derives the score and threshold jointly, more accurately estimates overall risk and has better operating characteristics compared with methods that derive the score first and the cutoff conditionally on the score especially for finite samples.

Evaluating the Population Impact on Racial/Ethnic Disparities in HIV in Adulthood of Intervening on Specific Targets: A Conceptual and Methodological Framework.

Reducing racial/ethnic disparities in human immunodeficiency virus (HIV) disease is a high priority. Reductions in HIV racial/ethnic disparities can potentially be achieved by intervening on important intermediate factors. The potential population impact of intervening on intermediates can be evaluated using observational data when certain conditions are met. However, using standard stratification-based approaches commonly employed in the observational HIV literature to estimate the potential population impact in this setting may yield results that do not accurately estimate quantities of interest. Here we describe a useful conceptual and methodological framework for using observational data to appropriately evaluate the impact on HIV racial/ethnic disparities of interventions. This framework reframes relevant scientific questions in terms of a controlled direct effect and estimates a corresponding proportion eliminated. We review methods and conditions sufficient for accurate estimation within the proposed framework. We use the framework to analyze data on 2,329 participants in the CFAR [Centers for AIDS Research] Network of Integrated Clinical Systems (2008-2014) to evaluate the potential impact of universal prescription of and ≥95% adherence to antiretroviral therapy on racial disparities in HIV virological suppression. We encourage the use of the described framework to appropriately evaluate the potential impact of targeted interventions in addressing HIV racial/ethnic disparities using observational data.

Modeling the causal effect of treatment initiation time on survival: Application to HIV/TB co-infection.

The timing of antiretroviral therapy (ART) initiation for HIV and tuberculosis (TB) co-infected patients needs to be considered carefully. CD4 cell count can be used to guide decision making about when to initiate ART. Evidence from recent randomized trials and observational studies generally supports early initiation but does not provide information about effects of initiation time on a continuous scale. In this article, we develop and apply a highly flexible structural proportional hazards model for characterizing the effect of treatment initiation time on a survival distribution. The model can be fitted using a weighted partial likelihood score function. Construction of both the score function and the weights must accommodate censoring of the treatment initiation time, the outcome, or both. The methods are applied to data on 4903 individuals with HIV/TB co-infection, derived from electronic health records in a large HIV care program in Kenya. We use a model formulation that flexibly captures the joint effects of ART initiation time and ART duration using natural cubic splines. The model is used to generate survival curves corresponding to specific treatment initiation times; and to identify optimal times for ART initiation for subgroups defined by CD4 count at time of TB diagnosis. Our findings potentially provide 'higher resolution' information about the relationship between ART timing and mortality, and about the differential effect of ART timing within CD4 subgroups.

A state transition framework for patient-level modeling of engagement and retention in HIV care using longitudinal cohort data.

The human immunodeficiency virus (HIV) care cascade is a conceptual model used to outline the benchmarks that reflects effectiveness of HIV care in the whole HIV care continuum. The models can be used to identify barriers contributing to poor outcomes along each benchmark in the cascade such as disengagement from care or death. Recently, the HIV care cascade has been widely applied to monitor progress towards HIV prevention and care goals in an attempt to develop strategies to improve health outcomes along the care continuum. Yet, there are challenges in quantifying successes and gaps in HIV care using the cascade models that are partly due to the lack of analytic approaches. The availability of large cohort data presents an opportunity to develop a coherent statistical framework for analysis of the HIV care cascade. Motivated by data from the Academic Model Providing Access to Healthcare, which has provided HIV care to nearly 200,000 individuals in Western Kenya since 2001, we developed a state transition framework that can characterize patient-level movements through the multiple stages of the HIV care cascade. We describe how to transform large observational data into an analyzable format. We then illustrate the state transition framework via multistate modeling to quantify dynamics in retention aspects of care. The proposed modeling approach identifies the transition probabilities of moving through each stage in the care cascade. In addition, this approach allows regression-based estimation to characterize effects of (time-varying) predictors of within and between state transitions such as retention, disengagement, re-entry into care, transfer-out, and mortality. Copyright © 2017 John Wiley & Sons, Ltd.

Resilience and HIV: a review of the definition and study of resilience.

We use a socioecological model of health to define resilience, review the definition and study of resilience among persons living with human immunodeficiency virus (PLWH) in the existing peer-reviewed literature, and discuss the strengths and limitations of how resilience is defined and studied in HIV research. We conducted a review of resilience research for HIV-related behaviors/outcomes of antiretroviral therapy (ART) adherence, clinic attendance, CD4 cell count, viral load, viral suppression, and/or immune functioning among PLWH. We performed searches using PubMed, PsycINFO and Google Scholar databases. The initial search generated 14,296 articles across the three databases, but based on our screening of these articles and inclusion criteria, n = 54 articles were included for review. The majority of HIV resilience research defines resilience only at the individual (i.e., psychological) level or studies individual and limited interpersonal resilience (e.g., social support). Furthermore, the preponderance of HIV resilience research uses general measures of resilience; these measures have not been developed with or tailored to the needs of PLWH. Our review suggests that a socioecological model of health approach can more fully represent the construct of resilience. Furthermore, measures specific to PLWH that capture individual, interpersonal, and neighborhood resilience are needed.