From 03d7d2556a9284dc0ac3e155b5c9a8d69b1b21ee Mon Sep 17 00:00:00 2001
From: xiangzhou
Date: Wed, 17 Aug 2016 15:00:35 -0400
Subject: version 0.95alpha

---
 src/gemma.cpp | 13 ++++++-------
 1 file changed, 6 insertions(+), 7 deletions(-)

(limited to 'src')

diff --git a/src/gemma.cpp b/src/gemma.cpp
index 54f1d3c..682835f 100644
--- a/src/gemma.cpp
+++ b/src/gemma.cpp
@@ -42,7 +42,6 @@
 #include "lmm_float.h"  //for LMM class, and functions CalcLambda, CalcPve, CalcVgVe
 #include "mvlmm_float.h"  //for MVLMM class
 #include "prdt_float.h"	//for PRDT class
-#include "varcov_float.h"  //for MVLMM class
 #include "mathfunc_float.h"	//for a few functions
 #else
 #include "io.h"
@@ -53,7 +52,6 @@
 #include "lmm.h"
 #include "mvlmm.h"
 #include "prdt.h"
-#include "varcov.h"
 #include "mathfunc.h"
 #endif
 
@@ -1349,6 +1347,7 @@ void GEMMA::BatchRun (PARAM &cPar)
 		gsl_matrix_free (G);
 	}
 
+	/*
 	//Compute the LDSC weights (not implemented yet)
 	if (cPar.a_mode==72) {
 		cout<<"Calculating Weights ... "<<endl;
@@ -1364,7 +1363,7 @@ void GEMMA::BatchRun (PARAM &cPar)
 
 		cVarcov.CopyToParam(cPar);
 	}
-
+	*/
 
 	//Compute the S matrix (and its variance), that is used for variance component estimation using summary statistics
 	if (cPar.a_mode==25 || cPar.a_mode==26) {
@@ -1472,7 +1471,7 @@ void GEMMA::BatchRun (PARAM &cPar)
 	  gsl_vector_free (s);
 	}
 
-
+	/*
     //Calculate SNP covariance
 	if (cPar.a_mode==71) {
 	  VARCOV cVarcov;
@@ -1486,7 +1485,7 @@ void GEMMA::BatchRun (PARAM &cPar)
 
 	  cVarcov.CopyToParam(cPar);
 	}
-
+	*/
 
 	//LM
 	if (cPar.a_mode==51 || cPar.a_mode==52 || cPar.a_mode==53 || cPar.a_mode==54) {  //Fit LM
@@ -2400,7 +2399,7 @@ void GEMMA::BatchRun (PARAM &cPar)
 	}
 
 
-
+	/*
 	//LDR
 	if (cPar.a_mode==14) {
 		gsl_vector *y=gsl_vector_alloc (cPar.ni_test);
@@ -2428,7 +2427,7 @@ void GEMMA::BatchRun (PARAM &cPar)
 		gsl_matrix_free (W);
 		gsl_matrix_free (G);
 	}
-
+	*/
 	cPar.time_total=(clock()-time_begin)/(double(CLOCKS_PER_SEC)*60.0);
 
 	return;
-- 
cgit v1.2.3