帮我看看这段代码好吗？proc nlmixed

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

送您一个全额奖学金名额~ !

立即领取

感谢您参与论坛问题回答

经管之家送您两个论坛币！

+2 论坛币

proc nlmixed data=clmfreq;
parms b0=0 b1=0 b2=0 b3=0;
linp=b0+b1*x1+ b2*x2+b3*x3;
mu=exp(linp);
model  count~Poisson(mu);
predict pdf('Poisson',count,mu) out=poiOut(keep= count pred);
run;
请问poiout数据集的pred是计数的预测值吗？
我将 poiout输出，然后代入估计的参数，发现估计值不是pred?这是为什么呢？

扫码加我 拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

分享0 收藏0 回帖

关键词：proc nlmixed nlmixed Mixed mix nlm count

不懂，帮顶

将得到的值，指数化，看等不等

ziyenano 发表于 2012-11-1 20:49
将得到的值，指数化，看等不等

我想问下 这个pred是不是预测值？
我得到了估计值，也指数化了，我用第一个值算exp(b0+b1*x1+ b2*x2+b3*x3)和第一个pred比，但是不等！

proc nlmixed有什么办法可以得到预测值？

Silence. 发表于 2012-11-1 21:25
我想问下 这个pred是不是预测值？
我得到了估计值，也指数化了，我用第一个值算exp(b0+b1*x1+ b2*x2+b3* ...

mu值是预测值；
predict pdf('Poisson',count,mu)
pred应该是：以mu为poison分布的均值，k=count时候的概率值；

data _null_;
a=pdf('poison',count值,exp(b0+b1*x1+ b2*x2+b3*x3)的值);
put a=;
run;
看a值是否和pred相等

ziyenano 发表于 2012-11-1 21:49
mu值是预测值；
predict pdf('Poisson',count,mu)
pred应该是：以mu为poison分布的均值，k=count时候的 ...

谢谢啦！

Silence. 发表于 2012-11-1 22:27
谢谢啦！

对否？因为nlmixed用的不多，poison回归一般用genmod来做的，这里只是根据语句大概猜的意思。

Silence. 发表于 2012-11-1 21:47
proc nlmixed有什么办法可以得到预测值？

proc nlmixed data=clmfreq;
parms b0=0 b1=0 b2=0 b3=0;
linp=b0+b1*x1+ b2*x2+b3*x3;
mu=exp(linp);
model  count~Poisson(mu);
predict pdf('Poisson',count,mu) out=poiOut(keep= count pred);
predict mu out=poiOut1(keep= count pred);
run;

Silence. 发表于 2012-11-1 22:27
谢谢啦！

是的 pred是那点预测的概率，而不是以预测的值。