从PROC GENMOD运行PROC NLMIXED

时间:2019-03-21 18:15:27

标签: sas cluster-analysis ordinal

我有一个数据集,其响应变量是聚簇序数数据(顺序为1、2、3)。

我使用SAS的以下代码拟合模型

PROC GENMOD DATA = data RORDER=data DESCENDING;
CLASS group (REF=‘1’) id/;
MODEL group = predictor/ dist=mult link=clogit;
REPEATED subject=id;
RUN;

我需要使用SAS NLMIXED运行我的模型,以验证我从PROC GENMOD获得的估计。这是我现在拥有的代码。

PROC NLMIXED DATA=DATA;
eta = i1 + predictor + u;
p1 = exp(eta1)/(1+exp(eta1))
ll = group*log(p1)
Model group ~ general(ll);
Random u ~ normal(0, sigma*sigma) subject = id;
Run;

尽管此模型运行,但这似乎非常错误。您能告诉我如何在SAS NLMIXED中进行编码吗?

0 个答案:

没有答案