Search PubMed⌕ Search

SEARCH · Search PubMed

Results for “Imputation”

Search indexed PubMed citations on genomics, clinical trials, systematic reviews and public health. Explore titles, authors and supplied subject terms, then open the PubMed record.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 343 records · Page 19Linked to original sources

Group comparisons involving missing data in clinical trials: a comparison of estimates and power (size) for some simple approaches.

When using 'intent-to-treat' approaches to compare outcomes between groups in clinical trials, analysts face a decision regarding how to account for missing observations. Most model-based approaches can be summarized as a process whereby the analyst makes assumptions about the distribution of the missing data in an attempt to obtain unbiased estimates that are based on functions of the observed data. Although pointed out by Rubin as often leading to biased estimates of variances, an alternative approach that continues to appear in the applied literature is to use fixed-value imputation of means for missing observations. The purpose of this paper is to provide illustrations of how several fixed-value mean imputation schemes can be formulated in terms of general linear models that characterize the means of distributions of missing observations in terms of the means of the distributions of observed data. We show that several fixed-value imputation strategies will result in estimated intervention effects that correspond to maximum likelihood estimates obtained under analogous assumptions. If the missing data process has been correctly characterized, hypothesis tests based on variances estimated using maximum likelihood techniques asymptotically have the correct size. In contrast, hypothesis tests performed using the uncorrected variance, obtained by applying standard complete data formula to singly imputed data, can provide either conservative or anticonservative results. Surprisingly, under several non-ignorable non-response scenarios, maximum likelihood based analyses can yield equivalent hypothesis tests to those obtained when analysing only the observed data.

Aged↗

Use of a single global assessment to reduce missing data in a clinical trial with follow-up at one year.

We conducted a randomized controlled trial (ISRCTN96537534) to assess the effects of acupuncture on migraine and chronic tension headache. Patients (n=401) completed a diary of headache severity four times a day for 4 weeks at baseline, immediately following a 3-month treatment period and 1 year after randomization. During the trial, it appeared that dropout might be higher than expected. We therefore obtained a rapid global assessment of headache from participants to aid imputation of missing data. Patients were contacted by telephone and asked to rate current and baseline headache on a 0-10 scale. Use of global assessment reduced the number of patients from whom we obtained no follow-up headache data from 69 (17%) to 24 (6%). Analysis of patients who provided both a diary and a global assessment demonstrated excellent properties of global assessment, with very similar results to the full diary. We therefore used the global assessment to help impute missing 1-year diary scores. Rapid global assessment can be easily implemented in any trial and aids imputation of missing data, though it should not be used instead of more intensive methods of assessment. Further research might usefully examine the value of global assessment for imputation of missing data in a variety of different settings.

Acupuncture↗

The long-term course of low-density lipoprotein cholesterol after initiation of statin treatment: retrospective database analysis over 3 years in health maintenance organization enrollees.

