Re: ˋ(′o‵")ˊ 再次求救~~~~

看板BIOSTAT作者 (小 美 )時間19年前 (2005/04/25 16:44), 編輯推噓2(200)
留言2則, 2人參與, 最新討論串3/3 (看更多)
: 都是在做 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
文章代碼(AID): #12RAtm4b (BIOSTAT)
文章代碼(AID): #12RAtm4b (BIOSTAT)