1

LIFEREGプロシージャでは、dist = gammaオプションを使用して一般化ガンマ分布を指定できます。このオプションは、3つのパラメータの一般化ガンマ分布に基づいて推定値を生成します。SASは、標準の2パラメーターのガンマ分布は利用できないと述べていますが、 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 に答える