Admissible and optimal sampling strategy for estimating finite population mean in randomized response surveys with multiple responses

ABSTRACT

We consider the problem of estimation of a finite population mean (or proportion) related to a sensitive character under a randomized response model when independent responses are obtained from each sampled individual as many times as he/she is selected in the sample and prove the admissibility of a sampling strategy in a class of comparable linear unbiased strategies. We prove that the admissible strategy is also optimal in this class under a super-population model.     

Analysis of mixed correlated semi-continuous and β-inflated Poisson responses

ABSTRACT

In this article, inflation at an arbitrary point β of a member of power series exponential family and mean-inflation as a cause of having semi-continuous distribution are discussed. Also, a joint modeling of such a semi-continuous response and β-inflated Poisson response is presented. Simultaneous effects of covariates on both responses, which have two-component mixture distributions, are investigated. To find the parameter estimates, the maximum likelihood approach is used. The proposed model is illustrated on some simulation studies and applied to a real survey dataset.     

A class of distributions with the quadratic mean residual quantile function

Abstract

The present paper introduces a new family of distributions with quadratic mean residual quantile function. Various distributional properties as well as reliability characteristics are discussed. Some characterizations of the class of distributions are presented. The estimation of parameters of the model using method of L-moments is studied. The practical application of the class of models is illustrated with a real life data set.     

An asymmetric multivariate weibull distribution

Abstract

A class of multivariate laws as an extension of univariate Weibull distribution is presented. A well known representation of the asymmetric univariate Laplace distribution is used as the starting point. This new family of distributions exhibits some similarities to the multivariate normal distribution. Properties of this class of distributions are explored including moments, correlations, densities and simulation algorithms. The distribution is applied to model bivariate exchange rate data. The fit of the proposed model seems remarkably good. Parameters are estimated and a bootstrap study performed to assess the accuracy of the estimators.     

A Marshall–Olkin family of heavy-tailed distributions which includes the lognormal one

Abstract

A new class of heavy-tailed distribution functions,, containing the lognormal distribution as a particular case is introduced. The class thus obtained depends on a set of three parameters, incorporating an additional distribution to the classical lognormal one. This new class of heavy-tailed distribution is presented as an alternative to other useful heavy-tailed distributions, such as the lognormal, Weibull, and Pareto distributions. The density and distribution functions of this new class are given by a closed expression which allows us to easily compute probabilities, quantiles, moments, and related measurements. Finally, some applications are shown as examples.     

On a bivariate Kumaraswamy type exponential distribution

ABSTRACT

This paper considers a class of absolutely continuous bivariate exponential distributions whose univariate margins are the ordinary exponential distributions. We study different mathematical properties of the proposed model. The estimation of the parameters by maximum likelihood is discussed. Application is made to a real data example to illustrate the flexibility of theproposed distribution for data analysis.     

A latent variable model for analyzing mixed longitudinal (k,l)-inflated count and ordinal responses

A random effects model for analyzing mixed longitudinal count and ordinal data is presented where the count response is inflated in two points (k and l) and an (k,l)-Inflated Power series distribution is used as its distribution. A full likelihood-based approach is used to obtain maximum likelihood estimates of parameters of the model. For data with non-ignorable missing values models with probit model for missing mechanism are used.The dependence between longitudinal sequences of responses and inflation parameters are investigated using a random effects approach. Also, to investigate the correlation between mixed ordinal and count responses of each individuals at each time, a shared random effect is used. In order to assess the performance of the model, a simulation study is performed for a case that the count response has (k,l)-Inflated Binomial distribution. Performance comparisons of count-ordinal random effect model, Zero-Inflated ordinal random effects model and (k,l)-Inflated ordinal random effects model are also given. The model is applied to a real social data set from the first two waves of the national longitudinal study of adolescent to adult health (Add Health study). In this data set, the joint responses are the number of days in a month that each individual smoked as the count response and the general health condition of each individual as the ordinal response. For the count response there is incidence of excess values of 0 and 30.     