OBJECTIVES: Our primary objective was to obtain robust estimates of the low-density lipoprotein cholesterol (LDL-C) decrease from baseline over a long period (ie, 3 years) after initiation of statin treatment in a usual-care setting. Our secondary objective was to investigate the predictors of the LDL-C time course. METHODS: We retrospectively analyzed the data for a sample of enrollees in a health maintenance organization (HMO) who started statin treatment between October 1, 1995, and December 31, 1998. Using the HMO's claims database, we examined the LDL-C change from baseline (as measured at the prescribing physicians' discretion) and computed mean estimates every 6 months up to 3 years. We investigated potential predictors of the LDL-C time course (ie, age, sex, baseline LDL-C, previous treatment, prescribing physician's specialty, and most recent treatment) with a mixed model applied to longitudinal data. This model enabled us to impute missing data for all enrollees still being followed, including those who had stopped treatment, and to discuss the robustness of our findings. We applied 2 methods of imputation, assuming either of the following: (1) data were missing at random but could be estimated from the parameters in the mixed model, or (2) LDL-C returned to the baseline value > or =15 days after treatment cessation. RESULTS: We examined data from 3193 individuals. In most cases, the statin used was fluvastatin or pravastatin. The observed mean (95% CI) LDL-C decrease from baseline widened progressively from 23.6% (23.0%-24.3%) at 6 months to 28.0% (27.1%-28.9%) at 18 months and 30.2% (28.7%-31.7%) at 36 months after treatment initiation. These results remained robust after the imputation of missing data, with mean LDL-C reduction consistently >20% at each 6-month time point during the 3 years after treatment initiation. Variations as a function of baseline characteristics were limited (demographics) or explicable by extraneous factors (baseline LDL-C). There were predictable variations as a function of the most recent treatment. CONCLUSIONS: This analysis indicates a long-term reduction in LDL-C among a sample of HMO enrollees who initiated statin treatment in a usual-care setting. The results were robust after imputation of missing data, with mean decrease from baseline consistently >20% over 3 years. However, given the retrospective design of our study and the absence of a control group, we cannot determine how much of the decrease was attributable to treatment.

Aged↗

Selection bias found in interpreting analyses with missing data for the prehospital index for trauma.

OBJECTIVE: To evaluate the effects of missing data on analyses of data from trauma databases, and to verify whether commonly used techniques for handling missing data work well in theses settings. STUDY DESIGN AND SETTING: Measures of trauma severity such as the Pre-Hospital Index (PHI) are used for triage and the evaluation of trauma care. As conditions of trauma patients can rapidly change over time, estimating the change in PHI from the arrival at the emergency room to hospital admission is important. We used both simulated and real data to investigate the estimation of PHI data when some data are missing. Techniques compared include complete case analysis, single imputation, and multiple imputation. RESULTS: It is well known that complete case analyses and single imputation methods often lead to highly misleading results that can be corrected by multiple imputation, an increasingly popular method for missing data situations. In practice, unverifiable assumptions may not hold, meaning that it may not be possible to draw definitive conclusions from any of the methods. CONCLUSION: Great care is required whenever missing data arises. This is especially true in trauma databases, which often have much missing data and where the data may not missing at random.

Database Management Systems↗

Multivariate modeling of missing data within and across assessment waves.

Missing data constitute a common but widely underappreciated problem in both cross-sectional and longitudinal research. Furthermore, both the gravity of the problems associated with missing data and the availability of the applicable solutions are greatly increased by the use of multivariate analysis. The most common approaches to dealing with missing data are reviewed, such as data deletion and data imputation, and their relative merits and limitations are discussed. One particular form of data imputation based on latent variable modeling, which we call Multivariate Imputation, is highlighted as holding great promise for dealing with missing data in the context of multivariate analysis. The recent theoretical extension of latent variable modeling to growth curve analysis also permitted us to extend the same kind of solution to the problem of missing data in longitudinal studies. Data simulations are used to compare the results of multivariate imputation to other common approaches to missing data.

Bias↗

A vision of how low-coverage sequence data should contribute to genetic evaluation in the future.

Low-coverage sequencing refers to sequencing DNA of individuals to a low depth of coverage (e.g., 0.5X) and imputing that sequence to a genomic sequence based on reference haplotypes from individuals sequenced to a high depth of coverage (e.g., ≥10X). It has been proposed as an alternative to genotyping by Single-nucleotide polymorphisms (SNP) arrays. At least one commercial product based on it is available for agricultural species. Concerns limiting adoption in its current form are: 1) the cost of storing the huge volume of data it generates and 2) whether that additional data will result in improved accuracy of genetic evaluation. This work envisions future implementation of low-coverage sequencing to reduce storage costs and enhance genetic evaluations by leveraging the additional information in the full sequence of the pangenome to account for more genetic variation. We propose addressing the storage issue by representing genomic sequence of an individual in a pair of haplotype arrays with each element pointing to an enumerated haplotype of the sequence within one of approximately 50,000 defined genome segments. Assuming 60 million genomic variants, the infrastructure required to translate the identifier of any enumerated haplotype into its genomic sequence would require less than 10 gigabytes of binary storage. Each haplotype array element would require 2 bytes, so the marginal binary storage required to represent the genomic sequence of an individual would be about 200 kilobytes (KB), similar to the genotypes from a SNP array with 200,000 markers. This assumes no pedigree and no ambiguity of the imputation, though the latter is unrealistic. Strategies to minimize, and when necessary, to manage and efficiently represent ambiguity are proposed. The genomic sequence of an individual could be stored in about 1 KB (binary) if both parents have unambiguous sequences stored as described above. The proposed system for representing the pangenome includes algorithms for read mapping and imputation intended to leverage all known genetic variation in the target population. It is also designed to use sequencing reads generated for imputing the genomic sequence of new individuals to identify unrecognized mutations, crossovers, and structural variants, thus continuously improving the genome representation, especially if widespread use of low-coverage sequencing in livestock industries is realized. This could make improved genetic merit and management of livestock feasible without computational burden.

