New corrections

Author: theogf

src/Predictions.jl

@@ -266,7 +266,7 @@ function studentpredictprobamc(model::GPModel,X_test;nSamples=100)
     st = TDist(model.ν)
     temp_array = zeros(Float64,nSamples)
     for i in 1:nTest
-        if cov_f[i] <= min_cov
+        if cov_f[i] <= 1e-3
             pyf =  LocationScale(m_f[i],1.0,st)
             for j in 1:nSamples
                 temp_array[j] = rand(pyf)

src/functions/PerturbativeCorrection.jl

@@ -54,7 +54,7 @@ function correctedpredict(model,X_test)
         f=function(x)
             return logit(x)*pdf(d,x)
         end
-        predic[i] = quadgk(f,meanfstar_corr[i]-10*covfstar_corr[i],meanfstar_corr[i]+10*covfstar_corr[i])[1]
+        predic[i] = expectation(logit,d)
     end
     return predic,meanfstar_corr,covfstar_corr
 end