问题已解决，应该是对的吧...自己来回复一下
后来是参考何宗武老师‘追踪资料分析_原理与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的问题，请问
: 是否有合适的套件呢？谢谢！