library(readr)
# 读取数据
mathylation <- read.csv("pca_principal_components.csv")
data <- read_tsv("TCGA-ACC.survival.tsv", na = c("NaN", "null", ""))
data <- read_tsv("TCGA-LGG.survival.tsv", na = c("NaN", "null", ""))
data = data[data$OS == 1, ]
# 更改列名
colnames(mathylation)[1] = 'sample'
# 合并数据
data = merge(data, mathylation, by = 'sample')
data = data[, -c(1, 2, 3)]
data[is.na(data)] <- 0
# 保存结果
write.csv(data, 'data.csv')