Question: Hypergeometric Test on Two Mutational Spectra
0
gravatar for Clint Valentine
5.5 years ago by
United States
Clint Valentine10 wrote:

I have two mutational spectra that I need to compare to see if they are significantly different from one another. Each column of the data below illustrates how many mutations were found at a given location in each sample. I am aware that Scipy has hypergeometric functions that will suit my needs but with my limited stats knowledge I am having a hard time distilling this raw information into the function to acquire the p-value.

Ultimately I am trying to write a script using Biopython/Scipy/Numpy to do what is outlined in the journal article Statistical Test for the Comparison of Samples from Mutational Spectra by W. Thomas Adams and Thomas R. Skopek in 1986.

Here is a sample data-set to work with that resembles the data I am looking at.

Pos.     Sample 1        Sample 2
2           0                3
3           0                0
6           0                1
8           0                0
12          2                5
15          1                1
26          1                0
34          2                0
47          0                2
77          4                4

 

Let me know if more details need to be provided!

ADD COMMENTlink modified 5.5 years ago • written 5.5 years ago by Clint Valentine10

To give an update on my progress. I am developing a script to handle the Adams & Skopek algorithm to answer my own question. I will post my answer once it's complete and working.

ADD REPLYlink written 5.5 years ago by Clint Valentine10
0
gravatar for Zev.Kronenberg
5.5 years ago by
United States
Zev.Kronenberg11k wrote:

Joint allele-frequency spectrum are cool.  Are you aware of the paper below?  The program associated with the paper (DADI) might have the test you are looking for.  If not the paper is still a good read.

 

http://www.plosgenetics.org/article/info%3Adoi%2F10.1371%2Fjournal.pgen.1000695

 

https://code.google.com/p/dadi/

ADD COMMENTlink written 5.5 years ago by Zev.Kronenberg11k

I'll give it a read and see if I can get anywhere. Thanks for the tip.

ADD REPLYlink written 5.5 years ago by Clint Valentine10
Please log in to add an answer.

Help
Access

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