A generalization of the logistic linear'model

Consider the logistic linear model, with some explanatory variables overlooked. Those explanatory variables may be quantitative or qualitative. In either case, the resulting true response variable is not a binomial or a beta-binomial but a sum of binomials. Hence, standard computer packages for logistic regression can be inappropriate even if an overdispersion factor is incorporated. Therefore, a discrete exponential family assumption is considered to broaden the class of sampling models. Likelihood and Bayesian analyses are discussed. Bayesian computation techniques such as Laplacian approximations and Markov chain simulations are used to compute posterior densities and moments. Approximate conditional distributions are derived and are shown to be accurate. The Markov chain simulations are performed effectively to calculate posterior moments by using the approximate conditional distributions. The methodology is applied to Keeler's hardness of winter wheat data for checking binomial assumptions and to Matsumura's Accounting exams data for detailed likelihood and Bayesian analyses.     

Admissibility of a variance estimator in finite population sampling under Warner's randomized response plan with multiple responses

ABSTRACT

Let P be the proportion of individuals in a finite population possessing a sensitive attribute. We consider the problem of unbiased estimation of (i) the variance of a linear unbiased estimator of P and (ii) the population variance P (1—P) for a given probability sampling design under Warner's (1965 Warner, S.L. (1965). Randomized response - A survey technique for eliminating evasive answer bias. J. Amer. Statist. Assoc. 60:63–69.[Taylor & Francis Online], [Web of Science ®] , [Google Scholar]) randomized response (RR) plan when independent responses are obtained from each sampled individual as many times as he/she is selected in the sample and prove the admissibility of a quadratic unbiased estimator for each.     

Poisson-Driven Stationary Markov Models

ABSTRACT

We propose a simple yet powerful method to construct strictly stationary Markovian models with given but arbitrary invariant distributions. The idea is based on a Poisson-type transform modulating the dependence structure in the model. An appealing feature of our approach is the possibility to control the underlying transition probabilities and, therefore, incorporate them within standard estimation methods. Given the resulting representation of the transition density, a Gibbs sampler algorithm based on the slice method is proposed and implemented. In the discrete-time case, special attention is placed to the class of generalized inverse Gaussian distributions. In the continuous case, we first provide a brief treatment of the class of gamma distributions, and then extend it to cover other invariant distributions, such as the generalized extreme value class. The proposed approach and estimation algorithm are illustrated with real financial datasets. Supplementary materials for this article are available online.     

Some New Results on the Multinomial Randomized Response Model

ABSTRACT

The randomized response technique is an effective survey method designed to elicit sensitive information while ensuring the privacy of the respondents. In this article, we present some new results on the randomization response model in situations wherein one or two response variables are assumed to follow a multinomial distribution. For a single sensitive question, we use the well-known Hopkins randomization device to derive estimates, both under the assumption of truthful and untruthful responses, and present a technique for making pairwise comparisons. When there are two sensitive questions of interest, we derive a Pearson product moment correlation estimator based on the multinomial model assumption. This estimator may be used to quantify the linear relationship between two variables when multinomial response data are observed according to a randomized-response protocol.     

A General Class of Distributions: Properties and Applications

ABSTRACT

In this article, we derive a general class of distributions and establish its relationship to χ² distribution. The proposed class includes normal, inverse Gaussian, lognormal, gamma, Rayleigh, and Maxwell distributions. Various statistical properties of the class are discussed. Some applications of the class are given.     

Sparsistent and constansistent estimation of the varying-coefficient model with a diverging number of predictors

ABSTRACT

The varying-coefficient model is a flexible class of approaches that extends simple linear relationships between covariates and responses. Two related problems concerning these models are selecting relevant variables and determining non-varying coefficients among those relevant ones. In this paper we study the sparsistency and constansistency of the regularized estimation approach when the number of predictors diverges with the sample size. Here, constansistency refers to the desired property that the non-zero, non-varying coefficients are identified with probability tending to one.     

