一牢屋、isTCGA是個什么參數(shù)族扰,到底什么時候加?
首先,對讀取gz格式還是txt格式?jīng)]有啥特別要求
image.png
二三热、但是讀取后支鸡,樣本名稱有顯著改變鼠渺,即sampl barcode有差異
isTCGA = T
image.png
isTCGA = F
及為默認值image.png
嘻嘻湾盗,在這里可能暫時看不出來有啥用,但你挑選標本就好用多了
這里報錯的原因就是捏萍,maf和你提供的樣本名不一致導致的
image.png
三太抓、這里提供一個篩選的標準流程
#---臨床信息獲取,并提取行名1——12位-----------------
clindata$sample <- sub_str(rownames(clindata),1,12)
table(clindata$group)
#---拿到想要的樣本名字-----------------------
ID1 <-clindata$sample[clindata$group=="high"]
ID2 <- clindata$sample[clindata$group=="low"]
ID <- clindata$sample
#-------提取特定樣本的突變信息-----------------
maf1 <- read.maf('All_maf.txt',isTCGA = T)%>%subsetMaf(tsb = ID1)
maf2 <- read.maf('All_maf.txt')%>%subsetMaf(tsb = ID2)