基尼系数简单算法_python简单计算器代码

基尼系数简单算法_python简单计算器代码用Python实现基尼系数的计算

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

Jetbrains全系列IDE稳定放心使用

参考资料:

(从文件读取数据到数组)https://blog.csdn.net/gaochen1412771148/article/details/94886847?utm_medium=distribute.pc_aggpage_search_result.none-task-blog-2~all~first_rank_v2~rank_v25-9-94886847.nonecase

(基尼系数计算方法-何史提的回答)https://www.zhihu.com/question/20219466/answer/25936162

简单的说明(吐槽):

最近对基尼系数有点感兴趣,于是就基于面向百度搜索的方法(?)写了一个python计算基尼系数的小文件(手动狗头)。

读入的文件格式设置为txt,然后它里面只有一列数据。

 

完整代码:

import os
import numpy as np
import pandas as pd

file = open("gini.txt", "r")    #以只读方式从文件夹中读取一列家庭年收入的数据文档
row = file.readlines()  #读取所有行

wealths = []    #定义一个wealths数组用于存放所有的数据

#以下开始将txt文件的数据逐个输入wealths数组
for line in row:
    line = list(line.strip().split('/n'))   #strip移除字符串头尾指定的字符(默认为空格或换行符),split指定分隔符对 字符串进行切片
    for i in line:
        wealths.append(int(i))  #逐个追加入wealths数组

#打印wealths数组,看一下有没有读入数据错误
for i in range(0,len(wealths)):
    print(wealths[i])

#基尼系数计算的函数
def gini_coef(wealths):
    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)
    yarray = cum_wealths / sum_wealths
    B = np.trapz(yarray, x=xarray)
    A = 0.5 - B
    return A / (A+B)

#打印基尼系数结果
print("\n基尼系数是%f"%gini_coef(wealths))

 

 测试代码:

首先是计算绝对平均的情况

我新建了一个“绝对平均”的txt文档,里面有100行数据,每一列都是1000,代表着100人的财产都是1000元钱。在这种绝对平均的情况下,基尼系数为0。

跑一下程序,发现结果确实为0。

基尼系数简单算法_python简单计算器代码

然后是计算绝对不平均的情况

新建一个“绝对不平均”的txt文档,里面有100行数据,前99列都为0,最后一列是100000,代表着99个人的财产都被第100个人独吞了(那99个人真悲催?) 。在这种绝对不平均的情况下,基尼系数为1。

跑了一下程序,却发现结果是0.99。

基尼系数简单算法_python简单计算器代码

呃,有没有大佬告诉一下我比较可靠的原因,现在我已经将原因归到(甩锅)电脑处理浮点数运算本来就不准确这个事实了。

最后我继续模拟了几组乱序的100个不同的数据,发现也能计算出结果,不过也还是不太知道结果的可靠性强不强。

所以如果有更加可靠的验证方法,我会继续更新的(应该吧,咕咕咕)。

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

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

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


相关推荐

  • 【贪玩巴斯】C/C++文件IO流操作的 seekp和seekg详解「建议收藏」

    【贪玩巴斯】C/C++文件IO流操作的 seekp和seekg详解「建议收藏」文件流对象有两个成员函数,分别是seekp和seekg。它们可以用于将读写位置移动到文件中的任何字节。一、如何来记这几个函数呢seek是寻找寻求的意思2.tell是告诉告知的意思那p即put放和输出的意思,在这里是保存到文件那g即get是获取,读入的意思,在这里是从文件读取所以seekp可用于将信息put(放入写入)到文件中2.seekg则可用于从文件中get(获取)信息。tellg()函数不需要带参数,它返回当前定位指针的位置,也代

    2022年6月1日
    43
  • jsp文件开头<!DOCTYPE HTML PUBLIC “-//W3C//DTD HTML 4.01 > 有什么用

    jsp文件开头<!DOCTYPE HTML PUBLIC “-//W3C//DTD HTML 4.01 > 有什么用http://www.w3.org/TR/html4/loose.dtd”> 声明位于文档中的最前面的位置,处于  标签之前。DOCTYPE是DocumentType(文档类型)的简写,此标签用来告知浏览器文档使用哪种 HTML 或 XHTML 规范。Eg:

    2022年7月14日
    31
  • FLAG_ACTIVITY_CLEAR_TOP的使用

    FLAG_ACTIVITY_CLEAR_TOP的使用本例使用FLAG_ACTIVITY_CLEAR_TOP退出整个应用程序:多activity中退出整个程序,例如从A->B->C->D,这时我需要从D直接退出程序。补充:finish()和system(0)都只能退出单个activity。我们知道Android的窗口类提供了历史栈,我们可以通过stack的原理来巧妙的实现,这里我们在D窗口打开A窗口时在Intent中直接加入标志Int

    2022年7月17日
    12
  • 安卓市场2016_鼓励大胆猜想

    安卓市场2016_鼓励大胆猜想时至今日,但凡中国的手机设计公司,要没有android手机项目,那都不好意思说自己是搞手机的。智能机替代功能机,是大势所趋,在新的一年里,结合去年一年所看所思,大胆做出一点今年的市场猜想,欢迎大家批评指教1.硬件性能瓶颈将不复存在       去年的低端android手机,基本上就是在“用户能接受多低的价格”与“用户能忍受多糟糕的体验”之间的危险博弈。那些个运营商所鼓吹的千元智能机,

    2022年9月23日
    4
  • Project interpreter not specified(eclipse+pydev)

    Project interpreter not specified(eclipse+pydev)

    2021年12月9日
    54
  • vc++菜鸟教程_java如何编写接口

    vc++菜鸟教程_java如何编写接口怎样编写自己的VCL控件    用过Delphi的朋友们,大概对Delphi的最喜欢Delphi的不是他的强类型的pascal语法,而是强大的VCL控件,本人就是一位VCL控件的爱好者。    VCL控件的开源,给我们带来了享之不尽的好处。不像以前的ole控件以及ActiveX,你完全可以重写Delphhi标准控件,而且网上这方面的资源很多。    关于如何编写VCL控件,和多

    2022年9月24日
    3

发表回复

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

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