blackbox is a function that takes a matrix of survey data in which individuals
place themselves on continuous scales across multiple issues, and locates those
citizens in a spatial model of voting. Mathematically, this function generalizes
the singular value of a matrix to cases in which there is missing data in the
matrix. Scales generated using perceptual data (i.e. scales of legislator locations
using liberal-conservative rankings by survey respondents) should instead use
the blackbox_transpose function in this package instead.
blackbox(data,missing=NULL,verbose=FALSE,dims=1,minscale)aldmck should print out detailed
output when scaling the data.blackbox.NNumber of respondents who provided a response to this stimulus.
cStimulus intercept.
w1Estimate of the stimulus weight on the first dimension. If viewing the
results for a higher dimension, higher dimension results will appear as w2, w3, etc.
R2The percent variance explained for the stimulus. This increases as
more dimensions are estimated.
c1Estimate of the individual intercept on the first dimension. If viewing the
results for a higher dimension, higher dimension results will appear as c2, c3, etc.
SSESum of squared errors.
SSE.explainedExplained sum of squared error.
percentPercentage of total variance explained.
SEStandard error of the estimate, with formula provided on pg. 973 of the article cited below.
singularSingluar value for the dimension.
Keith Poole, Jeffrey Lewis, Howard Rosenthal, James Lo, Royce Carroll (2016) ``Recovering a Basic Space from Issue Scales in R.'' Journal of Statistical Software. 69(7), 1--21. doi:10.18637/jss.v069.i07
Keith T. Poole (1998) ``Recovering a Basic Space From a Set of Issue Scales.'' American Journal of Political Science. 42(3), 954-993.
### Loads issue scales from the 1980 NES.
data(Issues1980)
Issues1980[Issues1980[,"abortion1"]==7,"abortion1"] <- 8 #missing recode
Issues1980[Issues1980[,"abortion2"]==7,"abortion2"] <- 8 #missing recode
### This command conducts estimates, which we instead load using data()
# Issues1980_bb <- blackbox(Issues1980,missing=c(0,8,9),verbose=FALSE,dims=3,minscale=8)
data(Issues1980_bb)
summary(Issues1980_bb)
Run the code above in your browser using DataLab