Un ami m'a posé la question suivante. Je n'ai pas pu l'aider mais j'espère que quelqu'un pourra me l'expliquer. Je n'ai pas pu trouver d'exemple similaire. Merci pour toute aide et explication.
Q: Les résultats de 100 expériences de lancer de pièces sont enregistrés comme 0 = "Queue" et 1 = "Tête". La sortie x est une chaîne de 0 et 1 de longueur 100. Et le nombre de fois où nous obtenons 1-0-0 en x est calculé et il est de 20 (ex: si x = (001001110100), 1-0-0 se produit 2 fois). Pensez-vous que c'est une pièce équitable?
probability
inference
bernoulli-distribution
Jimmy Dur
la source
la source
Réponses:
Résoudre le problème par simulation
Ma première tentative serait de simuler cela sur un ordinateur, ce qui peut retourner très rapidement de nombreuses pièces justes. Voici un exemple avec un million d'essais. L'événement « que le nombre de fois le motif « 1-0-0 » se produit dans pièces flips est de 20 ou plus » se produit à peu près une fois tous les trois mille essais, donc ce que vous avez observé est peu probable (pour une foire pièce de monnaie).X n = 100
Notez que l'histrogramme est pour la simulation et la ligne est le calcul exact expliqué plus loin ci-dessous.
Résoudre le problème avec un calcul exact
Pour une approche analytique, vous pouvez utiliser le fait que 'la probabilité d'observer 20 séquences ou plus' 1-0-0 'dans 100 tours de pièces est égale au 1 moins la probabilité qu'il faut plus de 100 tours pour faire 20 séquences' . Ceci est résolu dans les étapes suivantes:
Temps d'attente pour la probabilité de basculer «1-0-0»
La distribution, , du nombre de fois que vous devez retourner jusqu'à ce que vous obteniez exactement une séquence '1-0-0' peut être calculée comme suit:FN, x = 1( n )
Analysons les moyens d'arriver à «1-0-0» en tant que chaîne de Markov. Nous suivons les états décrits par le suffixe de la chaîne de flips: «1», «1-0» ou «1-0-0». Par exemple, si vous avez les huit flips 10101100 suivants, vous avez passé, dans l'ordre, les huit états suivants: "1", "1-0", "1", "1-0", "1", "1", '1-0', '1-0-0' et il a fallu huit flips pour atteindre '1-0-0'. Notez que vous n'avez pas la même probabilité d'atteindre l'état «1-0-0» à chaque flip. Ainsi, vous ne pouvez pas modéliser cela comme une distribution binomiale . Au lieu de cela, vous devez suivre un arbre de probabilités. L'état '1' peut aller dans '1' et '1-0', l'état '1-0' peut aller dans '1' et '1-0-0', et l'état «1-0-0» est un état absorbant. Vous pouvez l'écrire comme:
et la probabilité d'atteindre le modèle «1-0-0», après avoir obtenu un premier «1» (vous commencez par l'état «0», sans avoir encore renversé la tête), en flips, la moitié de la probabilité être dans l'état «1-0» dans flips:n n - 1
où est le ème nombre de Fibonnaci. La probabilité non conditionnelle est une sommeFje je
Temps d'attente pour la probabilité de retourner fois '1-0-0'k
Vous pouvez le calculer par convolution.
vous obtiendrez comme probabilité d'observer au moins 20 modèles `` 1-0-0 '' (selon l'hypothèse que la pièce est juste)
Voici le code R pour le calculer:
Calcul de pièces déloyales
Nous pouvons généraliser le calcul ci-dessus de la probabilité d'observer modèles dans flips, lorsque la probabilité de «1 = tête» est et que les flips sont indépendants.X n p
Nous utilisons maintenant une généralisation des nombres de Fibonacci:
les probabilités sont maintenant les suivantes:
et
Lorsque nous traçons cela, vous obtenez:
Ainsi, bien que la valeur de p soit petite pour une pièce de monnaie équitable 0,0003247, nous devons noter qu'elle n'est pas beaucoup meilleure (une seule commande) pour différentes pièces de monnaie déloyales. Le rapport de vraisemblance, ou facteur de Bayes , est d'environ 11 lorsque l'hypothèse nulle ( ) est comparée à l'hypothèse alternative . Cela signifie que le rapport de cotes postérieur n'est que dix fois plus élevé que le rapport de cotes précédent.p = 0,5 p = 0,33
Ainsi, si vous pensiez avant l'expérience que la pièce était improbable, alors vous devriez maintenant penser que la pièce est improbable.
Une pièce avec mais une injustice concernant les événements '1-0-0'ph e a ds= pt a i l s
On pourrait beaucoup plus facilement tester la probabilité d'une pièce équitable en comptant le nombre de têtes et de queues et utiliser une distribution binomiale pour modéliser ces observations et tester si l'observation est particulière ou non.
Cependant, il se peut que la pièce retourne, en moyenne, un nombre égal de têtes et de queues, mais n'est pas juste en ce qui concerne certains modèles. Par exemple, la pièce peut avoir une certaine corrélation pour les tours de pièces suivants (j'imagine un mécanisme avec des cavités à l'intérieur du métal de la pièce qui sont remplies de sable qui coulera comme un sablier vers l'extrémité opposée du tour de pièce précédent, qui charge la pièce pour tomber plus probablement du même côté que le côté précédent).
Soit le premier flip de pièces de monnaie à probabilité et à queues égales et les flips suivants sont avec probabilité du même côté que le flip précédent. Ensuite, une simulation similaire au début de ce post donnera les probabilités suivantes pour le nombre de fois que le modèle '1-0-0' dépasse 20:p
Vous pouvez voir qu'il est possible de le rendre légèrement plus susceptible d'observer le modèle `` 1-0-0 '' (quelque part autour de une pièce qui a une certaine corrélation négative), mais plus dramatique est que l'on peut le rendre beaucoup moins susceptibles d’observer le schéma «1-0-0». Pour un faible, vous obtenez plusieurs fois la queue après une tête, la première partie «1-0» du motif «1-0-0», mais vous n'obtenez pas si souvent deux queues d'affilée le «0-0» partie du motif. L'inverse est vrai pour les valeurs élevées.p = 0,45 p p
Utiliser les mathématiques dans les statistiques
Ce qui précède est très bien, mais ce n'est pas une réponse directe à la question
Pour répondre à cette question, on peut utiliser les mathématiques ci-dessus mais il faut vraiment d'abord très bien décrire la situation, les objectifs, la définition de l'équité, etc. Sans aucune connaissance du contexte et des circonstances, tout calcul ne sera qu'un exercice de mathématiques et non une réponse à la question explicite.
Une question ouverte est pourquoi et comment nous recherchons le modèle «1-0-0».
la source