下面是我利用proc nlmixed模块计算负二项模型：
title'NB model';
proc nlmixed data=one;
parms x1=0 x2=0 x3=0;
mu=exp(x1+x2*A1+x3*B1);
ll= lgamma (y+1/alpha)-lgamma(y+1)-lgamma(1/alpha)+y*log(alpha*mu)-(y+1/alpha)*log(1+alpha*mu);
model y~general (ll);
ods output parameterestimates=pe_nb;
run;

我想求出模型的pearson残差，画出预测的概率和实际的概率图，能否帮帮忙？谢谢了！