tracex {pscl} | R Documentation |
Produces a trace plot of the MCMC samples from the posterior density of
legislators' ideal
points.
tracex(object, legis=NULL, d=1, conf.int=0.95, showAll = FALSE, burnin=NULL)
object |
an object of class ideal . |
legis |
a vector of either the names of legislators (or
partial matches of the names as given in the
dimnames of object$x . |
d |
numeric, either a scalar or a vector of length two, the dimension(s) to be traced. |
conf.int |
numeric, the level of the confidence interval on the posterior mean to be plotted. |
showAll |
logical, if TRUE and length(d)==2 ,
display traces for all selected legislators' ideal points on
the one plot. |
burnin |
of the recorded MCMC samples, how many to discard as
burnin? Default is NULL , in which case the value of
burnin in the ideal object is used. |
Produces a trace plot showing the history of the MCMC
iterations for the ideal point of each of the legislators (partially)
named in legis
. For d=1
, each trace plot includes a
trace over iterations, the cumulative mean, a moving average, the
MCMC-based estimate of the mean of the posterior, and a confidence
interval (specified by conf.int
) around the mean of the
posterior (using the estimated quantiles
) of
the respective MCMC iterates). All of these values are calculated
discarding the initial burnin
iterations.
When d
is a vector of length two, a 2-dimensional trace plot is
displayed, with the d[1]
dimension on the horizontal axis, and
the d[2]
dimension on the vertical axis.
ideal
; pmatch
for matching
legislators' names.
data(s109) ## short run for demo purposes only id1 <- ideal(s109, meanzero=TRUE, maxiter=500, burnin=100,thin=10) tracex(id1,legis="KENNEDY") ## n.b., no such legislator has Horrendous Goblin tracex(id1,legis=c("KENNEDY","BOXER","KYL","Horrendous Goblin")) ## Not run: id2 <- ideal(s109, d=2, maxiter=5000, ## unidentified! burnin=0, thin=50) tracex(id2,d=1,legis=c("KENNEDY","BOXER","KYL","Horrendous Goblin")) tracex(id2,d=2,legis=c("KENNEDY","BOXER","KYL","Horrendous Goblin")) tracex(id2,d=1:2, legis=c("KENNEDY","BOXER","KYL","Horrendous Goblin")) ## partial matching tracex(id2,d=1:2, legis=c("KENN","BOX","BID","SNO","SPEC","MCCA","KYL", "Horrendous Goblin"), showAll=TRUE) ## End(Not run)