qRankParTop {BSTZINB}R Documentation

Bar plot for time-averaged log-q estimates over top ranking counties (descending order)

Description

Produce a descending order of bar plot for time-averaged log-q estimates over top ranking counties

Usage

qRankParTop(state.set,cname,bstfit,vn=12,
                   cex.title=18, cex.lab=18, cex.legend=18)

Arguments

state.set

character vector of set of states on which the the graphics is to be made

cname

character vector of the names of the counties

bstfit

the fitted data for BSTP, BSTNB or BSTZINB

vn

positive integer, number of sample counties to display

cex.title

Positive number to control the size of the text of the main title. Defaults to 18.

cex.lab

Positive number to control the size of the text in the axes labels. Defaults to 18.

cex.legend

Positive number to control the size of the text in the legend. Defaults to 18.

Value

bar graph

Examples

data(simdat)
y <- simdat$y
X <- cbind(simdat$V1,simdat$x)
data(county.adjacency)
data(USAcities)
IAcities <- subset(USAcities,state_id=="IA")
countyname <- unique(IAcities$county_name)
A <- get_adj_mat(county.adjacency,countyname,c("IA"))

res3 <- BSTZINB(y, X, A, LinearT=TRUE, nchain=3, niter=100, nburn=20, nthin=1)
qRankParTop(state.set=c("IA"),cname=countyname,bstfit=res3,vn=12,
             cex.title=18, cex.lab=12, cex.legend=12)



[Package BSTZINB version 1.0.1 Index]