rm(list = ls())
library(TCGAbiolinks)
f = "masked.Rdata"
if(!...<- GDCprepare(query)
colnames(seg)
seg = seg[,c(7,2:6)]
save(seg,file = f)
write.table(seg,...$Sample)) #人数
## [1] 85
#如果只需要tumor样本那么就:
#k = as.numeric(str_sub(seg$Sample,14,15))<10
#seg = seg[k...marker$freqcnv,] #只保留FALSE
marker =marker[,1:3]
colnames(marker) = c("Marker Name","Chromosome","Marker....*", "", rownames(cnvData)),代表把基因名中的点号后面的部分给去掉,这其实是给ensembel矩阵准备的。但如果是symbol为行名的矩阵,就要出错。
if(!