我想使用SAS进行对数回归中的标准似然比检验。我将有一个完整的逻辑模型,包含所有变量,名为A和嵌套逻辑模型B,它是从A中去掉一个变量得到的。
如果我想测试退出变量是否重要,我将执行模型A和B似然比测试。有没有一种简单的方法在SAS中使用程序执行此测试(本质上是卡方测试)?非常感谢你的帮助。
发布于 2011-10-20 09:56:28
如果您想执行完全模型与单变量丢弃模型的似然比测试,您可以使用带有type3选项的GENMOD过程。
脚本:
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
发布于 2011-10-13 00:05:31
我不是逻辑回归方面的专家,但我认为您正在尝试实现的目标可以使用PROC LOGISTIC完成,即在模型语句中使用"SELECTION=SCORE“选项。还有其他选择选项可用,例如逐步选择,但我认为得分匹配最接近您正在寻找的。我建议你好好读一读,因为还有一些相关的选项(BEST=,START= STOP=),你可能也会从中受益。
发布于 2015-11-23 17:28:05
我不确定是否有一条PROC语句可以专门执行LRT,但您可以计算嵌套模型的测试。
脚本
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;
https://stackoverflow.com/questions/7734495
复制相似问题