BayesSUR (version 1.1-2)

plot.BayesSUR: create a selection of plots for a "BayesSUR" class object

Description

Convenience function to create a selection of plots for a "BayesSUR" class object. They are plots of estimators, response graph, network, manhattan and MCMC diagnosis indexed by numbers 1:5.

Usage

# S3 method for BayesSUR
plot(x, which = c(1L:4L), ...)

Arguments

x

an object of class "BayesSUR".

which

if a subset of the plots is required, specify a subset of the numbers 1:5 which are plots of estimators, response graph, network, manhattan and MCMC diagnosis, respectively. Default is c(1L:4L) Only c(1,4,5) is valid for the HRR models.

...

other arguments

Examples

# NOT RUN {
data("example_eQTL", package = "BayesSUR")
hyperpar = list( a_w = 2 , b_w = 5 )

set.seed(9173)
fit <- BayesSUR(Y = example_eQTL[["blockList"]][[1]], 
                X = example_eQTL[["blockList"]][[2]],
                data = example_eQTL[["data"]], outFilePath = tempdir(),
                nIter = 100, burnin = 0, nChains = 2, gammaPrior = "hotspot",
                hyperpar = hyperpar, tmpFolder = "tmp/" )

## check output
# Show the interactive plots. Note that it needs at least 2000*(nbloc+1) iterations 
# for the diagnosis plots where nbloc=3 by default 
# }
# NOT RUN {
plot(fit)
# }
# NOT RUN {
# }