Presentations

Abstract: People living with HIV on antiretroviral therapy often have undetectable virus levels by standard assays, but latent" HIV still persists in viral reservoirs. Eliminating these reservoirs is the goal of HIV cure research. Dilution assays, including the quantitative viral outgrowth assay (QVOA) and more detailed Ultra Deep Sequencing Assay of the outgrowth virus (UDSA), are commonly used to estimate the reservoir size, i.e., the infectious units per million (IUPM) of HIV-persistent resting CD4+ T cells. This paper considers efficient statistical inference about the IUPM from combined dilution assay (QVOA) and deep viral sequencing (UDSA) data, even when some deep sequencing data are missing. Moreover, existing inference methods for the IUPM assumed that the assays are "perfect" (i.e., they have 100% sensitivity and specificity), which can be unrealistic in practice. The proposed methods accommodate assays with imperfect sensitivity and specificity, wells sequenced at multiple dilution levels, and include a novel bias-corrected estimator for small samples. The proposed methods are evaluated in a simulation study, applied to data from the University of North Carolina HIV Cure Center, and implemented in the open-source R package SLDeepAssay.

Abstract: Modeling symptom progression to identify informative subjects for a new Huntington’s disease clinical trial is problematic since time to diagnosis, a key covariate, can be heavily censored. Imputation is an appealing strategy where censored covariates are replaced with their conditional means, but existing methods saw over 200% bias under heavy censoring. Calculating these conditional means well requires estimating and then integrating over the survival function of the censored covariate from the censored value to infinity. To flexibly estimate the survival function, existing methods use the semiparametric Cox model with Breslow’s estimator. Then, for integration, the trapezoidal rule is used, but the trapezoidal rule is not designed for improper integrals and leads to bias. We propose calculating the conditional mean with adaptive quadrature instead, which can handle the improper integral. Yet, even with adaptive quadrature, the integrand (the survival function) is undefined beyond the observed data, so we identify the “Weibull extension” as the best method to extrapolate and then integrate. In simulation studies, we show that replacing the trapezoidal rule with adaptive quadrature and adopting the Weibull extension corrects the bias seen with existing methods. We further show how imputing with corrected conditional means helps to prioritize patients for future clinical trials.

Abstract: The allostatic load index (ALI) is an informative summary of whole-person health, drawing upon biomarkers to measure lifetime strain. Borrowing data from electronic health records (EHR) is a natural way to estimate whole-person health and identify at-risk patients on a large scale. However, these routinely collected data contain missingness and errors, and ignoring these data quality issues can lead to biased statistical results and incorrect clinical decisions. Validation of EHR data (e.g., through chart reviews) can provide better-quality data, but realistically, only a subset of patients' data can be validated. Thus, we devise an optimal study design to harness the error-prone surrogates from the EHR to target the most informative patient records for validation. Specifically, the optimal design seeks the best statistical precision to quantify the association between ALI and healthcare utilization in Poisson regression. Through simulations and an application to the EHR of an extensive academic learning health system, targeted partial validation is shown to be an effective and efficient way to correct data quality issues in EHR data before using them in research.

Abstract: Clinically meaningful variables are increasingly becoming available in observational databases like electronic health records (EHR). In particular,persons living with HIV engage in routine clinical care, generating especially large amounts of data in observational cohorts. However, these data can be error-prone, giving misleading results in statistical inference. Data auditing can help maintain data quality but is often unrealistic for entire databases (especially large ones like EHR). A cost-effective solution is the two-phase design: error-prone variables are observed for all patients during Phase I and that information is used to selectpatients for auditing during Phase II. However, even these partial audits can be expensive. To this end, we propose methods to promote the statistical efficiency of two-phase designs, ensuring the integrity of observational cohort data while maximizing our investment. To obtain efficient odds ratios with partially audited, error-prone data, we propose a semiparametric analysis approach that uses all information and accommodates many error mechanisms. The outcome and covariates can be error-prone, with correlated errors, and selection of Phase II records can depend on Phase I data in an arbitrary manner. We devise an EM algorithm to obtain estimators that are consistent, asymptotically normal, and asymptotically efficient. We demonstrate advantages of the proposed methods through extensive simulations and provide applications to a multi-national HIV cohort.

