tri.counts.method parameter in deconstructSigs
1
0
Entering edit mode
@jane-merlevede-5019
Last seen 5.5 years ago

Dear deconstructSigs users,

I am using deconstructSigs on WES data from hypermutated tumors:

whichSignatures = whichSignatures(tumor.ref = sigs.input, signatures.ref = signatures.cosmic, sample.id = "MyPatient", associated = c(),signatures.limit = NA, signature.cutoff = 0, contexts.needed = TRUE, tri.counts.method = "exome2genome")

In input, I gave the number of mutations found for the 96 trinucleotides, without any normalization.

I hesitate between three values "default", "exome" and "exome2genome" for tri.counts.method parameter.

 

In the vignette description from August 29, 2016 (https://cran.r-project.org/web/packages/deconstructSigs/deconstructSigs.pdf), one can read in the description of whichSignatures(): " For exome data, the ’exome2genome’ method is appropriate for the signatures included in this package.

On github, one can read: "For exome data, the default method is appropriate for the signatures included in this package."

So default or exome2genome? Is it related only to the dataset included in the package?

 

My WES data were generated from V5 Clinical Research Exome capture kit, ~50% of the mutations are non coding (mainly in intron). Do you recommend to exclude the non coding mutations because of their trinucleotide repartition differences?

 

Thank you in advance for your precisions.

deconstructsigs • 1.6k views
ADD COMMENT
0
Entering edit mode

I would contact the maintainer / developer directly. She is not registered in these forums, as far as I'm aware. Also, this package is on CRAN, not Bioconductor.

ADD REPLY
0
Entering edit mode
boostoo • 0
@boostoo-12289
Last seen 7.1 years ago

Dear,

I also have a similar question on that parameter. In my case we have WGS data so should we use 

tri.counts.method = "genome"
​instead of the "default"?

Thank you for your help,

Bastien 

 

ADD COMMENT
0
Entering edit mode

I would contact the maintainer / developer directly. She is not registered in these forums, as far as I'm aware. Also, this package is on CRAN, not Bioconductor.

ADD REPLY

Login before adding your answer.

Traffic: 711 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6