Wie durchführen, likelihood-ratio-test Logistische regression in SAS?

Ich möchte die standard-likelihood-ratio-test in logsitic regression mit SAS. Ich habe eine vollständige Logistische Modell mit allen Variablen namens A und einen verschachtelten Logistische Modell B, das abgeleitet ist von Drop-out-eine variable aus A.

Wenn ich will, um zu testen, ob das "drop-out" - variable signifikant ist oder nicht, werde ich führen Sie einen likelihood-ratio-test von Modell A und B. gibt es einen einfachen Weg, um diesen test durchzuführen (im wesentlichen ein chi-Quadrat-test) in SAS mit einem PROC? Vielen Dank für die Hilfe.

InformationsquelleAutor Steve | 2011-10-12
  1. 3

    Wenn Sie ausführen möchten likelihood-ratio-tests, die vollständige Modell v. s. eine variable gelöscht-Modell, können Sie mit der GENMOD-Prozedur mit den type3-option.

    Skript:

    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;

    Ausgabe:

         LR Statistics For Type 3 Analysis

                          Chi-
Source           DF     Square    Pr > ChiSq

y                 1      16.09        <.0001
z                 2       0.00        1.0000
    InformationsquelleAutor Triad sou.
  2. 1

    Ich bin kein Experte auf logistischer regression, aber ich denke, was Sie versuchen zu erreichen, kann man mit PROC LOGISTIC, mit der "AUSWAHL=SCORE" - option auf der MODELL-Anweisung. Es gibt auch andere AUSWAHL-Optionen zur Verfügung, wie SCHRITTWEISE, aber ich denke PUNKTZAHL entspricht am nächsten das, was du suchst. Ich würde vorschlagen, Lesen obwohl, da gibt es einige, die damit verbundenen Optionen (BESTEN= START= STOP=), die Sie nutzen könnten, auch.

    • Die AUSWAHL=PARTITUR ist eigentlich die beste Teilmenge automatische feature-Auswahl-config. Ist es nicht genau das, was ich Suche. Ich möchte wissen, ob es gibt einige in SAS PROC, die helfen können, mich führen likelihood-ratio-test der beiden bestehenden Modelle (eine verschachtelte anderen). Aber ich danke Ihnen sehr für den Beitrag. Es zu schätzen wissen.
    InformationsquelleAutor RWill
  3. 0

    Ich bin mir nicht sicher, über eine PROC-Anweisung kann explizit ausführen, LRT, aber Sie können berechnen Sie den test für genestete Modelle.

    Skript

    proc logistic 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 logistic 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;
    • Warum erstellen Sie die makro-Variablen? Sollten Sie kombinieren Sie einfach die beiden Datensätze direkt.
    • Ohne die makro-Variablen, wie identifizieren Sie den vollständigen und reduzierten Modell, um zu berechnen Ihre LR.
    • Verwenden Sie die RENAME= - dataset option.
    • Ich glaube, ich fand einen Tippfehler. Sollte es nicht proc logistic (nicht ligistic)?
    • Vielen Dank für den Hinweis auf den Tippfehler. Ich bearbeitet die post nun.
    InformationsquelleAutor Avinash Kumar
Schreibe einen Kommentar