У меня есть набор данных, где его переменной отклика являются кластерные порядковые данные (по порядку 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?