| 123456789101112 |
- #将文件转换为maf格式
- var.annovar.maf = annovarToMaf(annovar = "all_annovar2", Center = 'NA', refBuild = 'mm10', tsbCol = 'Tumor_Sample_Barcode', table = 'refGene',sep = "\t")
- write.table(x=var.annovar.maf,file="var_annovar_maf",quote= F,sep="\t",row.names=F)
- var_maf = read.maf(maf ="var_annovar_maf")
- #概览maf文件
- plotmafSummary(maf = var_maf, rmOutlier = TRUE, addStat = 'median')
|
