基尼系数计算方法_西南财经大学基尼系数

基尼系数计算方法_西南财经大学基尼系数上图图片来自维基百科·自由的百科全书最近在尝试做一些流量调控的事情,流量对于任何商业系统都是宝贵的资源,不能任由他人无限制获取,且需要通过一定手段加以调控,以期用户体验、平台生态、平台利益多方最大化。这个目标是宏大而美好的,但总要一步步慢慢来。先说如何平衡生态问题(当然指的是B端用户),生态问题就像国民收入均衡问题一样,常用的就是基尼系数,维基百科的中文链接要FQ才能打开,好伤心,英文的我…

大家好,又见面了,我是你们的朋友全栈君。如果您正在找激活码,请点击查看最新教程,关注关注公众号 “全栈程序员社区” 获取激活教程,可能之前旧版本教程已经失效.最新Idea2022.1教程亲测有效,一键激活。

Jetbrains全系列IDE稳定放心使用

上图
1244979-20171216201849889-1334466130.png

图片来自维基百科·自由的百科全书

最近在尝试做一些流量调控的事情,流量对于任何商业系统都是宝贵的资源,不能任由他人无限制获取,且需要通过一定手段加以调控,以期用户体验、平台生态、平台利益多方最大化。

这个目标是宏大而美好的,但总要一步步慢慢来。先说如何平衡生态问题(当然指的是B端用户),生态问题就像国民收入均衡问题一样,常用的就是基尼系数,维基百科的中文链接要FQ才能打开,好伤心,英文的我这里能打开。这里我并不要介绍概念和定义,只想把如何计算的方法贴上来。

最直观的计算公式是:

1244979-20171216203151124-1505368724.png

含义就是:把所有人(假设刚好 100 个人)的收入从小到大排序,然后从收入最少的开始累计,每计算一个人,横坐标为人数累计值占总人数比例,纵坐标为收入累计值占总收入比例,直到最后一个收入最大的人。

显然,图中横坐标和纵坐标都是 [ 0 -1 ] 之间。把图左下角和右上角连起来,表示人数累计占比恒等于收入累积占比,意味着收入完全均等。

好了,上面人人收入均等的线与实际收入曲线之间的面积就是 A, 实际收入曲线与 X 轴之间面积是 B

知道怎么算了,于是拿代码试试:

# coding=utf-8 # !/usr/bin/python from scipy.integrate import odeint import numpy as np from matplotlib import pyplot as pl # 解决matplotlib 中文显示问题 pl.rcParams['font.sans-serif'] = ['SimHei'] pl.rcParams['axes.unicode_minus'] = False fig, ax = pl.subplots() # 计算基尼系数的简单方法 def gini(): # 计算数组累计值,从 0 开始 wealths = [1.5, 2, 3.5, 10, 4.2, 2.1, 1.1, 2.2, 3.1, 5.1, 9.5, 9.7, 1.7, 2.3, 3.8, 1.7, 2.3, 5, 4.7, 2.3, 4.3, 12] cum_wealths = np.cumsum(sorted(np.append(wealths, 0))) # 取最后一个,也就是原数组的和 sum_wealths = cum_wealths[-1] # 人数的累积占比 xarray = np.array(range(0, len(cum_wealths))) / np.float(len(cum_wealths) - 1) # 均衡收入曲线 upper = xarray # 收入累积占比 yarray = cum_wealths / sum_wealths # 绘制基尼系数对应的洛伦兹曲线 ax.plot(xarray, yarray) ax.plot(xarray, upper) ax.set_xlabel(u'人数累积占比') ax.set_ylabel(u'收入累积占比') pl.show() # 计算曲线下面积的通用方法 B = np.trapz(yarray, x=xarray) # 总面积 0.5 A = 0.5 - B G = A / (A + B) print G 

得到下图

1244979-20171216200906999-1988143018.png


其他

当然还有很多其他的公式可以计算,比如:

1244979-20171216202119061-1910345819.png

参考

  1. https://en.wikipedia.org/wiki/Gini_coefficient 基尼系数定义
  2. https://www.zhihu.com/question/25404709 matplotlib 绘图中文乱码
  3. http://hyry.dip.jp/tech/book/page/scipynew/scipy-710-integrate.html python 科学计算

转载于:https://www.cnblogs.com/longwind09/p/8047539.html

版权声明:本文内容由互联网用户自发贡献,该文观点仅代表作者本人。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如发现本站有涉嫌侵权/违法违规的内容, 请联系我们举报,一经查实,本站将立刻删除。

发布者:全栈程序员-站长,转载请注明出处:https://javaforall.net/182480.html原文链接:https://javaforall.net

(0)
全栈程序员-站长的头像全栈程序员-站长


相关推荐

  • plc学习经验及摘录

    plc学习经验及摘录如题

    2022年10月10日
    3
  • 我的世界显示服务器领地指令,我的世界领地指令介绍 我的世界领地指令怎么设置…

    我的世界显示服务器领地指令,我的世界领地指令介绍 我的世界领地指令怎么设置…在我的世界这款经典有趣的建造类游戏中,为了让自己的领地不然其他玩家占用,我们可以设置一下领地。那我的世界领地怎么设置呢?下面是小编给大家分享的我的世界领地指令大全,大家赶紧来了解一下吧!一、我的世界设置领地:先用一块木头斧子左键敲击一方块设置点A,右键敲击一方块设置点B(可以输入“/resselectsize”查看所选区域的大小);之后输入“/rescreate123”(例)这样设置后,就…

    2022年9月23日
    2
  • html2canvas, JsPDF生成pdf

    html2canvas, JsPDF生成pdf创建 pdf js 引入依赖 importVuefro vue importhtml2c html2canvas importJsPDFf jspdf constPDF PDF install function Vue options targetDom 需要打印的 dom 对象 name pdf 的名字 callback 回调函数 Vue prototype create

    2025年7月12日
    3
  • centos7卸载MySQL_centos7下安装mysql

    centos7卸载MySQL_centos7下安装mysql**centos7下mysql卸载**1.关闭mysql服务servicemysqldstop2.删除mysql目录whereismysql——找到mysql目录rm-rfxxxxx——xxxxx代表mysql目录3.删除mysql相关配置文件my.cnf4.后记此方法针对于网络状况不好的服务器,不是yum或者rpm方法安装mysql的卸载同样有效…

    2022年9月30日
    2
  • 抽象类VS接口

    抽象类VS接口抽象类VS接口

    2022年4月24日
    33
  • html图片左右无缝循环滚动示例

    html图片左右无缝循环滚动示例首先,最外层的宽高要给个固定值,然后给个overflow:hidden;,放不下的图片隐藏。原理:(以左方向为例)先向左移动一个图片的宽度,此时第一个图片已经看不见了,这个时候将第一个标签元素克隆后追加到末尾,然后将第一个元素移除。达到循环滚动的效果。右移也是同理。下面是代码实现和demo<!DOCTYPEhtml><html><head><metacharse

    2022年7月18日
    14

发表回复

您的邮箱地址不会被公开。 必填项已用 * 标注

关注全栈程序员社区公众号