# NOT RUN {
library(Rfast)
library(mvtnorm)
library(plyr)
head(dataset1)
BP <- MVN_BayesianPosteriori(dataset1)
BP
BP_Gibbs <- MVN_GibbsSampler(5000, BP)
colMeans(BP_Gibbs)
colrange(BP_Gibbs)
result <- MVN_MCMC(BP, 5000, c(1), c(77.03))
result$Accept
# }
Run the code above in your browser using DataLab