SURE-tools 2.1.76__py3-none-any.whl → 2.1.78__py3-none-any.whl

SURE/PerturbFlow.py

@@ -881,25 +881,11 @@ class PerturbFlow(nn.Module):
         if self.loss_func == 'bernoulli':
             #counts = self.sigmoid(concentrate)
             counts = dist.Bernoulli(logits=concentrate).to_event(1).mean
-        elif self.loss_func == 'negbinomial':
-            rate = concentrate.exp()
-            theta = dist.DirichletMultinomial(total_count=1, concentration=rate).mean
-            
-            total_count = self.total_count
-            #total_count = pyro.param("inverse_dispersion")
-            #store = pyro.get_param_store()
-            #total_count = store['inverse_dispersion']
-            counts = dist.NegativeBinomial(total_count=total_count, probs=theta).to_event(1).mean
-        elif self.loss_func == 'poisson':
+        else:
             rate = concentrate.exp()
             theta = dist.DirichletMultinomial(total_count=1, concentration=rate).mean
             counts = theta * library_size
             #counts = dist.Poisson(rate=rate).to_event(1).mean
-        elif self.loss_func == 'multinomial':
-            rate = concentrate.exp()
-            theta = dist.DirichletMultinomial(total_count=1, concentration=rate).mean
-            #counts = dist.Multinomial(total_count=int(1e8), probs=theta).mean
-            counts = theta * library_size
         return counts
     
     def _count_sample(self,concentrate):
@@ -911,22 +897,17 @@ class PerturbFlow(nn.Module):
             counts = dist.Poisson(rate=counts).to_event(1).sample()
         return counts
     
-    def get_counts(self, zs, library_sizes = None,
+    def get_counts(self, zs, library_sizes,
                         batch_size: int = 1024, 
                         use_sampler: bool = False):
 
         zs = convert_to_tensor(zs, device=self.get_device())
         
-        if self.loss_func in ['multinomial','poisson']:
-            assert library_sizes is not None, 'Library sizes are required for multinomial!'
-            
-            if type(library_sizes) == list:
-                library_sizes = np.array(library_sizes).view(-1,1)
-            elif len(library_sizes.shape)==1:
-                library_sizes = library_sizes.view(-1,1)
-            ls = convert_to_tensor(library_sizes, device=self.get_device())
-        else:
-            ls = zs
+        if type(library_sizes) == list:
+            library_sizes = np.array(library_sizes).view(-1,1)
+        elif len(library_sizes.shape)==1:
+            library_sizes = library_sizes.view(-1,1)
+        ls = convert_to_tensor(library_sizes, device=self.get_device())
             
         dataset = CustomDataset2(zs,ls)
         dataloader = DataLoader(dataset, batch_size=batch_size, shuffle=False)
@@ -1084,7 +1065,8 @@ class PerturbFlow(nn.Module):
                 pbar.set_postfix({'loss': str_loss})
                 pbar.update(1)
                 
-        self.total_count = pyro.param('inverse_dispersion')
+        if self.loss_func == 'negbinomial':
+            self.total_count = pyro.param('inverse_dispersion')
 
     @classmethod
     def save_model(cls, model, file_path, compression=False):

{sure_tools-2.1.76.dist-info → sure_tools-2.1.78.dist-info}/METADATA

@@ -1,6 +1,6 @@
 Metadata-Version: 2.4
 Name: SURE-tools
-Version: 2.1.76
+Version: 2.1.78
 Summary: Succinct Representation of Single Cells
 Home-page: https://github.com/ZengFLab/SURE
 Author: Feng Zeng

{sure_tools-2.1.76.dist-info → sure_tools-2.1.78.dist-info}/RECORD

@@ -1,4 +1,4 @@
-SURE/PerturbFlow.py,sha256=sltfC-Cb5bFI_EJTMLqjksZiSRCHQvP89Q1NFId2fBg,54669
+SURE/PerturbFlow.py,sha256=BbpN3CSA-uSvIwVzHP1Vq7fuUxG3VC-x1KClCGa1euk,53695
 SURE/SURE.py,sha256=ko15a9BhvUqHviogZ0YCdTQjM-2zqkO9OvHZSpnGbg0,47458
 SURE/__init__.py,sha256=NOJI_K-eCqPgStXXvgl3wIEMp6d8saMTDYLJ7Ga9MqE,293
 SURE/assembly/__init__.py,sha256=jxZLURXKPzXe21LhrZ09LgZr33iqdjlQy4oSEj5gR2Q,172
@@ -17,9 +17,9 @@ SURE/utils/__init__.py,sha256=YF5jB-PAHJQ40OlcZ7BCZbsN2q1JKuPT6EppilRXQqM,680
 SURE/utils/custom_mlp.py,sha256=C0EXLGYsWkUQpEL49AyBFPSzKmasb2hdvtnJfxbF-YU,9282
 SURE/utils/queue.py,sha256=E_5PA5EWcBoGAZj8BkKQnkCK0p4C-4-xcTPqdIXaPXU,1892
 SURE/utils/utils.py,sha256=IUHjDDtYaAYllCWsZyIzqQwaLul6fJRvHRH4vIYcR-c,8462
-sure_tools-2.1.76.dist-info/licenses/LICENSE,sha256=TFHKwmrAViXQbSX5W-NDItkWFjm45HWOeUniDrqmnu0,1065
-sure_tools-2.1.76.dist-info/METADATA,sha256=h-xwe5IHtIZrIg8jlHXOxm2z9irVDcg1eFYnIoxr374,2678
-sure_tools-2.1.76.dist-info/WHEEL,sha256=_zCd3N1l69ArxyTb8rzEoP9TpbYXkqRFSNOD5OuxnTs,91
-sure_tools-2.1.76.dist-info/entry_points.txt,sha256=-nJI8rVe_qqrR0HmfAODzj-JNfEqCcSsyVh6okSqyHk,83
-sure_tools-2.1.76.dist-info/top_level.txt,sha256=BtFTebdiJeqra4r6mm-uEtwVRFLZ_IjYsQ7OnalrOvY,5
-sure_tools-2.1.76.dist-info/RECORD,,
+sure_tools-2.1.78.dist-info/licenses/LICENSE,sha256=TFHKwmrAViXQbSX5W-NDItkWFjm45HWOeUniDrqmnu0,1065
+sure_tools-2.1.78.dist-info/METADATA,sha256=39CFgVJlh26Bd1i8qkQcK5jQOiDGDz_JzoNKKVS3S5Q,2678
+sure_tools-2.1.78.dist-info/WHEEL,sha256=_zCd3N1l69ArxyTb8rzEoP9TpbYXkqRFSNOD5OuxnTs,91
+sure_tools-2.1.78.dist-info/entry_points.txt,sha256=-nJI8rVe_qqrR0HmfAODzj-JNfEqCcSsyVh6okSqyHk,83
+sure_tools-2.1.78.dist-info/top_level.txt,sha256=BtFTebdiJeqra4r6mm-uEtwVRFLZ_IjYsQ7OnalrOvY,5
+sure_tools-2.1.78.dist-info/RECORD,,