Skip to content

Commit

Permalink
Update intermediate data
Browse files Browse the repository at this point in the history
Some redundant example data was moved from TimiGP R package to this repo
  • Loading branch information
CSkylarL committed Sep 7, 2023
1 parent 103cf02 commit fff687d
Show file tree
Hide file tree
Showing 6 changed files with 109 additions and 12 deletions.
Binary file added Fig4/Immune3_COX_MP_SKCM06.rda
Binary file not shown.
Binary file added Fig4/Newman2015/Newman2015_COX_MP_SKCM06.rda
Binary file not shown.
60 changes: 53 additions & 7 deletions Fig4_Evaluation.R
Original file line number Diff line number Diff line change
Expand Up @@ -264,7 +264,7 @@ myinf <- list.files("~/Mypackage/MSofTimiGP/Fig4/",
pattern = "enrich_res.rda") %>%.[c(1,3, 7, 2, 5, 6, 4)]
# load data --------------------------------------------------------
DS <-
myinf %>% gsub(pattern = "/home/cli15/Mypackage/MSofTimiGP/Fig4/",
myinf %>% gsub(pattern = "~/Mypackage/MSofTimiGP/Fig4/",
replacement = "") %>%
gsub(pattern = "/enrich_res.rda",replacement = "") %>% strsplit("_") %>%
lapply("[[", 2) %>% unlist()
Expand Down Expand Up @@ -363,7 +363,7 @@ myinf <- list.files("~/Mypackage/MSofTimiGP/Fig4/",
pattern = "enrich_res.rda") %>%.[c(1,3, 7, 2, 5, 6, 4)]

DS <-
myinf %>% gsub(pattern = "/home/cli15/Mypackage/MSofTimiGP/Fig4/",
myinf %>% gsub(pattern = "~/Mypackage/MSofTimiGP/Fig4/",
replacement = "") %>%
gsub(pattern = "/enrich_res.rda",replacement = "") %>% strsplit("_") %>%
lapply("[[", 2) %>% unlist()
Expand Down Expand Up @@ -540,10 +540,32 @@ dev.off()

# [6] Charoentong2017_Bindea2013_Xu2018_Immune##################################
# The tutorial is example/example02_Charoentong2017_Bindea2013_Xu2018_Immune.R
# You can learn how to generate data(Immune3_COX_MP_SKCM06) there

# You can learn how to generate Immune3_COX_MP_SKCM06.rda there

#' COX regression Results from function TimiCOX with cell type meaker annotated by Charoentong2017_Bindea2013_Xu2018_Immune
#'
#' An intermediate result generated from function TimiCOX
#' that reveals the association between each marker pairs and favorable prognosis.
#'
#' @docType data
#'
#' @usage data(Immune3_COX_MP_SKCM06)
#'
#' @keywords intermediate result
#'
#' @format A list of 3 data frames
#' \describe{
#' \item{Row name}{Marker pair}
#' \item{HR}{Hazard.Ratio}
#' \item{PV}{P-Value}
#' \item{QV}{Adjust P-value}
#' }
#'
#'
#' @source intermediate result generated from function TimiCOX
# "Immune3_COX_MP_SKCM06"
rm(list=ls())
data(Immune3_COX_MP_SKCM06)
load("~/Mypackage/MSofTimiGP/Fig4/Immune3_COX_MP_SKCM06.rda")
cox_res <- Immune3_COX_MP_SKCM06

# a) TimiGP --------------------------------------------------------------------
Expand Down Expand Up @@ -791,12 +813,36 @@ dev.off()

#[7] Newman2015(Related to Figure S3) ==========================================
# The tutorial is example/example03_Newman2015_LM22.R
# You can learn how to generate data("Newman2015_COX_MP_SKCM06") there
# You can learn how to generate Newman2015_COX_MP_SKCM06.rda there


#' COX regression Results from function TimiCOX with cell type meaker annotated by Newman2015_LM22
#'
#' An intermediate result generated from function TimiCOX
#' that reveals the association between each marker pairs and favorable prognosis.
#'
#' @docType data
#'
#' @usage data(Newman2015_COX_MP_SKCM06)
#'
#' @keywords intermediate result
#'
#' @format A data frame with 115440 rows and 3 variables:
#' \describe{
#' \item{Row name}{Marker pair}
#' \item{HR}{Hazard.Ratio}
#' \item{PV}{P-Value}
#' \item{QV}{Adjust P-value}
#' }
#'
#'
#' @source intermediate result generated from function TimiCOX
# "Newman2015_COX_MP_SKCM06"

rm(list=ls())
outdir <- "~/Mypackage/MSofTimiGP/Fig4/Newman2015/"
dir.create(outdir)
data("Newman2015_COX_MP_SKCM06")
load(paste0(outdir,"Newman2015_COX_MP_SKCM06.rda"))
cox_res <- Newman2015_COX_MP_SKCM06

# a) Enrichment ----------------------------------------------------------------
Expand Down
Binary file not shown.
Binary file added Fig6/Zheng2021_Tcell/Tcell_COX_MP_SKCM06.rda
Binary file not shown.
61 changes: 56 additions & 5 deletions Fig6_Resolution.R
Original file line number Diff line number Diff line change
Expand Up @@ -10,11 +10,36 @@ library(TimiGP)

#[1] Tirosh2016_melanoma_TME(Related to Figure 6B,C)############################
# The tutorial is example/example04_Tirosh2016_melanoma_TME.R
# You can learn how to generate data("melanoma_TME_COX_MP_SKCM06") there
# You can learn how to generate melanoma_TME_COX_MP_SKCM06.rda there



#' COX regression Results from function TimiCOX with cell type meaker annotated by Tirosh2016_melanoma_TME
#'
#' An intermediate result generated from function TimiCOX
#' that reveals the association between each marker pairs and favorable prognosis.
#'
#' @docType data
#'
#' @usage data(melanoma_TME_COX_MP_SKCM06)
#'
#' @keywords intermediate result
#'
#' @format A data frame with 73536 rows and 3 variables:
#' \describe{
#' \item{Row name}{Marker pair}
#' \item{HR}{Hazard.Ratio}
#' \item{PV}{P-Value}
#' \item{QV}{Adjust P-value}
#' }
#'
#'
#' @source intermediate result generated from function TimiCOX
# "melanoma_TME_COX_MP_SKCM06"

rm(list=ls())
outdir <- "~/Mypackage/MSofTimiGP/Fig6/Tirosh2016_melanoma_TME/"
data("melanoma_TME_COX_MP_SKCM06")
load(paste0(outdir, "melanoma_TME_COX_MP_SKCM06.rda"))
cox_res <- melanoma_TME_COX_MP_SKCM06

# a) Enrichment ----------------------------------------------------------------
Expand Down Expand Up @@ -47,11 +72,37 @@ dev.off()

#[2] Zheng2021_Tcell(Related to Figure 6D,E)####################################
# The tutorial is example/example04_Zheng2021_Tcell.R
# You can learn how to generate data(Tcell_COX_MP_SKCM06) there
# You can learn how to generate Tcell_COX_MP_SKCM06.rda there




#' COX regression Results from function TimiCOX with cell type meaker annotated by Zheng2021_Tcell
#'
#' An intermediate result generated from function TimiCOX
#' that reveals the association between each marker pairs and favorable prognosis.
#'
#' @docType data
#'
#' @usage data(Tcell_COX_MP_SKCM06)
#'
#' @keywords intermediate result
#'
#' @format A data frame with 69006 rows and 3 variables:
#' \describe{
#' \item{Row name}{Marker pair}
#' \item{HR}{Hazard.Ratio}
#' \item{PV}{P-Value}
#' \item{QV}{Adjust P-value}
#' }
#'
#'
#' @source intermediate result generated from function TimiCOX
# "Tcell_COX_MP_SKCM06"

rm(list=ls())
outdir <- "~/Mypackage/MSofTimiGPt/Fig6/Zheng2021_Tcell/"
data(Tcell_COX_MP_SKCM06)
outdir <- "~/Mypackage/MSofTimiGP/Fig6/Zheng2021_Tcell/"
load(paste0(outdir, "Tcell_COX_MP_SKCM06.rda"))
cox_res <- Tcell_COX_MP_SKCM06

# a) Enrichment ----------------------------------------------------------------
Expand Down

0 comments on commit fff687d

Please sign in to comment.