我有一个逻辑回归模型:
proc genmod data=rats;
model cured/total = dose
/dist = binomial link=logit;
run;
我想要做的是根据固化/总量的给定比例计算或得出剂量。我意识到我可以根据它给出的模型来做代数,但我想知道是否有办法使用SAS。
答案 0 :(得分:1)
您正在寻找一个逆预测,可以使用带有INVERSECL选项的PROBIT程序获得:
proc probit data=rats plots=(predpplot ippplot lpredplot) ;
model cured/total = dose / dist=logistic inversecl;
run;
文档中有一个很好的例子: SAS/STAT 14.1 PROBIT Example 93.4 An epidemiology study
答案 1 :(得分:0)
是的,您可以使用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 documentation
现在,确实做了你想要的那种反向评分,只需制作一个从一些最小值到一些最大值的数据集,增量很小
答案 2 :(得分:0)
(我为第二个问题道歉"答案",我没有足够的评论点。)
Peter,PROBIT的优势在于可用于逆预测的区间估计(基准限制)。