A generalized ordered Probit model

Abstract

The ordered probit and logit models, based on the normal and logistic distributions, can yield biased and inconsistent estimators when the distributions are misspecified. A generalized ordered response model is introduced which can reduce the impact of distributional misspecification. An empirical exploration of various determinants of life satisfaction suggests possible benefits of allowing for diverse distributional characteristics. These improvements are confirmed using a Monte Carlo study to contrast the performance of the flexible parametric specifications to the probit and logit specifications.     

Test of parameter changes in a class of observation-driven models for count time series

Abstract

This paper investigates the parameter-change tests for a class of observation-driven models for count time series. We propose two cumulative sum (CUSUM) test procedures for detection of changes in model parameters. Under regularity conditions, the asymptotic null distributions of the test statistics are established. In addition, the integer-valued generalized autoregressive conditional heteroskedastic (INGARCH) processes with conditional negative binomial distributions are investigated. The developed techniques are examined through simulation studies and also are illustrated using an empirical example.     

A Bayesian random effects model for analyzing mixed negative binomial and continuous longitudinal responses

ABSTRACT

A general Bayesian random effects model for analyzing longitudinal mixed correlated continuous and negative binomial responses with and without missing data is presented. This Bayesian model, given some random effects, uses a normal distribution for the continuous response and a negative binomial distribution for the count response. A Markov Chain Monte Carlo sampling algorithm is described for estimating the posterior distribution of the parameters. This Bayesian model is illustrated by a simulation study. For sensitivity analysis to investigate the change of parameter estimates with respect to the perturbation from missing at random to not missing at random assumption, the use of posterior curvature is proposed. The model is applied to a medical data, obtained from an observational study on women, where the correlated responses are the negative binomial response of joint damage and continuous response of body mass index. The simultaneous effects of some covariates on both responses are also investigated.     

High breakdown point robust estimators with missing data

In this paper, we propose a new procedure to estimate the distribution of a variable y when there are missing data. To compensate the presence of missing responses, it is assumed that a covariate vector x is observed and that y and x are related by means of a semi-parametric regression model. Observed residuals are combined with predicted values to estimate the missing response distribution. Once the responses distribution is consistently estimated, we can estimate any parameter defined through a continuous functional T using a plug in procedure. We prove that the proposed estimators have high breakdown point.     

Scrambled randomized response methods for obtaining sensitive quantitative data

We study a multiplicative randomized response method for obtaining responses to sensitive questions when the answers are quantitative. The method involves the respondent multiplying his sensitive answer by a random number from a known distribution, and giving the product to the interviewer, who does not know the value of the random number and thus receives a scrambled response. Some particular distributions for the random scrambling number are proposed and studied, and ways of generating the scrambling numbers are discussed. Some modifications for increasing the efficiency of the method are proposed, and numerical results are given that show the scrambled response method is generally superior to the previously used method of randomizing questions.     

Analysis of mixed correlated bivariate zero-inflated count and (k,l)-inflated beta responses with application to social network datasets

This paper presents a new model that monitors the basic network formation mechanisms via the attributes through time. It considers the issue of joint modeling of longitudinal inflated (0, 1)-support continuous and inflated count response variables. For joint model of mentioned response variables, a correlated generalized linear mixed model is studied. The fraction response is inflated in two points k and l (k < l) and a k and l inflated beta distribution is introduced to use as its distribution. Also, the count response is inflated in zero and we use some members of zero-inflated power series distributions, hurdle-at-zero, members of zero-inflated double power series distributions and zero-inflated generalized Poisson distribution as our count response distribution. A full likelihood-based approach is used to yield maximum likelihood estimates of the model parameters and the model is applied to a real social network obtained from an observational study where the rate of the ith node’s responsiveness to the jth node and the number of arrows or edges with some specific characteristics from the ith node to the jth node are the correlated inflated (0, 1)-support continuous and inflated count response variables, respectively. The effect of the sender and receiver positions in an office environment on the responses are investigated simultaneously.