Abstract: Healthy foods are essential for a healthy life, but accessing healthy food can be more challenging for some people than others. With this disparity in food access comes disparities in well-being, leading to disproportionate rates of diseases in communities that face more challenges in accessing healthy food (i.e., low-access communities). Identifying low-access, high-risk communities for targeted interventions is a public health priority, but current methods to quantify food access rely on distance measures that are either computationally simple (the length of the shortest straight-line route) or accurate (the length of the shortest map-based route), but not both. We propose a hybrid statistical approach to combine these distance measures, allowing researchers to harness the computational ease of one with the accuracy of the other. The approach incorporates straight-line distances for all neighborhoods and map-based distances for just a subset, offering comparable estimates to the “gold standard” model using map-based distances for all neighborhoods and improved efficiency over the “complete case” model using map-based distances for just the subset. Through the adoption of a multiple imputation for measurement error framework, information from the straight-line distances can be used to fill in map-based measures for the remaining neighborhoods. Using data for Forsyth County, North Carolina, and its surrounding counties, we quantify and compare the associations between the prevalence of coronary heart disease and different measures of neighborhood-level food access. 

Abstract: Healthy foods are essential for a healthy life, but accessing healthy food can be more challenging for some people than others. With this disparity in food access comes disparities in well-being, leading to disproportionate rates of diseases in communities that face more challenges in accessing healthy food (i.e., low-access communities). Identifying low-access, high-risk communities for targeted interventions is a public health priority. Still, current methods to quantify food access rely on distance measures that are either computationally simple (the length of the straight-line route) or accurate (the length of the map-based route), but not both. Using data for Forsyth County, North Carolina, we first quantify and compare the associations between various health outcomes (coronary heart disease, diabetes, high blood pressure, and obesity) and neighborhood-level food access based on (i) straight-line distance versus (ii) map-based distance for all neighborhoods. We then propose a hybrid statistical approach to combine these methods, allowing researchers to harness the former’s computational ease with the latter’s accuracy. Specifically, we adopt a measurement error framework, incorporating food access based on straight-line distance for all neighborhoods (an “error-prone” covariate) and based on map-based distance for just a subset (an “error-free” covariate). By imputing map-based food access for the remaining neighborhoods, the hybrid model offers comparable estimates to the “gold standard” model with food access based on map-based distance for all neighborhoods and preserves the statistical efficiency of the whole study. We will also briefly discuss how open-source tools like the Google Maps API were used to collect data for this study, including grocery store addresses and geocoded locations. 

Abstract: Healthy foods are essential for a healthy life, but accessing healthy food can be more challenging for some people than others. With this disparity in food access comes disparities in well-being, leading to disproportionate rates of diseases in communitiesthat face more challenges in accessing healthy food (i.e., low-access communities). Identifying low-access, high-risk communities for targeted interventions is a public health priority, but current methods to quantify food access rely on distance measuresthat are either computationally simple (the length of the shortest straight-line route) or accurate (the length of the shortest map-based route), but not both. We propose a hybrid statistical approach to combine these distance measures, allowing researchers to harness the computational ease of one with the accuracy of the other. The approach incorporates straight-line distances for all neighborhoods and map-based distances for just a subset, offering comparable estimates to the “gold standard” model using map-based distances for all neighborhoods and improved efficiency over the “complete case” model using map-based distances for just the subset. Through the adoption of a multiple imputation for measurement error framework, information from the straight-line distances can be used to fill in map-based measures for the remaining neighborhoods. Using data for Forsyth County, North Carolina, and its surrounding counties, we quantify and compare the associations between various health outcomes (coronary heart disease, diabetes, high blood pressure, and obesity) and different measures of neighborhood-level food access.

