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

Description Usage Arguments Examples

View source: R/plot.BayesSUR.R

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

1
2
## S3 method for class '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

 1
 2
 3
 4
 5
 6
 7
 8
 9
10
11
12
13
14
data("example_eQTL", package = "BayesSUR")
hyperpar = list( a_w = 2 , b_w = 5 )

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 

plot(fit)

BayesSUR documentation built on Dec. 18, 2019, 1:42 a.m.