R Graphical Manual

Browse All

Last data update: 2014.03.03

R: Plot the summary of the bivariate model from Riley et al....

plot.riley

R Documentation

Plot the summary of the bivariate model from Riley et al. (2008).

Description

This function plots the summary sensitivity and false positive rate with their corresponding confidence regions.

Usage

## S3 method for class 'riley'
plot(x, plotsumm = TRUE, plotnumerics = TRUE, level = 0.95, 
main="", ylim = c(0,1), xlim = c(0,1), pch = 1, lty = 1, lwd = 1, 
cex.numerics=0.45, add=FALSE, ...)

Arguments

`x`	a `riley` object.
`plotsumm`	logical, should the plot draw the summary pair of sensitivity and false positive rate?
`plotnumerics`	logical, should the plot contain a summary table of sensitivity and false positive rate?
`level`	numeric, the level for calculations of confidence intervals
`main`	string, title of the plot
`ylim`	numeric of length 2, which section of the sensitivities to plot?
`xlim`	numeric of length 2, which section of the false positive rates to plot?
`pch`	integer, symbol for the pair of mean sensitivity and false positive rate
`lty`	integer, line type of confidence curve
`lwd`	integer, line width of the confidence curve
`cex.numerics`	numeric, text size
`add`	logical, should the confidence region be added to the current plot?
`...`	arguments to be passed on to other functions

Author(s)

Thomas Debray <thomas.debray@gmail.com>

References

Riley RD, Thompson JR, Abrams KR. An alternative model for bivariate random-effects meta-analysis when the within-study correlations are unknown. Biostatistics 2008; 9: 172–186.

Examples

data(Scheidler)

ds1 <- Scheidler[which(Scheidler$modality==1),]
ds2 <- Scheidler[which(Scheidler$modality==2),]
ds3 <- Scheidler[which(Scheidler$modality==3),]

#Perform the analyses
fit1 <- riley(ds1, type="test.accuracy")
fit2 <- riley(ds2, type="test.accuracy")
fit3 <- riley(ds3, type="test.accuracy")

plot(fit1,plotnumerics=FALSE,pch=0) #CT
plot(fit2,plotnumerics=FALSE,add=TRUE,pch=1) #LAG
plot(fit3,plotnumerics=FALSE,add=TRUE,pch=2) #MRI