计算SAS中的预测值

Question

我有一个逻辑回归模型：

 proc genmod data=rats;
     model cured/total = dose
           /dist = binomial link=logit;
 run;

我想做的是根据治愈/总剂量的给定比例计算或推导剂量。我意识到我可以根据它给出的模型来做代数，但我想知道是否有办法使用SAS来做。

Answer 1

您正在寻找一个反向预测，可以使用带有INVERSECL选项的PROBIT过程获得：

proc probit data=rats plots=(predpplot ippplot lpredplot)           ;
model cured/total = dose / dist=logistic inversecl;
run;

文档中有一个很好的例子：SAS/STAT 14.1 PROBIT示例93.4一项流行病学研究

Answer 2

是的，您可以在GENMOD中使用STORE命令，然后在新数据集上执行PROC PLM。简言之：

  proc genmod data=remiss descending;
        model remiss = smear blast / dist=binomial;
        store out=logmod;
  run;
  proc plm source=logmod;
        score data=new out=preds pred=pred lclm=lower uclm=upper / ilink;
  run;

有关详细信息，请参阅SAS文档

现在，做你想要的那种反向评分，只需制作一个数据集，以小增量从某个最小值到某个最大值

Answer 3

（我为第二个"答案"道歉，我没有足够的分数来评论。）

Peter，PROBIT的优点是可以获得用于反向预测的区间估计（基准极限）。