Пример #1
0
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)
}
// responsibility depends on the relative global probability of the hypotheses and the likelihoods of the observations under each hypothesis
// divide the global probabilities into "typical" data points and outliers
// divide the variant probabilities into each hypothesis (summing to 1)
// treat the 2 hypothesis case to start with
void CrossHypotheses::UpdateResponsibility(const vector<float > &hyp_prob, float typical_prob) {

    if (!success) {
        //cout << "alert: fail to splice still called" << endl;
        InitializeResponsibility();
    } else {
        //  vector<double> tmp_prob(3);
        tmp_prob_d[0] = (1.0f-typical_prob)*scaled_likelihood[0];   // i'm an outlier
        for (unsigned int i_hyp=1; i_hyp<scaled_likelihood.size(); i_hyp++)
            tmp_prob_d[i_hyp] = typical_prob * hyp_prob[i_hyp-1] * scaled_likelihood[i_hyp];

        double ll_denom = 0.0;
        for (unsigned int i_hyp=0; i_hyp<scaled_likelihood.size(); i_hyp++) {
            ll_denom += tmp_prob_d[i_hyp];
        }

        for (unsigned int i_hyp=0; i_hyp<responsibility.size(); i_hyp++)
            responsibility[i_hyp] = tmp_prob_d[i_hyp]/ll_denom;
    }
}