数学的プログラミング：データサイエンス能力を進化させるための大事な習慣(2/3)

１．数学的プログラミング：データサイエンス能力を進化させるための大事な習慣(2/3)

・ランダムダーツ投げのような検証は多数の試行が行われるという条件下でのみ有効
・他のモンテカルロ実験と同様に投げた回数が多いほど近似は良くなる
・科学的実験は単なるコーディングを超えた、仮説の探求と検証がなくては終わらない

２．数学的プログラミングの実例

このアイディアはどういう原理なのですか？
このアイディアは非常に単純です。沢山のダーツを投げた場合、ダーツが円の内側に刺さる確率は、正方形の板の面積に対する円の面積の比率になります。

基本的な数学の助けを借りれば、あなたはこの比率がpi / 4である事を示すことができます。 だから、円周率を得るためには、あなたはその数を4倍するだけです。

ここで重要なのは、(円の内側に刺さるダーツの割合が)真の確率と等しくなるように、たくさんのダーツを投げることです。

このランダムなダーツ投げによる検証は、多数の試行が行われるという条件下でのみ有効で、これは大数の法則、または頻度主義者による確率の定義(frequentist definition of probability)から来ています。

Python code
Pythonコードを説明するJupyterノートブックは、私のGithubレポジトリ「tirthajyoti/Stats-Maths-with-Python」にあります。 コピーまたはフォークしてください。手順は簡単です。

まず、ダーツをランダムに投げた際に刺さった位置をシミュレートする関数を作成します。

def throw_dart():
    """
    Simulates the randon throw of a dart. It can land anywhere in the square (uniformly randomly)
    """
    # Center point
    x,y = 0,0
    # Side of the square
    a = 2
    
    # Random final landing position of the dart between -a/2 and +a/2 around the center point
    position_x = x+a/2*(-1+2*random.random())
    position_y = y+a/2*(-1+2*random.random())

次に、刺さった座標が与えられた際にダーツが円の内側に刺さったかどうかを判定する関数を作成します。

def is_within_circle(x,y):
    """
    Given the landing coordinate of a dart, determines if it fell inside the circle
    """
    # Side of the square
    a = 2
    
    distance_from_center = sqrt(x**2+y**2)
    
    if distance_from_center < a/2:
        return True
    else:
        return False

最後に、ダーツを大量に投げた結果をシミュレートする関数を作成し、累積結果からpiの値を計算します。

def compute_pi_throwing_dart(n_throws):
    """
    Computes pi by throwing a bunch of darts at the square
    """
    n_throws = n_throws
    count_inside_circle=0
    for i in range(n_throws):
        r1,r2=throw_dart()
        if is_within_circle(r1,r2):
            count_inside_circle+=1
            
result = 4*(count_inside_circle/n_throws)

しかし、プログラミングをここで終えてはいけません。このプログラムによる近似値がどれほど優れているか、そしてランダムに投げた回数によってどのように変化するかをテストする必要があります。

他のモンテカルロ実験と同様に、投げた回数が多いほど近似は良くなると予想されます。

これがデータサイエンスと分析の中核です。期待される出力を表示する関数を書いて終わるだけでは不十分です。本質的な部分のプログラミングは完了しているかもしれませんが、科学的実験は単なるコーディングを超えた、仮説の探求と検証がなくては終わりません。

さて、この実験を繰り返して何が起こるか見てみましょう。

多数のランダムな投擲を繰り返すことで、より良い近似値が得られる事がわかります。

n = 5000000
sum=0
for i in range(20):
    p=compute_pi_throwing_dart(n)
    sum+=p
    print("Experiment number {} done. Computed value: {}".format(i+1,p))
print("-"*75)
pi_computed = round(sum/20,4)
print("Average value from 20 experiments:",pi_computed)

３．数学的プログラミング：データサイエンス能力を進化させるための大事な習慣(2/3)まとめ

１）www.kdnuggets.com
Mathematical programming —  Key Habit to Build Up for Advancing Data Science

２）medium.com
Essential Math for Data Science

３）github.com
tirthajyoti/Stats-Maths-with-Python