我对 R 编程非常陌生。我已经在 SAS 中实现了这个程序,以强制在逻辑回归模型中包含强制性变量。但我无法编写程序。下面是我用SAS编写的程序。当在程序中包含选项 include=1
时,它将自动包含模型中所需的变量。
我怎样才能在 R 中做到这一点,请告诉我。提前致谢
proc logistic data=Model_Data desc plots(only)=roc ;
model Target(event="1") = prop_sex Education Age / ctable pprob=0.51 pevent=0.40 include=1 ;
output out=out3 p=pred1;
run;
数据:
dput(Model_data)
structure(list(Age = c(68, 71, 84, 82, 72, 69, 62, 78, 64, 73,
77, 70, 96, 66, 72, 70, 78, 73, 65, 76, 72, 74, 82, 62, 78, 85,
77, 62, 75, 85, 85, 76, 63, 86, 75, 83, 75, 67, 78, 73, 69, 84,
63, 68, 64, 71, 64, 62, 71, 82, 72, 85, 72), prop_sex = c(1,
1, 2, 1, 2, 2, 1, 2, 2, 1, 2, 2, 2, 1, 2, 1, 2, 1, 1, 1, 2, 1,
1, 2, 2, 1, 2, 2, 2, 2, 1, 2, 2, 2, 2, 1, 1, 2, 2, 2, 2, 2, 2,
2, 2, 1, 1, 2, 1, 2, 2, 1, 2), Education = c(1, 2, 2, 1, 2, 1,
3, 3, 1, 3, 1, 2, 1, 2, 2, 999, 1, 999, 1, 999, 1, 1, 2, 1, 2,
2, 2, 1, 2, 1, 3, 1, 1, 2, 1, 1, 1, 1, 1, 2, 1, 3, 1, 1, 1, 1,
1, 1, 1, 1, 1, 2, 2), Target = c(1, 0, 1, 0, 0, 0, 1, 1, 1, 0,
0, 1, 0, 1, 0, 1, 1, 0, 0, 0, 0, 0, 1, 1, 1, 1, 1, 0, 0, 1, 0,
1, 0, 1, 1, 0, 1, 0, 1, 1, 1, 0, 1, 0, 0, 0, 0, 0, 0, 1, 0, 0,
0), cr_ci_1 = c(4, 4, 4, 4, 4, 5, 4, 4, 1, 1, 5, 3, 3, 1, 2,
5, 3, 5, 2, 4, 3, 4, 2, 1, 4, 3, 2, 2, 2, 5, 6, 4, 3, 4, 3, 5,
3, 2, 2, 2, 3, 4, 5, 4, 4, 5, 3, 4, 5, 3, 4, 1, 3), cr_ci_2 = c(3,
3, 3, 4, 3, 3, 3, 1, 4, 3, 3, 5, 3, 4, 2, 4, 2, 3, 5, 2, 2, 2,
2, 1, 1, 3, 4, 3, 2, 1, 4, 4, 3, 2, 1, 5, 1, 4, 1, 3, 2, 5, 4,
3, 4, 2, 3, 3, 4, 1, 3, 1, 1), cr_ci_3 = c(2, 2, 2, 2, 2, 2,
2, 1, 1, 1, 1, 1, 1, 1, 2, 1, 2, 1, 1, 2, 1, 2, 1, 1, 1, 1, 1,
1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 2, 1,
3, 1, 1, 1, 1), cr_ci_4 = c(2, 2, 3, 2, 2, 2, 2, 4, 5, 1, 1,
3, 1, 3, 2, 1, 4, 1, 2, 1, 1, 2, 1, 4, 4, 1, 1, 3, 1, 1, 3, 1,
4, 1, 4, 6, 1, 1, 3, 1, 1, 1, 1, 2, 1, 1, 2, 4, 1, 3, 3, 3, 1
), cr_ci_5 = c(4, 2, 3, 4, 4, 3, 2, 4, 1, 1, 1, 1, 1, 5, 2, 1,
2, 1, 1, 2, 1, 2, 2, 1, 1, 3, 4, 2, 3, 3, 1, 1, 1, 1, 1, 5, 1,
1, 2, 1, 1, 3, 1, 2, 3, 4, 2, 1, 4, 1, 1, 1, 1), cr_ci_6 = c(2,
4, 4, 2, 2, 2, 2, 1, 4, 1, 1, 1, 3, 1, 2, 1, 2, 1, 1, 1, 1, 2,
4, 1, 1, 1, 3, 3, 1, 1, 3, 1, 1, 1, 1, 5, 1, 1, 1, 1, 4, 3, 4,
3, 1, 4, 1, 3, 1, 1, 2, 1, 4), cr_ci_7 = c(2, 2, 4, 2, 2, 2,
2, 1, 1, 1, 1, 1, 1, 1, 2, 1, 2, 1, 1, 1, 1, 2, 1, 1, 1, 1, 1,
1, 1, 1, 1, 1, 1, 1, 1, 3, 1, 1, 1, 1, 2, 1, 1, 2, 1, 4, 2, 1,
1, 1, 2, 1, 1), cr_ci_8 = c(5, 3, 2, 3, 3, 3, 3, 1, 4, 2, 1,
1, 1, 1, 2, 3, 2, 2, 1, 4, 1, 2, 4, 1, 1, 1, 3, 1, 1, 1, 1, 1,
1, 1, 1, 1, 1, 2, 1, 1, 1, 3, 1, 2, 1, 1, 1, 1, 1, 1, 1, 1, 1
)), row.names = c(NA, -53L), class = c("tbl_df", "tbl", "data.frame"
))