bdribs.contour {bdribs}

contour plot - draws plot (optional) and returns a matrix/grid of posterior values

Description

Contour plot of posterior probabilities on a range of (y, E) values

Usage

bdribs.contour(ymax, pyrmax, eincr, tol, k, bg.rate, plt = TRUE, ...)

Arguments

ymax

maximum number of AESI event for which contour plot to be drawn

pyrmax

maximum risk exposure (in patient-year)

eincr

increment of patient-year exposures (default = 50)

tol

the maximum tolerance value of relative risk r (default =1)

k

allocation ratio (T:C)

bg.rate

estimated background rate (historical control rate) per patient-year (using inf.type=1)

plt

whether a contour plot to be drawn (default = TRUE)

...

to supply remaining parameters for bdribs call when supplied will override the default values

Value

returns contour plot matrix over the grid specified

Examples

## Sample calls
     #run 1: The contour plot
     
     bdribs.contour(ymax=15,pyrmax=2000,eincr=250,tol=1.5,k=2, bg.rate=0.0045)
     #run 2: Monitoring blinded AE over time using contour plot
     bdribs.contour(ymax=15,pyrmax=2000,eincr=250,tol=1.5,k=2, bg.rate=0.0045)
     obs.pyr=c(300,570,650,800, 1200, 1500)
     obs.y=c(2,4,5,6,10,12)
     points(obs.pyr, obs.y,type="p",pch=16, cex=1.4,col="maroon")
     if (length(obs.y)>1) points(c(0,obs.pyr), c(0,obs.y), type="s", lty=3, lwd=2,
     col="black")