Can'；SAS中使用零膨胀泊松回归模型的t分数测试集

我使用proc genmod运行了一个零膨胀泊松模型，我试图使用proc PLM对我的测试数据集进行评分，但它给了我以下错误：

proc genmod data = train2;
class region  / param=glm;
model response = var1 var2 var3 var4 var5
                    / dist=zip;
                    zeromodel;
output out = zeropoisson_output predicted= estprobzip;
store zero_poisson;
run;


proc plm source=zero_poisson;
  score data = test2 out= pred_zip;
run;

错误：此版本的PLM程序不支持为零充气模型评分

---

有什么办法可以解决这个问题吗

这将需要相当多的工作，但是您可以始终使用ODS输出选项来获取参数估计值，并从中解析数据。我从genmod上的SAS示例中获取了一些示例数据，并在下面演示了保存系数和解析系数的概念。输出是一个.sas文件，可以%包含在任何数据步骤中，以对验证样本进行评分

data drug;
   input drug$ x r n @@;
   datalines;
    A  .1   1  10   A  .23  2  12   A  .67  1   9
    B  .2   3  13   B  .3   4  15   B  .45  5  16   B  .78  5  13
    C  .04  0  10   C  .15  0  11   C  .56  1  12   C  .7   2  12
    D  .34  5  10   D  .6   5   9   D  .7   8  10
    E  .2  12  20   E  .34 15  20   E  .56 13  15   E  .8  17  20
    ;
run;

ods output ParameterEstimates = ZIP_COEFF_EST;
proc genmod data=drug;
      class drug;
      model r/n = x drug / dist = zip;
      zeromodel;
run;
ods output close;

data ZIP_COEFF_EST_Parsed;
    length equation $ 2500;
    set ZIP_COEFF_EST (rename=(estimate=coefficient)) end=last;
    where coefficient ne .;
    if upcase(Parameter) = "INTERCEPT" then do;
        equation = " = " || trim(left(put(coefficient,20.10)));
        output;
    end;
    else if LEVEL1 ne '' then do;
        equation = " + (" || trim(left(Parameter)) || " = '" || trim(left(LEVEL1)) || "') * (" || trim(left(put(coefficient,20.10))) || ")";
        output; 
    end;
    else do;
        equation = " + " || trim(left(Parameter)) || " * (" || trim(left(put(coefficient,20.10))) || ")";
        output;
    end;
    if last then do;
        equation=';';   
        output;
    end;
    keep equation;
run;

data _null_;
    set ZIP_COEFF_EST_Parsed;
    FILE  "C:/estimate_file.sas";;
    PUT equation;
run; 

---

具体来说，SAS/STAT的版本是什么？（从SAS 9.22开始，他们开始发布单独的版本号，从12.1开始，从SAS 9.3开始）。我使用的是SAS EG 4.3，但SAS Base是9.3-我不知道如何获得其他版本信息，或者是否有其他版本信息？您可以查看服务器（如果您在桌面上执行此操作，通常为本地服务器，否则选择服务器），右键单击，选择属性，然后选择“查看初始化日志”，其中将包含有关您的版本的信息，包括SAS/STAT版本。它看起来像SAS/STAT 13.1（基本9.4）是他们在PLM中添加了零膨胀法评分的地方，请参见并对比一下，在其左栏上没有9.4页的条目。啊……有什么想法我可以手动评分吗？如果是正常的logit/poisson模型，那么我可以开始使用可变系数手动评分，但使用零膨胀模型，我是l我不知道该怎么办？