View source: R/Bliss_Gibbs_Sampler.R
Bliss_Gibbs_Sampler | R Documentation |
A Gibbs Sampler algorithm to sample the posterior distribution of the Bliss model.
Bliss_Gibbs_Sampler(data, param, verbose = FALSE)
data |
a list containing:
|
param |
a list containing:
|
verbose |
write stuff if TRUE (optional). |
a list containing :
a matrix, the trace of the Gibbs Sampler.
a list containing parameters used to run the function.
# May take a while param_sim <- list(Q=1,n=25,p=50,grids_lim=list(c(0,1)),iter=1e4,K=2) data_sim <- sim(param_sim,verbose=FALSE) res_Bliss_Gibbs_Sampler <- Bliss_Gibbs_Sampler(data_sim,param_sim) theta_1 <- res_Bliss_Gibbs_Sampler$trace[1,] theta_1 # Resultat for few iterations param_sim <- list(Q=1,n=25,p=50,grids_lim=list(c(0,1)),iter=5e2,K=2) data_sim <- sim(param_sim,verbose=FALSE) res_Bliss_Gibbs_Sampler <- Bliss_Gibbs_Sampler(data_sim,param_sim) theta_1 <- res_Bliss_Gibbs_Sampler$trace[1,] theta_1
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.