作者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