about summary refs log tree commit diff
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,