Animals↗

A critical look at methods for handling missing covariates in epidemiologic regression analyses.

Epidemiologic studies often encounter missing covariate values. While simple methods such as stratification on missing-data status, conditional-mean imputation, and complete-subject analysis are commonly employed for handling this problem, several studies have shown that these methods can be biased under reasonable circumstances. The authors review these results in the context of logistic regression and present simulation experiments showing the limitations of the methods. The method based on missing-data indicators can exhibit severe bias even when the data are missing completely at random, and regression (conditional-mean) imputation can be inordinately sensitive to model misspecification. Even complete-subject analysis can outperform these methods. More sophisticated methods, such as maximum likelihood, multiple imputation, and weighted estimating equations, have been given extensive attention in the statistics literature. While these methods are superior to simple methods, they are not commonly used in epidemiology, no doubt due to their complexity and the lack of packaged software to apply these methods. The authors contrast the results of multiple imputation to simple methods in the analysis of a case-control study of endometrial cancer, and they find a meaningful difference in results for age at menarche. In general, the authors recommend that epidemiologists avoid using the missing-indicator method and use more sophisticated methods whenever a large proportion of data are missing.

Case-Control Studies↗

Predicting hospital mortality among injured children using a national trauma database.

PURPOSE: The purpose of this study was to develop a model that accurately predicts mortality among injured children based on components of the initial patient evaluation and that is generalizable to diverse acute care settings. Important predictive variables obtained in an emergency setting are frequently missing in even large national databases, limiting their effectiveness for developing predictions. In this study, a model predicting pediatric trauma mortality was developed using a national database and methods to handle missing data that may avoid biases that can occur restricting analyses to complete cases. METHODS: Records of pediatric patients included in the National Pediatric Trauma Registry (NPTR) between 1996 and 1999 were used as a training set in a logistic regression model to predict hospital mortality using vital signs, Glasgow Coma Scale (GCS) score, and intubation status. Multiple imputation was applied to handle missing data. The model was tested using independent data from the NPTR and National Trauma Data Bank (NTDB). RESULTS: Complete case analysis identified only GCS-eye and intubation status as predictors of mortality. A model based on complete case analysis had good discrimination (c-index = 0.784) and excellent calibration (Hosmer-Lemeshow c-statistic, 6.8) (p > 0.05). Using multiple imputation, three additional predictors of mortality (systolic blood pressure, pulse, and GCS-motor) were identified and improved model performance was observed. The model developed using multiple imputation had excellent discrimination (c-index, 0.947-0.973) in both test datasets. Calibration was better in the NPTR testing set than in the NTDB (Hosmer-Lemeshow c-statistic, 9.2 for NPTR [p > 0.05] and 258.2 for NTDB [p < 0.05]). At a probability cutoff that minimized misclassification in the training set, the false-negative and false-negative rates of the model were better than those obtained with either the Revised Trauma Score (RTS) or Pediatric Trauma Score using data from the NPTR testing set. Although the false-positive rates were lower with the RTS using data from the NTDB, the false-negative rates of the proposed model and the RTS were similar in this test dataset. CONCLUSIONS: Using multiple imputation to handle missing data, a model predicting pediatric trauma mortality was developed that compared favorably with existing trauma scores. Application of these methods may produce predictive trauma models that are more statistically reliable and applicable in clinical practice.

Child↗

Evaluating the performance of the Framingham risk equations in a population with diabetes.

AIMS: The Framingham risk equations are widely used to estimate risk of coronary heart disease (CHD). The purpose of this study was to evaluate the reliability of these equations in predicting CHD risk in people with diabetes and the reliability of using imputed mean HDL-cholesterol values. METHODS: Data describing the baseline characteristics of recognized CHD risk factors for 938 people aged 30-74 years were extracted from the Cardiff Diabetes Database. Data describing CHD events were available for up to 4 years following the baseline year (1996). Several mathematical techniques were used to assess the reliability of predictions provided by the Framingham equations in this population. RESULTS: Thirty-four percent of males and 25% of females who experienced CHD events had a predicted 10-year CHD risk >/= 30%. Seventy-five percent of males and 58% of females had a predicted 10-year CHD risk >/= 20%. Using imputed HDL-cholesterol values, 26% of males and 6% of females who later developed CHD events had a 10-year CHD risk >/= 30%. Using imputed HDL-cholesterol values, the CHD risk predicted by the Framingham equations consistently underestimated the actual risk of CHD events. However, refitting the Framingham risk equations to the Cardiff data resulted in only marginal improvements in discriminatory capabilities. CONCLUSIONS: The Framingham risk equations can be unreliable when applied to the diabetic population, tending to underestimate an individual's probability of progressing to CHD; the equations perform marginally better in women than in men. The use of imputed mean HDL-cholesterol values improved the reliability of the estimates of risk.

