[程式] genmod 多項式

看板Statistics作者時間12年前 (2012/07/04 16:04), 編輯推噓2(209)
留言11則, 2人參與, 最新討論串1/1
[軟體程式類別]: 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

07/04 16:46, , 1F
class id; <---改成這樣看看
07/04 16:46, 1F

07/05 11:21, , 2F
link=clogit 應該是執行cumulative logit model,所以
07/05 11:21, 2F

07/05 11:22, , 3F
beta應該只會有一個,而intercept會有兩個,此時您是將
07/05 11:22, 3F

07/05 11:24, , 4F
您的依變項當作ordinal的。
07/05 11:24, 4F

07/05 11:25, , 5F
若您的依變項不是Ordinal,則應該用baseline-category
07/05 11:25, 5F

07/05 11:26, , 6F
logit model,在proc logistic中link=glogit,可做這分析
07/05 11:26, 6F

07/05 11:27, , 7F
但是proc genmod中似乎沒有glogit這項功能,可能要另外
07/05 11:27, 7F

07/05 11:28, , 8F
找其他軟體做了
07/05 11:28, 8F

07/05 11:29, , 9F
最後,可參考作者為Agresti書名An introduction to
07/05 11:29, 9F

07/05 11:30, , 10F
categorical data analysis第六章Muticategory Logit
07/05 11:30, 10F

07/05 11:31, , 11F
Model裡面有一些解釋與範例
07/05 11:31, 11F
文章代碼(AID): #1Fy_aIQE (Statistics)