如何编写SAS代码估计下面这个logit模型的参数？

样本说明：

这是一项普通的社会调查，对于475个被访者，每人被访问两项，第一项是对于发生婚前性行为的态度，第二项是对于婚外性行为的态度（即婚后与配偶之外的第三者发生性行为的态度）。回答的态度以1、2、3、4四个水平给出，这四个水平有程度上的区别（ordinal）：1=always wrong，2=almost always wrong，3=wrong only sometimes，4=not wrong at all。

Table 10.5:

模型：

，t＝1，2， j＝1，…，I-1

Where and

说明：t=1表示对Premarital sex，t=2表示对Extramarital sex

表示对问题的答案

正确结果应该是： (SE=0.13), (df=2)

哪位SAS高手能写出代码？

如不能正确显示，请看附件。