作者shga (我爱慢慢跑)
看板Statistics
标题[问题] Logistic regression的adjust
时间Tue Sep 15 10:48:48 2015
如果是跟统计软体有关请重发文章,使用程式做为分类。
请详述问题内容,以利板友帮忙解答,过短文章依板规处置,请注意。
为避免版面混乱,请勿手动置底问题,擅用E做档案编辑
目前进行仿照某研究进行安眠焦虑药物与癌症的分析,使用case-control方式。
case(得癌)match 性别 ,年龄 ,就医日期(年月) 1:4 control。
请问logistic regression,照理说unadjust variable的pvalue
应该都会比adjust来的小。
但目前跑出结果为unadjust 的时候,该药物是不显着。
但adjust其他变数後,结果variable变显着。
确认一整夜,都无法得到了解为什麽,因此想要请问大家阿。
例如下列为unadjust
BZD 1 与 0 之间的关系 0.972 (0.944-1.002) p值是0.0687
下列是adjust的结果
BZD 1 与 0 之间的关系 0.895 (0.867-0.924) p值有显着??
下述是adjust 的sas语法提供参考。主要是调整CCI,sex,age,cci_score。
proc logistic data= cci.match_1;
class BZD(REF='0') / PARAM=REF;
strata match_id;
MODEL index(event='1')= BZD
MI CHF VASCUL CVD DEMENTIA COPD RHEUM ULCERS LIVER DIABET PARALYS RENAL
sex age cci_score ;
run;
下述是non adjust
proc logistic data= cci.match_1;
class BZD(REF='0') / PARAM=REF;
MODEL index(event='1')= BZD ;
run;
--
※ 发信站: 批踢踢实业坊(ptt.cc), 来自: 114.24.73.66
※ 文章网址: https://webptt.com/cn.aspx?n=bbs/Statistics/M.1442285331.A.E6D.html
1F:→ andrew43: 你前面说「无共变数调整的p值通常较小」并没有这回事。 09/15 11:10
2F:→ shga: 感谢楼上说明,另外还有不解的是。 09/15 11:31
3F:→ shga: 有遇过针对某变数unadjust没有显着,但adjust变成显着? 09/15 11:31
4F:→ andrew43: 当共变数有影响力且在组间不相近时不就都会这样? 09/15 19:51
5F:→ recorriendo: 试着思考:adjust和unadjust下得到系数代表的意义是? 09/15 19:51
6F:→ yhliu: 不调整相关共变数时不显着, 调整後显着, 这应不罕见. 09/16 00:36
7F:→ yhliu: 模型 Y = f(X,error) 时 X 不显着而 Y = f(X,Z,error) 模型 09/16 00:37
8F:→ yhliu: 下 X 效应显着, 表示 X 的效应在不校正 Z 之效应时被 Z 所 09/16 00:38
9F:→ yhliu: 抑制. 如果 X,Z 对 Y 的效应是同向的, 但 X, Z 有负相关; 09/16 00:38
10F:→ yhliu: 或 X, Z 对 Y 之效应反向, 而 x, z 是正相关, 那麽不校正 Z 09/16 00:39
11F:→ yhliu: 之效应, 则 X 之效应被抑制, 甚至有时候被扭曲, Simpson's 09/16 00:40
12F:→ yhliu: paradox 就是说这种现象: 因为模型未考虑 Z 而致 X 与 Y 的 09/16 00:41
13F:→ yhliu: 关联方向被扭曲. 09/16 00:42