X-Git-Url: https://git.donarmstrong.com/?a=blobdiff_plain;f=linearalgebra.h;h=7453f4ebd8b86ed20f530c2e1ba7989808361fb2;hb=6b32d112bb60e9f7eb6d4407a4eed4c49b67bced;hp=5819ea85d8deefba15ca1c45ccf03efd21c9321a;hpb=372fb21ea66ced432b109225851a1b80ef0491a3;p=mothur.git

diff --git a/linearalgebra.h b/linearalgebra.h
index 5819ea8..7453f4e 100644
--- a/linearalgebra.h
+++ b/linearalgebra.h
@@ -20,9 +20,12 @@ public:
 	~LinearAlgebra() {}
 	
 	vector<vector<double> > matrix_mult(vector<vector<double> >, vector<vector<double> >);
+    vector<vector<double> >transpose(vector<vector<double> >);
 	void recenter(double, vector<vector<double> >, vector<vector<double> >&);
-	int tred2(vector<vector<double> >&, vector<double>&, vector<double>&);
+	//eigenvectors
+    int tred2(vector<vector<double> >&, vector<double>&, vector<double>&);
 	int qtli(vector<double>&, vector<double>&, vector<vector<double> >&);
+    
 	vector< vector<double> > calculateEuclidianDistance(vector<vector<double> >&, int); //pass in axes and number of dimensions
 	vector< vector<double> > calculateEuclidianDistance(vector<vector<double> >&); //pass in axes
 	vector<vector<double> > getObservedEuclideanDistance(vector<vector<double> >&);
@@ -43,8 +46,10 @@ public:
     vector<double> solveEquations(vector<vector<double> >, vector<double>);
     vector<float> solveEquations(vector<vector<float> >, vector<float>);
     vector<vector<double> > getInverse(vector<vector<double> >);
-    
-    
+    double choose(double, double);
+    double normalvariate(double mu, double sigma);
+    vector< vector<double> > lda(vector< vector<double> >& a, vector<string> groups, vector< vector<double> >& means, bool&); //Linear discriminant analysis - a is [features][valuesFromGroups] groups indicates which group each sampling comes from. For example if groups = early, late, mid, early, early. a[0][0] = value for feature0 from groupEarly.
+    int svd(vector< vector<double> >& a, vector<double>& w, vector< vector<double> >& v); //Singular value decomposition
 private:
 	MothurOut* m;
 	
@@ -57,6 +62,7 @@ private:
     double gcf(double&, const double, const double, double&);
     double erfcc(double);
     double gammp(const double, const double);
+    double pnorm(double x);
     
     double ran0(int&); //for testing 
     double ran1(int&); //for testing
@@ -70,7 +76,6 @@ private:
     
     void ludcmp(vector<vector<float> >&, vector<int>&, float&);
     void lubksb(vector<vector<float> >&, vector<int>&, vector<float>&);
-
     
 };