/*二因素回归*/
data ex;
input y x1 x2 x3 @@;
cards;
1 0 20 1 1 0 23 1
1 0 32 1 1 0 38 1
1 1 25 1 1 0 20 0
1 0 24 0 1 0 28 0
1 0 30 0 1 0 32 0
1 0 38 0 1 1 26 0
1 1 29 0 1 1 34 0
1 1 33 1 1 1 38 0
0 0 26 1 0 0 29 1
1 1 40 0 0 0 22 1
0 0 34 1 0 0 30 1
0 0 38 1 0 0 37 1
0 1 24 1 0 1 25 1
0 1 29 1 0 1 32 1
0 1 34 1 0 1 37 1
0 1 40 1 0 1 40 1
0 0 33 0 0 0 36 0
0 1 24 0 0 1 34 0
0 1 32 0 0 1 36 0
0 1 38 0 0 0 39 0
;
proc logistic data=ex;
model y= x2/link=probit
sls=0.05
scale=none aggregate stb;/*R检验——模型拟合优度*/
run;
proc logistic data=ex;
model y= x1 x3 /scale=none aggregate;
run;
--------------------dong
/*二因素回归*/
data ex;
input y x1 x2 x3 @@;
cards;
1 0 20 1 1 0 23 1
1 0 32 1 1 0 38 1
1 1 25 1 1 0 20 0
1 0 24 0 1 0 28 0
1 0 30 0 1 0 32 0
1 0 38 0 1 1 26 0
1 1 29 0 1 1 34 0
1 1 33 1 1 1 38 0
0 0 26 1 0 0 29 1
1 1 40 0 0 0 22 1
0 0 34 1 0 0 30 1
0 0 38 1 0 0 37 1
0 1 24 1 0 1 25 1
0 1 29 1 0 1 32 1
0 1 34 1 0 1 37 1
0 1 40 1 0 1 40 1
0 0 33 0 0 0 36 0
0 1 24 0 0 1 34 0
0 1 32 0 0 1 36 0
0 1 38 0 0 0 39 0
;
proc logistic des;
model y=x1 x2 x3/scale=none aggregate;
run;
proc logistic descending;
freq count;
model effect=sex degrees/scale=none aggregate;
run;
