This is a test version of Biostars. For the public version, visit https://www.biostars.org.
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

assembly chip-seq sequence next-gen rna-seq

Is this homework?

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!

0 answers

No answers yet.

Log in to answer this question.