请选择 进入手机版 | 继续访问电脑版
楼主: moonstone
35277 28

[学习分享] 一般线性模型(GLM)的相关思考   [推广有奖]

讲师

74%

还不是VIP/贵宾

-

威望
0
论坛币
10435 个
通用积分
335.8255
学术水平
160 点
热心指数
169 点
信用等级
124 点
经验
262166 点
帖子
237
精华
1
在线时间
519 小时
注册时间
2007-4-27
最后登录
2024-2-1

moonstone 发表于 2013-11-19 21:12:31 |显示全部楼层 |坛友微信交流群

+2 论坛币
k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群
赵安豆老师微信:zhaoandou666

经管之家联合CDA

送您一个全额奖学金名额~ !

感谢您参与论坛问题回答

经管之家送您两个论坛币!

+2 论坛币
最近遇到一个问题,如果因变量为一个连续变量(如胰岛素水平),主要考察的变量为分组变量(如正常血糖组,前糖尿病组,糖尿病组三组),现在的目的是想看调整多种变量(包括多个连续性变量和分类变量)后,胰岛素水平是否一致。


一开始的思路想到的是采用GLM进行协方差分析来解决。

但是有觉得似乎不是很对,因为经典的协方差分析通常只考虑一个连续变量(协变量)和一个分组变量,同时协变量和分组变量只有不存在交互的时候(经典协方差分析的前提)才能使用协方差分析。
针对我目前的问题,如果想调整多个连续性变量和分类变量,这种方法能否再叫协方差分析?如果可以认为是协方差的思想,用不用检验协方差分析的前提(如协变量与分组变量之间的交互)?多个连续性变量和分类变量存在时,该前提应该怎么检验?

通过跟别人交流之后,有一句话非常受用:线性模型其实最重要的不在于用的方差分析还是协方差分析,而主要是检验残差是否符合线性的几个条件。

受到该启发后,认真再复习GLM的相关资料,得到更加重要的总结如下(来自高惠璇SAS/STAT软件使用手册,实际是SAS8.2的User's guide的中文版,但是目前SAS 9.2,9.3的User's guide关于GLM模型的介绍中已经删去了这么经典的总结,实在可惜,倒让人看不到GLM的真正长处了):
如果X1-X3,Y1-Y2为连续性变量,Y3为分类变量,a-c为分类变量,time为时间变量,目前我们熟悉的模型可以简单概括如下:
(1) y1= x1    简单回归
(2) y1= x1 x2 x3  多重回归(multiple regression)
(3) y1 y2=x1 x2   多元回归(multivariate regression)
(4) y1= a   单因素方差分析
(5) y1= a b   (析因设计的)主效应分析
(6) y1= a b  a*b (析因设计的)主效应加交互项分析
(7) y1= a x1   协方差分析
(8) y3= a   单因素logistic回归
(9) y3= a b c x1 x2 x3   多因素logistic回归
(10) y3(time) =a   单因素cox回归
(11) y3(time) = a b c x1 x2 x3  多因素cox回归

1-7采用SAS的一般线性模型GLM都能实现,而1-9采用SAS的广义线性模型GENMOD都能实现,具体验证详见后面举例

再次回到开始的问题:掌握上述的基本思路后,因为因变量为连续变量,所以采用线性模型肯定是对的。如果因变量可以认为是正态的,那采用一般线性模型是合理的。所以现在的关键问题是:如果调整多个变量(包含分类和连续变量)后看不同分组间因变量(连续变量)是否仍有差异时,能否再称为协方差分析?我目前认为应该是可以的,但是事实上我们遇到这种情况后,并不再去强调它是协方差的思想,而只是回到线性模型分析的最初的起点,也即是检查残差是否符合线性的基本条件即可。

但是现实中,我们在使用GLM解决前面遇到的类似问题时,只是简单地用了,而很多时候我们都没有认真去检验残差是否符合这个条件,这可能是我们滥用GLM的表现之一,因为我们更多只关注模型的参数是否有意义,而不去关心对结果“无关紧要”的前提条件。

再次思考一个问题,上述列举的1-7模型,在GLM中并没有特定的选项指定是哪一种模型,而采用一种表达方式。由此,可以进一步深入概括一句话,GLM模型,对于上述列举的1-7模型并没有本质区别,唯一的区别只是模型中自变量的属性和数量不同。但是我们对1-7模型的叫法却不相同。而其原因是我们对事物的认识是一个由浅到深的过程,之前我们认为他们是不同的7件事情,随着认识的加深,发现原来这些问题可以用一个方式表达出来。而SAS的GENMOD则更能说明这一问题。现在还没有一个模型能把上述模型1-11用一种表达方式表示,但是COX回归在抛开基线生存函数之后剩下的部分也是线性模型,所以说不定哪天真的能够把上述所有模型用一种表达方式表示出来。到时候更应该相信人们对事物的认识绝对是一个由浅到深的过程啦。

