問題已解決，應該是對的吧...自己來回覆一下
後來是參考何宗武老師『追蹤資料分析_原理與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的問題，請問
: 是否有合適的套件呢？謝謝！