Python program to ascertain the paternity index in a single mutation case
PYTHON PROGRAMMING CODE TO ASSERTAIN THESE:
person genotype
mother ab
child bc
alleged father de
X = P(man without C will contribute C)
X = P(contributed gene will mutate) * P(mutated gene will be a C)
m = observed rate of mutations/meiosis for the locus
P(mutated gene will be a C) ie. Frequency of C allele = C
X = C * m
For the numerator (X) Probability = 2ab x 2de x 0.5 x m x C
In order to explain this evidence the denominator must calculate the probability that the paternal allele is C and a random man would have a genotype inconsistent with paternity at this locus Y = P(paternal allele is C and random man has no C allele) = P(paternal gene is C) * P(random man has no C allele)
P(paternal allele will be a C) ie. Frequency of C allele = C
P(random man has no C allele) = probability of exclusion
Y = C * A
For denominator (Y) Probability = 2ab x 2de x 0.5 x C x A
Hence the final result should be:
PI = X/Y
PI = (2ab * 2de * 0.5 * C * m ) / (2ab * 2de * 0.5 * C * A)
PI = m/A
N:B where PI is Paternity Index
• 1,914 views
•
link
0 answers
No answers yet.
Log in to answer this question.
Is this homework?
It is a practice work of research
You won't get much practice if someone just gives you the code. Try to do it by yourself and ask a question if you get stuck. Also, seems like the problem is mostly just programming so most likely this would be better suited to StackOverflow.
Hello leequest77!
We believe that this post does not fit the main topic of this site.
Really old post that went nowhere
For this reason we have closed your question. This allows us to keep the site focused on the topics that the community can help with.
If you disagree please tell us why in a reply below, we'll be happy to talk about it.
Cheers!