Adult↗

Drotrecogin alfa (activated) (recombinant human activated protein C) reduces host coagulopathy response in patients with severe sepsis.

Drotrecogin alfa (activated) improved survival in patients with severe sepsis in PROWESS, a double-blind, study of 1690 adult patients randomized to drotrecogin alfa (activated) at 24 microg/kg/h (N=850) or placebo (N=840) infused for 96 hours. Pharmacodynamic effects of drotrecogin alfa (activated) were assessed with 15 prospectively defined systemic biomarkers of hemostasis, inflammation and endothelial injury. The last-observation-carried-forward (LOCF) method of imputation for missing observations was the prospectively defined statistical method. The results were also analyzed with only the observed values without imputation for missing data (repeated measures analysis). With both statistical methods, drotrecogin alfa (activated)-treated patients demonstrated antithrombotic (reduced markers of thrombin generation and accelerated normalization of anticoagulant factor, protein C and fibrinolytic factors) and anticoagulant (prolonged PT and APTT) effects compared with placebo. A profibrinolytic (reduction in plasminogen activator inhibitor-1) effect was significant only with the LOCF imputation method in observed case and percent change from baseline analyses. An anti-inflammatory (reduction in interleukin-6) effect was significant only with the LOCF imputation method in change from baseline and percent change from baseline analyses. Drotrecogin alfa (activated) is a new and promising agent for treatment of patients with severe sepsis. The extensive analysis of systemic biomarkers confirms the previously published antithrombotic effects. However, the present results using different statistical methods do not provide a strong basis for systemic anti-inflammatory or pro-fibrinolytic effects. These latter two effects may occur at the local or cellular level. The systemic biomarkers reported here might not be the most appropriate approach to demonstrate these potential effects of drotrecogin alfa (activated).

Anti-Inflammatory Agents, Non-Steroidal↗

Methods for handling missing data in palliative care research.

Missing data is a common problem in palliative care research due to the special characteristics (deteriorating condition, fatigue and cachexia) of the population. Using data from a palliative study, we illustrate the problems that missing data can cause and show some approaches for dealing with it. Reasons for missing data and ways to deal with missing data (including complete case analysis, imputation and modelling procedures) are explored. Possible mechanisms behind the missing data are: missing completely at random, missing at random or missing not at random. In the example study, data are shown to be missing at random. Imputation of missing data is commonly used (including last value carried forward, regression procedures and simple mean). Imputation affects subsequent summary statistics and analyses, and can have a substantial impact on estimated group means and standard deviations. The choice of imputation method should be carried out with caution and the effects reported.

Adult↗

Quality determination and the repair of poor quality spots in array experiments.

BACKGROUND: A common feature of microarray experiments is the occurrence of missing gene expression data. These missing values occur for a variety of reasons, in particular, because of the filtering of poor quality spots and the removal of undefined values when a logarithmic transformation is applied to negative background-corrected intensities. The efficiency and power of an analysis performed can be substantially reduced by having an incomplete matrix of gene intensities. Additionally, most statistical methods require a complete intensity matrix. Furthermore, biases may be introduced into analyses through missing information on some genes. Thus methods for appropriately replacing (imputing) missing data and/or weighting poor quality spots are required. RESULTS: We present a likelihood-based method for imputing missing data or weighting poor quality spots that requires a number of biological or technical replicates. This likelihood-based approach assumes that the data for a given spot arising from each channel of a two-dye (two-channel) cDNA microarray comparison experiment independently come from a three-component mixture distribution--the parameters of which are estimated through use of a constrained E-M algorithm. Posterior probabilities of belonging to each component of the mixture distributions are calculated and used to decide whether imputation is required. These posterior probabilities may also be used to construct quality weights that can down-weight poor quality spots in any analysis performed afterwards. The approach is illustrated using data obtained from an experiment to observe gene expression changes with 24 hr paclitaxel (Taxol) treatment on a human cervical cancer derived cell line (HeLa). CONCLUSION: As the quality of microarray experiments affect downstream processes, it is important to have a reliable and automatic method of identifying poor quality spots and arrays. We propose a method of identifying poor quality spots, and suggest a method of repairing the arrays by either imputation or assigning quality weights to the spots. This repaired data set would be less biased and can be analysed using any of the appropriate statistical methods found in the microarray literature.

