# 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 = 50, nChains = 2, gammaPrior = "hotspot",
hyperpar = hyperpar, tmpFolder = "tmp/" )
## check output
# show the graph relationship between responses
Gy <- getEstimator(fit, estimator = "Gy")
plot(Gy)
# }
Run the code above in your browser using DataCamp Workspace