Question: Questions about SomaticSignatures package
gravatar for Asma rabe
2.8 years ago by
Asma rabe290
Asma rabe290 wrote:


I would like to thank the authors of SomaticSignatures  for the helpful package and for support on bioconductor forum.

Although the vignette is clear, I suggest that you add an example on how to start from text file in this format:

Chr Stt end Ref Var SampleID study

This will be much more easy for many readers to start processing there own data



sca_vr = VRanges(seqnames=Data$Chr,ranges=IRanges(start=Data$Start,width=1),ref=Data$Ref,alt=Data$Var,study=Data$study)



other steps are as exactly in vignette.


In case of using fast to get mutation context, I could not find mutational context function

sca_vr = VRanges(seqnames=Data$Chr,ranges=IRanges(start=Data$Start,width=1),ref=Data$Ref,alt=Data$Var,study=Data$study)

fa_A = FastaFile("human_hg19.ucsc.fa")

vr_A = mutationalContext(sca_vr, fa_A)
Error: could not find function "mutationalContext"

When i tried mutationContext , i got this error

vr_A = mutationContext(sca_vr,fa_A)
Error in (function (classes, fdef, mtable)  : 
  unable to find an inherited method for function ‘getSeq’ for signature ‘"FastaFile"’

Any idea??


In mutational signatures plot y-axis represents the contribution of each motif in the signatures. How this contribution is estimated?

Thank you very much.

ADD COMMENTlink modified 2.7 years ago • written 2.8 years ago by Asma rabe290

Hi Asma,

Have you figure out how to input your own data with VRanges?

I also faced with same problem. If you already success, please let me know.






ADD REPLYlink written 2.2 years ago by shengfeng0
Please log in to add an answer.


Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.2.0
Traffic: 305 users visited in the last hour