位置: IT常识 - 正文

Python统计学11——分位数回归(python 统计学包)

编辑:rootadmin
Python统计学11——分位数回归

推荐整理分享Python统计学11——分位数回归(python 统计学包),希望有所帮助,仅作参考,欢迎阅读内容。

文章相关热门搜索词:python中的统计函数,python基本统计值计算解析,python 统计学,python 统计学包,python统计图教程,python中的统计函数,用python动手学统计学,python 统计学,内容如对您有帮助,希望把文章链接给更多的朋友!

分位数回归也是数理统计里面经典的模型,他相对于在最小二乘模型上进行了改进,虽然本身还是线性的参数模型,但对损失函数进行了改进。我们都知道最小二乘的损失函数是均方误差最小,分位数的损失函数是:

可以看到分位数损失函数会对高估的值和低估的值给予一个不同的权重,这样就可以做到‘’分位‘’。

该模型对于存在异方差的数据有很好的的效果。能准确计算出5%~95%的置信区间

具体看代码理解:

导入包,加载自带的案例数据

import numpy as npimport pandas as pdimport statsmodels.api as smimport statsmodels.formula.api as smfimport matplotlib.pyplot as pltdata = sm.datasets.engel.load_pandas().datadata.head()

 

 X是收入,y是食物支出,很经典的发散数据,因为不同收入区间的家庭的食物支出比例不一样,随着X增大,Y的波动也增大。存在异方差。

q=0.5时候的分位数回归

mod = smf.quantreg("foodexp ~ income", data)res = mod.fit(q=0.5)print(res.summary())

Python统计学11——分位数回归(python 统计学包)

 当q不一样是回归出来的系数是不一样的。我们计算0.05,0.15,0.25.....0.95分位数出来的回归系数,还有最小二乘的回归系数。

quantiles = np.arange(0.05, 0.96, 0.1)def fit_model(q): res = mod.fit(q=q) return [q, res.params["Intercept"], res.params["income"]] + res.conf_int().loc["income"].tolist()models = [fit_model(x) for x in quantiles]models = pd.DataFrame(models, columns=["q", "a", "b", "lb", "ub"])ols = smf.ols("foodexp ~ income", data).fit()ols_ci = ols.conf_int().loc["income"].tolist()ols = dict(a=ols.params["Intercept"], b=ols.params["income"], lb=ols_ci[0], ub=ols_ci[1])print(models)print(ols)

画图对比x = np.arange(data.income.min(), data.income.max(), 50)get_y = lambda a, b: a + b * xfig, ax = plt.subplots(figsize=(6, 4))for i in range(models.shape[0]): y = get_y(models.a[i], models.b[i]) ax.plot(x, y, linestyle="dotted", color="grey")y = get_y(ols["a"], ols["b"])ax.plot(x, y, color="red", label="OLS")ax.scatter(data.income, data.foodexp, alpha=0.2)ax.set_xlim((200, 3000))ax.set_ylim((200, 2000))legend = ax.legend()ax.set_xlabel("Income", fontsize=16)ax.set_ylabel("Food expenditure", fontsize=16)

 透明蓝色散点为样本点。蓝色虚线为不同分位数上的回归方程。

可以看出的几个结论:

'''粮食支出随着收入的增加而增加

粮食支出的分散度随着收入的增加而增加

最小二乘估计值与低收入观测值的拟合度相当差(即OLS线越过大多数低收入家庭)'''

而且5%~95%的回归方程区间涵盖了所有的真实样本点,置信区间很准确。

 画出回归系数随着分位数的变化图

n = models.shape[0]plt.plot(models.q, models.b, color="black", label="Quantile Reg.")plt.plot(models.q, models.ub, linestyle="dotted", color="black")plt.plot(models.q, models.lb, linestyle="dotted", color="black")plt.plot(models.q, [ols["b"]] * n, color="red", label="OLS")plt.plot(models.q, [ols["lb"]] * n, linestyle="dotted", color="red")plt.plot(models.q, [ols["ub"]] * n, linestyle="dotted", color="red")plt.ylabel(r"$\beta_{income}$")plt.xlabel("Quantiles of the conditional food expenditure distribution")plt.legend()plt.show()

 

