Rの出力を扱う (1)
lm の出力で与えられる p 値は、デフォルトでは両側検定の p 値です。双子の研究では、双子の IQ スコアが正の関連をもつという片側の科学的仮説に従うほうが、より妥当に思えるかもしれません。p 値は「観測されたデータ以上に極端な結果が得られる確率」なので、両側検定の p 値は片側検定の結果の 2 倍になります。つまり、R の両側の出力から片側の p 値を得るには、p 値を 2 で割ります。
Foster 対 Biological の線形回帰モデル model はスクリプト内に用意されています。
この演習はコースの一部です
R による線形回帰の推測
演習の手順
- モデルから Biological の係数を取得します。
- モデルを tidy します。
termが"Biological"の行でフィルタします。
- 片側検定の p 値を表す列
one_sided_p_valueを追加します。
実践的なインタラクティブ演習
このサンプルコードを完成させて、この演習に挑戦してみましょう。
model <- lm(Foster ~ Biological, data = twins)
# Get the Biological model coefficient
biological_term <- model %>%
# Tidy the model
___ %>%
# Filter for term equal to "Biological"
___
biological_term %>%
# Add a column of one-sided p-values
___(one_sided_p_value = ___)