############################# PARP2_summary ###############################################
data(PARP2_summary)
Breast <- PARP2_summary$Breast
Thyroid <- PARP2_summary$Thyroid
Betah <- list(Breast$beta, Thyroid$beta)
Sigmah <- list(diag(Breast$se), diag(Thyroid$se))
genename <- "PARP2"
snpnames <- Breast$snp
K <- 2
m <- 6
set.seed(123)
RES <- HS(Betah, Sigmah,
kappa0 = 0.5, kappastar0 = 0.5, sigma20 = 1, s20 = 1,
m = m, K = K, niter = 800, burnin = 400, nthin = 1, nchains = 1,
a1 = 0.1, a2 = 0.1, d1 = 0.1, d2 = 0.1, c1 = 1, c2 = 1, e2 = 1, snpnames, genename
)
summaryHS(RES)
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.