R/tempExchange.R

In eifer4/stochasticSampling: Stochastic partitioning sampler.

# tempExchange <- function(temps, N, X, Y, comp, priors) {
#   R <- length(temps)
#   newTemps <- sample(temps, R)
#   ldens <- numeric(R)
#   
#   h <- priors$h
#   h0 <- priors$h0
#   beta0 <- priors$beta0
#   alpha0 <- priors$alpha0
#   nu <- prios$alpha0
#   sigma02 <- priors$sigma02
#   
#   m <- n <- NULL
#   
#   for(r in 1:R){
#     for(k in 1:length(comp[[r]])){
#       m <- length(comp[[r]][[k]]$X)
#       n <- length(comp[[r]][[k]]$Y)
#       dens[r] <- dens[r] + min(0,(1/temps[r] - 1/newTemps[r]) * logCompDens(m,n, h, h0, beta0, alpha0, nu, sigma02, N, X, Y))
#     } 
#   }
#   U <- unif(1)
#   if(U<=exp(dens)){
#     temps <- newTemps
#   }
#   
#   return(temps)
# }