aboutsummaryrefslogtreecommitdiff
path: root/src/lmm.h
diff options
context:
space:
mode:
authorPjotr Prins2017-10-06 08:41:17 +0000
committerPjotr Prins2017-10-06 08:41:17 +0000
commit2add397847701d8f939eab55bacddb54fdb8c641 (patch)
treef793f52557dd9ef071aae77a34ea5b000466caa5 /src/lmm.h
parentdf161be507ac0ad1d67a6528ebc664acec89fc9c (diff)
downloadpangemma-2add397847701d8f939eab55bacddb54fdb8c641.tar.gz
LMM: lifted out file access
Diffstat (limited to 'src/lmm.h')
-rw-r--r--src/lmm.h12
1 files changed, 9 insertions, 3 deletions
diff --git a/src/lmm.h b/src/lmm.h
index 72ce523..bd303b0 100644
--- a/src/lmm.h
+++ b/src/lmm.h
@@ -23,6 +23,7 @@
#include "gsl/gsl_vector.h"
#include "io.h"
#include "param.h"
+#include <functional>
using namespace std;
@@ -89,13 +90,18 @@ public:
void AnalyzeGene(const gsl_matrix *U, const gsl_vector *eval,
const gsl_matrix *UtW, const gsl_vector *Utx,
const gsl_matrix *W, const gsl_vector *x);
- void AnalyzePlink(const gsl_matrix *U, const gsl_vector *eval,
- const gsl_matrix *UtW, const gsl_vector *Uty,
- const gsl_matrix *W, const gsl_vector *y);
+ void Analyze(std::function< string(void) >& fetch_line,
+ const gsl_matrix *U, const gsl_vector *eval,
+ const gsl_matrix *UtW, const gsl_vector *Uty,
+ const gsl_matrix *W, const gsl_vector *y,
+ const set<string> gwasnps);
void AnalyzeBimbam(const gsl_matrix *U, const gsl_vector *eval,
const gsl_matrix *UtW, const gsl_vector *Uty,
const gsl_matrix *W, const gsl_vector *y,
const set<string> gwasnps);
+ void AnalyzePlink(const gsl_matrix *U, const gsl_vector *eval,
+ const gsl_matrix *UtW, const gsl_vector *Uty,
+ const gsl_matrix *W, const gsl_vector *y);
void AnalyzePlinkGXE(const gsl_matrix *U, const gsl_vector *eval,
const gsl_matrix *UtW, const gsl_vector *Uty,
const gsl_matrix *W, const gsl_vector *y,