SASのロジスティック回帰に関する尤度比検定の実行方法は？

Question

私はSASを使用してLogistic回帰で標準尤度比検定を実行したいと思います。 Aという名前のすべての変数と、Aから1つの変数を削除して得られるネストロジスティックモデルBを含む完全なロジスティックモデルを持ちます。

ドロップアウト変数が有意であるかどうかをテストしたい場合モデルAとBの尤度比テストを実行します.SASではPROCを使用してこのテスト（本質的にカイ二乗検定）を実行する簡単な方法はありますか？助けてくれてありがとうございました。

Answer 1

フルモデルv.s.の尤度比テストを実行する場合は、 1つの変数をドロップしたモデルでは、GENMODプロシージャをtype3オプションとともに使用できます。

スクリプト：

data d1; 
do z = 0 to 2; 
do y = 0 to 1; 
do x = 0 to 1; 
    input n @@; 
    output; 
end; end; end; 
cards; 
100 200 300 400 
50 100 150 200 
50 100 150 200 
; 
proc genmod data = d1; 
    class y z; 
    freq n; 
    model x = y z/error = bin link = logit type3; 
run; 

出力：

 LR Statistics For Type 3 Analysis 

          Chi- 
Source   DF  Square Pr > ChiSq 

y     1  16.09  <.0001 
z     2  0.00  1.0000 

Answer 2

私はロジスティック回帰について専門家はいませんが、あなたが達成しようとしていることは、MODELステートメントの "SELECTION = SCORE"オプションを使ってPROC LOGISTICで行うことができると思います。 STEPWISEのような他のSELECTIONオプションもありますが、SCOREはあなたが探しているものに最も近いものと考えています。私はそれを読むことをお勧めします。あなたが恩恵を受けるかもしれないいくつかの関連するオプション（BEST =、START = STOP =）があるので。

Answer 3

私は、具体的LRTを実行することはできますが、ネストされたモデルのためのテストを計算することができPROCステートメントについてはよく分かりません。

スクリプト

proc ligistic data = full_model; 
model dependent_var = independent_var(s); 
ods output GlobalTests = GlobalTests_full; 
run; 

data _null_; 
set GlobalTests_full; 
if test = "Likelihood Ratio" then do; 
    call symput("ChiSq_full", ChiSq); 
    call symput("DF_full", DF); 
    end; 
run; 

proc ligistic data = reduced_model; 
model dependent_var = independent_var(s); 
ods output GlobalTests = GlobalTests_reduced; 
run; 

data _null_; 
set GlobalTests_reduced; 
if test = "Likelihood Ratio" then do; 
    call symput("ChiSq_reduced", ChiSq); 
    call symput("DF_reduced", DF); 
    end; 
run; 

data LRT_result; 
LR = &ChiSq_full - &ChiSq_reduced; 
DF = &DF_full - &DF_reduced; 
p = 1 - probchi(ChiSq,DF); 
run;