plot.singletable: Plot Method for 'singletable' objects

Description Usage Arguments Details References See Also Examples

View source: R/singleplot.r

Description

Produces various plots for single table analysis.

Usage

1
2
3
    ## S3 method for class 'singletable'
plot(x,type=type,file=NULL,select=c(1,2)
	,xlab=NULL,ylab=NULL,addline=NULL,xlim=NULL,ylim=NULL,...)

Arguments

x

an object inheriting from class singletable.

type

a chracter string specifying the type of plots to produce. Options are sidebyside and overlap. See details

file

a character string specifying the filename as which the plots are saved. Default is NULL to view on screen. See details.

select

a numeric value or vector specifying which distribution should be plotted. select=1 is posterior distribution. select=2 is prior distribution. select=c(1,2) is both. Default is c(1,2). This argument is only used when type="sidebyside".

xlab

a character string specifying the x-axis label in the plot. Default is the name of the measure of association

ylab

a character string specifying the x-axis label in the plot. Default is "Density"

addline

a numeric value specifying the x-value for a vertical reference line at x=addline. Default is NULL

xlim, ylim

a numeric vectors of length 2 specifying the lower and upper limits of the axes

...

Other arguments can be passed to plot function

Details

If type="sidebyside", the posterior distribution of measure and the prior distribution are drawn side by side in two plots. If type="overlap", the posterior distribution of measure and the prior distribution are overlaid in one plot.

If file=NULL, the plots will be displayed on screen. Or else, the plots will be saved as "./mmeta/codefile.pdf", where "./" denotes current working directory.

References

Luo, S., Chen, Y., Su, X., Chu, H., (2014). mmeta: An R Package for Multivariate Meta-Analysis. Journal of Statistical Software, 56(11), 1-26.

Chen, Y., Luo, S., (2011a). A Few Remarks on "Statistical Distribution of the Difference of Two Proportions' by Nadarajah and Kotz, Statistics in Medicine 2007; 26(18):3518-3523" . Statistics in Medicine, 30(15), 1913-1915.

Chen, Y., Chu, H., Luo, S., Nie, L., and Chen, S. (2014a). Bayesian analysis on meta-analysis of case-control studies accounting for within-study correlation. Statistical Methods in Medical Research, doi: 10.1177/0962280211430889. In press.

Chen, Y., Luo, S., Chu, H., Su, X., and Nie, L. (2014b). An empirical Bayes method for multivariate meta-analysis with an application in clinical trials. Communication in Statistics: Theory and Methods. In press.

Chen, Y., Luo, S., Chu, H., Wei, P. (2013). Bayesian inference on risk differences: an application to multivariate meta-analysis of adverse events in clinical trials. Statistics in Biopharmaceutical Research, 5(2), 142-155.

See Also

singletable

Examples

 1
 2
 3
 4
 5
 6
 7
 8
 9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
# Inference under Jeffreys prior distribution
#single.OR.Jeffreys <- singletable(a1=0.5, b1=0.5, a2=0.5,
#                                  b2=0.5, y1=40, n1=96, y2=49, n2=109,
#                                  model="Independent",
#                                  measure="OR", method="exact")
#summary(single.OR.Jeffreys)

# Inference under Laplace prior distribution
#single.OR.Laplace <- singletable(a1=1, b1=1, a2=1, b2=1,
#                                 y1=40, n1=96, y2=49, n2=109,
#                                 model="Independent", measure="OR",
#                                 method="exact")
# Inference under Sarmanov prior distribution with positive correlation
#single.OR.Sar1 <- singletable(a1=0.5, b1=0.5, a2=0.5, b2=0.5,
#                              rho=0.5, y1=40, n1=96, y2=49, n2=109,
#                              model="Sarmanov",
#                              measure="OR", method="exact")
# Inference under Sarmanov prior distribution with negative correlation
#single.OR.Sar2 <- singletable(a1=0.5, b1=0.5, a2=0.5, b2=0.5,
#                              rho=-0.5, y1=40, n1=96, y2=49, n2=109,
#                              model="Sarmanov",
#                              measure="OR", method="exact")
# generate a 2X2 panel plot
#par(mfrow=c(2,2))
#plot(single.OR.Jeffreys, type="overlap", xlim=c(0.5, 2),
#    main="Jefferys Prior")
#plot(single.OR.Laplace, type="overlap", xlim=c(0.5, 2),
 #    main="Laplace Prior")
#plot(single.OR.Sar1, type="overlap", xlim=c(0.5, 2),
#     main=expression(paste("Sarmanov Prior ",rho," = 0.5")))
#plot(single.OR.Sar2, type="overlap", xlim=c(0.5, 2),
#     main=expression(paste("Sarmanov Prior ",rho," = -0.5")))

Example output



mmeta documentation built on May 2, 2019, 6:47 a.m.