#上图画出了回归系数随着分位数的变化而变化,OLS系数是恒定的,分位数回归的系数随着分位数变大而变大

#在大多数情况下,分位数回归点估计值位于OLS置信区间之外,这表明收入对食品支出的影响在整个分布区间内可能不是恒定的

本文链接地址:https://www.jiuchutong.com/zhishi/299819.html 转载请保留说明!

上一篇:前端项目性能优化方案有哪些(前端项目性能优势怎么写)

下一篇:第四季新星计划即将开启,博客之星取消拉票你怎么看?(新星计划片头好可爱啊)

  • 影响微博营销效果的因素有哪些?(影响微博营销效益的因素)

    影响微博营销效果的因素有哪些?(影响微博营销效益的因素)

  • vivo手机滚动字幕怎么弄(vivo手机滚动字幕怎么设置)

    vivo手机滚动字幕怎么弄(vivo手机滚动字幕怎么设置)

  • 苹果猪头特效在哪(苹果猪头特效怎么变声)

    苹果猪头特效在哪(苹果猪头特效怎么变声)

  • 微信消息不同步到电脑上的原因是什么(为什么手机和电脑微信消息不同步)

    微信消息不同步到电脑上的原因是什么(为什么手机和电脑微信消息不同步)

  • 微信运动怎么删除(微信运动怎么删除关注的人)

    微信运动怎么删除(微信运动怎么删除关注的人)

  • 苹果11微信没有长按显示(苹果11微信没有提示消息 设置都开启了)

    苹果11微信没有长按显示(苹果11微信没有提示消息 设置都开启了)

  • 抖音怎么设置让别人下载自己的作品(抖音怎么设置让别人看不到我在线)

    抖音怎么设置让别人下载自己的作品(抖音怎么设置让别人看不到我在线)

  • 网络无法访问互联网怎么办(如何访问外网)

    网络无法访问互联网怎么办(如何访问外网)

  • 爱奇艺聊天室可以回看吗(爱奇艺聊天室可以投屏吗)

    爱奇艺聊天室可以回看吗(爱奇艺聊天室可以投屏吗)

  • qq标识小海螺是什么(qq海螺是什么意思)

    qq标识小海螺是什么(qq海螺是什么意思)

  • 苹果x来消息灯光闪怎么设置(苹果x消息提醒闪光灯开了不会闪怎么设置)

    苹果x来消息灯光闪怎么设置(苹果x消息提醒闪光灯开了不会闪怎么设置)

  • 苹果录音发生异常无法分享(苹果录音发生异常无法分享微信怎么办)

    苹果录音发生异常无法分享(苹果录音发生异常无法分享微信怎么办)

  • 退群后分享的二维码还有效吗(退群后分享的二维码是啥)

    退群后分享的二维码还有效吗(退群后分享的二维码是啥)

  • 拼多多免拼券怎么获得(拼多多免拼券怎么用)

    拼多多免拼券怎么获得(拼多多免拼券怎么用)

  • 苹果红包来了提醒功能(iphone红包来了)

    苹果红包来了提醒功能(iphone红包来了)

  • 华为畅享9plus带不带nfc(华为畅享9plus带红外线功能吗)

    华为畅享9plus带不带nfc(华为畅享9plus带红外线功能吗)

  • 拼多多怎么手动赠送免拼卡(拼多多怎么手动打印快递单)

    拼多多怎么手动赠送免拼卡(拼多多怎么手动打印快递单)

  • 华为bkk一al00是什么型号(bkl-al00华为什么型号)

    华为bkk一al00是什么型号(bkl-al00华为什么型号)

  • vivox9怎么备份(vivox9怎么备份微信数据)

    vivox9怎么备份(vivox9怎么备份微信数据)

  • 小红书怎么找商品链接(小红书怎么找商务合作)

    小红书怎么找商品链接(小红书怎么找商务合作)

  • 微机的开机顺序是什么(微机的开机顺序类似答案)

    微机的开机顺序是什么(微机的开机顺序类似答案)

  • 闲鱼小法庭能二次申诉吗(闲鱼小法庭能二次开庭吗)

    闲鱼小法庭能二次申诉吗(闲鱼小法庭能二次开庭吗)

  • 华为mate20手机怎么打开使用密码保险箱(华为mate20手机怎么强制关机)

    华为mate20手机怎么打开使用密码保险箱(华为mate20手机怎么强制关机)

  • Chat GPT实用案例——VUE+Chat GPT实现聊天功能教程

    Chat GPT实用案例——VUE+Chat GPT实现聊天功能教程

  • 递延所得税负债怎么算
  • 销售赠送赠品会计处理
  • 背书行为有效的是指什么
  • 资产负债表的货币资金根据什么填
  • 其他综合收益影响所有者权益吗
  • 开办费计入管理费用后是否要转出
  • 法定盈余公积金的提取比例一般按照()
  • 处置固定资产增值税和所得税报税能一致吗
  • 固定资产发票未到怎么提折旧
  • 实收资本借方增加还是减少
  • 存款取款业务记入那个日记账?
  • 含税和不含税哪个对卖方有利
  • 哪种情形
  • 企业所得税核定征收最新政策
  • 股权案例分析2020
  • 发票验旧后可领新发票吗
  • 增值税要转到营业税金及附加吗
  • 购买土地建厂房还要交土地使用费吗
  • 限额领料单是一种多次使用的
  • 事业基金的主要账务处理
  • 应用程序错误0x000000
  • thinkphp5控制器
  • php循环语句
  • php调用sql
  • 如何修复错误 OC3 INVT CHF100a
  • 企业所得税计税依据
  • 2019年下半年中小学教师资格考试综合素质试题
  • RuntimeError: CUDA out of memory 已解决
  • php语句和php变量都是区分大小写的
  • 贴现会不会损害所有者权益
  • 基于springboot的毕设
  • 上月开了红字发票当月如何申报
  • HTTP 错误 500.19- Internal Server Error 错误解决方法
  • 劳务公司成本核算方法选择什么
  • 房地产预缴所得税计算公式
  • tcpdump 监听指定端口
  • 金融机构存放的保证金存款
  • 解决什么情况
  • java微信公众号开发源码
  • 固定资金的概念及其特点
  • 民间非盈利组织会计要素组成
  • 打印银行电子回单在打印机上怎么设置A5
  • Ficstar Software � Web Data Extraction Tools
  • 缴纳税收的凭证怎么打印
  • 罚款可以直接在12123
  • 出口免税申报流程视频
  • 收到退还的增值税
  • 企业期末预收账款怎么算
  • 取得收入未开具发票
  • 建筑工程人工费调整
  • 差旅费记入应付职工薪酬吗
  • 分配现金股利如何做会计分录长投
  • 车辆按揭贷款需要什么
  • 怎么处理固定资产残值
  • 以前年度损益调整借贷方向
  • 学校收到政府补助的会计分录
  • 一般纳税人只交社保无发票往来需要做账吗
  • 销售折让的会计账务处理分录
  • 公转私钱最好的方法
  • 备查账目
  • mysql中具体到删某一个数据
  • windows mysql my.cnf
  • xp win10 打印机
  • linux中samba服务器的设计方案
  • os x yosemite10.10.5怎么样?os x yosemite10.10.5公测版下载地址
  • linux怎么用u盘传输文件
  • 怎么解决xp不能安装软件
  • win10怎么把中文系统改成英文
  • win7微软账户
  • css行与行之间的间距怎么调
  • 安卓字库ic
  • python的理解
  • 怎么做一个扇形图
  • 举例讲解生产可能性曲线
  • jquery的validate前端表单验证
  • 安卓开发常见问题
  • 电子发票未验真是假发票吗
  • 合伙人转让出资的法律规定
  • 广告发布合同属于什么合同
  • 关于新时代新征程
  • 免责声明:网站部分图片文字素材来源于网络,如有侵权,请及时告知,我们会第一时间删除,谢谢! 邮箱:opceo@qq.com

    鄂ICP备2023003026号

    网站地图: 企业信息 工商信息 财税知识 网络常识 编程技术

    友情链接: 武汉网站建设