SAS 与 Stata：伽玛分布函数

Question

我正在尝试将 Stata 程序逆向工程到 SAS 中，但我遇到了 Stata 中的 gammap 函数以及它与 SAS 中的关联的问题。从 Stata 文档看来，gammap 函数返回了累积伽马分布。

测试数据：

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 结果一致。