MIMCA: multiple imputation for categorical variables with multiple correspondence analysis

We propose a multiple imputation method to deal with incomplete categorical data. This method imputes the missing entries using the principal component method dedicated to categorical data: multiple correspondence analysis (MCA). The uncertainty concerning the parameters of the imputation model is reflected using a non-parametric bootstrap. Multiple imputation using MCA (MIMCA) requires estimating a small number of parameters due to the dimensionality reduction property of MCA. It allows the user to impute a large range of data sets. In particular, a high number of categories per variable, a high number of variables or a small number of individuals are not an issue for MIMCA. Through a simulation study based on real data sets, the method is assessed and compared to the reference methods (multiple imputation using the loglinear model, multiple imputation by logistic regressions) as well to the latest works on the topic (multiple imputation by random forests or by the Dirichlet process mixture of products of multinomial distributions model). The proposed method provides a good point estimate of the parameters of the analysis model considered, such as the coefficients of a main effects logistic regression model, and a reliable estimate of the variability of the estimators. In addition, MIMCA has the great advantage that it is substantially less time consuming on data sets of high dimensions than the other multiple imputation methods.     

Asymptotic behaviour of binned kernel density estimators for locally non-stationary random fields

We investigate the asymptotic behaviour of binned kernel density estimators for dependent and locally non-stationary random fields converging to stationary random fields. We focus on the study of the bias and the asymptotic normality of the estimators. A simulation experiment conducted shows that both the kernel density estimator and the binned kernel density estimator have the same behavior and both estimate accurately the true density when the number of fields increases. We apply our results to the 2002 incidence rates of tuberculosis in the departments of France.     

An Ecological Model for Quantitative Risk Assessment for Schistosomiasis: The Case of a Patchy Environment in the Coastal Tropical Area of Northeastern Brazil

We developed a stochastic model for quantitative risk assessment for the Schistosoma mansoni (SM) parasite, which causes an endemic disease of public concern. The model provides answers in a useful format for public health decisions, uses data and expert opinion, and can be applied to any landscape where the snail Biomphalaria glabrata is the main intermediate host (South and Central America, the Caribbean, and Africa). It incorporates several realistic and case‐specific features: stage‐structured parasite populations, periodic praziquantel (PZQ) drug treatment for humans, density dependence, extreme events (prolonged rainfall), site‐specific sanitation quality, environmental stochasticity, monthly rainfall variation, uncertainty in parameters, and spatial dynamics. We parameterize the model through a real‐world application in the district of Porto de Galinhas (PG), one of the main touristic destinations in Brazil, where previous studies identified four parasite populations within the metapopulation. The results provide a good approximation of the dynamics of the system and are in agreement with our field observations, i.e., the lack of basic infrastructure (sanitation level and health programs) makes PG a suitable habitat for the persistence and growth of a parasite metapopulation. We quantify the risk of SM metapopulation explosion and quasi‐extinction and the time to metapopulation explosion and quasi‐extinction. We evaluate the sensitivity of the results under varying scenarios of future periodic PZQ treatment (based on the Brazilian Ministry of Health's plan) and sanitation quality. We conclude that the plan might be useful to slow SM metapopulation growth but not to control it. Additional investments in better sanitation are necessary.     

Risk Assessment for Hospital-Acquired Diseases: A Risk-Theory Approach

We introduce a new approach to hospital-acquired disease risk assessment from public health databases. In a spirit similar to actuarial risk theory, we define an adjustment coefficient that can quantify the risk associated with a hospital department, allowing comparisons of similar departments. The adjustment coefficient characterizes the tail of the distribution of the total patient length of stay in a department before the first disease event occurs. We show that this coefficient is the solution of a Lundberg-like equation, and we provide a nonparametric estimation procedure for this measure, based on a Cramér-Lundberg approximation for the tail of the distribution. Using simulations, we provide evidence of the robustness of the approximation to various individual risk models. In addition, we illustrate the relevance of this approach by evaluating the risk associated with a standard patient safety indicator in 20 hospitals of southeastern France.     

An Empirical Process View of Inverse Regression

A common approach taken in high‐dimensional regression analysis is sliced inverse regression, which separates the range of the response variable into non‐overlapping regions, called ‘slices’. Asymptotic results are usually shown assuming that the slices are fixed, while in practice, estimators are computed with random slices containing the same number of observations. Based on empirical process theory, we present a unified theoretical framework to study these techniques, and revisit popular inverse regression estimators. Furthermore, we introduce a bootstrap methodology that reproduces the laws of Cramér–von Mises test statistics of interest to model dimension, effects of specified covariates and whether or not a sliced inverse regression estimator is appropriate. Finally, we investigate the accuracy of different bootstrap procedures by means of simulations.     

The Comparison of Performances of Widely Used Cointegration Tests

The aim of this study is to compare performances of commonly cointegration tests used in literature in terms of their empirical power and type I error probabilty for various sample sizes. As a result of the study, it has been found that some tests are not appropriate in testing cointegration in terms of empirical power and type I error probability. As a result of simulation study, λ_max test for any values of ρ and sample sizes have been found most appropriate test in conclusion.     

New Shrinkage Parameters for the Liu-type Logistic Estimators

The binary logistic regression is a widely used statistical method when the dependent variable has two categories. In most of the situations of logistic regression, independent variables are collinear which is called the multicollinearity problem. It is known that multicollinearity affects the variance of maximum likelihood estimator (MLE) negatively. Therefore, this article introduces new shrinkage parameters for the Liu-type estimators in the Liu (2003) in the logistic regression model defined by Huang (2012) in order to decrease the variance and overcome the problem of multicollinearity. A Monte Carlo study is designed to show the goodness of the proposed estimators over MLE in the sense of mean squared error (MSE) and mean absolute error (MAE). Moreover, a real data case is given to demonstrate the advantages of the new shrinkage parameters.     

Welfare analysis of a tax reform for Germany: a comparison of the unitary and collective models of household labour supply

We compare results of a tax reform analysis obtained with the collective and unitary models of household behaviour. We simulate real world micro-data by means of a collective approach, using a compound procedure of estimation and calibration based on the 1998 wave of the German socio-economic panel. We estimate a unitary model on this ‘collective’ data set. Investigating a move from joint to individual taxation on the basis of both models, we obtain important discrepancies between predicted adjustments to labour supply and distortions in the welfare analysis of the reform on the basis of unitary estimates.