ランダム化による推測とtによる推測の比較

技術的な条件（次の章を参照）が満たされると、ランダム化検定とt分布に基づく検定からの推測は同等の結論になります。手法が異なるため、まったく同じ数値にはなりませんが、p値や信頼区間は十分に近いはずです。

この演習はコースの一部です

R による線形回帰の推測

演習の手順

観測された傾き obs_slope の絶対値を計算します。
perm_slope に、各置換レプリケートでの傾きの絶対値の列を追加します。傾きの列名は stat です。
summarize() の呼び出しで、p値を「観測された傾きの絶対値と同程度以上に極端な傾きの絶対値」の割合として計算します。

実践的なインタラクティブ演習

このサンプルコードを完成させて、この演習に挑戦してみましょう。

# The slope in the observed data and each permutation replicate
obs_slope
perm_slope

# Calculate the absolute value of the observed slope
abs_obs_slope <- ___

# Find the p-value
perm_slope %>%
  # Add a column for the absolute value of stat
  ___ %>%
  summarize(
    # Calculate prop'n permuted values at least as extreme as observed
    p_value = ___
  )

コードを編集して実行