附:GLM与GENMOD在协方差分析结果的比较
   data drugtest;
      input Drug $ PreTreatment PostTreatment @@;
      datalines;
   A 11  6   A  8  0   A  5  2   A 14  8   A 19 11
   A  6  4   A 10 13   A  6  1   A 11  8   A  3  0
   D  6  0   D  6  2   D  7  3   D  8  1   D 18 18
   D  8  4   D 19 14   D  8  9   D  5  1   D 15  9
   F 16 13   F 13 10   F 11 18   F  9  5   F 21 23
   F 16 12   F 12  5   F 12 16   F  7  1   F 12 20
   ;

   proc glm;
      class Drug;
      model PostTreatment = Drug PreTreatment / solution;
      lsmeans Drug / stderr pdiff cov out=adjmeans;
   run;

   proc genmod data=drugtest;
      class Drug;
      model PostTreatment = Drug PreTreatment / dist=NOR link=ID obstats type1;
   run;

   proc print data=adjmeans;
   run;




二维码

扫码加我 拉你入群

请注明:姓名-公司-职位

以便审核进群资格,未注明则拒绝

关键词:线性模型 Multivariate multivariat regression logistic回归 模型

已有 5 人评分经验 论坛币 学术水平 热心指数 信用等级 收起 理由
dxystata + 40 鼓励积极发帖讨论
李会超 + 100 精彩帖子
Tigflanker + 5 + 4 + 4 + 4 观点有启发
jingju11 + 5 + 5 + 4 鼓励积极发帖讨论
qqface123 + 1 + 1 + 1 观点有启发

总评分: 经验 + 140  论坛币 + 5  学术水平 + 10  热心指数 + 10  信用等级 + 9   查看全部评分

本帖被以下文库推荐

luisluan 发表于 2013-11-20 01:44:53 来自手机 |显示全部楼层 |坛友微信交流群
学习了,总结的真好,特别是那句关于前提的提醒

使用道具

nkunku 发表于 2013-11-20 08:49:34 |显示全部楼层 |坛友微信交流群

使用道具

jingju11 发表于 2013-11-22 11:29:17 |显示全部楼层 |坛友微信交流群
很有意思。京剧

使用道具

vincentwbb 发表于 2013-11-24 09:30:13 |显示全部楼层 |坛友微信交流群
学习了。

使用道具

是风啊 发表于 2013-11-24 12:49:26 |显示全部楼层 |坛友微信交流群
学习了,,“我们更多只关注模型的参数是否有意义,而不去关心对结果“无关紧要”的前提条件
”这句话是该提醒了很多人啊,个人感觉不要只知道应用,还要知道原理,所以说理论性的东西还是不可缺的,,谢谢楼主了哈,,

使用道具

思考的txt 发表于 2013-12-12 11:21:16 |显示全部楼层 |坛友微信交流群
好东西,收藏了,谢谢

使用道具

suzhzh 发表于 2013-12-14 12:28:02 |显示全部楼层 |坛友微信交流群

使用道具

anyachao 发表于 2014-1-18 16:18:57 |显示全部楼层 |坛友微信交流群
深入浅出,简练的表达,传达了复杂的内容。
貌似目前好多著者、师者都缺少这种思考与总结,有的书、课件就是靠东拼西凑来的。

使用道具

荒漠之驴 发表于 2014-4-26 11:42:47 |显示全部楼层 |坛友微信交流群
楼主,你好 关于你开始提到的那个问题“最近遇到一个问题,如果因变量为一个连续变量(如胰岛素水平),主要考察的变量为分组变量(如正常血糖组,前糖尿病组,糖尿病组三组),现在的目的是想看调整多种变量(包括多个连续性变量和分类变量)后,胰岛素水平是否一致。”glm的代码是怎样的呢 我最近也在探讨类似的问题,希望得到指点哈

使用道具

您需要登录后才可以回帖 登录 | 我要注册

本版微信群
加好友,备注cda
拉您进交流群

京ICP备16021002-2号 京B2-20170662号 京公网安备 11010802022788号 论坛法律顾问:王进律师 知识产权保护声明   免责及隐私声明

GMT+8, 2024-3-29 17:04