Remove ttest from pySpikeDB. Better to just use python libraries.

baubie · Apr 1, 2012 · b92d07e · b92d07e
1 parent a2a7796
commit b92d07e
Show file tree

Hide file tree

Showing 3 changed files with 0 additions and 36 deletions.
diff --git a/src/pySpikeDB.cpp b/src/pySpikeDB.cpp
@@ -724,40 +724,6 @@ double pySpikeDB::stddev(boost::python::list &v)
 	return sqrt(r/(rN-1));
 }
 
-bp::object pySpikeDB::ttest(bp::list &a, bp::list &b, bool eqvar)
-{
-	bp::dict r;
-	double t = 0;
-	double p = 0;
-	double df = 0;
-	double ma = mean(a);
-	double mb = mean(b);
-	double sa = stddev(a);
-	double sb = stddev(b);
-	double na = bp::len(a);
-	double nb = bp::len(b);
-	df = na + nb - 2;
-
-	if (eqvar) 
-	{
-		if (bp::len(a) == bp::len(b))
-		{
-			double S = sqrt(0.5*(sa*sa+sb*sb));
-			t = (ma-mb)/(S*sqrt(2.0/na));
-		} else {
-			double S = sqrt(((na-1)*sa*sa+(nb-1)*sb*sb)/(na+nb-2));
-			t = (ma-mb)/(S*sqrt(1.0/na + 1.0/nb));
-		}
-	} else {
-		// Unequal variance not yet implemented
-	}
-	r["t"] = t>0?t:-t; // Absolute value
-	r["p1"] = p;
-	r["p2"] = p;
-	r["df"] = df;
-	return r;
-}
-
 pySpikeDB::Option* pySpikeDB::findOptionByName(const Glib::ustring &name) {
 
 	std::vector<boost::shared_ptr<Option> >::iterator it;

diff --git a/src/pySpikeDB.h b/src/pySpikeDB.h
@@ -95,7 +95,6 @@ class pySpikeDB {
 
 		double stddev(boost::python::list &v);
 		double mean(boost::python::list &v);
-		boost::python::object ttest(boost::python::list &a, boost::python::list &b, bool eqvar);
 
 		void print(const std::string &s);
 

diff --git a/src/uiAnalysis.cpp b/src/uiAnalysis.cpp
@@ -375,7 +375,6 @@ if (!uiAnalysis::setupPython)
 		.def("getFilesSingleCell", &pySpikeDB::getFilesSingleCell)
 		.def("mean", &pySpikeDB::mean)
 		.def("stddev", &pySpikeDB::stddev)
-		.def("ttest", &pySpikeDB::ttest)
 		.def("filterSpikesAbs", &pySpikeDB::filterSpikesAbs)
 		.def("filterSpikesRel", &pySpikeDB::filterSpikesRel)
 		.def("plotXLabel", &pySpikeDB::plotXLabel)