Fixed wrong global in EstimParam

diyabc · Dec 31, 2019 · 80d32f5 · 80d32f5
1 parent 17ceab7
commit 80d32f5
Showing 1 changed file with 2 additions and 1 deletion.
diff --git a/src/EstimParam.cpp b/src/EstimParam.cpp
@@ -230,7 +230,8 @@ EstimParamResults EstimParam_fun(Reftable &myread,
         size_t j = p.second;
         std::vector<double> quants = forestQuantiles(obs,preds[5][j],probs);
         auto reality = y(p.first);
-        auto diff = expectation - reality;
+        auto diff = preds[0][0][p.first] - reality;
+        // auto diff = expectation - reality;
         auto sqdiff = diff * diff;
         auto CI = quants[2] - quants[0];
         MSE += sqdiff;