Vous suivez un programme de biostatistiques ou un programme universitaire de santé publique/épidémiologie? Je peux calculer la probabilité de manière exacte mais je me demande si c’est ce qui est attendu vu que c’est un QCM. Si vous suivez un programme de mathématiques, c’est très certainement le cas. Si c’est un QCM pour le concours de PACES en fin d’année, cela m’étonnerait, ils s’attendent probablement à ce que vous preniez des raccourcis et répondiez à toute vitesse.
Le calcul exact va être super long car beaucoup les probabilités du problème sont liés. Dans le cas exact, l’information sur les fils change les probabilités sur les parents. Pour la suite je vais commencer le calcul exact de ces probabilités, vous allez voir que c’est très long. Je ne suis pas certain que cela vous sera utile (cela dépend du cursus que vous suivez et à quel point on s’attend à ce que vous fassiez un calcul exact).
Tout d’abord que sait-on? Le test ne créé pas de faux positif mais créé des faux négatifs. Il va falloir calculer la probabilité que la mère de A soit un faux négatif. Mais dans le calcul exact, cette probabilité n’est pas de 1/49 (comme pour une personne testée négative dont on ne sait rien de plus) puisqu’on sait que A1 et A2 ne sont pas homozygotes (sinon ils n’auraient pas atteint l’âge adulte) et on sait même qu’ils sont soit faux négatifs hétérozygotes, soit vrais négatifs. Toutes ces informations vont changer la probabilité que leur mère testé négative soit hétérozygote et rendre le calcul exact très pénible. Et comme on sait que le père est hétérozygote, cela change aussi la probabilité que les fils soient faux négatifs. Je suppose aussi que si les parents de A1 et A2 avaient eu un fils homozygotes qui n’avait pas atteint l’âge adulte, nous le saurions.
Pour faire le calcul exact, il va falloir procéder méthodiquement et poser des égalités sans essayer de faire de racourcis: on note F_A le père de A, M_A la mère de A, mêmes notations pour les parents de B. On note X\in\mathcal{H} l’évènement X est hétérozygote, et on note X\in\mathcal{T} l’évènement X testé positif.
On utilise le symbole \notin pour la négation de ces évènements. On souhaite d’abord calculer
P\big(A_1\in\mathcal{H}\mid A_2\notin\mathcal{T}\cap A_1\notin\mathcal{T}\cap A_2\notin\mathcal{T}\cap M_A\notin\mathcal{T}\cap P_A\in\mathcal{T}\big)
et
P\big(A_2\in\mathcal{H}\mid A_2\notin\mathcal{T}\cap A_1\notin\mathcal{T}\cap A_2\notin\mathcal{T}\cap M_A\notin\mathcal{T}\cap P_A\in\mathcal{T}\big)
Par symétrie, ces deux valeurs sont égales.
Il est clair que le risque de se tromper est grand si on essaie d’utiliser son intuition.
Il y a une question à laquelle je ne peux pas répondre car y répondre demande des connaissances sur le fonctionnement des tests, ce que je n’ai pas en tant que mathématicien: la probabilité d’être faux négatif est-elle indépendante de la mutation? IE, le test détecte-t-il 50% des mutations parmi celles qui causent la maladie ou detecte-il 50% des hétérozygotes indépendamment de la mutation qu’ils portent. Si le test a détecté qu’un parent est hétérozygote, et qu’un fils a hérité du gène détecté sur le parent, la probabilité que ce fils soit faux-négatif est-elle toujours de 1/2 ? Je vais supposer que c’est le cas. IE, je suppose que la probabilité qu’un test sur une personne donne un résultat faux négatif est complètement indépendant de l’état (test faux négatifs, hétérozygotie) des autres personnes.
Que sait-on?
- Pour tout individu X P(X\in\mathcal{H}\mid X\in\mathcal{T})=1 (pas de faux positifs).
- P(A_2\in\mathcal{H})=1/25.
- P(A_2\in\mathcal{H}\mid M_A\in\mathcal{H}\cap P_A\in\mathcal{H})=2/3 car on sait que A_2 n’est pas homozygote.
- P(A_2\in\mathcal{H}\mid M_A\notin\mathcal{H}\cap P_A\in\mathcal{H})=1/2 car on sait que A_2 n’est pas homozygote.
Faisons, nos calculs:
P\big(A_2\in\mathcal{H}\mid A_1\notin\mathcal{T}\cap A_2\notin\mathcal{T}\cap M_A\notin\mathcal{T}\cap F_A\in\mathcal{T}\big)\
=\frac{P\big(A_2\in\mathcal{H}\cap A_1\notin\mathcal{T}\cap A_2\notin\mathcal{T}\cap M_A\notin\mathcal{T}\cap F_A\in\mathcal{T}\big)}{P\big( A_1\notin\mathcal{T}\cap A_2\notin\mathcal{T}\cap M_A\notin\mathcal{T}\cap F_A\in\mathcal{T}\big)}\
=\frac{P\big(A_2\notin\mathcal{T}\mid A_1\notin\mathcal{T}\cap A_2\in\mathcal{H}\cap M_A\notin\mathcal{T}\cap F_A\in\mathcal{T}\big)P\big( A_1\notin\mathcal{T}\cap A_2\in\mathcal{H}\cap M_A\notin\mathcal{T}\cap F_A\in\mathcal{T}\big)}{P\big( A_1\notin\mathcal{T}\cap A_2\notin\mathcal{T}\cap M_A\notin\mathcal{T}\cap F_A\in\mathcal{T}\big)}\
=\frac{1/2\cdot P\big( A_1\notin\mathcal{T}\cap A_2\in\mathcal{H}\cap M_A\notin\mathcal{T}\cap F_A\in\mathcal{T}\big)}{P\big(A_1\notin\mathcal{T}\cap A_2\notin\mathcal{T}\cap M_A\notin\mathcal{T}\cap F_A\in\mathcal{T}\big)}\
=\frac{1/4\cdot P\big(A_1\notin\mathcal{T}\cap A_2\in\mathcal{H}\cap M_A\notin\mathcal{T}\cap F_A\in\mathcal{H}\big)}{1/2\cdot P\big( A_1\notin\mathcal{T}\cap A_2\notin\mathcal{T}\cap M_A\notin\mathcal{T}\cap F_A\in\mathcal{H}\big)}\
=\frac{1}{2}\cdot\frac{a+b+c+d}{(a+b+c+d)/2+w+x+y+z}
avec
a=P\big(A_1\notin\mathcal{H}\cap A_2\in\mathcal{H}\cap M_A\notin\mathcal{H}\cap F_A\in\mathcal{H}\big)\
b=P\big(A_1\in\mathcal{H}\cap A_2\in\mathcal{H}\cap M_A\notin\mathcal{H}\cap F_A\in\mathcal{H}\big)/2\
c=P\big(A_1\notin\mathcal{H}\cap A_2\in\mathcal{H}\cap M_A\in\mathcal{H}\cap F_A\in\mathcal{H}\big)/2\
d=P\big(A_1\in\mathcal{H}\cap A_2\in\mathcal{H}\cap M_A\in\mathcal{H}\cap F_A\in\mathcal{H}\big)/4
et
w=P\big(A_1\notin\mathcal{H}\cap A_2\notin\mathcal{H}\cap M_A\notin\mathcal{H}\cap F_A\in\mathcal{H}\big)\
x=P\big(A_1\in\mathcal{H}\cap A_2\notin\mathcal{H}\cap M_A\notin\mathcal{H}\cap F_A\in\mathcal{H}\big)/2\
y=P\big(A_1\notin\mathcal{H}\cap A_2\notin\mathcal{H}\cap M_A\in\mathcal{H}\cap F_A\in\mathcal{H}\big)/2\
z=P\big(A_1\in\mathcal{H}\cap A_2\notin\mathcal{H}\cap M_A\in\mathcal{H}\cap F_A\in\mathcal{H}\big)/4
Je ne sais pas si cela vaut le coup que je continue. Il y a encore du travail pour calculer ces 8 valeurs Mais je pense que c’est faisable. N’oubliez pas qu’implicitement, tout est conditionné par le fait qu’il n’y a pas d’homozygotes dans la famille (les formules auraient été trop longues à écrire si on avait explicité ce conditionnement).
Une fois calculé a,b,c,d,w, x,y,z, il faudra faire de même pour la famille B. Et après ce sera facile de calculer la probabilité pour leur enfant d’être homozygote.