DianaPerri2 {RRTCS}R Documentation

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 nn (the sample size)

mu

vector with the means of the scramble variables WW and UU

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

Details

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

To estimate Yˉ\bar{Y} a sample of respondents is selected according to simple random sampling with replacement. The transformed variable is

ri=ziβμWμU(1β)μWr_i=\frac{z_i-\beta\mu_W\mu_U}{(1-\beta)\mu_W}

where μW,μU\mu_W,\mu_U are the means of W,UW,U scramble variables, respectively.

The estimated variance in this model is

V^(Yˉ^R)=sz2n(1β)2μW2\widehat{V}(\widehat{\bar{Y}}_R)=\frac{s_z^2}{n(1-\beta)^2\mu_W^2}

where sz2=i=1n(zizˉ)2n1s_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

V^(Yˉ^R)=sz2n(1β)2μW2(1nN)\widehat{V}(\widehat{\bar{Y}}_R)=\frac{s_z^2}{n(1-\beta)^2\mu_W^2}\left(1-\frac{n}{N}\right)

Value

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

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

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")

[Package RRTCS version 0.0.4 Index]