PROC Lifereg - 保存一個固定參數

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