LOGISTICプロシジャの結果でパラメータ推定値の符号が逆転する
[OS] ALL
[リリース] SAS 8.2以降
[キーワード] LOGISTIC, DESCENDING, REF=, ORDER=
[質問]パラメータ推定値の符号が正負逆となるのですが、どうしてでしょうか。また、応答変数のどの水準に対して、モデルが推定されているのでしょうか。
[回答]応答変数において、目的の水準に対するモデルではないことが考えられます。また、水準が3つ以上である場合には、水準の順序が意図されているとおりではないかもしれません。LOGISTICプロシジャを実行すると、ログにおけるNOTEメッセージ、およびアウトプットにおける‘Response Profile’にて、2値応答変数(binary response)の場合には、どの水準に対するモデルであるかが表示されます。また、順序応答変数(ordinal response)の場合には水準の順序、名義応答変数(nominal response)の場合には参照の水準が表示されます。2値応答変数、または順序応答変数の場合、モデルにおける水準、または水準の順序の変更は、パラメータ推定値の符号に影響します。 SAS 8.2以降では、MODELステートメントにおける応答変数のオプションとして、応答変数の水準、水準の順序、参照水準を指定することができます。これらのオプションを用いることによって、LOGISTICプロシジャにて意図されているモデルを推定することができます。
2値応答変数に対するモデル
EVENT=オプションを用いて、目的の応答水準を指定することができます。 proc logistic; model y(event='1') = <your model effects>; run;
順序応答変数に対するモデル
算出される結果は、応答変数の水準が一様に昇順であるか、降順であるかに影響されます。
たとえば、応答変数の水準が‘lo’,‘med’,‘hi’であるとします。
入力データセットにおいて、水準が昇順、または降順となっている場合、(たとえば、‘lo’が‘med’の前、‘med’が‘hi’の前となっている場合)、ORDER=DATAオプションを指定することによって、適切な並び順とすることができます。 proc logistic; model y(order=data) = <your model effects>; run;
この場合、LOGISTICプロシジャでは、‘Ordered Value’が最小である‘lo’に対するモデルを推定します。 proc logistic; model y(order=data descending) = <your model effects>; run;
名義応答変数に対するモデル
MODELステートメントにて、LINK=GLOGITオプションを指定することにより、一般化ロジットモデル(generalized logit model)を推定することができます。REF=オプションにて常に参照水準を指定することをお勧めします。 proc logistic; model y(ref='1') = <your model effects>; run; 上記の指定によって、以下の一般化ロジットに対するモデルを推定することになります。
参照水準1は、両方のロジットにおける分母となります。 上記の内容は、米国SAS Institute社のWebページに記載されているFAQをもとに作成しています。 |