PredictiveThrower.cpp File Reference

#include "PredictiveThrower.h"
#include "Parameters/ParameterHandlerGeneric.h"
#include "TH3.h"

Include dependency graph for PredictiveThrower.cpp:

Go to the source code of this file.

Functions
bool	CheckBounds (const std::vector< const M3::float_t * > &BoundValuePointer, const std::vector< std::pair< double, double >> &ParamBounds)
void	AccumulateWAICToy (const double neg_LLH_temp, double &mean_llh, double &mean_llh_squared, double &sum_exp_llh)
void	AccumulateWAICBin (double &mean_llh, double &mean_llh_squared, double &sum_exp_llh, const unsigned int Ntoys, double &lppd, double &p_WAIC)

Function Documentation

AccumulateWAICBin()

void AccumulateWAICBin	(	double &	mean_llh,
		double &	mean_llh_squared,
		double &	sum_exp_llh,
		const unsigned int	Ntoys,
		double &	lppd,
		double &	p_WAIC
	)

Definition at line 1932 of file PredictiveThrower.cpp.

                                                                               {
// ****************
  // Compute the mean log-likelihood and the squared mean
  mean_llh /= Ntoys;
  mean_llh_squared /= Ntoys;
  sum_exp_llh /= Ntoys;
  sum_exp_llh = std::log(sum_exp_llh);
 
  // Log pointwise predictive density based on Eq. 4 in Gelman2014
  lppd += sum_exp_llh;
 
  // Compute the effective number of parameters for WAIC
  p_WAIC += mean_llh_squared - (mean_llh * mean_llh);
}

AccumulateWAICToy()

void AccumulateWAICToy	(	const double	neg_LLH_temp,
		double &	mean_llh,
		double &	mean_llh_squared,
		double &	sum_exp_llh
	)

Definition at line 1917 of file PredictiveThrower.cpp.

                                            {
// ****************
  // Negate the negative log-likelihood to get the actual log-likelihood
  double LLH_temp = -neg_LLH_temp;
 
  mean_llh += LLH_temp;
  mean_llh_squared += LLH_temp * LLH_temp;
  sum_exp_llh += std::exp(LLH_temp);
}

CheckBounds()

bool CheckBounds	(	const std::vector< const M3::float_t * > &	BoundValuePointer,
		const std::vector< std::pair< double, double >> &	ParamBounds
	)

Definition at line 460 of file PredictiveThrower.cpp.

                                                                       {
// *************************
  for (size_t i = 0; i < BoundValuePointer.size(); ++i) {
    const double val = *(BoundValuePointer[i]);
    const double minVal = ParamBounds[i].first;
    const double maxVal = ParamBounds[i].second;
 
    if (val < minVal || val > maxVal)
      return false; // out of bounds
  }
  return true; // all values are within bounds
}