Question

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

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

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

Answer 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

PROC Lifereg - 固定一个参数

1 个答案: