Partially parametric techniques for multiple imputation

Correlated pseudorandom variables with prescribed marginal distribution functions sometimes are required in simulation such as in Monte Carlo studies. In this paper, we present a general procedure and a simple but effective numerical approach to generating correlated random variables sampling sequence with prescribed marginal probability distribution functions and correlation coefficient matrix based on linear transformation-nonlinear transformation with Choesky factor. Some simulation results are reported. Simulation results show that the collections of random numbers generated by the presented procedure have desired correlations and pass the Kolmogorov-Smirnov non-parametric hypothesis test of specified marginal distribution. Some restrictions on the application of this method are discussed.

Missing Data and Imputation: A Practical Illustration in a Prognostic Study on Low Back Pain

Missing data are nearly always a problem in research, and missing values represent a serious threat to the validity of inferences drawn from findings. Increasingly, social science researchers are turning to multiple imputation to handle missing data. Multiple imputation, in which missing values are replaced by values repeatedly drawn from conditional probability distributions, is an appropriate method for handling missing data when values are not missing completely at random. However, use of this method requires developing an imputation model from the observed data. This is typically a rigorous and time-consuming process. To encourage wider adoption of multiple imputation in social work research, a simple framework for designing imputation models is presented. The framework and its ability to generate unbiased estimates are demonstrated in a simulation study. KEY WORDS: missing data; multiple imputation; nonresponse ********** Missing data are ubiquitous in social research, and missingness or nonresponse can represent a threat to the validity of inferences because of undue effects on efficiency, power, and parameter bias (Shadish, Cook, & Campbell, 2002). Social work researchers are now addressing missing data in a more rigorous manner. Recently, Saunders et al. (2006) and Choi, Golder, Gillmore, and Morrison (2005) described important data imputation methods and dispelled misunderstandings regarding popular imputation methods, such as mean substitution. Recent advances in analytic methods, such as multiple imputation (MI), are taking hold in social work research. With MI, missing values are replaced with values repeatedly drawn from simulated conditional probability distributions (Schafer, 1997), thus creating multiple versions of the data set. Each version of the data set is analyzed according to the data analysis model, and the multiple results are combined into point estimates (Rubin, 1996). A critical task in MI is to devise an imputation model (Allison, 2002) or missing data model (Graham, Olchowski, & Gilreath, 2007), which involves specifying the measures that are putatively associated with the missing values. Although this process adds additional steps, the specification of an imputation model and the creation of multiple data sets can produce less-biased estimates in the presence of missing data across a wide variety of data analysis techniques (Schafer, 1997). Besides MI, there are many other methods for addressing missing data (Schafer, 1999; Schafer & Graham, 2002). An equally rigorous method known as direct or full information maximum likelihood (FIML) estimation can produce unbiased estimates and correct standard errors in the presence of missing data. When the number of imputations is sufficiently large, identical missing data models will produce the same estimates under MI and FIML (Graham et al., 2007). Unlike MI, FIML is limited to maximum likelihood analytic techniques and the missing data model must be included in the analysis model. Although we focus on MI, the steps we describe for developing an imputation model are equally appropriate for use in the missing data model for FIML. On the basis of the MI literature, this article describes a framework for developing an imputation model for use with any free or commercial software package that performs MI. BRIEF REVIEW OF MISSING DATA CONCEPTS Generally, both MI and a broad range of missing data issues have received ample attention in the applied literature (Graham et al., 2007). We briefly discuss the three types of distributions that describe the randomness of nonresponse given that this property has consequences for the development of an imputation model. For a discussion of general missing data concepts that are not critical to understanding our discussion of imputation model development, we refer readers to Schafer and Graham (2002). Distribution of Nonresponse The probability distribution of nonresponse--more frequently referred to as the missing data or nonresponse mechanism (Rubin, 1976)--is both an important factor in the decision to impute with MI and a context for the development of an imputation model. …

BackgroundMissing data are a common problem in prospective studies with a long follow-up, and the volume, pattern and reasons for missing data may be relevant when estimating the cost of illness. We aimed to evaluate the effects of different methods for dealing with missing longitudinal cost data and for costing caregiver time on total societal costs in Alzheimer’s disease (AD).MethodsGERAS is an 18-month observational study of costs associated with AD. Total societal costs included patient health and social care costs, and caregiver health and informal care costs. Missing data were classified as missing completely at random (MCAR), missing at random (MAR) or missing not at random (MNAR). Simulation datasets were generated from baseline data with 10–40 % missing total cost data for each missing data mechanism. Datasets were also simulated to reflect the missing cost data pattern at 18 months using MAR and MNAR assumptions. Naïve and multiple imputation (MI) methods were applied to each dataset and results compared with complete GERAS 18-month cost data. Opportunity and replacement cost approaches were used for caregiver time, which was costed with and without supervision included and with time for working caregivers only being costed.ResultsTotal costs were available for 99.4 % of 1497 patients at baseline. For MCAR datasets, naïve methods performed as well as MI methods. For MAR, MI methods performed better than naïve methods. All imputation approaches were poor for MNAR data. For all approaches, percentage bias increased with missing data volume. For datasets reflecting 18-month patterns, a combination of imputation methods provided more accurate cost estimates (e.g. bias: −1 % vs −6 % for single MI method), although different approaches to costing caregiver time had a greater impact on estimated costs (29–43 % increase over base case estimate).ConclusionsMethods used to impute missing cost data in AD will impact on accuracy of cost estimates although varying approaches to costing informal caregiver time has the greatest impact on total costs. Tailoring imputation methods to the reason for missing data will further our understanding of the best analytical approach for studies involving cost outcomes.Electronic supplementary materialThe online version of this article (doi:10.1186/s12874-016-0188-1) contains supplementary material, which is available to authorized users.

ObjectiveTo compare the effects of accounting for different missing data mechanisms in an investigation of the role of lifecourse socio-economic position (SEP) on later-life verbal cognitive ability.DesignTwo UK prospective cohort...

The occurrence of missing data due to protocol deviations is inevitable in clinical trials. When missing data exist, analyses rely on assumptions about the behaviour of the individuals after dropping out. As a result, sensitivity analysis, which is now advocated by regulatory bodies, should be performed to explore the robustness of the inference to those assumptions. These assumptions should be relevant to the estimand of the study and be practically accessible by all parties. The aim of this document is twofold: to assess the statistical validity of a new method for sensitivity analysis, and apply this method to a published Alzheimer’s study. At the beginning of the thesis, a description of the Alzheimer’s study and issues with missing data encountered therein, take place. This study was mainly set up to investigate the effect rosiglitazone, as an adjunct therapy in Alzheimer’s patients. Two different doses of the drug were compared to placebo. The study suffered from a moderate degree of missing data in each treatment arm. The thesis proceeds with a critique on the per-protocol and intention-to-treat estimands, and revisits their meaning when missing data occur. Two new estimands are introduced, which are particularly amenable to studies with missing data. They are termed de-jure and de-facto. Following that, the main methods for dealing with missing data are introduced, with a particular emphasis on multiple imputation, and how it can easily incorporate missing not at random (MNAR) analyses. A thorough presentation of the new methodology is given. This is built around a set of assumptions, that reflect possible distributional behaviours of the subjects after protocol deviation. The assumptions are Randomised-treatment Missing at Random (MAR), Jump to Reference, Last Mean Carried Forward, Copy Increments in Reference, and Copy Reference. Estimation and inference is achieved via multiple imputation, and it is shown how the predictive distribution of the imputation model can be constructed from parameters borrowed from an MAR model and manipulated in a pattern-mixture model approach, to obtain the five assumptions for the 4 unobserved component of groups of individuals. A number of simulations whose aim is to explore the statistical properties of the new method, are carried out. The simulated datasets, which are based on the parameters from the Alzheimer’s study, focus on the estimator bias, the size and power of the methods, the bias of the variance estimator, and coverage. The results obtained from the simulations show the method has sensible properties; no bias for the estimator was detectable and the sizes and power of the methods agreed closely with their theoretical equivalents. The main result however, pertains to Rubin’s variance estimator, which proves to appropriately reflect the loss of information from missing data. It is therefore argued, it is the right estimator to use in this setting. The results from the application of the proposed method on the Alzheimer’s dataset are presented in tables. Inferences from the sensitivity analysis assumptions were consistent with those from the original MAR analyses. The comparison between high dose rosiglitazone and placebo did not show any evidence in favour of the treatment under any sensitivity assumption. Some evidence of treatment difference existed when the low dose treatment was compared to placebo. This finding though, should be interpreted with caution, as the differences were obtained from analyses not subjected to the rigorous inferential process that was used in the Harringtons study. It was further argued, this finding might have been due to chance, and it was not replicated in a different study.

Background:The National Trauma Data Bank (NTDB) is plagued by the problem of missing physiological data. The Glasgow Coma Scale score, Respiratory Rate and Systolic Blood Pressure are an essential part of risk adjustment strategies for trauma system evaluation and clinical research. Missing data on these variables may compromise the feasibility and the validity of trauma group comparisons.Aims:To evaluate the validity of Multiple Imputation (MI) for completing missing physiological data in the National Trauma Data Bank (NTDB), by assessing the impact of MI on 1) frequency distributions, 2) associations with mortality, and 3) risk adjustment.Methods:Analyses were based on 170,956 NTDB observations with complete physiological data (observed data set). Missing physiological data were artificially imposed on this data set and then imputed using MI (MI data set). To assess the impact of MI on risk adjustment, 100 pairs of hospitals were randomly selected with replacement and compared using adjusted Odds Ratios (OR) of mortality. OR generated by the observed data set were then compared to those generated by the MI data set.Results:Frequency distributions and associations with mortality were preserved following MI. The median absolute difference between adjusted OR of mortality generated by the observed data set and by the MI data set was 3.6% (inter-quartile range: 2.4%-6.1%).Conclusions:This study suggests that, provided it is implemented with care, MI of missing physiological data in the NTDB leads to valid frequency distributions, preserves associations with mortality, and does not compromise risk adjustment in inter-hospital comparisons of mortality.

ObjectivesResearchers are concerned whether multiple imputation (MI) or complete case analysis should be used when a large proportion of data are missing. We aimed to provide guidance for drawing conclusions from data with a large proportion of missingness. Study Design and SettingVia simulations, we investigated how the proportion of missing data, the fraction of missing information (FMI), and availability of auxiliary variables affected MI performance. Outcome data were missing completely at random or missing at random (MAR). ResultsProvided sufficient auxiliary information was available; MI was beneficial in terms of bias and never detrimental in terms of efficiency. Models with similar FMI values, but differing proportions of missing data, also had similar precision for effect estimates. In the absence of bias, the FMI was a better guide to the efficiency gains using MI than the proportion of missing data. ConclusionWe provide evidence that for MAR data, valid MI reduces bias even when the proportion of missingness is large. We advise researchers to use FMI to guide choice of auxiliary variables for efficiency gain in imputation analyses, and that sensitivity analyses including different imputation models may be needed if the number of complete cases is small.

Multiple imputation is a common technique for dealing with missing values and is mostly applied in regression settings. Its application in cluster analysis problems, where the main objective is to classify individuals into homogenous groups, involves several difficulties which are not well characterized in the current literature. In this paper, we propose a framework for applying multiple imputation to cluster analysis when the original data contain missing values. The proposed framework incorporates the selection of the final number of clusters and a variable reduction procedure, which may be needed in data sets where the ratio of the number of persons to the number of variables is small. We suggest some ways to report how the uncertainty due to multiple imputation of missing data affects the cluster analysis outcomes-namely the final number of clusters, the results of a variable selection procedure (if applied), and the assignment of individuals to clusters. The proposed framework is illustrated with data from the Phenotype and Course of Chronic Obstructive Pulmonary Disease (PAC-COPD) Study (Spain, 2004-2008), which aimed to classify patients with chronic obstructive pulmonary disease into different disease subtypes.

On the weak convergence of the empirical conditional copula under a simplifying assumption

e15003 Background: Stella Trial is a phase III, multinational, double-blind and randomized study to confirm clinical similarity between MB02 and EU-bevacizumab in patients with stage IIIB/IV no squamous NSCLC. 627 subjects with newly diagnosed or recurrent stage IIIb/IV NSCLC were randomized 1:1 to receive either MB02 or EU-bevacizumab plus chemotherapy (paclitaxel and carboplatin) every 3-week cycle for six cycles (week 18) followed by MB02/bevacizumab in blinded monotherapy until disease progression, treatment intolerance, death, patient withdrawal or end of study (w52). An approach targeting the MI of the continuous sum of target diameters (mm), and subsequent categorization of Overall Response (OR) was performed as sensitivity analysis for the assessment of the primary endpoint. Multiple Imputation (MI) has emerged as a credible method to assess the effects of missing data (MD), an inescapable problem with a potential ability to undermine research results' strength and validity in a clinical study, providing the user with a valuable toolset to sufficiently account for the varying types of MD and appropriately adjust the assumptions. Methods: The primary endpoint of the study was Objective Response Rate (ORR) at week 18 per an Independent Radiological Committee (IRC) in the Intention to treat set (ITT). The ORR was analysed with a Cochran-Mantel-Haenszel model, including the stratification factors of sex, smoking status, disease diagnosis and disease stage, comparing the stratified estimates risk difference (RD). In addition, these results were analysed implementing a pattern-mixture MI process accounting for missing at random (MAR) and missing not at random (MNAR) data based on the sum of target diameters in subjects without tumour response data or falling into the categories NonCR/NonPD or non-evaluable (NE). Results: The clinical equivalence of MB02 with EU-bevacizumab is demonstrated by the efficacy data provided through the primary endpoint and the MI sensitive analysis applied. The MI (MAR and MNAR) represents an added value that supports the biosimilarity of MB02 and EU-bevacizumab. The results for the primary analysis of the RD in ORR at w18 in the ITT set were entirely contained within the boundaries of predefined margins (±12): -4.02 [95% CI: -11.76 to 3.71]). Under MI, the ORR RD showed similarity at 95% CIs (-1.92; 95% CI: -10.02 to 6.19) and using the multiple imputations for subjects without tumour response data (missing, NonCR/NonPD or NE), the ORR RD was -2.22 with 95% CI of (-10.54 to 6.10) at w18 in the ITT. Conclusions: The clinical equivalence of MB02 with EU-bevacizumab is demonstrated by the efficacy data provided through the primary endpoint and the MI sensitive analysis applied. The MI (MAR and MNAR) represent an added value that supports the biosimilarity of MB02 and EU-bevacizumab. Clinical trial information: NCT03296163.

IntroductionWe assessed the added value of multiple imputation (MI) of missing values in longitudinal datasets with carotid intima-media thickness (CIMT) as primary endpoint subsequently analysed with linear mixed effects (LME)...

In this study, we consider the nonparametric quantile regression model with the covariates Missing at Random (MAR). Multiple imputation is becoming an increasingly popular approach for analyzing missing data, which combined with quantile regression is not well-developed. We propose an effective and accurate two-stage multiple imputation method for the model based on the quantile regression, which consists of initial imputation in the first stage and multiple imputation in the second stage. The estimation procedure makes full use of the entire dataset to achieve increased efficiency and we show the proposed two-stage multiple imputation estimator to be asymptotically normal. In simulation study, we compare the performance of the proposed imputation estimator with Complete Case (CC) estimator and other imputation estimators, e.g., the regression imputation estimator and k-Nearest-Neighbor imputation estimator. We conclude that the proposed estimator is robust to the initial imputation and illustrates more desirable performance than other comparative methods. We also apply the proposed multiple imputation method to an AIDS clinical trial data set to show its practical application.

Karmakar and Simonovic (2008)describe the methodology of assigning appropriate marginal distributions for three flood characteristics. It is found that the gamma distribution is best fitted for peak flow (P), and a nonparametric distribution from the orthonormal series method best fits to volume (V) and duration (D), based on the root mean square error, Akaike information criterion and Bayesian information criteria. In addition, the chi‐square test is performed to check the significance of fitness. In this paper, a methodology is developed to derive bivariate joint distributions of the flood characteristics using the concept of copulas, considering a set of parametric and nonparametric marginal distributions forP, VandDto mathematically model the correlated structure among them. In the conventional method of flood frequency analysis, the marginal distribution functions of peak flow, volume and duration are assumed to follow some specific parametric distribution function. The concept of copulas relaxes the restriction of traditional flood frequency analysis by selecting marginals from different families of probability distribution functions for flood characteristics. The present study performs a better selection of marginal distribution functions for flood characteristics by parametric and nonparametric estimation procedures, and demonstrates how the concept of copulas may be used for establishing a joint distribution function with mixed marginal distributions. The results obtained are useful for hydrologic design and planning purposes. The methodology is demonstrated with 70 years of stream flow data of Red River at Grand Forks of North Dakota, USA.

BackgroundThe occurrence and timing of mycobacterial culture conversion is used as a proxy for tuberculosis treatment response. When researchers serially sample sputum during tuberculosis studies, contamination or missed visits leads to missing data points. Traditionally, this is managed by ignoring missing data or simple carry-forward techniques. Statistically advanced multiple imputation methods potentially decrease bias and retain sample size and statistical power.MethodsWe analyzed data from 261 participants who provided weekly sputa for the first 12 weeks of tuberculosis treatment. We compared methods for handling missing data points in a longitudinal study with a time-to-event outcome. Our primary outcome was time to culture conversion, defined as two consecutive weeks with no Mycobacterium tuberculosis growth. Methods used to address missing data included: 1) available case analysis, 2) last observation carried forward, and 3) multiple imputation by fully conditional specification. For each method, we calculated the proportion culture converted and used survival analysis to estimate Kaplan-Meier curves, hazard ratios, and restricted mean survival times. We compared methods based on point estimates, confidence intervals, and conclusions to specific research questions.ResultsThe three missing data methods lead to differences in the number of participants achieving conversion; 78 (32.8%) participants converted with available case analysis, 154 (64.7%) converted with last observation carried forward, and 184 (77.1%) converted with multiple imputation. Multiple imputation resulted in smaller point estimates than simple approaches with narrower confidence intervals. The adjusted hazard ratio for smear negative participants was 3.4 (95% CI 2.3, 5.1) using multiple imputation compared to 5.2 (95% CI 3.1, 8.7) using last observation carried forward and 5.0 (95% CI 2.4, 10.6) using available case analysis.ConclusionWe showed that accounting for missing sputum data through multiple imputation, a statistically valid approach under certain conditions, can lead to different conclusions than naïve methods. Careful consideration for how to handle missing data must be taken and be pre-specified prior to analysis. We used data from a TB study to demonstrate these concepts, however, the methods we described are broadly applicable to longitudinal missing data. We provide valuable statistical guidance and code for researchers to appropriately handle missing data in longitudinal studies.