Abstract: Healthy foods are essential for a healthy life, but accessing healthy food can be more challenging for some people than others. With this disparity in food access comes disparities in well-being, leading to disproportionate rates of diseases in communities that face more challenges in accessing healthy food (i.e., low-access communities). Identifying low-access, high-risk communities for targeted interventions is a public health priority, but current methods to quantify food access rely on distance t are either computationally simple (the length of the shortest straight-line route) or accurate (the length of the shortest map-based route), but not both. We propose a hybrid statistical approach to combine these distance measures, allowing researchers to harness the computational ease of one with the accuracy of the other. The approach incorporates straight-line distances for all neighborhoods and map-based distances for just a subset, offering comparable estimates to the “gold standard” model using map-based distances for all neighborhoods and improved efficiency over the “complete case” model using map-based distances for just the subset. Through the adoption of a multiple imputation for measurement error framework, information from the straight-line distances can be used to fill in map-based measures for the remaining neighborhoods. Using data for Forsyth County, NorthCarolina, and its surrounding counties, we quantify and compare the associations between health outcomes (coronary heart disease, diabetes, high blood pressure, and obesity)and different measures of neighborhood-level food access.

Abstract: The growing availability of observational databases like electronic health records (EHR) provides unprecedented opportunities for secondary use of such data in biomedical research. However, these data can be error-prone and must be validated before use. It is usually unrealistic to validate the whole database due to resource constraints. A cost-effective alternative is to implement a two-phase design that validates a subset of patient records that are enriched for information about the research question of interest. Herein, we consider odds ratio estimation under differential outcome and exposure misclassification. We propose optimal designs that minimize the variance of the maximum likelihood odds ratio estimator. We develop a novel adaptive grid search algorithm that can locate the optimal design in a computationally feasible and numerically accurate manner.  Because the optimal design requires specification of unknown parameters at the outset and thus is unattainable without prior information, we introduce a multi-wave sampling strategy to approximate it. We demonstrate the proposed designs' efficiency gains through extensive simulations and two large observational studies.

Abstract: People living with HIV on antiretroviral therapy often have undetectable virus levels by standard assays, but “latent" HIV still persists in viral reservoirs. Eliminating these reservoirs is the goal of HIV cure research. The quantitative viral outgrowth assay (QVOA) is commonly used to estimate the reservoir size, i.e., the infectious units per million (IUPM) of HIV-persistent resting CD4+ T cells. A new variation of the QVOA, the Ultra Deep Sequencing Assay of the outgrowth virus (UDSA), was recently developed that further quantifies the number of viral lineages within a subset of infected wells. Performing the UDSA on a subset of wells provides additional information that can improve IUPM estimation. This paper considers statistical inference about the IUPM from combined dilution assay (QVOA) and deep viral sequencing (UDSA) data, even when some deep sequencing data are missing. The proposed methods accommodate assays with wells sequenced at multiple dilution levels and include a novel bias-corrected estimator for small samples. The proposed methods are evaluated in a simulation study, applied to data from the University of North Carolina HIV Cure Center, and implemented in the open-source R package SLDeepAssay.

Abstract: Modeling symptom progression to prioritize subjects for a new Huntington's disease clinical trial is problematic since key covariate time to diagnosis can be censored. Imputation is an appealing strategy where censored covariates are replaced with their conditional means, but existing methods saw over 100% bias. Calculating these conditional means well requires estimating and integrating over the survival function of the censored covariate from the censored value to infinity. To flexibly estimate the survival function, existing methods use the Cox model with Breslow's estimator. Then, for integration, the trapezoidal rule, which is not designed for indefinite integrals, is used. This leads to bias. We propose a calculation that handles the indefinite integral with adaptive quadrature. Yet, even with adaptive quadrature, the integrand (the survival function) is undefined beyond the observed data. We identify the best method to extrapolate. In simulations, we show that replacing the trapezoidal rule with adaptive quadrature (plus extrapolation) corrects the bias. We further show how imputing with corrected conditional means helps prioritize patients for clinical trials.

