r - ロジスティック回帰におけるカテゴリ予測変数の参照グループの変更

Question

レベル0〜6のカテゴリ予測変数を使用してロジスティック回帰を実行しています。デフォルトでは、Rはレベル0を参照グループと見なします。

レベルの名前を変更せずに、たとえばレベル0の代わりにレベル3を参照として使用するようにRに指示するにはどうすればよいですか？

Answer 1

この関数を使用Cして、データフレームのコントラストを定義します。

データフレームがDFであり、因子変数がfctである場合、

DF$fct <- C(DF$fct, contr.treatment, base=3)

（未テスト）。

Answer 2

これは、Frank Harrellのパッケージ（多くの追加機能を可能にする）を使用すると非常に簡単です。
たとえば、「df」という名前のデータフレームの場合

library(Hmisc) 
library(rms)
dd=datadist(df)
options(datadist='dd')
(m=lrm(y ~ catvar, data=df)) #uses the mode as the reference group
summary(m, catvar=3) #using level=3 as the reference group
summary(m, catvar=0) #using level=0 as the reference group

（そして、いくつかのサンプルデータは良かったでしょう-あなたはいつでもdput小さなデータセットを出力するために使うことができます）

Answer 3

このC関数はすでに提案されています。また、、、、なども見contrastsてrelevelくださいreorder。

Answer 4

再レベル化関数を使用できます：dataframe $ x1 <-relevel（datafrmae $ x1、 "ここに参照カテゴリを入力"）