TCGA 临床数据的下载与整理
<h1>背景</h1><p>前面我们下载了新版tcga的rnaseq 数据,同时下载了样本的临床数据，此次就可以对临床数据进行整理了。</p>
<h1>一、临床数据整理</h1>
<pre><code>rm(list = ls())####魔幻操作，一键清空~
getwd()
setwd('C:/Users/Administrator/project_gdc2')
#先把下载的临床数据文件解压
tar_file <- "clinical.cart.2024-12-22.tar.gz"# 导入tar.gz文件
dir.create('clinical_data')
extract_dir <- "./clinical_data" #解压后存放位置
untar(tar_file, exdir = extract_dir) #解压tar.gz文件
# 提取临床数据整理生存分析需要的数据
setwd("C:/Users/Administrator/project_gdc2/clinical_data") ##设置路径
library(readr)
library(dplyr)
#将下载好的metadata.json文件放入clinical文件夹
# install.packages("jsonlite")
library(jsonlite)
json <- jsonlite::fromJSON("../metadata.cart.2024-12-22.json") #读取JSON文件
entity_submitter_id <- sapply(json$associated_entities, function(x) unlist(x[, 1]))
case_id <- sapply(json$associated_entities, function(x) unlist(x[, 3]))
sample_case <- t(rbind(entity_submitter_id, case_id))
clinical <- read_tsv('clinical.tsv') #读取tsv文件
clinical <- as.data.frame(clinical[!duplicated(clinical$case_id),]) #去除重复的sample
# 371
str(sample_case)# 查看sample_case的结构
str(clinical) # 查看clinical的结构
sample_case <- as.data.frame(sample_case)
## 将sample_case$case_id和clinical$case_id转化为字符串格式，便于后续操作。
sample_case$case_id <- as.character(sample_case$case_id)
clinical$case_id <- as.character(clinical$case_id)
matrix <- merge(sample_case,clinical,by="case_id",all.x=T)
#424
colnames(clinical)
demo <- c("case_submitter_id","age_at_index","ethnicity","gender","race",
"vital_status","days_to_death","days_to_last_follow_up",
"ajcc_pathologic_stage","ajcc_pathologic_t","ajcc_pathologic_m",
"ajcc_pathologic_n","treatment_type")
matrix = matrix[,demo] #筛选需要的临床信息
dim(matrix)# 424 13
colnames(matrix) <- c("ID","Age","Ethnicity","Gender","Race",
"Status","days_to_death","days_to_last_follow_up",
"Stage","T","M","N","Treatment")
#排除结局为"Not Reported"的Sample,保留Alive和Dead的数据
matrix = matrix
# 422
# 把matrix数值列转换为数值型，便于记录生存信息
matrix$days_to_last_follow_up <- as.numeric(matrix$days_to_last_follow_up)
matrix$days_to_death <- as.numeric(matrix$days_to_death)
matrix$Age <- as.numeric(matrix$Age)
# 去除NA，替换为0
matrix$days_to_last_follow_up = 0
matrix$days_to_death = 0
matrix$Age = 0
matrix$days <- ifelse(matrix$Status=='Alive',matrix$days_to_last_follow_up,matrix$days_to_death)
## 添加生存分析需要的信息：存活状态、月、年
matrix$OS <- ifelse(matrix$Status == "Alive", 0, 1)
matrix$month=round(matrix$days/30,0) #以month为单位，小数不保留
matrix$OS.time <- floor(matrix$month/12)
floor(12.1)#12 整数年，也可以不取整。
head(matrix)
# 保存临床信息
fwrite(matrix,"../RawData/03.LIHC_clin.txt") # txt格式
write.csv(matrix, "../RawData/csv/03.LIHC_clin.csv", row.names = F) # csv格式
</code></pre>
<p>结果如下：<br />
<img src="https://roim-picx-bpc.pages.dev/rest/bXxgSlK.png" alt="" /><br />
后续可以进一步进行生存分析及可视化等操作。</p>
<p>自38行单独提取出肿瘤组织 的数据</p>
<pre><code>matrix <- merge(sample_case,clinical,by="case_id",all.x=T)
#424
# 去除癌旁单独保存一个原始的数据格式
library(stringr)
group_list <- ifelse(as.numeric(str_sub(matrix$entity_submitter_id,14,15))<10,'tumor','normal')
table(group_list)
matrix_T <- matrix
fwrite(matrix_T,"../RawData/03.LIHC_clin_T.txt") # 肿瘤的txt格式
</code></pre>
<p>tcga癌旁是配对了癌组织，是一个人。所以一般只选择癌组织数据进行生存分析。</p>
页:
[1]