Algorithms↗

Integrative missing value estimation for microarray data.

BACKGROUND: Missing value estimation is an important preprocessing step in microarray analysis. Although several methods have been developed to solve this problem, their performance is unsatisfactory for datasets with high rates of missing data, high measurement noise, or limited numbers of samples. In fact, more than 80% of the time-series datasets in Stanford Microarray Database contain less than eight samples. RESULTS: We present the integrative Missing Value Estimation method (iMISS) by incorporating information from multiple reference microarray datasets to improve missing value estimation. For each gene with missing data, we derive a consistent neighbor-gene list by taking reference data sets into consideration. To determine whether the given reference data sets are sufficiently informative for integration, we use a submatrix imputation approach. Our experiments showed that iMISS can significantly and consistently improve the accuracy of the state-of-the-art Local Least Square (LLS) imputation algorithm by up to 15% improvement in our benchmark tests. CONCLUSION: We demonstrated that the order-statistics-based integrative imputation algorithms can achieve significant improvements over the state-of-the-art missing value estimation approaches such as LLS and is especially good for imputing microarray datasets with a limited number of samples, high rates of missing data, or very noisy measurements. With the rapid accumulation of microarray datasets, the performance of our approach can be further improved by incorporating larger and more appropriate reference datasets.

Algorithms↗

Source of bias in prenatal care utilization indices: implications for evaluating the Medicaid expansion.

BACKGROUND: Recent expansions in eligibility for coverage of prenatal care services by the Medicaid program reflect national initiatives to improve pregnancy outcomes. This study investigates the potential impact that completeness of reporting of prenatal care and gestational age variables and strategies to impute missing data may have on evaluations of the Medicaid expansion. METHODS: This study, examining 15 years of vital record data from a single state and comparing 1 year of data from four mid-Atlantic states, selected single live births to resident mothers for analyses. The "day 15" and the "preceding case" methods were used to impute missing gestational age data. RESULTS: Considerable temporal and geographic variation was detected in completeness of reporting of variables used to construct prenatal care indices. After imputing values for cases with missing data, the proportion of cases for which adequacy of prenatal care utilization could not be determined ranged from 3% to 24% among the states investigated. For those cases where gestational age data could be imputed, the distribution of prenatal care utilization was not markedly disparate from those cases with complete reporting of gestational age. CONCLUSIONS: The results indicate that variations in reporting, decisions regarding the treatment of missing data, and the choice of the denominator can alter prenatal care utilization percentages and have implications for evaluations of the impact of the recent Medicaid expansion on prenatal care utilization.

Bias↗

A prospective study of calf factors affecting age, body size, and body condition score at first calving of holstein dairy heifers.

Data were collected prospectively on parameters related to first calving on 18 farms located in Northeastern Pennsylvania. This project was designed to study possible residual effects of calf management practices and events occurring during the first 16 wk of life on age, BW, skeletal growth, and body condition score at first calving. Multiple imputation method for handling missing data was incorporated in these analyses. This method has the advantage over ad hoc single imputations because the appropriate error structure is maintained. Much similarity was found between the multiple imputation method and a traditional mixed model analysis, except that some estimates from the multiple imputation method seemed more logical in their effects on the parameter measured. Factors related to increased age at first calving were increased difficulty of delivery, antibiotic treatment of sick calves, increased amount of milk or milk replacer fed before weaning, reduced quality of forage fed to weaned calves, maximum humidity, mean daily temperature, and maximum ammonia levels in calf housing areas. Body weight at calving tended to increase with parity of the dam, increased amount of grain fed to calves, increased ammonia levels, and increased mean temperature of the calf housing area. Body condition score at calving tended to be positively influenced by delivery score at first calving, dam parity, and milk or milk replacer dry matter intake. Withers height at calving was positively affected by treatment of animals with antibiotics and increased mean temperature in the calf area. This study demonstrated that nutrition, housing, and management factors that affect health and growth of calves have long-term effects on the animal at least through first calving.

Aging↗

Costs, benefits and effectiveness of worksite physical activity counseling from the employer's perspective.

