楼主: 有福有德
5497 5

[学习分享] 累积logit模型 [推广有奖]

教师

院士

6%

还不是VIP/贵宾

-

威望
4
论坛币
56500 个
通用积分
8.4735
学术水平
344 点
热心指数
203 点
信用等级
274 点
经验
37400 点
帖子
1138
精华
10
在线时间
945 小时
注册时间
2010-10-14
最后登录
2024-3-21

+2 论坛币
k人 参与回答

经管之家送您一份

应届毕业生专属福利!

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

经管之家联合CDA

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

感谢您参与论坛问题回答

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

+2 论坛币
  1. formula:以因变量为3个水平数为例(p1、p2、p3):
  2. logit(p1/1-p1)=-a_1+b_1*x_1+~~+b_p*x_p,
  3. logit((p1+p2)/1-(p1+p2))=-a_2+b_1*x_1+~~+b_p*x_p
  4. notes:系数不变,截距在变。

  5. 计算概率:
  6. p1=exp(-a_1+b_1*x_1+~~+b_p*x_p)/(1+exp(-a_1+b_1*x_1+~~+b_p*x_p)),
  7. p2=exp(-a_2+b_1*x_1+~~+b_p*x_p)/(1+exp(-a_2+b_1*x_1+~~+b_p*x_p))-p1,
  8. p3=1-p1-p2;
复制代码
    累积logit模型的建模过程有点特别,看一下formula:是3个取值的因变量对应的概率p1、p2、p3
    等号左边计算的是累积概率,右边的回归截距是不同的,但所有的自变量及其系数都是相同的,下面是计算概率的公式。
    下面是SAS中经常使用的语句:

  1. proc genmod data=temp rorder=data;
  2. class sex area;
  3. model purchase=area consume
  4. /dist=multinomial
  5. link=cumlogit
  6. aggregate;
  7. output out=a predicted=p1;
  8. run;
复制代码
           说明:1)predicted:计算Pr(Y<=_LEVEL_)的概率,从而可以得到p1、p2、p3的概率。
                       2)dist=multinomial和link=cumlogit指定分布形式和连接函数,这个固定的,不需要改;
                       3)aggregate:报告数据的离散(或集中)信息,广义线性模型均需要提供数据离散信息(devience/df是否等于1),如果有问题可以通过scale调整,也可以贝叶斯的方法处理。

          SAS logistic语句也可以实现这种方法,不过genmod使用的更频繁。
         可以看出该模型对数据有很严格的要求(系数不变,截距在变),经常在使用这种模型的时候,也会将这种分类视为无序分类从而使用其他模型处理。例如无序logistic模型等。

二维码

扫码加我 拉你入群

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

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

关键词:logit模型 logit Log Multinomial Predicted 模型 计算 因变量

已有 2 人评分经验 论坛币 学术水平 热心指数 信用等级 收起 理由
hamsik11 + 1 + 1 + 1 观点有启发
webgu + 60 + 10 + 1 + 1 + 1 分析的有道理

总评分: 经验 + 60  论坛币 + 10  学术水平 + 2  热心指数 + 2  信用等级 + 2   查看全部评分

本帖被以下文库推荐

所有模型都是错的
沙发
stata18 发表于 2013-4-4 09:41:08 |只看作者 |坛友微信交流群
学习

使用道具

藤椅
fxwh123 发表于 2013-4-9 11:37:17 |只看作者 |坛友微信交流群
感谢楼主分享~~

使用道具

板凳
1道疤 发表于 2013-4-10 10:11:02 |只看作者 |坛友微信交流群
学习学习

使用道具

报纸
liyuzhong7856 发表于 2014-12-28 15:11:16 |只看作者 |坛友微信交流群
多谢楼主。请问能实现累积logit的SAS模块都有哪些选择?有没有相关学习资料可以推荐一下,谢谢!

使用道具

地板
mxn2019 发表于 2021-2-1 10:15:09 |只看作者 |坛友微信交流群
请问累计概率模型的系数如何解读? 现知道系数求出了风险发生比率(or=1.63)

使用道具

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

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

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

GMT+8, 2024-4-25 08:45