merge q2info

CNIC-Proteomics · Feb 19, 2024 · e3375ec · e3375ec
1 parent ef3211b
commit e3375ec
Show file tree

Hide file tree

Showing 2 changed files with 8 additions and 0 deletions.
diff --git a/4_qTableReport/qReportMaker.py b/4_qTableReport/qReportMaker.py
@@ -374,6 +374,11 @@ def qReportDesign(config, quan, qTableD, contrast):
         for i,j in qTableD.columns
         ])
 
+    # Add info contained in q2info
+    if os.path.isfile(config['q2info']):
+        q2info = pd.read_csv(config['q2info'], sep='\t')
+        q2info.columns = pd.MultiIndex.from_tuples([qTableD.columns[0] if n==0 else (i,'','') for n,i in enumerate(q2info.columns)])
+        qTableD = pd.merge(q2info, qTableD, how='right', on=[qTableD.columns[0]])
 
     plotted_q = [os.path.splitext(i)[0] for i in os.listdir(config['plotFolder'])]
     qTableD[qTableD.columns[0]] = \

diff --git a/4_qTableReport/qReportMaker.yaml b/4_qTableReport/qReportMaker.yaml
@@ -8,6 +8,9 @@ infile: D:\CNIC\Scripts\ReportStats\test\FDR_LIMMA_NM_qfq_table_pgmfreq.tsv
 # Output file
 outfolder: D:\CNIC\Scripts\ReportStats\test\
 
+# Path to file relating protein Uniprot ID (first col) to other information (e.g. Category)
+q2info: S:\U_Proteomica\UNIDAD\NextCloud\1_PTM_Analysis\Heteroplasmia_Heart\ReportStats-v0.6\myMitocarta.tsv
+
 # FDR threshold used
 qvThr:
     - 0.01