Могу ли я использовать шаг Proc внутри другого шага proc?
Я пишу nlmixed
процедуры, и для функции правдоподобия я хочу использовать значения, сгенерированные с iml
процедура. Поэтому мне интересно, есть ли способ использования proc iml
внутри proc nlmixed
,
proc nlmixed data = xxx;
parms b0=0 b1=0;
mu = exp(b0 + b1*Age);
ll = log(((mu**y)*exp(-mu))/gamma(y+1));
model y~ general(ll);
run;
proc iml;
v = {5,6,7,8,9,10,11,12,13,14};
z = j(10,1,.);
do i = 1 to 6;
z[i] = ((v[i]-5)/5)*((mu**v[i])*exp(-mu))/gamma(v[i]+1);
end;
ll=log(sum(z));
quit;
Идея состоит в том, чтобы:
- использовать му из
nlmixed
внутриproc iml
- но оба шага должны быть внутри
nlmixed