楼主: 序幕、
3042 2

[有偿编程] 如何拟合泊松回归? [推广有奖]

  • 0关注
  • 0粉丝

初中生

66%

还不是VIP/贵宾

-

威望
0
论坛币
8 个
通用积分
0
学术水平
0 点
热心指数
0 点
信用等级
0 点
经验
143 点
帖子
7
精华
0
在线时间
13 小时
注册时间
2016-9-27
最后登录
2017-6-19

+2 论坛币
k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群
赵安豆老师微信:zhaoandou666

经管之家联合CDA

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

感谢您参与论坛问题回答

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

+2 论坛币
  1. data insure;
  2.    input n c car$ age;
  3.    ln = log(n);
  4.    datalines;
  5. 500   42  small  1
  6. 1200  37  medium 1
  7. 100    1  large  1
  8. 400  101  small  2
  9. 500   73  medium 2
  10. 300   14  large  2
  11. ;
  12. run;
复制代码

上述为sashelp中的例子,请问如何用proc genmod过程拟合出该泊松回归并计算出car变量的OR值(基线为small)?
二维码

扫码加我 拉你入群

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

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

关键词:泊松回归 SASHELP Medium genmod Small 如何

沙发
序幕、 发表于 2017-6-19 22:43:19 |只看作者 |坛友微信交流群
有没有大神来解答一下~~

使用道具

藤椅
地下爆菊 发表于 2017-6-20 16:48:56 |只看作者 |坛友微信交流群
c = 0.767 - 3.36/lambda
beta = PI/sqrt(3.0*lambda)
alpha = beta*lambda
k = log(c) - lambda - log(beta)

forever
{
        u = random()
        x = (alpha - log((1.0 - u)/u))/beta
        n = floor(x + 0.5)
        if (n < 0)
                continue
        v = random()
        y = alpha - beta*x
        lhs = y + log(v/(1.0 + exp(y))^2)
        rhs = k + n*log(lambda) - log(n!)
        if (lhs <= rhs)
                return n
}

使用道具

您需要登录后才可以回帖 登录 | 我要注册

本版微信群
加好友,备注cda
拉您进交流群

京ICP备16021002-2号 京B2-20170662号 京公网安备 11010802022788号 论坛法律顾问:王进律师 知识产权保护声明   免责及隐私声明

GMT+8, 2024-4-20 22:45