test/inverse_temperatures_choosing_example.R

In kkdey/tmcmcR: tmcmcR : a R/C++ package for Transformation based Markov Chain Monte Carlo

## Select inverse temperatures
library(tmcmcR)
unipdf <- function(x) return (dnorm(x,0,1,log=TRUE))
d <- 1
nsamples <- 5000;

base <- 0

temp <- tmcmcR:::tmcmc_metrop(unipdf,base=base, scale=1,nsamples=5000,burn_in = NULL)$chain;

temp1 <- tmcmcR:::.rand_generate(unipdf, method="TMCMC")

inv_temp_fixed_tmcmc <- tmcmcR:::select_inverse_temp(unipdf, minbeta=0.05, L_iter =50, sim_method="TMCMC", inv_temp_scheme="fixed")
inv_temp_rand_tmcmc <- tmcmcR:::select_inverse_temp(unipdf, minbeta=0.05, L_iter =50, sim_method="TMCMC", inv_temp_scheme="randomized")
inv_temp_fixed_rwmh <- tmcmcR:::select_inverse_temp(unipdf, minbeta=0.05, L_iter =50, sim_method="RWMH", inv_temp_scheme="fixed")
inv_temp_rand_rwmh <- tmcmcR:::select_inverse_temp(unipdf, minbeta=0.05, L_iter =50, sim_method="RWMH", inv_temp_scheme="randomized")

print(inv_temp_fixed_tmcmc)
print(inv_temp_rand_tmcmc)
print(inv_temp_fixed_rwmh)
print(inv_temp_rand_rwmh)