1

在 LIFEREG 过程中，您可以使用选项指定广义 gamma 分布，该dist = gamma选项基于三参数广义 gamma 分布生成估计值。SAS 声明标准的两个参数 gamma 分布不可用，但如果可以根据http://en.wikipedia.org/wiki/Generalized_gamma_distribution将 Shape 参数固定为等于 1 就可以了。

是否可以在 LIFEREG 中修复特定参数的值，或者是否有类似 NLMIXED 的设置可能有效。作为参考，我将使用的完整代码如下所示：

proc lifereg data=work.data;
    model t*event(0) = X  / D= Gamma;   
run;
4

1 回答 1

1

您可以在数据步骤中对 2 参数伽马分布进行 MLE。片段：

s = log(meanvar) - meanlogvar;
k = (3 - s + sqrt( (s - 3)**2 + 24 * s )) / (12 * s);

do j=1 to &iterations until( abs(k - ki) < &condition );
  ki = k;
  k = ki - ( (log(ki) - digamma(ki) - s) / ((1/ki) - trigamma(ki)) );
end;
theta = meanvar / k;

请参阅：http ://en.wikipedia.org/wiki/Gamma_distribution#Maximum_likelihood_estimation

于 2013-02-08T10:10:37.810 回答