Covariate measurement error adjustment for matched case-control studies

We propose a conditional scores procedure for obtaining bias-corrected estimates of log odds ratios from matched case-control data in which one or more covariates are subject to measurement error. The approach involves conditioning on sufficient statistics for the unobservable true covariates that are treated as fixed unknown parameters. For the case of Gaussian nondifferential measurement error, we derive a set of unbiased score equations that can then be solved to estimate the log odds ratio parameters of interest. The procedure successfully removes the bias in naive estimates, and standard error estimates are obtained by resampling methods. We present an example of the procedure applied to data from a matched case-control study of prostate cancer and serum hormone levels, and we compare its performance to that of regression calibration procedures.     

Analysis of case-control age-at-onset data using a modified case-cohort method

Case-control designs are widely used in rare disease studies. In a typical case-control study, data are collected from a sample of all available subjects who have experienced a disease (cases) and a sub-sample of subjects who have not experienced the disease (controls) in a study cohort. Cases are oversampled in case-control studies. Logistic regression is a common tool to estimate the relative risks of the disease with respect to a set of covariates. Very often in such a study, information of ages-at-onset of the disease for all cases and ages at survey of controls are known. Standard logistic regression analysis using age as a covariate is based on a dichotomous outcome and does not efficiently use such age-at-onset (time-to-event) information. We propose to analyze age-at-onset data using a modified case-cohort method by treating the control group as an approximation of a subcohort assuming rare events. We investigate the asymptotic bias of this approximation and show that the asymptotic bias of the proposed estimator is small when the disease rate is low. We evaluate the finite sample performance of the proposed method through a simulation study and illustrate the method using a breast cancer case-control data set.     

Comparison of racial differences in childhood cancer risk in case-control studies and population-based cancer registries

Introduction: Although selection bias in case-control studies has been studied extensively, little is known about selection of cases and controls among various ethnic groups. This study compares racial differences in childhood cancer rates as estimated by case-control studies with various design features. It also compares estimates of racial distribution among cases as reported by case-control studies to those observed for an ideal case series with complete ascertainment of cases for these studies or in population-based cancer registries in corresponding geographic regions and calendar periods. Methods: Peer-reviewed publications on childhood leukemia and brain tumors from North America, published between 1980 and 2007, were reviewed. Incidence data by race/ethnicity were compiled from research publications, federal cancer statistics, and cancer registries. Meta-analysis was conducted to assess racial/ethnic differences by study characteristics. Racial distributions of cases from published case-control studies were compared to those of a presumably noncensored case distribution (i.e. include both participating and non-participating cases in a case-control study) or cases recorded by cancer registries. Results: In interview-based case-control studies of childhood cancer, the proportion of Whites compared to non-Whites tended to be higher among controls than among cases; however, the opposite was true for record-based case-control studies. Additionally, the proportion of Whites tended to be higher among the participating cases in the published case-control studies compared to the proportion of Whites among the non-participating cases or in cancer registries. Conclusions: Investigators need to consider differential participation by racial group as a potential source of bias in the interpretation of case-control study results.     

Ascertainment bias in rate ratio estimation from case-sibling control studies of variable age-at-onset diseases

Motivated by a Finnish case-control study of early onset diabetes in which diabetic children are matched to sibling controls, we investigate ascertainment bias of the usual rate ratio estimator from case-control data under simplex complete ascertainment of families during a fixed interval of time. Analytic results indicate that the assumptions necessary for valid estimation are that the disease is rare and the factors under study are exchangeable--essentially that the covariate distribution does not depend on calendar time or birth order. Further, we found that the rare disease assumption could be dropped by restricting to cases that were diagnosed during the enrollment period of the study or including all cases but eliminating the proband as a control for non-enrollment-period cases. An important consequence of this work is that standard family-based case-control studies are subject to ascertainment bias if exchangeability of the covariates under investigation does not hold.     

On the design of synthetic case-control studies

A design is proposed for "case-control within cohort" studies. In this design, controls are sampled without replacement from failure-free members of the cohort at each distinct failure time. Upon selection, a subject ceases to be eligible for control selection at later failure times. Also, if a subject failing at time t had been selected as a control at t' less than t, then the matched controls at t are selected to have also been at risk at t'. In these circumstances correlation exists between score statistic contributions at t and t'. An estimator is developed for this correlation. A small simulation study compares the design just described to other possible synthetic case-control designs.     

Improving sandwich variance estimation for marginal Cox analysis of cluster randomized trials

