始める無料で始める

binomial family 引数

lm() で線形回帰を実行する場合と、glm() でロジスティック回帰を実行する場合の大きな違いは、glm()family 引数を binomial に設定する必要があることです。binomial() は関数で、回帰の計算方法を glm() に伝える複数の関数を要素とするリストを返します。中でも重要なのが linkinvlinkfun で、これは実数全体(マイナス無限大から無限大)と確率(0 から 1)との間で変換するために使われます。

値のベクトル x と、確率のベクトル p が用意されています。

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

Rで学ぶ中級回帰分析

コースを見る

演習の手順

  • binomial() 関数の構造を確認します。要素として関数が 2 つ含まれていることに注目してください。 binomial()$linkinv binomial()$linkfun です。
  • x に対して binomial()$linkinv() を呼び出し、linkinv_x に代入します。
  • all.equal() を使って、linkinv_xxplogis() が同じ結果になることを確認します。
  • p に対して binomial()$linkfun() を呼び出し、linkfun_p に代入します。
  • all.equal() を使って、linkfun_ppqlogis() が同じ結果になることを確認します。

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

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

# Look at the structure of binomial() function
___

# Call the link inverse on x
linkinv_x <- ___

# Check linkinv_x and plogis() of x give same results 
___

# Call the link fun on p
linkfun_p <- ___

# Check linkfun_p and qlogis() of p give same results  
___
コードを編集して実行