added correction for m*MMD change in quadratic time MMD

shogun-toolbox · Aug 11, 2012 · d3c42f6 · d3c42f6
1 parent 2f96a8a
commit d3c42f6
Showing 1 changed file with 3 additions and 3 deletions.
diff --git a/examples/undocumented/libshogun/statistics_quadratic_time_mmd.cpp b/examples/undocumented/libshogun/statistics_quadratic_time_mmd.cpp
@@ -35,14 +35,14 @@ void test_quadratic_mmd_fixed()
 
 	CQuadraticTimeMMD* mmd=new CQuadraticTimeMMD(kernel, features, n);
 
-	/* unbiased statistic */
+	/* unbiased statistic, shogun return m*MMD */
 	mmd->set_statistic_type(UNBIASED);
-	float64_t difference=CMath::abs(mmd->compute_statistic()-0.051325806508381);
+	float64_t difference=CMath::abs(mmd->compute_statistic()/n-0.051325806508381);
 	ASSERT(difference<=10E-16);
 
 	/* biased statistic */
 	mmd->set_statistic_type(BIASED);
-	difference=CMath::abs(mmd->compute_statistic()-1.017107688196714);
+	difference=CMath::abs(mmd->compute_statistic()/n-1.017107688196714);
 	ASSERT(difference<=10E-16);
 
 	SG_UNREF(mmd);