Re: ˋ(′o‵")ˊ 再次求救~~~~
: 都是在做 cumulative logit model
: model 應該是一樣的
: 只是response超過兩類 若加以covariates也超過兩類
: 不同procedure的預設表示法可能不同
: 也就是它們可能採用不一樣的 reference group 之類的
我來自我解答了 ^^a
謝謝學姊的提醒
經過了測試之後發現真的是reference group 設定的問題
以性別為例 gender=0,1
proc logistic的預設是
beta_gender0 + beta_gender1 = 0
( 即 gender0 = 1, beta_gender1 = -1 )
然後會列出 beta_gender0 的估計係數
proc genmod則是跟我們習慣使用的dummy設法一樣
令beta_gender1 = 0
( 即 gender0 = 1, beta_gender1 = 0 )
然後列出 beta_gender0 的估計係數
所以如果以這個例子來說的話
這兩個模型的beta_gender0剛好會是兩倍的關係
( 但是level>2的話 就不一定囉 )
☆☆ 讓 genmod 和 logistic 可以得到相同的參數估計值的方法 ☆☆
proc logistic ;
class gender / PARAM = REFERENCE ;
model Y = X ;
run;
----------------------------------------------------
真的笨哪~~~
一個指令搞兩天 Orz""
--
※ 發信站: 批踢踢實業坊(ptt.cc)
◆ From: 140.109.74.44
推
140.112.122.14 04/25, , 1F
140.112.122.14 04/25, 1F
推
140.112.5.71 04/25, , 2F
140.112.5.71 04/25, 2F
討論串 (同標題文章)
本文引述了以下文章的的內容:
完整討論串 (本文為第 3 之 3 篇):