1

次のコードを検討してください。

x <- read.table("http://data.princeton.edu/wws509/datasets/cuse.dat",
                header=TRUE)[,1:2]

fit <- glm(education ~ age, family="binomial", data=x)

summary(fit)

年齢には 4 つのレベルがあります: "<25" "25-29" "30-39" "40-49"

結果は次のとおりです。

ここに画像の説明を入力

したがって、デフォルトでは、レベルの 1 つが参照レベルとして使用されます。4 つのレベルすべて + インターセプトの glm 出力係数を取得する方法はありますか (つまり、基準レベルはありません)。SAS のようなソフトウェア パッケージはデフォルトでこれを行うので、これに対するオプションがあるかどうか疑問に思っていました。

ありがとう！

4

2 に答える 2

1

このmodel.matrix関数を使用して、age 変数の因子をバイナリ変数に変換します。

この回答を参照してください。

編集：ここに例があります：

x <- read.table("http://data.princeton.edu/wws509/datasets/cuse.dat",
                header=TRUE)[,1:2]
binary_variables <- model.matrix(~ x$age -1, x)
fit <- glm(x$education ~ binary_variables, family="binomial")
summary(fit)
于 2014-04-28T17:35:51.193 に答える