大神们!帮帮我吧!
我的问题是这样,我要对多重应答题(多选题)进行分析,多选题中有6个选项,用GEE拟合了六个常数项,分别代表六个选项,因为六个选项选中的平均概率都不同。然后分别对各自变量做哑变量,结果如下。
女性生殖权利影响因素分析 SAS程序
因变量为“您认为哪些属于女性生殖健康权利?”1. 信息权, 2. 服务可及性权利,3知情选择权,4. 安全权,5. 隐私权,6.尊严权。
自变量有:
l 地区(a),分为4类:2. 内蒙古,6. 重庆,7.广州,8.河南;
l 年龄(age),分为4类:1. <25,2. 25-34,3. 35-44,4. 45-49;
l 户籍(hr),分为两类:1. 城市,2. 农村;
l 您是否做过性别鉴定(fst),分为两类:1. 是,2. 否;
l 生育孩子数 (g2aa),分为三类:0, 1, 2
l 您对生殖健康健康服务满意度如何(sat),分为三类:1. 很满意, 2. 满意,3. 不满意
根据文献拟合的模型:
Logit(E(Pij)=β01z1+β02z2+β03z3+β04z4+β05z5+β06z6+β1az2_8+β2az6_8+β3az7_8+β4agez1_4+β5agez2_4+β6agez3_4+β7hrz1_2+β8fstz1_2+β9g2aaz0_2+β10g2aaz1_2+β11satz1_3+β12satz2_3
sas代码如下:
proc genmod data=jixianguding_3 desc;class id a;
model outcome=int1 int2 int3 int4 int5 int6 az2_8 az6_8 az7_8 fdz1_4 fdz2_4 fdz3_4 a3z1_2 a4az1_4 a4az2_4 a4az3_4 a5az1_4 a5az2_4 a5az3_4 a7z1_5 a7z2_5
a7z3_5 a7z4_5 b2aaz0_2 b2aaz1_2 b14z1_2 b15z1_2 c4az1_2 c14z1_3 c14z2_3 int1*b2aaz0_2 int2*b2aaz0_2 int3*b2aaz0_2 int4*b2aaz0_2 int5*b2aaz0_2 int6*b2aaz0_2/d=binomial link=logit noint;
repeated subject=id(a)/TYPE=un corrw;run;
分析结果中发现最后一项“int6*b2aaz0_2”估计值全为0,考虑是因为以该项作为基准来估计其他五项。但是,我也需要“int6*b2aaz0_2”的结果啊!这怎么办?!