Niinleslie / MesKit

A tool kit for dissecting cancer evolution from multi-region derived tumor biopsies via somatic mutations
GNU General Public License v3.0
35 stars 9 forks source link

@Niinleslie readMaf function has issues with return object with right number of tumor samples per patient #255

Closed avaghasia1 closed 2 years ago

avaghasia1 commented 2 years ago

Hi All, I am having issues reading in MAFtools generated merged file from multiple patients. The readMaf gives error even-though I have more than 2 samples per patient. Second, thing I notice is it removes one of tumor sample and never gets stored. Is there any filtering on mutations while reading in maf files? I tried to debug but couldn't figure it out. I think issue is happening post this line of code in readMaf function. Can you please help find solution for this? maf <- Maf(data = data.table::setDT(data), sample.info = as.data.frame(sample.info), nonSyn.vc = nonSyn.vc, ref.build = refBuild) maf_patient_list[[patient]] <- maf