From 85788254a88f16938bacd8c844bf8717a172d522 Mon Sep 17 00:00:00 2001
From: Achilleas <anp4r@neverland>
Date: Tue, 24 Jul 2018 09:20:53 -0400
Subject: [PATCH] Use same segpi for cross validation iterations

---
 src/SNPs.cpp | 4 ++++
 1 file changed, 4 insertions(+)

diff --git a/src/SNPs.cpp b/src/SNPs.cpp
index 4e45b9f..7dc8ad3 100644
--- a/src/SNPs.cpp
+++ b/src/SNPs.cpp
@@ -742,10 +742,14 @@ double SNPs::cross10(bool penalize){
 	// L = 1/10* \sum_i L*(i)
 	// where L*(i) is the likelihood of data in group i after optimizing model without it
 	//
+        double initial_segpi = segpi; // Store local segpi
 	double toreturn =0;
 	vector< set<int> > split10 = make_cross10();
 	vector<double> Lstar;
 	for (vector<set<int> >::iterator it = split10.begin(); it != split10.end(); it++){
+                // Ensure the use of local segpi and not the estimated one from the previous
+                // cross validation
+                segpi = initial_segpi;  
 		GSL_xv_optim(*it, penalize);
 		double tmpl = 0;
 		for (set<int>::iterator it2 = it->begin(); it2 != it->end(); it2++) tmpl += llk(*it2);