系列文章

  1. R做面板数据分析:R语言 面板数据分析 plm包实现(一) ——LSDV和固定效应模型

  2. 如果想看随机效应模型怎么做,参见这篇文章
    R语言 面板数据分析 plm包实现(二)——随机效应模型

  3. 如果想看如何判断面板数据适用随机效应模型还是固定效应模型,参见这篇文章:
    R语言 面板数据分析 plm包实现(三)——面板数据与面板模型的检验

  4. 使用随机效应模型,且一些时间或个体存在数据缺失,应当使用Swamy Arora估计,如何用R语言来实现,参见这篇文章:
    R语言 面板数据如何做Swamy Arora估计。

  5. 如何从统计年鉴获知某地接受转移支付额度

目录

  • 系列文章
  • 1.模型描述
    • 数据导入
  • 2.假定γt =0,直接估计随机 “Individual effects” 模型
  • 3.假定ui =0,直接估计随机“Time effects”模型;

首先是plm包安装和数据导入部分,参见文章:
《R语言 面板数据分析 plm包实现(一) ——LSDV和固定效应模型》
使用随机效应模型,且一些时间或个体存在数据缺失,应当使用Swamy Arora估计,如何用R语言来实现,参见这篇文章:
R语言 面板数据如何做Swamy Arora估计。

1.模型描述

有数据集:Ex1_1.dta
数据样式:
点击下载
其中FN代表公司,总共有三家;YR代表年份;I是总投资;F是企业实际价值;C是企业实际资本存量。

更多解释:

数据导入

这个数据集是stata的数据集,因此在Rstudio中你可以选择文件–>导入数据集(import dataset)–>导入stata文件,即可完成导入工作

此外,我好像在其它地方也看见过此数据集,如果你无法下载,可以在其它地方寻找数据集(我印象里是在某个面板相关的R程序包里自带的数据集)。

很多童鞋反映数据集获取困难,我把这个数据集上传到github的一个项目里了(免费),注意,只有一个文件是数据集。如果有帮到你,请给文章点个赞哦~

2.假定γt =0,直接估计随机 “Individual effects” 模型

# R codes
rankData<-pdata.frame(Ex1_1,index=c("FN","YR")) #index里是个体和时间,转化为面板数据
random <- plm(I~ F + C,data=rankData,effect = "individual", model="random")
summary(random)

输出结果:

