R語言計算tmb值

tmb值與免疫檢查點抑制劑療效相關，而TCGA數據庫中的tmb值可以通過TCGAmutations包來計算

# =================================================================
#   
#   
# 
#   
#   
# =================================================================



rm(list=ls())

setwd("D:\\SCIwork\\F25\\results\\tmb")
# BiocManager::install("PoisonAlien/TCGAmutations")
# TCGAmutations包整合了TCGA中全部樣本的maf文件
# devtools::install_github(repo = "PoisonAlien/TCGAmutations")
library(TCGAmutations)

tmp=as.data.frame(tcga_available())

dt <- TCGAmutations::tcga_load(study = "KIRC")

dt <- dt@data

dt1 <- as.data.frame( table(dt$Tumor_Sample_Barcode))

names(dt1) <- c('Barcode', 'Freq')

dt1$tmb <- dt1$Freq/38

names(dt1)
 
 write.csv(dt1, file = 'KIRC_TMB.csv')