diff options
author | Peter Carbonetto | 2017-08-02 13:07:23 -0500 |
---|---|---|
committer | GitHub | 2017-08-02 13:07:23 -0500 |
commit | adaf9557f776ca274b51e921af0542ef1b84eb61 (patch) | |
tree | 8f7581bec20eb815848d02f926cb949ed038008b /src/lm.h | |
parent | 84360c191f418bf8682b35e0c8235fcc3bd19a06 (diff) | |
parent | 1e49e7a27d0f4b811a87c64db1e875779766f6b0 (diff) | |
download | pangemma-adaf9557f776ca274b51e921af0542ef1b84eb61.tar.gz |
Merge pull request #62 from genenetwork/spacing
Spacing and coding style.
Diffstat (limited to 'src/lm.h')
-rw-r--r-- | src/lm.h | 92 |
1 files changed, 46 insertions, 46 deletions
@@ -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 |