R/updateParam.R

In eifer4/stochasticSampling: Stochastic partitioning sampler.

updateParam <- function(comp, priors, invariantData, compCount, temp=1) {
  nu <- priors$nu
  sigma02 <- priors$sigma02
  
  h <- priors$h
  h0 <- priors$h0
  
  beta0 <- priors$beta0
  alpha0 <- priors$alpha0
  rho <- priors$rho

  move1out   <- move1(comp=comp, priors=priors, invariantData=invariantData, 
                      compCount = compCount, temp= temp)
  move2out   <- move2(comp=move1out$comp, priors = priors, invariantData = invariantData,
                      compCount = move1out$compCount, temp = temp)
  
  
  accept <- c(move1out$accept, move2out$accept)
  
  
  return(list(comp=move2out$comp, accept=accept, compCount = move2out$compCount))
}