各位大牛：本人工作中遇到一模型，具体内容见附件文献中所示，有数据、有例子，已经有用R和STATA软件实现的，但我一直希望能用一套SAS宏搞定，无奈小弟才疏学浅，特向各位大牛有偿请教，经费都好商量。
希望：通过宏程序，设置不同的参数，控制模型，能输出所有参数估计结果和预测值，并能教会小弟如何使用。
模型简介：
被测（y）为个体主观报告值，为等级变量（1,2,3,4,5），其报告值由潜在变量（y*）和一套切点（t0，t1，。。。t5）值决定。y*，y以及切点值之间有如下关系：
y*=xβ+e， e~N(0,1)，x为个体特征向量，β为带估计系数；
y=1，if t0<y*<=t1；
y=2，if t1<y*<=t2；
y=3，if t2<y*<=t3；
y=4，if t3<y*<=t4；
y=5，if t4<y*<t5；
由于每个个体的切点不一样，因此有：
t1=xγ；
t(k)=t(k-1)+exp(xγ)；，k=2,3,4,5且 -∞=t0<t1<t2<t3<t4<t5=+∞
为了估计以上参数，引入若干情景，即让同一个体对虚拟个体的同一被测问题进行报告，因此，对于每一个虚拟个体，有（z和z*）；
zj*=θj+ej， ej~N(0,σj2)，θj为第j虚拟个体的真实被测，z为个体依据同一套切点（可理解为个人评价标准）对该虚拟个体的报告值，z和z*以及切点之间有如下关系；
z=1，if t0<z*<=t1；
z=2，if t1<z*<=t2；
z=3，if t2<z*<=t3；
z=4，if t3<z*<=t4；
z=5，if t4<z*<t5；
通过构建联合似然函数对以上模型求解，获得关心的参数（y*,β,t1,t2,t3,t4,t5,γ等）
模型大致如此~~



小弟联系方式：QQ：10217081，邮箱：oliverwist@163.com
期待大牛的出现！