SAS PROC NLMIXED PARMS - PullRequest
Новая
       25

SAS PROC NLMIXED PARMS

0 голосов
/ 23 марта 2019

Я пытаюсь подогнать нелинейные смешанные модели, используя PROC NLMIXED, чтобы увидеть, как оценки сравниваются между PROC GENMOD и PROC NLMIXED.Ниже показан мой код PROC GENMOD, в котором результат представляет собой коррелированные порядковые данные (1,2,3). 

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

Затем я запускаю нелинейные смешанные модели, используя следующий код PROC NLMIXED. 

PROC NLMIXED DATA=data.gee;
PARMS B0=0 B1=0 SD=1 THRES1=1 THRES2=1;
Z=B0+B1*pred+U;
IF (group=1) THEN P=1/(1+EXP(-(0-Z)));
ELSE IF (group=2) THEN P=(1/(1+EXP(-(THRES1-Z))))-(1/(1+EXP(-(0-Z))));
ELSE IF (group=3) THEN P=1-(1/(1+EXP(-(THRES1+THRES2-Z))));
LL=LOG(P);
MODEL group ~ general(LL);
RANDOM U ~ normal(0, SD*SD) subject=id;
RUN;

Результаты (бета-оценки, p-значение) двух методов PROC слишком сильно различаются, и я считаю, что мой код PROC NLMIXED неверен.Видите ли вы какие-либо немедленные ошибки с моим кодом PROC NLMIXED?Кроме того, как мы можем получить начальные значения бета-версий и порогов в заявлении «PARMS»?В настоящее время я использую случайные 0, 1 значения.

...