我有一个数据集,其响应变量是聚簇序数数据(顺序为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中进行编码吗?