Entering edit mode
2.8 years ago
melissachua90
▴
70
How do I find the most common mutations in a MAF dataset? Should I use pfamDomains
from maftools
?
library(maftools)
d <- merge_mafs(lapply(Sys.glob("mafs/Patient*.maf"), read.maf))
d.pfam = pfamDomains(maf=d, AACol="Protein_Change", top=15)
d.pfam$proteinSummary[,1:7, with=FALSE]