SAS / 最大似然
< 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 语言并使用内置的优化器之一。