全部版块 我的主页
论坛 数据科学与人工智能 数据分析与数据科学 SAS专版
2801 11
2012-11-01
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?这是为什么呢?

二维码

扫码加我 拉你入群

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

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

全部回复
2012-11-1 20:38:48
不懂,帮顶
二维码

扫码加我 拉你入群

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

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

2012-11-1 20:49:34
将得到的值,指数化,看等不等
二维码

扫码加我 拉你入群

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

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

2012-11-1 21:25:37
ziyenano 发表于 2012-11-1 20:49
将得到的值,指数化,看等不等
我想问下 这个pred是不是预测值?
我得到了估计值,也指数化了,我用第一个值算exp(b0+b1*x1+ b2*x2+b3*x3)和第一个pred比,但是不等!
二维码

扫码加我 拉你入群

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

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

2012-11-1 21:47:01
proc nlmixed有什么办法可以得到预测值?
二维码

扫码加我 拉你入群

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

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

2012-11-1 21:49:35
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相等
二维码

扫码加我 拉你入群

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

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

点击查看更多内容…
相关推荐
栏目导航
热门文章
推荐文章

说点什么

分享

扫码加好友,拉您进群
各岗位、行业、专业交流群