Abstract: Clinical trials to test experimental treatments for Huntington's disease are expensive, so it is prudent to enroll subjects whose symptoms may be most impacted by the treatment during follow-up. However, modeling how symptoms progress to identify such subjects is problematic since time to diagnosis, a key predictor, can be censored. Imputation is an appealing strategy where censored predictors are replaced with their conditional means, the calculation of which requires estimating and integrating over its conditional survival function from the censored value to infinity. However, despite efforts to make conditional mean imputation as flexible as possible, it still makes restrictive assumptions about the censored predictor (such as proportional hazards) that may not hold in practice. We develop a suite of extensions to conditional mean imputation to encourage its applicability to a wide range of clinical settings. We adopt new estimators for the conditional survival function to offer more efficient and robust inference and propose an improved conditional mean calculation. We discuss in simulations when each version of conditional mean imputation is most appropriate and evaluate our methods as we model symptom progression from Huntington's disease data. Our imputation suite is implemented in the open-source R package, imputeCensRd.

Abstract: Ever since the discovery of therapies that target the genetic root cause of Huntington disease, researchers have worked to test if these therapies can slow or halt the disease symptoms. A first step towards achieving this is modeling how symptoms progress to know when the best time is to initiate a therapy. Because symptoms are most detectable before and after a clinical diagnosis, modeling how symptoms progress has been problematic since the time to clinical diagnosis is often censored (i.e., for patients who have not yet been diagnosed). This creates a pressing statistical challenge for modeling how symptoms (the outcome) change before and after time to clinical diagnosis (a censored predictor). Strategies to tackle this challenge include fitting a generalized linear model with a censored covariate using maximum likelihood estimation. Still, implementation of these models can be taxing because each new setting (i.e., different outcome models and distributions for the censored predictor) requires a new algorithm to be derived. To this end, we have created the glmCensRd package, which includes generalized linear model fitting functions for a multitude of outcome and (censored) predictor specifications and various random and non-random censoring types. The glmCensRd package makes fiKng generalized linear models in R as accessible with censored predictors as without. 

Abstract: 

Abstract: Since discovering therapies that target the genetic root of Huntington's disease, researchers have investigated whether these therapies can slow or halt symptoms and, if so when is best to initiate treatment. Symptoms are most detectable before and after clinical diagnosis, but modeling their progression is problematic since the time to clinical diagnosis is often censored. This creates a pressing statistical challenge: modeling how symptoms (the outcome) change across time to clinical diagnosis (a censored predictor). Conditional mean imputation is an appealing strategy, replacing censored times of clinical diagnoses. However, despite efforts to make conditional mean imputation flexible, it still makes restrictive assumptions (like proportional hazards) that may be unrealistic. We develop a suite of extensions to conditional mean imputation by incorporating estimators for the conditional distribution of the censored predictor to offer more efficient and robust inference. We discuss in simulations when each version of conditional mean imputation is most appropriate and evaluate our methods in Huntington's disease data.

Abstract: Since discovering therapies that target the genetic root cause of Huntington's disease, researchers have been investigating whether these therapies can slow or halt the disease symptoms. To identify the best time to initiate treatment, a first step is modeling how symptoms progress in prospective cohorts. Symptoms are most detectable before and after a clinical diagnosis, but modeling how symptoms progress has been problematic since the time to clinical diagnosis is often censored (i.e., for patients thus far undiagnosed). In contrast to traditional survival analysis, this creates a pressing statistical challenge as we model how symptoms (the outcome) change relative to clinical diagnosis (a censored covariate). Conditional mean imputation is an appealing strategy since it is an intuitive way to “fill in the gaps” in our data by estimating all censored times to clinical diagnoses. Depending on how we calculate the conditional means, imputation can offer both efficiency and robustness in statistical inference.  However, despite efforts to make conditional mean imputation flexible, it still makes restrictive assumptions (such as proportional hazards) that may not hold in practice. We set out to develop a suite of extensions that could offer even more robustness and efficiency in a wide range of clinical settings. Along the way, we discovered that the published conditional mean formula itself needed some improvements.  Through extensive simulations, we quantify the bias introduced by the published formula and show that our proposed method does a good job of correcting for it. We also apply our methods to Huntington's disease study data to evaluate how imputing with biased conditional means can impact clinical decision-making.

