SAS 与 Stata：伽玛分布函数

Question

我正在尝试将 Stata 程序逆向工程为 SAS，但我在 Stata 中的 gammap 函数及其与 SAS 中的关联方面遇到了问题。从 Stata 文档看来，gammap 函数 returns 是累积伽玛分布。

测试数据：

PSCORE  PALPHA  PBETA
0.032352097 21.4639 0.002864125
0.030794526 21.4639 0.002864125
0.032952468 21.4639 0.002864125
0.041141297 21.4639 0.002864125
0.033376449 21.4639 0.002864125
0.032352097 5.7865  0.005516187
0.030794526 5.7865  0.005516187
0.032952468 5.7865  0.005516187
0.041141297 5.7865  0.005516187
0.033376449 5.7865  0.005516187

所以在 Stata 程序中，我有：

RESULT = gammap(PALPHA,PSCORE/PBETA)

我在 SAS 中将其翻译为 CDF 函数的一部分：

RESULT = CDF('GAMMA',PALPHA,PSCORE/PBETA);

然而，结果并不匹配，所以我的理论是错误的。

统计数据：

0.0045025
0.0025791
0.0055082
0.0474779
0.0063245
0.5680494
0.5214182
0.5854155
0.7804033
0.5974566

SAS：

0.99394
0.99605
0.99291
0.95694
0.99209
0.54212
0.58959
0.52384
0.29454
0.51097

有人可以提供有关正确的关联 SAS 代码可能是什么以及我在使用 CDF 函数时出错的地方的见解吗？我不应该使用 CDF 函数吗？

Answer 1

SAS CDF 函数采用 2 个 Gamma 分布参数。

data have;
input PSCORE  PALPHA  PBETA;

datalines;
0.032352097 21.4639 0.002864125
0.030794526 21.4639 0.002864125
0.032952468 21.4639 0.002864125
0.041141297 21.4639 0.002864125
0.033376449 21.4639 0.002864125
0.032352097 5.7865  0.005516187
0.030794526 5.7865  0.005516187
0.032952468 5.7865  0.005516187
0.041141297 5.7865  0.005516187
0.033376449 5.7865  0.005516187
;

data want;
set have;
result = cdf('gamma',pscore,palpha,pbeta);
run;

WANT 的结果似乎与您的 STATA 结果一致。