楼主: 序幕、
375 2

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

  • 0关注
  • 0粉丝

初中生

66%

还不是VIP/贵宾

-

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

序幕、 发表于 2017-6-19 15:42:53 |显示全部楼层
  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)?

支持楼主:购买VIP购买贵宾 购买后,论坛将把您花费的资金全部奖励给楼主,以表示您对TA发好贴的支持
 
载入中......
stata SPSS
序幕、 发表于 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
}
回复

使用道具 举报

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

GMT+8, 2017-11-22 05:24