X-Git-Url: https://git.donarmstrong.com/?p=mothur.git;a=blobdiff_plain;f=linearalgebra.h;h=e0933ee52564915d35bcf6e9f91f73c2cea80d2d;hp=58a60dae8f04d9eae12fff9f8d729cc2119d4c92;hb=a8e2df1b96a57f5f29576b08361b86a96a8eff4f;hpb=f0a594f6676ef5a52d1f122b6de70de2fda08c81

diff --git a/linearalgebra.h b/linearalgebra.h
index 58a60da..e0933ee 100644
--- a/linearalgebra.h
+++ b/linearalgebra.h
@@ -29,11 +29,47 @@ public:
 	double calcPearson(vector<vector<double> >&, vector<vector<double> >&);
 	double calcSpearman(vector<vector<double> >&, vector<vector<double> >&);
 	double calcKendall(vector<vector<double> >&, vector<vector<double> >&);
+	
+	double calcPearson(vector<double>&, vector<double>&, double&);
+	double calcSpearman(vector<double>&, vector<double>&, double&);
+	double calcKendall(vector<double>&, vector<double>&, double&);
+    
+	double calcSpearmanSig(double, double, double, double); //length, f^3 - f where f is the number of ties in x, f^3 - f where f is the number of ties in y, sum of squared diffs in ranks. - designed to find the sif of one score.
+    double calcPearsonSig(double, double); //length, coeff.
+    double calcKendallSig(double, double); //length, coeff.
+    
+    vector<double> solveEquations(vector<vector<double> >, vector<double>);
+    vector<float> solveEquations(vector<vector<float> >, vector<float>);
+    vector<vector<double> > getInverse(vector<vector<double> >);
 
+    
 private:
 	MothurOut* m;
 	
 	double pythag(double, double);
+    double betacf(const double, const double, const double);
+    double betai(const double, const double, const double);
+    double gammln(const double);
+    double gammp(const double, const double);
+    double gammq(const double, const double);
+    double gser(double&, const double, const double, double&);
+    double gcf(double&, const double, const double, double&);
+    double erfcc(double);
+    
+    double ran0(int&); //for testing 
+    double ran1(int&); //for testing
+    double ran2(int&); //for testing
+    double ran3(int&); //for testing
+    double ran4(int&); //for testing
+    void psdes(unsigned long &, unsigned long &); //for testing
+    
+    void ludcmp(vector<vector<double> >&, vector<int>&, double&);
+    void lubksb(vector<vector<double> >&, vector<int>&, vector<double>&);
+    
+    void ludcmp(vector<vector<float> >&, vector<int>&, float&);
+    void lubksb(vector<vector<float> >&, vector<int>&, vector<float>&);
+
+    
 };
 
 #endif