heatmap.simulations: Multiple simulations for a heatmap plot

In bardwell/BARD:

Description

Multiple simulations for a heatmap plot

Usage

heatmap.simulations(res, no.draws = 1000)

Arguments

res	A bard result
no.draws	Number of draws to take from posterior distribution

Value

A list of matrices

Examples

set.seed(1)
k_N <- 10
p_N <- 0.1
k_A <- 15
p_A <- 0.3
pi_N <- 0.8
affected_dim <- 8
lower <- 0.3
upper <- 0.7
simparams <- c(k_N, p_N, k_A, p_A, pi_N, affected_dim, lower, upper)
sim <- sim.data(T = 1000, N = 200, simparams)

bardparams <- c(k_N, p_N, k_A, p_A, pi_N, affected_dim)
mu_seq <- c(seq(-0.7, -0.3, by=0.05), seq(0.3, 0.7, by=0.05))
res <- Rbard(sim$data, bardparams, mu_seq, alpha = 1e-4)

hs = heatmap.simulations(res, no.draws = 100)

bardwell/BARD documentation built on Dec. 3, 2019, 12:39 a.m.