Learn R Programming

RRTCS (version 0.0.4)

DianaPerri2: Diana-Perri-2 model

Description

Computes the randomized response estimation, its variance estimation and its confidence interval through the Diana-Perri-2 model. The function can also return the transformed variable. The Diana-Perri-2 model was proposed by Diana and Perri (2010, page 1879).

Usage

DianaPerri2(z,mu,beta,pi,type=c("total","mean"),cl,N=NULL,method="srswr")

Arguments

z

vector of the observed variable; its length is equal to \(n\) (the sample size)

mu

vector with the means of the scramble variables \(W\) and \(U\)

beta

the constant of weighting

pi

vector of the first-order inclusion probabilities

type

the estimator type: total or mean

cl

confidence level

N

size of the population. By default it is NULL

method

method used to draw the sample: srswr or srswor. By default it is srswr

Value

Point and confidence estimates of the sensitive characteristics using the Diana-Perri-2 model. The transformed variable is also reported, if required.

Details

In the Diana-Perri-2 model, each respondent is asked to report the scrambled response \(z_i=W(\beta U+(1-\beta)y_i)\) where \(\beta \in [0,1)\) is a suitable constant controlled by the researcher and \(W,U\) are scramble variables whose distribution is assumed to be known.

To estimate \(\bar{Y}\) a sample of respondents is selected according to simple random sampling with replacement. The transformed variable is $$r_i=\frac{z_i-\beta\mu_W\mu_U}{(1-\beta)\mu_W}$$ where \(\mu_W,\mu_U\) are the means of \(W,U\) scramble variables, respectively.

The estimated variance in this model is $$\widehat{V}(\widehat{\bar{Y}}_R)=\frac{s_z^2}{n(1-\beta)^2\mu_W^2}$$ where \(s_z^2=\sum_{i=1}^n\frac{(z_i-\bar{z})^2}{n-1}\).

If the sample is selected by simple random sampling without replacement, the estimated variance is $$\widehat{V}(\widehat{\bar{Y}}_R)=\frac{s_z^2}{n(1-\beta)^2\mu_W^2}\left(1-\frac{n}{N}\right)$$

References

Diana, G., Perri, P.F. (2010). New scrambled response models for estimating the mean of a sensitive quantitative character. Journal of Applied Statistics 37 (11), 1875-1890.

See Also

DianaPerri2Data

DianaPerri1

ResamplingVariance

Examples

Run this code
# NOT RUN {
N=100000
data(DianaPerri2Data)
dat=with(DianaPerri2Data,data.frame(z,Pi))
beta=0.8
mu=c(50/48,5/3)
cl=0.95
DianaPerri2(dat$z,mu,beta,dat$Pi,"mean",cl,N,"srswor")
# }

Run the code above in your browser using DataLab