ShixiangWang / sigminer

🌲 An easy-to-use and scalable toolkit for genomic alteration signature (a.k.a. mutational signature) analysis and visualization in R https://shixiangwang.github.io/sigminer/reference/index.html
https://shixiangwang.github.io/sigminer/
Other
147 stars 19 forks source link

Error in get_sig_similarity() function #389

Closed aditisk closed 3 years ago

aditisk commented 3 years ago

Hi @ShixiangWang,

I am using some of my scripts that worked previously with v2.1.1 and getting some error in the get_sig_similarity() function. The error is pasted below, I tried a few different databases but that doesn't seem to fix the issue.

all_sim <- get_sig_similarity(all_sig) Error in get_sig_similarity(all_sig) : The following components cannot be found in reference! A[C>A]A A[C>A]C A[C>A]G A[C>A]T A[C>G]A A[C>G]C A[C>G]G A[C>G]T A[C>T]A A[C>T]C A[C>T]G A[C>T]T A[T>A]A A[T>A]C A[T>A]G A[T>A]T A[T>C]A A[T>C]C A[T>C]G A[T>C]T A[T>G]A A[T>G]C A[T>G]G A[T>G]T C[C>A]A C[C>A]C C[C>A]G C[C>A]T C[C>G]A C[C>G]C C[C>G]G C[C>G]T C[C>T]A C[C>T]C C[C>T]G C[C>T]T C[T>A]A C[T>A]C C[T>A]G C[T>A]T C[T>C]A C[T>C]C C[T>C]G C[T>C]T C[T>G]A C[T>G]C C[T>G]G C[T>G]T G[C>A]A G[C>A]C G[C>A]G G[C>A]T G[C>G]A G[C>G]C G[C>G]G G[C>G]T G[C>T]A G[C>T]C G[C>T]G G[C>T]T G[T>A]A G[T>A]C G[T>A]G G[T>A]T G[T>C]A G[T>C]C G[T>C]G G[T>C]T G[T>G]A G[T>G]C G[T>G]G G[T>G]T T[C>A]A T[C>A]C T[C>A]G T[C>A]T T[C>G]A T[C>G]C T[C>G]G T[C>G]T T[C>T]A T[C>T]C T[C>T]G T[C>T]T T[T>A]A T[T>A]C T[T>A]G T[T>A]T T[T>C]A T[T>C]C T[T>C]G T[T>C]T T[T>G]A T[T>G]C T[T>G]G T[T>G]T NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA

I used sig_auto_extract() to get "all_sig" after making a tally from my MAF files. Please let me know how to proceed. Thanks a lot.

ShixiangWang commented 3 years ago

thanks for reporting. Could you provide your data so i can check?

@aditisk

aditisk commented 3 years ago

Hi @ShixiangWang, just sent you an email on shixiang1994wang@gmail.com to continue this discussion privately since these data are unpublished. Thanks for your help.

aditisk commented 3 years ago

Here is the solution in case anyone encounters the same error:

get_sig_similarity(all_sig, sig_db = “TSB")