Abstract: The growing availability of observational databases like electronic health records (EHR) provides unprecedented opportunities for secondary use of such data in biomedical research. However, these data can be error-prone and need to be validated before use. It is usually unrealistic to validate the whole database due to resource constraints. A cost-effective alternative is to implement a two-phase design that validates a subset of patient records that are enriched for information about the research question of interest. In this talk, I will discuss proper statistical approaches to analyze such two-phase studies, which can efficiently use the information in the unvalidated data in Phase I and address the potential biased validation sample selection in Phase II. I will demonstrate the advantages of the proposed methods over existing ones through extensive simulations and an application to an ongoing HIV observational study.

Abstract: Since discovering therapies that target the genetic root of Huntington's disease, researchers have investigated whether these therapies can slow or halt symptoms and, if so, when the best time is to initiate treatment. Symptoms are most detectable before and after clinical diagnosis, but modeling their progression is problematic since the time to clinical diagnosis is often censored. This creates a pressing statistical challenge: modeling how symptoms (the outcome) change across time to clinical diagnosis (a censored predictor). Conditional mean imputation is an appealing strategy since it is a simple way to estimate all times of clinical diagnoses that are censored. However, despite efforts to make conditional mean imputation as flexible as possible, it still makes restrictive assumptions about the censored predictor (such as proportional hazards) that may not hold in practice. We develop a suite of extensions to conditional mean imputation to encourage its applicability to a wide range of clinical settings. We incorporate additional estimators for the conditional distribution of the censored predictor to offer more efficient and robust inference and propose an improved conditional mean calculation. Through extensive simulations, we discuss when each version of conditional mean imputation is most appropriate.

Abstract: Since discovering therapies that target the genetic root of Huntington's disease, researchers have investigated whether these therapies can slow or halt symptoms and, if so, when the best time is to initiate treatment. Symptoms are most detectable before and after clinical diagnosis, but modeling their progression is problematic since the time to clinical diagnosis is often censored. This creates a pressing statistical challenge: modeling how symptoms (the outcome) change across time to clinical diagnosis (a censored predictor). Conditional mean imputation is an appealing strategy since it is a simple way to estimate all times of clinical diagnoses that are censored. However, despite efforts to make conditional mean imputation as flexible as possible, it still makes restrictive assumptions about the censored predictor (such as proportional hazards) that may not hold in practice. We develop a suite of extensions to conditional mean imputation to encourage its applicability to a wide range of clinical settings. We incorporate additional estimators for the conditional distribution of the censored predictor to offer more efficient and robust inference and propose an improved conditional mean calculation. Through extensive simulations, we discuss when each version of conditional mean imputation is most appropriate.

Abstract: The growing availability of observational databases like electronic health records (EHR) provides unprecedented opportunities for secondary use of such data in biomedical research. However, these data can be error-prone and need to be validated before use. It is usually unrealistic to validate the whole database due to resource constraints. A cost-effective alternative is to implement a two-phase design that validates a subset of patient records that are enriched for information about the research question of interest. Herein, we consider odds ratio estimation under differential outcome and exposure misclassification. We propose optimal designs that minimize the variance of the maximum likelihood odds ratio estimator. We develop a novel adaptive grid search algorithm that can locate the optimal design in computationally feasible and accurate manner. Because the optimal design requires specification of unknown parameters at the outset, and thus is unattainable without prior information, we introduce a multi-wave sampling strategy to approximate it in practice. We demonstrate the efficiency gains of the proposed designs over existing ones through extensive simulations and two large observational studies. We provide an R package and Shiny app to facilitate the use of the optimal designs.

