| 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') | |