Closed ebrintn closed 1 week ago
Should be working now - @afarineshp can you check to see if you are getting the error you were getting before?
Fixed all bugs - see iss #37. Closing.
Talked in meeting about needing to add denominators from stationary distributions.
These values are really weird when implemented. I think we need to meet and talk more.
Not all codon distributions will be the same. Change to have users specify codon stationary distributions instead.
Updated with most recent push to include mutation matrices!
Counted each value of dN/dS and output the specific substitutions. Users can now see the denominators. I also fixed a bug in the way that fixations are determined. This may change the way that we were seeing substitutions before.
I definately think that it needs to be multiplied not divided by the selection coefficient we calculated
I started thinking about this again and I realize where I was wrong last time - the denominator of the dN/dS values should be individual. I will fix today
Talked with @jasondk and this is correct except we need to multiply by 3 instead of divide.
Use fixation counting commands to calculate dN and dS