Cluster randomized trials (CRTs) frequently recruit a small number of clusters, therefore necessitating the application of small-sample corrections for valid inference. A recent systematic review indicated that CRTs reporting right-censored, time-to-event outcomes are not uncommon and that the marginal Cox proportional hazards model is one of the common approaches used for primary analysis. While small-sample corrections have been studied under marginal models with continuous, binary, and count outcomes, no prior research has been devoted to the development and evaluation of bias-corrected sandwich variance estimators when clustered time-to-event outcomes are analyzed by the marginal Cox model. To improve current practice, we propose nine bias-corrected sandwich variance estimators for the analysis of CRTs using the marginal Cox model and report on a simulation study to evaluate their small-sample properties. Our results indicate that the optimal choice of bias-corrected sandwich variance estimator for CRTs with survival outcomes can depend on the variability of cluster sizes and can also slightly differ whether it is evaluated according to relative bias or type I error rate. Finally, we illustrate the new variance estimators in a real-world CRT where the conclusion about intervention effectiveness differs depending on the use of small-sample bias corrections. The proposed sandwich variance estimators are implemented in an R package CoxBcv .     

Evaluating the use of friend or family controls in epidemiologic case-control studies

BackgroundTraditional methodologies for identifying and recruiting controls in epidemiologic case-control studies, such as random digit dialing or neighborhood walk, suffer from declining response rates. Here, we revisit the feasibility and comparability of using alternative sources of controls, specifically friend and family controls.MethodsWe recruited from a recently completed case-control study of non-Hodgkin lymphoma (NHL) among women in Los Angeles County where controls from the parent study were ascertained by neighborhood walk. We calculated participation rates and compared questionnaire responses between the friend/family controls and the original matched controls from the parent study.ResultsOf the 182 NHL case patients contacted, 111 (61%) agreed to participate in our feasibility study. 70 (63%) provided contact information for potential friend and/or family member controls. We were able to successfully contact and recruit a friend/family member for 92% of the case patients. This represented 46 friend controls and 54 family controls. Family controls significantly differed from original matched controls by sex and household income. Other characteristics were similar between friend controls and the original study’s neighborhood controls.ConclusionThe apparent comparability of neighborhood controls to friend and family controls among respondents in this study suggests that these alternative methods of control identification can serve as a complementary source of eligible controls in epidemiologic case-control studies.     

Impact of random and systematic recall errors and selection bias in case--control studies on mobile phone use and brain tumors in adolescents (CEFALO study)

Whether the use of mobile phones is a risk factor for brain tumors in adolescents is currently being studied. Case--control studies investigating this possible relationship are prone to recall error and selection bias. We assessed the potential impact of random and systematic recall error and selection bias on odds ratios (ORs) by performing simulations based on real data from an ongoing case--control study of mobile phones and brain tumor risk in children and adolescents (CEFALO study). Simulations were conducted for two mobile phone exposure categories: regular and heavy use. Our choice of levels of recall error was guided by a validation study that compared objective network operator data with the self-reported amount of mobile phone use in CEFALO. In our validation study, cases overestimated their number of calls by 9% on average and controls by 34%. Cases also overestimated their duration of calls by 52% on average and controls by 163%. The participation rates in CEFALO were 83% for cases and 71% for controls. In a variety of scenarios, the combined impact of recall error and selection bias on the estimated ORs was complex. These simulations are useful for the interpretation of previous case-control studies on brain tumor and mobile phone use in adults as well as for the interpretation of future studies on adolescents.     

Inverse sampling of controls in a matched case control study

A method of inverse sampling of controls in a matched case-control study is described in which, for each case, controls are sampled until a discordant set is achieved. For a binary exposure, inverse sampling is used to determine the number of controls for each case. When most individuals in a population have the same exposure, standard case-control sampling may result in many case-control sets being concordant with respect to exposure and thus uninformative in the conditional logistic analysis. The method using inverse control sampling is proposed as a solution to this problem in situations when it is practically feasible. In many circumstances, inverse control sampling is found to offer improved statistical efficiency relative to a comparable study with a fixed number of controls per case.     

Attributable risk estimation from matched case-control data

A methodology is proposed for obtaining summary estimators, variances, and confidence intervals for attributable risk measures from data obtained through a case-control study design where one or more controls have been matched to each case. The sampling design for obtaining these data is conceptualized as a simple random sample of cases being equivalent to a simple random sample of matched sets. By combining information across the strata determined by the matched sets, this approach provides all of the benefits associated with the Mantel-Haenszel procedure for the estimators of attributable risk among the exposed and population attributable risk. Asymptotic variances are derived under the assumption that the frequencies of the unique response patterns follow the multinomial distribution. Simulation results indicate that these methods fare very well with respect to bias and coverage probability.     

Sample size for individually matched case-control studies

