SAS / 最大似然

• PROC NLMIXED 允許您指定似然

proc nlmixed data = recidiv ;
parms cst = 0 b_male = 0 b_black = 0 b_release =0 b_age = 0;
xb = cst + b_male*male + b_black*black + b_release*release + b_age*age__days_1000_ ; 
if c = 1 then prob = exp(-exp(xb)*y);
if c = 0 then prob = exp(xb)*exp(-exp(xb)*y);
ll= log(prob);
model y ~ general(ll);
run;

• parms 向引數提供初始值
• xb 是一種中間步驟，提供協變數列表
• 然後我們定義似然並指定模型

• 還可以使用 IML 語言並使用內建的最佳化器之一。