Re: [問題] Ordered Probit for panel data

看板R_Language作者 (阿速)時間7年前 (2016/12/23 15:18), 編輯推噓0(000)
留言0則, 0人參與, 最新討論串2/2 (看更多)
問題已解決,應該是對的吧...自己來回覆一下 後來是參考何宗武老師『追蹤資料分析_原理與R程式實務』書裡面的 附錄一,採用pglm套件裡的pglm函數對panel data跑Ordered Probit model, 原本用rms套件的lrm函數是錯的(未考慮資料為panel的形式),下面是我使用 的方式。 pglm(Y ~ X, mydata, family = ordinal("probit"), R = 5, print.level = 3, method = 'bfgs', index = 'country', model = "random") 其中, Y是被解釋變數,依照原資料切成五個百分位數來分成5個rank,也就是說 1 < 2 < 3 < 4 < 5 X就是我的解釋變數,mydata就是放我的資料,index = 'country'是因為 我的觀察值是各個國家2007-2015年的資料,國家那一欄名稱即"country" 比較特別的是R = 5,據我所知是跟數值計算需要給個初始值有關,查了 help文件預設是20,但會run不出來,用5是因為文件裡的範例是用5,這 部分我不懂數值計算,只能等待有識之士來詳解了 print.level則是拿掉也能跑,純粹照抄文件範例 第一次在本版po文 & 不太會用ptt,若排版不好閱讀還請見諒,謝謝大家! ※ 引述《v303ofsky (阿速)》之銘言: : 我碩士論文的資料是Panel data,然後需要跑Ordered Probit model, : 我原本是用rms套件的lrm,但是這樣似乎沒有考慮到panel的問題,請問 : 是否有合適的套件呢?謝謝! -- ※ 發信站: 批踢踢實業坊(ptt.cc), 來自: 122.117.94.232 ※ 文章網址: https://www.ptt.cc/bbs/R_Language/M.1482477492.A.B13.html
文章代碼(AID): #1ONC-qiJ (R_Language)
文章代碼(AID): #1ONC-qiJ (R_Language)