The standard formulas used to calculate sample size for an individually matched case-control study assume a constant probability of exposure throughout the pool of possible controls. We propose new formulas that allow for heterogeneity in the probability of exposure among controls in different matched sets. Since matching factors are suspected of being confounders, they are expected to divide the total population into subgroups with different proportions exposed. Thus, the assumption of homogeneity of exposure among controls, made by the currently used formulas, is inconsistent with the assumptions used to design a matched study. The proposed formulas avoid this inconsistency. We present an example to illustrate how heterogeneity can affect the required sample size.     

A flexible approach to measurement error correction in case-control studies

SUMMARY: We investigate the use of prospective likelihood methods to analyze retrospective case-control data where some of the covariates are measured with error. We show that prospective methods can be applied and the case-control sampling scheme can be ignored if one adequately models the distribution of the error-prone covariates in the case-control sampling scheme. Indeed, subject to this, the prospective likelihood methods result in consistent estimates and information standard errors are asymptotically correct. However, the distribution of such covariates is not the same in the population and under case-control sampling, dictating the need to model the distribution flexibly. In this article, we illustrate the general principle by modeling the distribution of the continuous error-prone covariates using the skewnormal distribution. The performance of the method is evaluated through simulation studies, which show satisfactory results in terms of bias and coverage. Finally, the method is applied to the analysis of two data sets which refer, respectively, to a cholesterol study and a study on breast cancer.     

Inference methods for the conditional logistic regression model with longitudinal data

This paper considers inference methods for case-control logistic regression in longitudinal setups. The motivation is provided by an analysis of plains bison spatial location as a function of habitat heterogeneity. The sampling is done according to a longitudinal matched case-control design in which, at certain time points, exactly one case, the actual location of an animal, is matched to a number of controls, the alternative locations that could have been reached. We develop inference methods for the conditional logistic regression model in this setup, which can be formulated within a generalized estimating equation (GEE) framework. This permits the use of statistical techniques developed for GEE-based inference, such as robust variance estimators and model selection criteria adapted for non-independent data. The performance of the methods is investigated in a simulation study and illustrated with the bison data analysis.     

Detection of bias in harvest-based estimates of chronic wasting disease prevalence in mule deer

Diseased animals may exhibit behavioral shifts that increase or decrease their probability of being randomly sampled. In harvest-based sampling approaches, animal movements, changes in habitat utilization, changes in breeding behaviors during harvest periods, or differential susceptibility to harvest via behaviors like hiding or decreased sensitivity to stimuli may result in a non-random sample that biases prevalence estimates. We present a method that can be used to determine whether bias exists in prevalence estimates from harvest samples. Using data from harvested mule deer (Odocoileus hemionus) sampled in northcentral Colorado (USA) during fall hunting seasons 1996-98 and Akaike's information criterion (AIC) model selection, we detected within-yr trends indicating potential bias in harvest-based prevalence estimates for chronic wasting disease (CWD). The proportion of CWD-positive deer harvested slightly increased through time within a yr. We speculate that differential susceptibility to harvest or breeding season movements may explain the positive trend in proportion of CWD-positive deer harvested during fall hunting seasons. Detection of bias may provide information about temporal patterns of a disease, suggest biological hypotheses that could further understanding of a disease, or provide wildlife managers with information about when diseased animals are more or less likely to be harvested. Although AIC model selection can be useful for detecting bias in data, it has limited utility in determining underlying causes of bias. In cases where bias is detected in data using such model selection methods, then design-based methods (i.e., experimental manipulation) may be necessary to assign causality.     

Method for Evaluating Multiple Mediators: Mediating Effects of Smoking and COPD on the Association between the CHRNA5-A3 Variant and Lung Cancer Risk

A mediation model explores the direct and indirect effects between an independent variable and a dependent variable by including other variables (or mediators). Mediation analysis has recently been used to dissect the direct and indirect effects of genetic variants on complex diseases using case-control studies. However, bias could arise in the estimations of the genetic variant-mediator association because the presence or absence of the mediator in the study samples is not sampled following the principles of case-control study design. In this case, the mediation analysis using data from case-control studies might lead to biased estimates of coefficients and indirect effects. In this article, we investigated a multiple-mediation model involving a three-path mediating effect through two mediators using case-control study data. We propose an approach to correct bias in coefficients and provide accurate estimates of the specific indirect effects. Our approach can also be used when the original case-control study is frequency matched on one of the mediators. We employed bootstrapping to assess the significance of indirect effects. We conducted simulation studies to investigate the performance of the proposed approach, and showed that it provides more accurate estimates of the indirect effects as well as the percent mediated than standard regressions. We then applied this approach to study the mediating effects of both smoking and chronic obstructive pulmonary disease (COPD) on the association between the CHRNA5-A3 gene locus and lung cancer risk using data from a lung cancer case-control study. The results showed that the genetic variant influences lung cancer risk indirectly through all three different pathways. The percent of genetic association mediated was 18.3% through smoking alone, 30.2% through COPD alone, and 20.6% through the path including both smoking and COPD, and the total genetic variant-lung cancer association explained by the two mediators was 69.1%.     

