X-Git-Url: https://git.donarmstrong.com/?a=blobdiff_plain;f=linearalgebra.cpp;h=e19861912d9c7fc78f9b2cbbb87ecd76e9388ba9;hb=f0a594f6676ef5a52d1f122b6de70de2fda08c81;hp=6970c6429481f94cb5d9b7a26bd990407410ff8d;hpb=65f2be9d8f9b4cf72b4aafd800fdd11ad26435fd;p=mothur.git

diff --git a/linearalgebra.cpp b/linearalgebra.cpp
index 6970c64..e198619 100644
--- a/linearalgebra.cpp
+++ b/linearalgebra.cpp
@@ -696,4 +696,30 @@ double LinearAlgebra::calcKendall(vector< vector<double> >& euclidDists, vector<
 
 /*********************************************************************************************************************************/
 
+vector<vector<double> > LinearAlgebra::getObservedEuclideanDistance(vector<vector<double> >& relAbundData){
 
+	int numSamples = relAbundData.size();
+	int numOTUs = relAbundData[0].size();
+	
+	vector<vector<double> > dMatrix(numSamples);
+	for(int i=0;i<numSamples;i++){
+		dMatrix[i].resize(numSamples);
+	}
+	
+	for(int i=0;i<numSamples;i++){
+		for(int j=0;j<numSamples;j++){
+			
+			double d = 0;
+			for(int k=0;k<numOTUs;k++){
+				d += pow((relAbundData[i][k] - relAbundData[j][k]), 2.0000);
+			}
+			dMatrix[i][j] = pow(d, 0.50000);
+			dMatrix[j][i] = dMatrix[i][j];
+			
+		}
+	}
+	return dMatrix;
+	
+}
+
+/*********************************************************************************************************************************/