diff options
| author | xiangzhou | 2015-07-11 15:36:03 -0400 |
|---|---|---|
| committer | xiangzhou | 2015-07-11 15:36:03 -0400 |
| commit | 3ab77854a52ac016b9e2b824858f5f0c695d4170 (patch) | |
| tree | 040c78f03f5bf22df7c4b2b8263702a842249adf /src/param.cpp | |
| parent | 2e82184182d8bbaa2cc3fa8088058bf3de2048bd (diff) | |
| download | pangemma-3ab77854a52ac016b9e2b824858f5f0c695d4170.tar.gz | |
add GxE test
Diffstat (limited to 'src/param.cpp')
| -rw-r--r-- | src/param.cpp | 15 |
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;} } |