Abstract: Clinically meaningful variables are increasingly becoming available in observational databases like electronic health records (EHR). However, these data can be error-prone, giving misleading results in statistical inference. Data auditing can help maintain data quality but is often unrealistic for entire databases (especially large ones like EHR). A cost-effective solution is the two-phase design: error-prone variables are observed for all patients during Phase I and that information is used to select patients for auditing during Phase II. However, even these partial audits can be expensive. To this end, we propose methods to promote the statistical efficiency of two-phase designs, ensuring the integrity of observational cohort data while maximizing our investment. First, given the resource constraints imposed upon data audits, targeting the most informative patients is paramount for efficient statistical inference. Using the asymptotic variance of the maximum likelihood estimator, we compute the most efficient design under complex outcome and exposure misclassification. Since the optimal design depends on unknown parameters, we propose a multi-wave design to approximate it in practice. We demonstrate the superior efficiency of the optimal designs through extensive simulations and illustrate their implementation in observational HIV studies. Then, to obtain efficient odds ratios with partially audited, error-prone data, we propose a semiparametric analysis approach that uses all information and accommodates many error mechanisms. The outcome and covariates can be error-prone, with correlated errors, and selection of Phase II records can depend on Phase I data in an arbitrary manner. We devise an EM algorithm to obtain estimators that are consistent, asymptotically normal, and asymptotically efficient. We demonstrate advantages of the proposed methods through extensive simulations and provide applications to a multi-national HIV cohort.

Abstract: Persons living with HIV engage in clinical care often, so observational HIV research cohorts generate especially large amounts of routine clinical data. Increasingly, these data are being used in biomedical research, but available information can be error-prone, and biased statistical estimates can mislead results. The Caribbean, Central, and South America network for HIV epidemiology is one such cohort; fortunately, data audits have been conducted. Risk of AIDS defining event after initiating antiretroviral therapy is of clinical interest, expected to be associated with CD4 lab value and AIDS status. Error-prone values for 5109 patients were in the research database, and validated data were available (substantiated by clinical source documents) on only 117 patients. Instead of naive (unaudited) or complete case (audited) analysis, we propose a novel semi-parametric likelihood method using all available information (unau- dited and audited) to obtain unbiased, efficient odds ratios with error-prone outcome and covariates. Point estimates were farther from the null than the naive analysis, directionality agreed with the complete case analysis but had narrower confidence intervals.

Abstract: Persons living with HIV engage in clinical care often, so observational HIV research cohorts generate especially large amounts of routine clinical data. Increasingly, these data are being used in biomedical research, but available information can be error prone and biased statistical estimates can mislead results. The Caribbean, Central, and South America network for HIV epidemiology is one such cohort; fortunately, data audits have been conducted. Risk of AIDS defining event after initiating antiretroviral therapy is of clinical interest, expected to be associated with CD4 lab value and AIDS status. Error-prone values for 5109 patients were in the research database, and validated data were available (substantiated by clinical source documents) on only 117 patients. Instead of naive (unaudited) or complete case (audited) analysis, we propose a novel semiparametric likelihood method using all available information (unaudited and audited) to obtain unbiased, efficient odds ratios with error prone outcome and covariates. Point estimates were farther from the null than the naive analysis, directionality agreed with the complete case analysis, but had narrower confidence intervals.  

Abstract: Persons living with HIV engage in routine clinical care, generating large amounts of data in observational cohorts. These data are usually error-prone, and directly using them in biomedical research can yield misleading results. A cost-effective solution is the two-phase design, under which the error-prone variables are observed for all patients during Phase I and that information is used to select patients for data auditing during Phase II. For example, the Caribbean, Central, and South America Network for HIV Epidemiology (CCASAnet) selected random samples from each site for data auditing. We consider efficient odds ratio estimation with partially audited, error-prone data and propose a semiparametric approach that accommodates a number of error mechanisms. We allow both the outcome and covariates to be error-prone and errors to be correlated, and selection of Phase II can depend on Phase I data in an arbitrary manner. We devise an EM algorithm to obtain estimators that are consistent, asymptotically normal, and asymptotically efficient. We demonstrate the superiority of the proposed methods over existing ones through extensive simulations and provide applications to CCASAnet.

