Update Fig1_Prognostic_bias.R

CSkylarL · Aug 16, 2023 · 103cf02 · 103cf02
1 parent f2a2642
commit 103cf02
Showing 1 changed file with 10 additions and 1 deletion.
diff --git a/Fig1_Prognostic_bias.R b/Fig1_Prognostic_bias.R
@@ -80,7 +80,16 @@ cox_img_res <- data.frame(HR=hr, PV=pval, QV=QV)
 row.names(cox_img_res) <- row.names(mygrp)
 
 # b) example(Related to Figure 1D and S2)----------------------------------------
-good.gene <- rownames(resdata)[which(resdata$group == "Favorable")]
+tmp <- cox_img_res
+tmp$group <- as.factor(
+  ifelse(
+    tmp$QV < 0.05,
+    ifelse(tmp$HR < 1, "Favorable", "Unfavorable"),
+    "NS"
+  )
+)
+
+good.gene <- rownames(tmp)[which(tmp$group == "Favorable")]
 good.exp <- TimiPrePropress(marker = good.gene,rna = SKCM06rna,
                             cohort = rownames(info),
                             log = T,GMNorm = T)