Entering edit mode
4.3 years ago
kittys
•
0
Hi I am trying to analyse double base substitution (DBS) signatures using the deconstruct sigs R package. I have run the example code below with success:
sigs.input = mut.to.sigs.input(mut.ref = sample.mut.ref, sample.id = "Sample",
chr = "chr", pos = "pos", ref = "ref", alt = "alt", bsg =
BSgenome.Hsapiens.UCSC.hg19, sig.type = 'DBS', dbs_table = dbs_possible)
However when I try to run this on my own data, I get a completely empty output. My data is in the same format with the same column names and data types per column. No errors come up just the output is empty. The only other thing I have to change is the reference genome as my data uses hg38.
When I change the sig.type to SBS with my data, this works fine. So it is something specifically wrong with the DBS part.
Any suggestions would be really helpful!
Thanks, Kitty
