void CrossHypotheses::UpdateRelevantLikelihoods() {
for (unsigned int i_hyp=0; i_hyp<basic_likelihoods.size(); i_hyp++) {
for (unsigned int t_flow=0; t_flow<test_flow.size(); t_flow++) {
int j_flow = test_flow[t_flow];
basic_likelihoods[i_hyp][j_flow] = TDistOddN(residuals[i_hyp][j_flow],sigma_estimate[i_hyp][j_flow],skew_estimate,heavy_tailed); // pure observational likelihood depends on residual + current estimated sigma under each hypothesis
}
}
ComputeLogLikelihoods(); // automatically over relevant likelihoods
}
void CrossHypotheses::InitializeDerivedQualities() {
InitializeResponsibility(); // depends on hypotheses
InitializeTmpVariables();
// in theory don't need to compute any but test flows
ComputeBasicResiduals(); // predicted and measured
InitializeSigma(); // depends on predicted
ComputeBasicLikelihoods(); // depends on residuals and sigma
// compute log-likelihoods
ComputeLogLikelihoods(); // depends on test flow(s)
}
void CrossHypotheses::InitializeDerivedQualities() {
InitializeResponsibility(); // depends on hypotheses
// in theory don't need to compute any but test flows
SetModPredictions(); // make sure that mod-predictions=predictions
ComputeBasicResiduals(); // predicted and measured
InitializeSigma(); // depends on predicted
my_t.SetV(heavy_tailed);
ComputeBasicLikelihoods(); // depends on residuals and sigma
// compute log-likelihoods
ComputeLogLikelihoods(); // depends on test flow(s)
}
void CrossHypotheses::UpdateRelevantLikelihoods() {
ComputeBasicLikelihoods();
ComputeLogLikelihoods(); // automatically over relevant likelihoods
}
