Fixed issue #31: GSEA error

setup.py

@@ -14,7 +14,7 @@
 
 
 setup(name='wishbone',
-      version='0.4',
+      version='0.4.1',
       description='Wishbone algorithm for identifying bifurcating trajectories from single-cell data',
       author='Manu Setty',
       author_email='[email protected]',
@@ -24,7 +24,7 @@
           'numpy>=1.10.0',
           'pandas>=0.18.0',
           'scipy>=0.14.0',
-          'tsne==0.1.3',
+          'bhtsne',
           'matplotlib',
           'seaborn',
           'sklearn',

src/wishbone/__init__.py

@@ -3,4 +3,4 @@
 from . import wishbone_gui
 from . import autocomplete_entry
 
-__version__ = "0.4"
+__version__ = "0.4.1"

src/wishbone/wb.py

@@ -24,7 +24,7 @@
     warnings.simplefilter('ignore')  # catch experimental ipython widget warning
     import seaborn as sns
 
-from tsne import bh_sne
+import bhtsne
 from scipy.sparse import csr_matrix, find
 from scipy.sparse.linalg import eigs
 from numpy.linalg import norm
@@ -417,8 +417,9 @@ def run_tsne(self, n_components=15, perplexity=30):
         if data.shape[0] < 100 and perplexity > perplexity_limit:
             print('Reducing perplexity to %d since there are <100 cells in the dataset. ' % perplexity_limit)
             perplexity = perplexity_limit
-        self.tsne = pd.DataFrame(bh_sne(data, perplexity=perplexity),
-                                 index=self.data.index, columns=['x', 'y'])
+        self.tsne = pd.DataFrame(bhtsne.tsne(data, perplexity=perplexity),
+                         index=self.data.index, columns=['x', 'y'])
+
 
     def plot_tsne(self, fig=None, ax=None, title='tSNE projection'):
         """Plot tSNE projections of the data
@@ -782,8 +783,8 @@ def plot_gene_component_correlations(
 
     @staticmethod
     def _gmt_options():
-        mouse_options = os.listdir(os.path.expanduser('~/.seqc/tools/mouse'))
-        human_options = os.listdir(os.path.expanduser('~/.seqc/tools/human'))
+        mouse_options = os.listdir(os.path.expanduser('~/.wishbone/tools/mouse'))
+        human_options = os.listdir(os.path.expanduser('~/.wishbone/tools/human'))
         print('Available GSEA .gmt files:\n\nmouse:\n{m}\n\nhuman:\n{h}\n'.format(
                 m='\n'.join(mouse_options),
                 h='\n'.join(human_options)))
@@ -863,7 +864,7 @@ def run_gsea(self, output_stem, gmt_file=None,
         else:
             if not len(gmt_file) == 2:
                 raise ValueError('gmt_file should be a tuple of (organism, filename).')
-            gmt_file = os.path.expanduser('~/.seqc/tools/{}/{}').format(*gmt_file)
+            gmt_file = os.path.expanduser('~/.wishbone/tools/{}/{}').format(*gmt_file)
 
         if components is None:
             components = self.diffusion_map_correlations.columns