Dagum基尼系数分解R语言实现

最新推荐文章于 2024-03-31 09:15:17 发布

皖渝

最新推荐文章于 2024-03-31 09:15:17 发布

阅读量3.7k

点赞数

分类专栏：学术相关文章标签：数据分析 r语言数据挖掘

本文链接：https://blog.csdn.net/shine4869/article/details/123900953

版权

一、原理介绍

在这里插入图片描述

二、示例数据集

列名	含义
A列	分组列1，2，3，本次案例1，2，3为中国东、中、西部
B列	year=2008年某指标数据
C列	year=2009年某指标数据
…	…
N 列	year=2020年某指标数据

在这里插入图片描述

三、运行结果

在这里插入图片描述

四、R语言代码

rm(list = ls())
da<- read.csv('Dagum分解sample_data.csv',header=F)
group <- da[,1] # 读取分组数据
R  <-  unique(group)# 计算分组标识
nR  <-  length(R)# 分组个数
N <- dim(da)[1]
Tt <- dim(da)[2]

data0 <- da
out <- matrix(0,nrow=Tt-1,ncol=4) # 输出分解结果，每年保存一行
out_Gw <- matrix(0,nrow=Tt-1,ncol=nR) # 输出组内差异
out_Gnb <- matrix(0,nrow=Tt-1,ncol=nR*(nR-1)/2) # 输出组间差异
gR <- matrix(0,nrow=Tt-1,ncol=nR)

Dagum_Gini <- function(Data,group){
   
  if(dim(as.matrix(Data))[2] != 1)
    Data  <-  t(Data)
  if(dim(as.matrix(group))[2] != 1)
    group  <-  t(group)
  if(dim(as.matrix(Data))[2