Oneway (individual) effect Random Effect Model (Swamy-Arora's transformation)Call:
plm(formula = I ~ F + C, data = rankData, effect = "individual", model = "random")Balanced Panel: n = 10, T = 20, N = 200Effects:var std.dev share
idiosyncratic 2784.46   52.77 0.282
individual    7089.80   84.20 0.718
theta: 0.8612Residuals:Min.   1st Qu.    Median   3rd Qu.      Max.
-177.6063  -19.7350    4.6851   19.5105  252.8743 Coefficients:Estimate Std. Error z-value Pr(>|z|)
(Intercept) -57.834415  28.898935 -2.0013  0.04536 *
F             0.109781   0.010493 10.4627  < 2e-16 ***
C             0.308113   0.017180 17.9339  < 2e-16 ***
---
Signif. codes:  0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1Total Sum of Squares:    2381400
Residual Sum of Squares: 548900
R-Squared:      0.7695
Adj. R-Squared: 0.76716
Chisq: 657.674 on 2 DF, p-value: < 2.22e-16

即广义最小二乘法GLS,结果解读课参见上一篇文章,此处略去。

3.假定ui =0,直接估计随机“Time effects”模型;

# R codes
rankData<-pdata.frame(Ex1_1,index=c("FN","YR")) #index里是个体和时间,转化为面板数据
random2 <- plm(I~ F + C,data=rankData,effect = "time", model="random")
summary(random2)

结果:

Oneway (time) effect Random Effect Model (Swamy-Arora's transformation)Call:
plm(formula = I ~ F + C, data = rankData, effect = "time", model = "random")Balanced Panel: n = 10, T = 20, N = 200Effects:var std.dev share
idiosyncratic 9623.4    98.1     1
time             0.0     0.0     0
theta: 0Residuals:Min.   1st Qu.    Median   3rd Qu.      Max.
-291.6757  -30.0137    5.3033   34.8293  369.4464 Coefficients:Estimate  Std. Error z-value  Pr(>|z|)
(Intercept) -42.7143694   9.5116760 -4.4907 7.098e-06 ***
F             0.1155622   0.0058357 19.8026 < 2.2e-16 ***
C             0.2306785   0.0254758  9.0548 < 2.2e-16 ***
---
Signif. codes:  0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1Total Sum of Squares:    9359900
Residual Sum of Squares: 1755900
R-Squared:      0.81241
Adj. R-Squared: 0.8105
Chisq: 853.151 on 2 DF, p-value: < 2.22e-16

下一期预告:如何做检验,来确定该使用哪一种模型。

R语言 面板数据分析 plm包实现(二)——随机效应模型相关推荐

  1. R语言 面板数据分析 plm包实现(一) ——LSDV和固定效应模型

    系列文章 R做面板数据分析:R语言 面板数据分析 plm包实现(一) --LSDV和固定效应模型 如果想看随机效应模型怎么做,参见这篇文章 R语言 面板数据分析 plm包实现(二)--随机效应模型 如 ...

  2. mcem r语言代码_处理潜变量或随机效应的高维积分时,我们做些什么之MCEM

    本文面向对象是心理和教育等社科统计的初入门者. 广义线性潜变量模型(Generalized Linear Latent Variable Models,在心理学和教育学常用的是连接函数为probit或 ...

  3. r语言 面板数据回归_工具方法 | “名牌包”:面板、时间序列模型常用R语言包...

    计量经济学是数学.统计技术和经济分析的综合,即运用数学.统计方法和相关经济理论,通过计量模型来揭示经济数量关系和规律.R语言包,已经实现了现代计量经济学的很多统计分析功能,下面从面板数据模型和时间序列 ...

  4. R语言与数据分析练习:使用ARIMA模型预测网站访问量

    R语言与数据分析练习:使用ARIMA模型预测网站访问量 使用ARIMA模型预测网站访问量 一.实验背景: 随着流量的增大,某网站的数据信息量也在以一定的幅度增长 基于该网站2016年9月~2017年2 ...

  5. R语言与数据分析练习:创建和使用R语言数据集数据的导入导出

    R语言与数据分析练习:创建和使用R语言数据集&数据的导入导出 实验一 创建和使用R语言数据集 一.实验目的: 了解R语言中的数据结构. 熟练掌握他们的创建方法,和函数中一些参数的使用. 对创建 ...

  6. 《R语言游戏数据分析与挖掘》一2.2 数据对象

    2.2 数据对象 R拥有许多用于存储数据的对象类型,包括向量.矩阵.数组.数据框和列表.它们在存储数据的类型.创建方式.结构复杂度,以及用于定位和访问其中个别元素的标记等方面均有所不同.多样化的数据对 ...

  7. R语言入门第三集 实验二:基本数据处理

    R语言入门第三集 实验二:基本数据处理 一.资源 [R语言]R语言数据处理--东北大学大数据班R实训第二次作业 二.答案更新纠正 2.11.从df中选取date . item_id . cate_id ...

  8. R语言数据科学程序包:Tidyverse介绍

    R语言数据科学程序包:Tidyverse介绍 1. R语言简介 2. 数据科学简介 3. Tidyverse简介 1. R语言简介 R语言是用于统计计算和绘图的免费软件.它可以在Windows, Un ...

  9. 2014全球软件技术峰会WOT:R语言金融数据分析

    跨界知识聚会系列文章,"知识是用来分享和传承的",各种会议.论坛.沙龙都是分享知识的绝佳场所.我也有幸作为演讲嘉宾参加了一些国内的大型会议,向大家展示我所做的一些成果.从听众到演讲 ...

  10. 【R语言与数据分析实战】绘图

    目录 1.散点图 2.图像选项 2.1 坐标轴名称 2.2 图形标题 2.3 点的类型 2.4 点的大小 2.5 颜色 2.6 坐标轴的取值范围 2.7 图像类型 2.8 线型 2.9 图形排列 2. ...

最新文章

  1. 积微论坛报告视频+PPT:用微生物组时序数据重现生物膜装配动态过程
  2. golang中的strings.EqualFold
  3. 20万用户同时访问一个热点Key,如何优化缓存架构?
  4. 通向架构师的道路(第十一天)之Axis2 Web Service(二)
  5. border-collapse 关于继承问题
  6. docker拉取的pytorch-gpu版找不到cuda和cudnn的位置
  7. java五子棋人机对战_实现简单的人机对战五子棋(实践)
  8. LETTers比赛第三场 --1002 Ignatius and the Princess III解题报告
  9. I2C总线协议原理介绍
  10. Hibernate 二级缓存和查询缓存
  11. minHJ函数matlab黄金分割法,黄金分割法求极小值
  12. android开发学习:打电话和发短信
  13. Linux中ps命令详解
  14. 树的高度单位换算 c语言编辑,单位换算公式大全
  15. PMP之假设日志Assumption Log
  16. Python 发送邮件 email 模块、smtplib 模块
  17. WMS与企业ERP/MES数据交换方法
  18. Java 17最新免费商业版本
  19. Python语言程序设计基础 第二版(嵩天著)课后答案第三章(程序练习题)
  20. 电商工具箱针对补单平台商品主图获取,照妖镜验号,插旗,打标,淘口令等API接口提供

热门文章

  1. QCC3040---FS(file system)中添加文件方法
  2. 报错Ordinal parameter not bound
  3. 数字视频编解码基础知识大全
  4. 过程现场总线解决方案的全球与中国市场2022-2028年:技术、参与者、趋势、市场规模及占有率研究报告
  5. cosh和acosh--双曲余弦和反双曲余弦函数
  6. 应用未安装!安装包似乎已经损坏
  7. java xml特殊字符_mybatis xml中特殊字符处理及特殊符号
  8. AutoPatch java_EBS开发技术之Patch安装
  9. linux配置路由器命令手册,Linux配置路由器
  10. Vim下快速删除字符串