subsetMaf (and depending functions?) not working on data imported with annovarToMaf about maftools

[PoisonAlien]

Hi,

Thanks for reporting this !
This was a mistake from my side, should have been more careful. You were right about adding hugo_symbol to barcodes. I have fixed it. Thanks again.

P.S: that after subsetting in the above example you may not be be able to draw heat maps or rainfall plots since you're left with only one sample.

[sribi]

Hi,
Thanks a lot for fixing this so promptly! subsetMaf() and oncoplot() now work as expected. However, I am still having issues with 'rainfallPlot()':

#import
> var.annovar.maf <- annovarToMaf(annovar = var.annovar, Center = 'CSI-NUS', refBuild = 'hg19', tsbCol = 'Tumor_Sample_Barcode', table = 'ensGene', header = TRUE, MAFobj=TRUE)

#plot
> rainfallPlot(maf=var.annovar.maf, detectChangePoints = TRUE, fontSize = 12, pointSize = 0.6)
Error in seg.spl[[1]] : subscript out of bounds
In addition: Warning message:
In is.factor(x) : NAs introduced by coercion

#check if more than one sample 
> getSampleSummary(var.annovar.maf)
   Tumor_Sample_Barcode Frame_Shift_Del Frame_Shift_Ins Missense_Mutation total
1:               fake_4               1               0                 1     2
2:               fake_5               0               1                 1     2
3:               fake_7               0               1                 1     2
4:               fake_1               0               0                 1     1
5:               fake_2               0               0                 1     1
6:               fake_3               0               0                 1     1
7:               fake_6               0               0                 1     1

I was assuming the issue might be due to the low number of mutations in the test dataset, but I'm running into the same issue with my own, larger dataset with plenty of mutations...

Thanks again,
sribi

[PoisonAlien]

Hi,
This was an issue with empty chromosomes. i.e, chromosomes with no variants resulted in error. I have fixed it. I would really appreciate if you could try installing it again and let me know if it works.

Target level branch, since I haven't merged it with master yet.

devtools::install_github(repo = "PoisonAlien/maftools", ref = 'devel')

Thanks again. This was helpful.

[sribi]

Hi,

The bug in rainfallPlot() , as described above, seems to be fixed in your latest push to devel!

Thanks for fixing it so quickly.