作者aerin5168 ()
看板Statistics
标题[程式] genmod 多项式
时间Wed Jul 4 16:04:32 2012
[软体程式类别]: SAS
[程式问题]: 用gee跑多项式logistic reg
[软体熟悉度]: 低(1~3个月)
[问题叙述]:
我的资料是重覆测量的,原本依变项是2类的类别变项,但现在改为3类,
所以改用multinomial,但还是要用gee去跑
原本的2类是proc genmod是用 dist=binomial link=logit type=un,去跑出来的,
现在改为3类的程式,改为 dist=multinomial link=clogit type=ind,
是可以跑出来,不过只有一条回归呈现,
但应该要有2条才对啊??
不知道是不有哪个地方写错了~
拜托大家了,已经试了很久找不到答案！
以下呈现依变项为3类类别的程式
[程式范例]:
DATA two; SET one;
PROC GENMOD DESCENDING;
weight count;
CLASS id age gender ;
MODEL new= age gender/dist=multinomial link=clogit ;
REPEATED SUBJECT=id / TYPE=ind ;
ESTIMATE "60-64岁" age_1 -1 1 0 0 0 /EXP;
ESTIMATE "65-69岁" age_1 -1 0 1 0 0 /EXP;
ESTIMATE "70-74岁" age_1 -1 0 0 1 0 /EXP;
ESTIMATE "75岁以上" age_1 -1 0 0 0 1 /EXP;
ESTIMATE "女性" gender -1 1 /EXP;
run;
-----------------------------------------------------------------------------
--
※ 发信站: 批踢踢实业坊(ptt.cc)
◆ From: 111.240.136.80
1F：推 tew:class id; <---改成这样看看 07/04 16:46
2F：推 sadricky:link=clogit 应该是执行cumulative logit model，所以 07/05 11:21
3F：→ sadricky:beta应该只会有一个，而intercept会有两个，此时您是将 07/05 11:22
4F：→ sadricky:您的依变项当作ordinal的。 07/05 11:24
5F：→ sadricky:若您的依变项不是Ordinal，则应该用baseline-category 07/05 11:25
6F：→ sadricky:logit model，在proc logistic中link=glogit，可做这分析 07/05 11:26
7F：→ sadricky:但是proc genmod中似乎没有glogit这项功能，可能要另外 07/05 11:27
8F：→ sadricky:找其他软体做了 07/05 11:28
9F：→ sadricky:最後，可参考作者为Agresti书名An introduction to 07/05 11:29
10F：→ sadricky:categorical data analysis第六章Muticategory Logit 07/05 11:30
11F：→ sadricky:Model里面有一些解释与范例 07/05 11:31