maftools icon indicating copy to clipboard operation
maftools copied to clipboard

Subset genes in the initial steps

Open camcecc opened this issue 2 years ago • 5 comments

Is there a way to subset genes and remove some of them (namely) from the maf file? Or, the opposite, only keep the ones we want to investigate more and would like to see more deeply (i.e. in the oncoprint)? Thanks

camcecc avatar Jul 14 '22 12:07 camcecc

Hi,

Check out filterMaf function. I hope this helps.

#remove genes of intrest
filterMaf(maf = laml, genes =c("TTN", "AHNAK2"))

PoisonAlien avatar Jul 14 '22 12:07 PoisonAlien

Wow that was quick! Thank you so much.

How about selecting genes from the ratio mutation / length? I know this would be MutSig material, but is there a way (maybe simplified) to obtain this data from maftools?

camcecc avatar Jul 14 '22 13:07 camcecc

Hi,

This is a very specific case. But if you know what genes, samples, or genomic regions, you are particularly interested in, you could use subsetMaf.

PoisonAlien avatar Jul 14 '22 13:07 PoisonAlien

I see. Is there some point in the pipeline where I can find the genes length, so I can manually do the ratio and then subset? Thank you again!

camcecc avatar Jul 14 '22 13:07 camcecc

Unfortunately it is not possible :|

PoisonAlien avatar Jul 14 '22 13:07 PoisonAlien

This issue is stale because it has been open for 60 days with no activity.

github-actions[bot] avatar Aug 07 '23 02:08 github-actions[bot]

This issue was closed because it has been inactive for 14 days since being marked as stale.

github-actions[bot] avatar Aug 22 '23 01:08 github-actions[bot]