about summary refs log tree commit diff
path: root/src/param.cpp
diff options
context:
space:
mode:
authorxiangzhou2015-07-11 15:36:03 -0400
committerxiangzhou2015-07-11 15:36:03 -0400
commit3ab77854a52ac016b9e2b824858f5f0c695d4170 (patch)
tree040c78f03f5bf22df7c4b2b8263702a842249adf /src/param.cpp
parent2e82184182d8bbaa2cc3fa8088058bf3de2048bd (diff)
downloadpangemma-3ab77854a52ac016b9e2b824858f5f0c695d4170.tar.gz
add GxE test
Diffstat (limited to 'src/param.cpp')
-rw-r--r--src/param.cpp15
1 files changed, 7 insertions, 8 deletions
diff --git a/src/param.cpp b/src/param.cpp
index c4b234a..33b7b48 100644
--- a/src/param.cpp
+++ b/src/param.cpp
@@ -169,7 +169,6 @@ void PARAM::ReadFiles (void)
 		//post-process covariates and phenotypes, obtain ni_test, save all useful covariates
 		ProcessCvtPhen();
 
-
 		//obtain covariate matrix
 		gsl_matrix *W=gsl_matrix_alloc (ni_test, n_cvt);
 		CopyCvt (W);
@@ -242,12 +241,12 @@ void PARAM::ReadFiles (void)
 			indicator_idv.push_back(k);
 		}
 
-        //post-process covariates and phenotypes, obtain ni_test, save all useful covariates
-        ProcessCvtPhen();
+		//post-process covariates and phenotypes, obtain ni_test, save all useful covariates
+		ProcessCvtPhen();
 
-        //obtain covariate matrix
-        gsl_matrix *W=gsl_matrix_alloc (ni_test, n_cvt);
-        CopyCvt (W);
+		//obtain covariate matrix
+		gsl_matrix *W=gsl_matrix_alloc (ni_test, n_cvt);
+		CopyCvt (W);
 
 		if (ReadFile_gene (file_gene, vec_read, snpInfo, ng_total)==false) {error=true;}
 	}
@@ -381,9 +380,9 @@ void PARAM::CheckParam (void)
 	}
 
 	if (!file_oxford.empty()) {
-		str=file_bfile+".bgen";
+		str=file_oxford+".bgen";
 		if (stat(str.c_str(),&fileInfo)==-1) {cout<<"error! fail to open .bgen file: "<<str<<endl; error=true;}
-		str=file_bfile+".sample";
+		str=file_oxford+".sample";
 		if (stat(str.c_str(),&fileInfo)==-1) {cout<<"error! fail to open .sample file: "<<str<<endl; error=true;}
 	}