about summary refs log tree commit diff
path: root/src/lm.h
diff options
context:
space:
mode:
authorPeter Carbonetto2017-08-07 13:23:44 -0500
committerPeter Carbonetto2017-08-07 13:23:44 -0500
commitecc5c0b5b11cbc2501b60d6d4e9a9b3e8f270d02 (patch)
treeb27cd975b19b85e75255cacf1c3d771266c56a10 /src/lm.h
parent8be2ce433c015d8c6776d270537d645c7e761bdf (diff)
parent7360d14216400b8f12fbfda03ac2f4827b102711 (diff)
downloadpangemma-ecc5c0b5b11cbc2501b60d6d4e9a9b3e8f270d02.tar.gz
Merge branch 'master' of github.com:genetics-statistics/GEMMA
Diffstat (limited to 'src/lm.h')
-rw-r--r--src/lm.h92
1 files changed, 46 insertions, 46 deletions
diff --git a/src/lm.h b/src/lm.h
index cf428f0..cb22d3b 100644
--- a/src/lm.h
+++ b/src/lm.h
@@ -19,61 +19,61 @@
 #ifndef __LM_H__
 #define __LM_H__
 
-#include "gsl/gsl_vector.h"
 #include "gsl/gsl_matrix.h"
-#include "param.h"
+#include "gsl/gsl_vector.h"
 #include "io.h"
+#include "param.h"
 
 using namespace std;
 
 class LM {
 
 public:
-	// IO-related parameters.
-	int a_mode;	// Analysis mode: 50+1/2/3/4 for Frequentist tests.
-	size_t d_pace;	// Display pace.
-
-	string file_bfile;
-	string file_geno;
-	string file_oxford;
-	string file_out;
-	string path_out;
-
-	string file_gene;
-
-	// Summary statistics.
-	size_t ni_total, ni_test;  // Number of individuals.
-	size_t ns_total, ns_test;  // Number of SNPs.
-	size_t ng_total, ng_test;  // Number of genes.
-	size_t n_cvt;
-	double time_opt;	   // Time spent.
-
-        // Indicator for individuals (phenotypes): 0 missing, 1
-        // available for analysis.
-	vector<int> indicator_idv;
-
-        // Sequence indicator for SNPs: 0 ignored because of (a) maf,
-        // (b) miss, (c) non-poly; 1 available for analysis.
-	vector<int> indicator_snp;
-
-	vector<SNPINFO> snpInfo;  // Record SNP information.
-
-	// Not included in PARAM.
-	vector<SUMSTAT> sumStat;  // Output SNPSummary Data.
-
-	// Main functions.
-	void CopyFromParam (PARAM &cPar);
-	void CopyToParam (PARAM &cPar);
-	void AnalyzeGene (const gsl_matrix *W, const gsl_vector *x);
-	void AnalyzePlink (const gsl_matrix *W, const gsl_vector *y);
-	void AnalyzeBimbam (const gsl_matrix *W, const gsl_vector *y);
-	// WJA added.
-	void Analyzebgen (const gsl_matrix *W, const gsl_vector *y);
-
-	void WriteFiles ();
+  // IO-related parameters.
+  int a_mode;    // Analysis mode: 50+1/2/3/4 for Frequentist tests.
+  size_t d_pace; // Display pace.
+
+  string file_bfile;
+  string file_geno;
+  string file_oxford;
+  string file_out;
+  string path_out;
+
+  string file_gene;
+
+  // Summary statistics.
+  size_t ni_total, ni_test; // Number of individuals.
+  size_t ns_total, ns_test; // Number of SNPs.
+  size_t ng_total, ng_test; // Number of genes.
+  size_t n_cvt;
+  double time_opt; // Time spent.
+
+  // Indicator for individuals (phenotypes): 0 missing, 1
+  // available for analysis.
+  vector<int> indicator_idv;
+
+  // Sequence indicator for SNPs: 0 ignored because of (a) maf,
+  // (b) miss, (c) non-poly; 1 available for analysis.
+  vector<int> indicator_snp;
+
+  vector<SNPINFO> snpInfo; // Record SNP information.
+
+  // Not included in PARAM.
+  vector<SUMSTAT> sumStat; // Output SNPSummary Data.
+
+  // Main functions.
+  void CopyFromParam(PARAM &cPar);
+  void CopyToParam(PARAM &cPar);
+  void AnalyzeGene(const gsl_matrix *W, const gsl_vector *x);
+  void AnalyzePlink(const gsl_matrix *W, const gsl_vector *y);
+  void AnalyzeBimbam(const gsl_matrix *W, const gsl_vector *y);
+  // WJA added.
+  void Analyzebgen(const gsl_matrix *W, const gsl_vector *y);
+
+  void WriteFiles();
 };
 
-void MatrixCalcLmLR (const gsl_matrix *X, const gsl_vector *y,
-		     vector<pair<size_t, double> > &pos_loglr);
+void MatrixCalcLmLR(const gsl_matrix *X, const gsl_vector *y,
+                    vector<pair<size_t, double>> &pos_loglr);
 
 #endif