作者smallmei (小 美 )
看板BIOSTAT
標題Re: ˋ(′o‵")ˊ 再次求救~~~~
時間Mon Apr 25 17:00:34 2005
: 都是在做 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
1F:推 FrancisWang: 難怪一個碩士要唸... 140.112.122.14 04/25
2F:推 bugu:好有毅力啊~給你拍拍手... 140.112.5.71 04/25