Foi-me feita a seguinte pergunta por um amigo. Não pude ajudá-la, mas espero que alguém possa me explicar. Não encontrei nenhum exemplo semelhante. Obrigado por qualquer ajuda e explicação.
P: Os resultados de 100 experimentos de sorteio são registrados como 0 = "Cauda" e 1 = "Cabeça". A saída x é uma sequência de 0 e 1 de comprimento 100. E o número de vezes que obtemos 1-0-0 em x é calculado e é 20 (ex: se x = (001001110100), 1-0-0 ocorre 2 vezes). Você acha que essa é uma moeda justa?
probability
inference
bernoulli-distribution
Jimmy Dur
fonte
fonte
Respostas:
Resolvendo o problema por simulação
Minha primeira tentativa seria simular isso em um computador, que pode jogar muitas moedas justas muito rapidamente. Abaixo está um exemplo com um milhão de tentativas. O evento 'que o número de vezes o padrão '1-0-0' ocorre em coin flips é de 20 ou mais' ocorre aproximadamente uma vez a cada três mil ensaios, de modo que você observou não é muito provável (para uma feira moeda).X n = 100
Observe que o histrograma é para a simulação e a linha é o cálculo exato explicado mais adiante.
Resolvendo o problema com um cálculo exato
Para uma abordagem analítica, você pode usar o fato de que 'a probabilidade de observar 20 ou mais seqüências' 1-0-0 'em 100 lançamentos de moedas é igual a 1 menos a probabilidade de que são necessários mais de 100 lançamentos para fazer 20 sequências' . Isso é resolvido nas seguintes etapas:
Tempo de espera pela probabilidade de virar '1-0-0'
A distribuição, , do número de vezes que você precisa inverter até obter exatamente uma sequência '1-0-0' pode ser calculada da seguinte forma:fN, x = 1( N )
Vamos analisar as maneiras de chegar ao '1-0-0' como uma cadeia de Markov. Seguimos os estados descritos pelo sufixo da sequência de lançamentos: '1', '1-0' ou '1-0-0'. Por exemplo, se você tiver os oito lançamentos a seguir 10101100, passou por oito estados seguintes: '1', '1-0', '1', '1-0', '1', '1', '1-0', '1-0-0' e foram necessários oito movimentos para chegar a '1-0-0'. Observe que você não tem a mesma probabilidade de atingir o estado '1-0-0' em cada jogada. Portanto, você não pode modelar isso como uma distribuição binomial . Em vez disso, você deve seguir uma árvore de probabilidades. O estado '1' pode entrar em '1' e '1-0', o estado '1-0' pode entrar em '1' e '1-0-0', e o estado '1-0-0' é um estado absorvente. Você pode anotá-lo como:
e a probabilidade de atingir o padrão '1-0-0', depois de rolar o primeiro '1' (você começa com o estado '0', ainda não tendo virado a cabeça), dentro de flips é metade da probabilidade estar no estado '1-0' dentro de invertidos:n n - 1
onde é o ésimo número da Fibonnaci. A probabilidade não condicional é uma somaFEu Eu
Tempo de espera pela probabilidade de virar vezes '1-0-0'k
Isso você pode calcular por uma convolução.
você terá a probabilidade de observar 20 ou mais padrões '1-0-0' (com base na hipótese de que a moeda é justa)
Aqui está o código R para calculá-lo:
Computando para moedas injustas
Podemos generalizar o cálculo acima da probabilidade de observar padrões em , quando a probabilidade de '1 = cabeça' é e os movimentos são independentes.x n p
Agora usamos uma generalização dos números de Fibonacci:
as probabilidades são agora como:
e
Quando traçamos isso, você obtém:
Portanto, embora o valor-p seja pequeno para uma moeda justa 0,0003247, devemos observar que não é muito melhor (apenas uma única ordem) para diferentes moedas injustas. A razão de verossimilhança, ou fator de Bayes , é de cerca de 11 quando a hipótese nula ( ) é comparada com a hipótese alternativa . Isso significa que o odds ratio posterior é apenas dez vezes maior que o odds ratio anterior.p = 0,5 p = 0,33
Portanto, se você pensou antes do experimento que a moeda era improvável injusta, agora você ainda deve pensar que a moeda é improvável injusta.
Uma moeda com mas injustiça em relação a ocorrências '1-0-0'ph um e um ds= pt a i l s
Poder-se-ia testar muito mais facilmente a probabilidade de uma moeda justa contando o número de caras e coroas e usar uma distribuição binomial para modelar essas observações e testar se a observação é particular ou não.
No entanto, pode ser que a moeda esteja lançando, em média, um número igual de cara e coroa, mas não é justo em relação a certos padrões. Por exemplo, a moeda pode ter alguma correlação para lançamentos sucessivos de moedas (imagino algum mecanismo com cavidades no interior do metal da moeda que são preenchidas com areia que flui como uma ampulheta em direção à extremidade oposta do lançamento anterior da moeda, que está carregando a moeda cair mais provavelmente no mesmo lado que o lado anterior).
Seja a primeira moeda lançada com igual probabilidade de cara e coroa e os lançamentos subsequentes terão a probabilidade do mesmo lado que o lançamento anterior. Em seguida, uma simulação semelhante ao início deste post fornecerá as seguintes probabilidades para o número de vezes que o padrão '1-0-0' excede 20:p
Você pode ver que é possível aumentar ligeiramente a probabilidade de observar o padrão '1-0-0' (algo em torno de uma moeda que tem alguma correlação negativa), mas mais dramático é que se pode torná-lo muito menos provável que supere o padrão '1-0-0'. Para baixo, você obtém muitas vezes as caudas após as cabeças, a primeira parte '1-0' do padrão '1-0-0', mas você não recebe tantas vezes duas caudas seguidas no '0-0' parte do padrão. O oposto é verdadeiro para os altos valores de .p = 0,45 p p
Usando a matemática em estatística
O acima exposto está bem, mas não é uma resposta direta à pergunta
Para responder a essa pergunta, pode-se usar a matemática acima, mas deve-se primeiro descrever muito bem a situação, os objetivos, a definição de justiça, etc. a pergunta explícita.
Uma questão em aberto é por que e como estamos procurando o padrão '1-0-0'.
fonte