Description Usage Arguments Details Value References Examples
Function to calculate predictions and uncertainties of predictions from estimates from multivariate regression analysis of survey data with the item count technique, using predicted responses to list experiments as predictors in outcome regressions.
1 2 3 4 5 
object 
Object of class inheriting from "ictreg.joint" 
newdata 
An optional data frame containing data that will be used to make predictions from. If omitted, the data used to fit the regression are used. 
newdata.diff 
An optional data frame used to compare predictions with predictions from the data in the provided newdata data frame. 
se.fit 
A switch indicating if standard errors are required. 
interval 
Type of interval calculation. 
level 
Significance level for confidence intervals. 
avg 
A switch indicating if the mean prediction and associated statistics across all obserations in the dataframe will be returned instead of predictions for each observation. 
sensitive.value 
Userspecified value for the sensitive item. 
sensitive.diff 
A switch indicating if the difference in predictions when the sensitive item = 1 and when the sensitive item = 0 is calculated. 
return.draws 
A switch indicating if the draws from the simulations used to generate predictions will be returned. 
predict.sensitive 
A switch indicating whether predictions from the sensitive item model are generated. 
... 
further arguments to be passed to or from other methods. 
predict.ictreg.joint
produces predicted values, obtained by
evaluating the regression function in the frame newdata (which defaults to
model.frame(object)). By using sensitive.value, users must set the value of
z – the latent response to the sensitive item – to be either zero or one,
depending on the prediction that the user requires.
Two additional types of mean prediction are also available. The first, if a newdata.diff data frame is provided by the user, calculates the mean predicted values across two datasets, as well as the mean difference in predicted value. Standard errors and confidence intervals are also added. For newdata.diff predictions, sensitive.value must be set to 1 or 0, not "both" (and sensitive.diff must also be set to FALSE). Users may also set the logical sensitive.diff to TRUE and sensitive.value to "both", which will output the mean predicted values across all observations for z = 0 as well as z = 1, in addition to the mean difference in predicted value. Standard errors and confidence intervals are also added. For difference predictions (sensitive.diff and newdata.diff), the option avg must be set to TRUE.
Users can also use the predict.sensitive = TRUE option to generate predictions of responses to the sensitive item, with standard errors and confidence intervals.
NOTE: In order to generate predictions from userprovided data frames
(newdata and newdata.diff), users MUST run models using ictreg.joint
on data that does not contain any missingness. Further, the data frames
provided to predict.ictreg.joint
must also not contain any
missingness.
predict.ictreg.joint
produces a vector of predictions or a
matrix of predictions and bounds with column names fit, lwr, and upr if
interval is set. If sensitive.value = "both", predict.ictreg.joint
will produce a list, where the first element corresponds to when the
sensitive item = 0 and the second element corresponds to when the sensitive
item = 1. If sensitive.diff = TRUE, the third element in the list
corresponds to the difference (sensitive = 0 subtracted from sensitive = 1).
If se.fit is TRUE, a list with the following components is returned:
fit 
vector or matrix as above. 
se.fit 
standard error of prediction(s) 
If return.draws is TRUE, the list includes
draws.predict 
A matrix of draws from a multivariate normal
distribution with mean equal to the vector of estimated coefficients from
the outcome regression model, and sigma equal to the variancecovariance
matrix of the outcome regression model. Rows are observations; colums are
10,000 draws. If sensitive.value = both, will be a list of two elements
where each element is a matrix as described; the first matrix will be for
when the sensitive item = 0, the second matrix will be for when the
sensitive item = 1. If newdata.diff is provided, 
draws.mean 
The

sens.diff 
If
sensitive.diff = TRUE, a vector of 10,000 draws generated from subtracting
the first item in 
If predict.sensitive = TRUE, the list also includes
fitsens 
a vector of predictions and bounds with column names fit, lwr, and upr if interval is set, generated from the sensitive item model. 
draws.predict.sens 
A matrix of draws from a multivariate normal distribution with mean equal to the vector of estimated coefficients from the sensitive item model, and sigma equal to the variancecovariance matrix of the sensitive item model. Rows are observations; colums are 10,000 draws (only returned if return.draws is TRUE). If newdata.diff is provided, this will be a list of two matrices as described. The first will correspond to newdata, and the second to newdata.diff. 
draws.mean.sens 
The

Imai, Kosuke, Bethany Park, and Kenneth F. Greene. (2014) “Using the Predicted Responses from List Experiments as Explanatory Variables in Regression Models.” available at http://imai.princeton.edu/research/files/listExp.pdf
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58  data(mexico)
loyal < mexico[mexico$mex.loyal == 1,]
notloyal < mexico[mexico$mex.loyal == 0,]
## Not run:
## Logistic outcome regression
## (effect of voteselling on turnout)
## This replicates Table 4 in Imai et al. 2014
loyalreg < ictreg.joint(formula = mex.y.all ~ mex.male + mex.age + mex.age2 + mex.education +
mex.interest + mex.married +
mex.wealth + mex.urban + mex.havepropoganda + mex.concurrent, data = loyal,
treat = "mex.t", outcome = "mex.votecard", J = 3, constrained = TRUE,
outcome.reg = "logistic", maxIter = 1000)
## Linear outcome regression
## (effect of voteselling on candidate approval)
## This replicates Table 5 in Imai et al. 2014
approvalreg < ictreg.joint(formula = mex.y.all ~ mex.male + mex.age + mex.age2 +
mex.education +
mex.interest + mex.married +
mex.urban +
mex.cleanelections + mex.cleanelectionsmiss +
mex.havepropoganda +
mex.wealth + mex.northregion +
mex.centralregion + mex.metro + mex.pidpriw2 +
mex.pidpanw2 + mex.pidprdw2,
data = mexico, treat = "mex.t", outcome = "mex.epnapprove",
J = 3, constrained = TRUE,
outcome.reg = "linear", maxIter = 1000)
summary(approvalreg)
## Generate predicted probability of turnout, averaged over the whole sample,
## for vote sellers (z = 1), nonvote sellers (z = 0), and the difference
## between vote sellers and nonvote sellers, in the sample of party supporters.
## This replicates the results in the righthand panel of Figure 2 in Imai et al. 2014
loyalpred < predict.ictreg.joint(loyalreg, se.fit = TRUE, interval = "confidence",
level = 0.95, avg = TRUE,
sensitive.value = "both",
sensitive.diff = TRUE, return.draws = TRUE,
predict.sensitive = TRUE)
loyalpred$fit
## View predicted probability of vote selling, in the sample of party supporters.
## This replicates the results in the lefthand panel of Figure 2 in Imai et al. 2014
loyalpred$fitsens
## End(Not run)

Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.