Bliss_Gibbs_Sampler {bliss} | R Documentation |
Bliss_Gibbs_Sampler
Description
A Gibbs Sampler algorithm to sample the posterior distribution of the Bliss model.
Usage
Bliss_Gibbs_Sampler(data, param, verbose = FALSE)
Arguments
data |
a list containing:
|
param |
a list containing:
|
verbose |
write stuff if TRUE (optional). |
Value
a list containing :
- trace
a matrix, the trace of the Gibbs Sampler.
- param
a list containing parameters used to run the function.
Examples
# 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