OBJECTIVES: This study evaluated the impact of worksite physical activity counseling using cost-benefit and cost-effectiveness analyses. METHODS: Civil servants (N = 299) were randomly assigned to an intervention (N = 131) or control (N = 168) group for 9 months. The intervention costs were compared with the monetary benefits gained from reduced sick leave. In addition, the intervention costs minus the monetary benefits from sick leave reduction were compared with the effects (percentage meeting the public health recommendation for moderate-intensity physical activity, energy expenditure, cardiorespiratory fitness, and upper extremity symptoms). RESULTS: The intervention costs were EUR 430 per participant, and the benefits were EUR 125 due to sick leave during the intervention period, for net total costs of EUR 305 for the intervention. During the same 9-month period the year after the intervention, the benefits from sick leave reduction were EUR 635. No statistically significant differences in costs and benefits were found between the groups. As to the cost-effectiveness, improvement in energy expenditure and cardiorespiratory fitness was observed at higher costs. The point estimates of the cost-effectiveness ratios were EUR 5.2 (without imputation of effect data) and EUR 2.7 (with imputation of effect data) per extra kilocalorie of energy expenditure per day and EUR 235 (without imputation of effect data) and EUR 45.9 (with imputation of effect data) per beat per minute of decrease in submaximal heart rate. CONCLUSIONS: This study does not provide a financial reason for implementing worksite counseling intervention on physical activity on the short-term. However, positive effects were shown for energy expenditure and cardiorespiratory fitness.

Adult↗

Equivalence of the acute phase reactants C-reactive protein, plasma viscosity, and Westergren erythrocyte sedimentation rate when used to calculate American College of Rheumatology 20% improvement criteria or the Disease Activity Score in patients with early rheumatoid arthritis. Western Consortium of Practicing Rheumatologists.

OBJECTIVE: In an additive cohort of patients with early rheumatoid arthritis (RA), to determine the effect of substituting one acute phase reactant for another on the number of patients satisfying the American College of Rheumatology (ACR) 20% preliminary criteria for improvement, and on calculated Disease Activity Scores (DAS). METHODS: A total of 251 patients with 6.4 months average disease duration had detailed clinical assessments at entry and 6, 12, and 24 months in a multicenter prospective longterm observational study. Matched erythrocyte sedimentation rate (ESR), C-reactive protein (CRP), and plasma viscosity (PV) assays were done at 366 time points. Disease modifying antirheumatic drugs were not started until after the baseline evaluation. RESULTS: After 6, 12, and 24 months, 50%, 53%, and 57% of patients were responders, as defined by the ACR 20% improvement criteria. The difference in response rates when ESR, CRP, or PV was used as the acute phase reactant ranged from 0.4% at 12 months to 3% at 24 months. Percentile distributions of the 366 matched CRP, ESR, and PV values were used to prepare a nomogram that can be used to calculate the other acute phase reactant values if the value of one is known. When the nomogram was used to impute ESR values from observed PV or CRP values, average DAS scores calculated with the actual ESR values were not different from average DAS scores calculated from the imputed ESR values. CONCLUSION: ESR, CRP, and PV are equally useful in calculating ACR 20% response rates in patients with active early RA. A nomogram can be used to impute ESR values from CRP or PV values; use of the imputed ESR values is as accurate as use of the actual ESR values to calculate average DAS.

Acute-Phase Reaction↗

A preprocessing method for improving data mining techniques. Application to a large medical diabetes database.

The Knowledge Discovery in Databases (KDD) methodology seems to be attractive on the analyze of large clinical databases. In the KDD process, the preprocessing step (data cleaning and handling of missing values) is paramount since it conditions the quality of the results obtained by data mining procedures and represents about 80% of the whole project time. The aims of the present study were to analyze this step and provide tools to handle inconsistent data and missing values. We have broken down the process into 3 main stages: data cleaning--explanatory study of missing values--choice of the procedure used for handling missing values. The data cleaning stage was based on a system of logical rules to correct mistakes and on cluster analysis to discard the poorly filled files. The missing-data mechanism was analyzed by means of multivariate statistical procedures. Two methods to deal with missing values were compared: imputation by the most common value (mode) and imputation using decision trees. This study was performed on a large medical diabetes database (23,601 patients) including numerous missing values. A system of logical rules allowed to correct mistakes on essential parameters (for example, the type of diabetes). Cluster analysis allowed to identify 10% of poorly filled files. After multivariate analysis, the missing-data mechanism could be considered as random. For variables with low number of missing values (< 10%) and categories (< 4), imputation using decision trees provided better results than imputation by mode.

Data Interpretation, Statistical↗