R/betadraw_rmse.R

In scalablebayesm: Distributed Markov Chain Monte Carlo for Bayesian Inference in Marketing

betadraw_rmse = function(betadraw, beta_true) {
  rmse = c()
  for (i in 1:dim(betadraw)[1]) #iterate over n
  {
    squared_error = (beta_true[i, ] - betadraw[i,,])^2 # 5x80
    sum_squared_error = colSums(squared_error)
    rmse[i] = sqrt(mean(sum_squared_error))
  }
  return(mean(rmse))
}