A Gaussian copula approach for the analysis of secondary phenotypes in case-control genetic association studies

In many case-control genetic association studies, a set of correlated secondary phenotypes that may share common genetic factors with disease status are collected. Examination of these secondary phenotypes can yield valuable insights about the disease etiology and supplement the main studies. However, due to unequal sampling probabilities between cases and controls, standard regression analysis that assesses the effect of SNPs (single nucleotide polymorphisms) on secondary phenotypes using cases only, controls only, or combined samples of cases and controls can yield inflated type I error rates when the test SNP is associated with the disease. To solve this issue, we propose a Gaussian copula-based approach that efficiently models the dependence between disease status and secondary phenotypes. Through simulations, we show that our method yields correct type I error rates for the analysis of secondary phenotypes under a wide range of situations. To illustrate the effectiveness of our method in the analysis of real data, we applied our method to a genome-wide association study on high-density lipoprotein cholesterol (HDL-C), where "cases" are defined as individuals with extremely high HDL-C level and "controls" are defined as those with low HDL-C level. We treated 4 quantitative traits with varying degrees of correlation with HDL-C as secondary phenotypes and tested for association with SNPs in LIPG, a gene that is well known to be associated with HDL-C. We show that when the correlation between the primary and secondary phenotypes is >0.2, the P values from case-control combined unadjusted analysis are much more significant than methods that aim to correct for ascertainment bias. Our results suggest that to avoid false-positive associations, it is important to appropriately model secondary phenotypes in case-control genetic association studies.     

Design Aspects and Odds Ratio Estimation in Multicenter Case-Control Studies

In this paper we propose a method to be used in the planning stage of a case-control study. An allocation rule for controls in multicenter case-control studies is proposed which would assure a simple, efficient and unbiased estimation of the odds ratio in the pooled data. It is shown that the efficiency of the design increases with increasing correlation between study center and risk factor. Sources of bias and their implications for relative risk estimation are discussed. The method is demonstrated with data from a case-control study.     

The bias due to incomplete matching

Observational studies comparing groups of treated and control units are often used to estimate the effects caused by treatments. Matching is a method for sampling a large reservoir of potential controls to produce a control group of modest size that is ostensibly similar to the treated group. In practice, there is a trade-off between the desires to find matches for all treated units and to obtain matched treated-control pairs that are extremely similar to each other. We derive expressions for the bias in the average matched pair difference due to the failure to match all treated units--incomplete matching, and the failure to obtain exact matches--inexact matching. A practical example shows that the bias due to incomplete matching can be severe, and moreover, can be avoided entirely by using an appropriate multivariate nearest available matching algorithm, which, in the example, leaves only a small residual bias due to inexact matching.     

Two Methods for Engaging with the Community in Setting Priorities for Child Health Research: Who Engages?

ObjectiveThe aims of this study were to assess participatory methods for obtaining community views on child health research.BackgroundCommunity participation in research is recognised as an important part of the research process; however, there has been inconsistency in its implementation and application in Australia. The Western Australian Telethon Kids Institute Participation Program employs a range of methods for fostering active involvement of community members in its research. These include public discussion forums, called Community Conversations. While participation levels are good, the attendees represent only a sub-section of the Western Australian population. Therefore, we conducted a telephone survey of randomly selected households to evaluate its effectiveness in eliciting views from a broader cross-section of the community about our research agenda and community participation in research, and whether the participants would be representative of the general population. We also conducted two Conversations, comparing the survey as a recruitment tool and normal methods using the Participation Program.ResultsWhile the telephone survey was a good method for eliciting community views about research, there were marked differences in the profile of study participants compared to the general population (e.g. 78% vs 50% females). With a 26% response rate, the telephone survey was also more expensive than a Community Conversation. The cold calling approach proved an unsuccessful recruitment method, with only two out of a possible 816 telephone respondents attending a Conversation.ConclusionWhile the results showed that both of the methods produced useful input for our research program, we could not conclude that either method gained input that was representative of the entire community. The Conversations were relatively low-cost and provided more in-depth information about one subject, whereas the telephone survey provided information across a greater range of subjects, and allowed more quantitative analysis.