tracex {pscl}R Documentation

trace plot of MCMC iterates, posterior density of legislators' ideal points

Description

Produces a trace plot of the MCMC samples from the posterior density of legislators' ideal points.

Usage

tracex(object, legis=NULL, d=1, conf.int=0.95,
       showAll = FALSE, burnin=NULL)

Arguments

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.

Details

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.

See Also

ideal; pmatch for matching legislators' names.

Examples

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)

[Package pscl version 0.61 Index]