如何编写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:
Premarital sex | Extramarital sex | total | |||
1 | 2 | 3 | 4 | ||
1 | 144 | 2 | 0 | 0 | 146 |
2 | 33 | 4 | 2 | 0 | 39 |
3 | 84 | 14 | 6 | 1 | 105 |
4 | 126 | 29 | 25 | 5 | 185 |
Total | 387 | 49 | 33 | 6 | 475 |
模型:
Where
说明:t=1表示对Premarital sex,t=2表示对Extramarital sex
正确结果应该是:
哪位SAS高手能写出代码?
如不能正确显示,请看附件。