Python - Répéter une expérience de retournement de pièces biaisée

Question

J'essaie d'écrire une fonction qui simule X nombre d'expériences de retournement de pièces biaisées (H = 0.6, T = 0.4), qui consiste en N retournements de pièces, pour répondre à la question " Quelle est la valeur attendue du nombre de groupes après les retournements de pièces N? " Le groupe, par définition, est la séquence maximale des mêmes valeurs dans un ordre consécutif. Par exemple: ['H', 'H', 'H', 'H', 'H', 'H', 'T', 'H', 'H', 'H'] a 3 groupes et ['T', 'H', 'H', 'H', 'T', 'T', 'H', 'H', 'H', 'H'] a 4 groupes.

# Libraries 
import random 
from itertools import groupby 
from itertools import chain 

# Function for biased coin 
def flip(p): 
    return 'H' if random.random() < p else 'T' 

# Number of coin flips 
N = 10 
flips = [flip(0.6) for i in range(N)] 
print (len(list(groupby(flips)))) 

# Function to simulate X iterations of N coin flips 
def simulate(X, N): 
    Outcome = [] # Empty list to store each experiment's result 
    # X Number of coin simulations 
    for i in range(X): 
     # Coin simulation of N flips 
     flips = [flip(0.6) for j in range(N)] 
     # Append experiment's result into Outcome list 
     Outcome.append(len(list(groupby(flips)))) 
    # Expected Value of the number of groups 
    sum(Outcome)/X 

Une idée pourquoi cela ne fonctionne pas? J'ai l'erreur suivante TypeError: unsupported operand type(s) for /: 'list' and 'int'

Answer 1

Vous divisez une liste par un nombre. Et vous devez utiliser return dans votre fonction.

def simulate(X, N): 

    Outcome = [] # Empty list to store each experiment's result 
    # X Number of coin simulations 
    for i in range(X): 
     # Coin simulation of N flips 
     flips = [flip(0.6) for j in range(N)] 
     # Append experiment's result into Outcome list 
     Outcome.append(len(list(groupby(flips)))) 
    X = float(X) 
    # Expected Value of the number of groups 
    return sum(Outcome)/X 

Answer 2

# Libraries 
import random 
from itertools import groupby 
from itertools import chain 

# Function for biased coin 
def flip(p): 
    return 'H' if random.random() < p else 'T' 

# Number of coin flips 
N = 10 
flips = [flip(0.6) for i in range(N)] 
print (len(list(groupby(flips)))) 

# Function to simulate X iterations of N coin flips 
def simulate(X, N): 
    Outcome = [] 
    for i in range(X): 
     flips = [flip(0.6) for j in range(N)] 
     Outcome.append(len(list(groupby(flips)))) 
    print(sum(Outcome)/(X)) 

Answer 3

Vous essayez de prendre la moyenne des résultats, donc au lieu de les stocker dans une liste, vous voulez les additionner et diviser par le nombre d'essais.

def simulate(x, n): 
    outcome = 0.0 
    for i in range(x): 
     flips = [flip(0.6) for j in range(n)] 
     outcome += len(list(groupby(flips))) 
    return outcome/x