fix ase in case no DNA data is provided

CompEpigen · Nov 20, 2024 · 14b904a · 14b904a
1 parent bcbdf46
commit 14b904a
Show file tree

Hide file tree

Showing 4 changed files with 5 additions and 4 deletions.
diff --git a/figeno/ase.py b/figeno/ase.py
@@ -32,7 +32,8 @@ def read_SNPs_RNA(ase_file,gene=None,chr=None,start=None,end=None,exons=None,min
     if "refCount_DNA" in df.columns: variables_selected+=["refCount_DNA","altCount_DNA"]
     df = df[variables_selected]
     df["VAF"] = [min(df.loc[i,"refCount"],df.loc[i,"altCount"]) / (df.loc[i,"refCount"]+df.loc[i,"altCount"]) for i in df.index]
-    df["VAF_DNA"] = [min(df.loc[i,"refCount_DNA"],df.loc[i,"altCount_DNA"]) / (df.loc[i,"refCount_DNA"]+df.loc[i,"altCount_DNA"]) for i in df.index]
+    if "refCount_DNA" in df.columns:
+        df["VAF_DNA"] = [min(df.loc[i,"refCount_DNA"],df.loc[i,"altCount_DNA"]) / (df.loc[i,"refCount_DNA"]+df.loc[i,"altCount_DNA"]) for i in df.index]
     in_exons=[]
     if exons is not None:
         for i in df.index:

diff --git a/figeno/cli/cli.py b/figeno/cli/cli.py
@@ -2,7 +2,7 @@
 
 from figeno.cli import gui, init,make
 
-__version__ = "1.6.3"
+__version__ = "1.6.4"
 
 def main():
     parser = ArgumentParser("figeno",formatter_class=ArgumentDefaultsHelpFormatter)

diff --git a/figeno/gui/package.json b/figeno/gui/package.json
@@ -1,6 +1,6 @@
 {
   "name": "figeno",
-  "version": "1.6.3",
+  "version": "1.6.4",
   "private": true,
   "homepage": "./",
   "dependencies": {

diff --git a/pyproject.toml b/pyproject.toml
@@ -9,7 +9,7 @@ packages = ["figeno", "figeno.data", "figeno.cli", "figeno.gui"]
 
 [project]
 name = 'figeno'
-version = "1.6.3"
+version = "1.6.4"
 description = 'Package for generating genomics figures.'
 readme = 'README.md'
 authors = [