用R语言处理来自《中国工业经济》的stata数据

2022-05-28

学习来自《中国工业经济》的一篇文章(宋敏 et al., 2021)，提供了不少有价值的数据，尝试用R语言进行进一步加工。

本文主要工作如下：

构建了一个处理GDP面板数据的函数，可以计算各城市的年度GDP增长率。
构建了一个处理A股上市公司面板数据的函数，解决开头数字为0的问题，如果是4位，就在前面填充00（002中小板），如果是3位，就在前面填充000（000深A）

城市GDP面板数据

城市GDP面板数据
- 数据来源：(宋敏 et al., 2021)
- 注意：数据中的year，实际上是出报告的year，其实对应的是上一年度的year

数据处理过程

data=import("D:/博士/参考文献/一手资料/中国工业经济stata程序源代码/宋敏,2021–金融科技与企业全要素生产率——“赋能”和信贷配给的视角/数据/L1_GDP.dta")
# export(data,"地级市GDP面板数据2011_2020.xlsx")
# 计算地区GDP增长率
## 增长率 = （本年度 - 上一年度）/上一年度 
## 第一年的增长率为NA
growth_rate <- function(vec) {
(vec - dplyr::lag(vec)) / dplyr::lag(vec)
}

a <- data %>% 
group_by(city2) %>% 
arrange(year, .group = TRUE) %>%
# 假设要计算增长率的列名是GDP
mutate(growth_rate = growth_rate(GDP))

b <- a %>% drop_na()

export(list(data,a, b), sheet=c("data","all", "drop_na"), file="地级市GDP面板数据_含增长率2011_2020.xlsx")

上市公司面板数据

文章回归用的上市公司面板数据

文章回归用到的上市公司面板数据
- 数据来源：(宋敏 et al., 2021)

文章回归用到的上市公司面板数据处理过程

data=import("D:/博士/参考文献/一手资料/中国工业经济stata程序源代码/宋敏,2021–金融科技与企业全要素生产率——“赋能”和信贷配给的视角/数据/Baseline_BigLimited.dta")
# head(data)


#处理证券代码

# nchar((data[[1]]))
# 构造一个处理函数，如果是4位，就在前面填充00（002中小板），如果是3位，就在前面填充000（000深A）
id_error <- function(id) {
  if(nchar(id) == 4) {
    paste("00",id,sep = "")
    } else 
      if(nchar(id) == 3) {
        paste("000",id,sep = "")
        } else
        if(nchar(id) == 2) {
          paste("0000",id,sep = "")
          } else
            if(nchar(id) == 1) {
              paste("00000",id,sep = "")
              } else
                id
}

#将获得的id利用map函数处理成一个数据框（一开始获得的是一个list）
id <- purrr::map(data[[1]],~ id_error(.)) %>% as.data.frame() %>% t()


df <- data %>% 
  mutate(id) %>% 
  mutate(realyear=df$year-1,.before = year) %>% 
  select(id,everything())

export(df,"宋敏.2021回归用到的上市公式面板数据.xlsx")

原始的上市公司面板数据

数据来源：(宋敏 et al., 2021)

处理过程:

data=import("D:/博士/参考文献/一手资料/中国工业经济stata程序源代码/宋敏,2021–金融科技与企业全要素生产率——“赋能”和信贷配给的视角/数据/Chanel3.dta")
# head(data)


#处理证券代码

# nchar((data[[1]]))
# 构造一个处理函数，如果是4位，就在前面填充00（002中小板），如果是3位，就在前面填充000（000深A）
id_error <- function(id) {
  if(nchar(id) == 4) {
    paste("00",id,sep = "")
    } else 
      if(nchar(id) == 3) {
        paste("000",id,sep = "")
        } else
        if(nchar(id) == 2) {
          paste("0000",id,sep = "")
          } else
            if(nchar(id) == 1) {
              paste("00000",id,sep = "")
              } else
                id
}

#将获得的id利用map函数处理成一个数据框（一开始获得的是一个list）
id <- purrr::map(data[[1]],~ id_error(.)) %>% as.data.frame() %>% t()


df <- data %>% 
  mutate(id) %>% 
  select(id,everything())

export(df,"上市公式面板数据.xlsx")

参考文献

宋敏, 周鹏, & 司海涛. (2021). 金融科技与企业全要素生产率——“赋能”和信贷配给的视角. 中国工业经济, 04, 138–155. https://doi.org/10.19581/j.cnki.ciejournal.2021.04.006