適合度檢定

適合度檢定的虛無假設，是針對每個比例給出一組具體的參數值。在你的分析中，等價的假設是：Benford 定律適用於市層級總得票數的首位數分佈。你可以寫成：

$$ H_0: p_1 = .30, p_2 = .18, \ldots, p_9 = .05 $$

其中 $p_1$ 是 Benford 長條圖中第一個長條的高度。對立假設是其中至少有一個比例不同；也就是首位數分佈「不」遵循 Benford 定律。

在這個練習中，你會用模擬來建立虛無分配，觀察在這些計數確實遵循 Benford 定律時，可能出現的卡方統計量分佈。

本練習屬於課程

R 中的類別資料推論

檢視課程

練習說明

把 p_benford 印到畫面上檢視。
以 iran 開始，使用 chisq_stat 計算卡方統計量。注意你必須在資料框中指定作為應變數的變數，以及要與之比較的機率向量。
在 point 虛無假設下，將比例向量 p 設為 p_benford，透過模擬產生 500 筆 Chisq 統計量樣本來建立虛無分配。將得到的統計量儲存為 null。

動手互動練習

試著完成這個範例程式碼，體驗一下這個練習。

# Inspect p_benford
p_benford

# Compute observed stat
chi_obs_stat <- ___
  chisq_stat(response = ___, p = ___)

# Form null distribution
null <- ___
  # Specify the response
  ___
  # Set up the null hypothesis
  ___
  # Generate 500 reps
  ___
  # Calculate statistics
  ___

編輯並執行程式碼