Abstract: While electronic health records were first intended to support clinical care, financial billing, and insurance claims, these databases are now used for clinical investigations aimed at preventing disease, improving patient care, and informing policymaking. However, both responses and predictors of interest can be captured with errors and their discrepancies correlated. Odds ratios and their standard errors estimated via logistic regression using error-prone data will be biased. A cost-effective solution to a complete data audit is a two-phase design. During Phase I error-prone variables are observed for all subjects, and this information then used to select a Phase II validation subsample. Previous approaches to outcome misclassification using two-phase design data are limited to error-prone categorical predictors and make distributional assumptions about the errors. We propose a semiparametric approach to two-phase designs with a misclassified, binary outcome and categorical or continuous error-prone predictors, allowing for dependent errors and arbitrary second-phase selection. The proposed method is robust because it yields consistent estimates without making assumptions about the predictors’ error mechanisms. An EM algorithm was devised to maximize the likelihood function. The resulting estimators possess desired statistical properties. Performance is compared to existing approaches through extensive simulation studies and illustrated in an observational HIV study.

Abstract: Modern electronic health records systems routinely collect data on variables of clinical interest. In epidemiology, logistic regression can be used to analyze the association between a binary response (e.g. disease status) and predictors of interest (e.g. age or exposure status). However, in observational databases responses and predictors can be differentially misclassified, and their misclassification can be dependent. A cost-effective solution to a complete data audit is the two-phase design. During the first phase, error-prone outcome and covariates are observed for all subjects and this information is then used to select a validation subsample for accurate measurements of these variables in the second phase. Previous measurement error corrections focused on error-prone covariates only or rely on a validation subsample that is simple or stratified. Herein, we propose a semiparametric approach to general two-phase measurement error problems with a binary outcome, allowing for correlated errors in the outcome and covariates and arbitrary second-phase selection. We devise a computationally efficient and numerically stable EM algorithm to maximize the nonparametric likelihood function. The resulting estimators possess desired statistical properties. We compare the proposed methods to existing approaches through extensive simulation studies, and we illustrate their use in an observational HIV study.

Abstract: In modern electronic health records systems, both the outcome and covariates of interest can be error-prone and their errors often correlated. A cost-effective solution is the two-phase design, under which the error-prone outcome and covariates are observed for all subjects during the first phase and that information is used to select a validation subsample for accurate measurements of these variables in the second phase. Previous research on two-phase measurement error problems largely focused on scenarios where there are errors in covariates only or the validation sample is a simple random sample. Herein, we propose a semiparametric approach to general two-phase measurement error problems with a continuous outcome, allowing for correlated errors in the outcome and covariates and arbitrary second-phase selection. We devise a computationally efficient and numerically stable EM algorithm to maximize the nonparametric likelihood function. The resulting estimators possess desired statistical properties. We compare the proposed methods to existing approaches through extensive simulation studies, and we illustrate their use in an observational HIV study.

Abstract: The CDC publishes the CHSI dataset, providing county-level key health indicators for the United States. The current analysis focuses on the density of primary care physicians per 100,000 people as a gauge of healthcare access in a subset of 878 counties in Tennessee. Access to primary care may be measured by travel distance or travel time. We investigate the implications of travel time and two distance measures: travel distance and travel time using the Google Maps API and shortest distance ("as the crow flies") using the Haversine formula. These measures were used to create empirically estimated (classical) semivariograms. We utilized ordinary kriging methods (OK) to interpolate the density of primary care physicians in the areas between the geographic centroids of the counties (three analyses total). We will compare kriged results between the three methods, illustrate results with maps, and explore the relationships between travel distance/time and Haversine distance.