【因果推断python】16_工具变量2

目录

出生季度和教育对工资的影响

第一阶段


出生季度和教育对工资的影响

到目前为止,我们一直将这些工具视为一些神奇的变量 Z,它们具有仅通过干预变量影响结果的神奇特性。老实说,好的工具变量来之不易,我们不妨将它们视为奇迹。让我们说它不适合胆小的人。有传言说,芝加哥经济学院的酷孩子们谈论他们是如何在酒吧里想出这种或那种工具变量的。

不过,我们确实有一些有趣的工具示例,可以让事情变得更具体一些。我们将再次尝试估计教育对工资的影响。为此,我们将使用该人的出生季度作为工具 Z

这个想法利用了美国强制出勤法。通常,他们声明孩子必须在他们入学当年的 1 月 1 日之前满 6 岁。因此,年初出生的孩子入学年龄较大。强制出勤法还要求学生在学校上学直到他们年满 16 岁,届时他们在法律上被允许退学。结果是,与年初出生的人相比,年末出生的人平均受教育年限更长。如果我们接受出生季度与能力因素无关,即它不会混淆教育对工资的影响,我们可以将其用作工具。换句话说,我们需要相信出生季度对工资没有影响,除了对教育的影响。如果你不相信占星术,这是一个非常有说服力的论点。

g = gr.Digraph()g.edge("能力", "教育")
g.edge("能力", "工资")
g.edge("教育", "工资")
g.edge("出生季度", "教育")
g

为了进行这种分析,我们可以使用来自三次十年一次的人口普查的数据,这些数据与Angrist和Krueger 在他们关于IV的文章中使用的数据相同。该数据集包含有关我们的结果变量,即工资取对数的结果,以及我们的干预变量,即受教育年限的信息。它还包含我们的工具变量,即出生季度,以及其他控制变量的数据,例如出生年份和出生状态。

data = pd.read_csv("./data/ak91.csv")
data.head()

第一阶段

在我们使用出生季度作为工具变量之前,我们需要确保它是有效的。这意味着我们支持工具变量的两个假设:

  1. Cov(Z\text{,}T)\neq0。这是说我们应该有一个强大的第一阶段,或者工具变量确实会影响干预变量。
  2. Y\perp Z|T.这是排除限制,声明工具变量Z仅通过干预T影响结果Y。

幸运的是,第一个假设是可以验证的。我们从数据中可以看出Cov(Z,T)不为零。在我们的例子中,如果出生的季度确实是一个工具变量,就像我们所说的那样,我们应该期望在一年中最后一个季度出生的人比年初出生的人有更多的受教育时间。在运行任何统计测试来验证这一点之前,让我们绘制我们的数据并亲眼看到它。

group_data = (data.groupby(["year_of_birth", "quarter_of_birth"])[["log_wage", "years_of_schooling"]].mean().reset_index().assign(time_of_birth = lambda d: d["year_of_birth"] + (d["quarter_of_birth"])/4))plt.figure(figsize=(15,6))
plt.plot(group_data["time_of_birth"], group_data["years_of_schooling"], zorder=-1)
for q in range(1, 5):x = group_data.query(f"quarter_of_birth=={q}")["time_of_birth"]y = group_data.query(f"quarter_of_birth=={q}")["years_of_schooling"]plt.scatter(x, y, marker="s", s=200, c=f"C{q}")plt.scatter(x, y, marker=f"${q}$", s=100, c=f"white")plt.title("Years of Education by Quarter of Birth (first stage)")
plt.xlabel("Year of Birth")
plt.ylabel("Years of Schooling");

值得注意的是,在一年中的四分之一之后,学校教育的年份有一个季节性的模式。从图形上我们可以看到,一年中第一季度出生的人的受教育程度几乎总是低于最后一个季度出生的人(毕竟,一旦我们控制了出生年份,那些晚年出生的人通常受教育程度更高)。

为了更严格一点,我们可以将第一阶段作为线性回归运行。我们首先将出生季度转换为虚拟变量:

factor_data = data.assign(**{f"q{int(q)}": (data["quarter_of_birth"] == q).astype(int)for q in data["quarter_of_birth"].unique()})factor_data.head()

为简单起见,现在只使用最后一个季度,即Q4,作为工具。我们将使用干预变量:受教育年限,对工具变量:出生季度,进行回归。这将向我们展示出生在哪个季度是否确实像我们在上图中看到的那样对教育时间产生了积极影响。我们还需要在这里控制出生年份,我们将添加出生状态作为额外的控制。

first_stage = smf.ols("years_of_schooling ~ C(year_of_birth) + C(state_of_birth) + q4", data=factor_data).fit()print("q4 parameter estimate:, ", first_stage.params["q4"])
print("q4 p-value:, ", first_stage.pvalues["q4"])q4 parameter estimate:,  0.10085809272785906
q4 p-value:,  5.464829416638474e-15

看起来,在一年中最后一个季度出生的人平均比在一年中其他季度出生的人多受教育0.1年。p 值接近于零。这结束了关于出生在哪个季度是否导致更多或更少受教育年限的案例。

本文来自互联网用户投稿,该文观点仅代表作者本人,不代表本站立场。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如若转载,请注明出处:http://www.mzph.cn/diannao/23089.shtml

如若内容造成侵权/违法违规/事实不符,请联系多彩编程网进行投诉反馈email:809451989@qq.com,一经查实,立即删除!

相关文章

Windows11系统 和Android 调试桥(Android Debug Bridge,ADB)工具安装,app抓取日志内容

文章目录 目录 文章目录 安装流程 小结 概要安装流程技术细节小结 概要 Android调试桥(ADB)是一种多功能命令行工具,它允许开发者与连接到计算机上的Android设备进行通信和控制。ADB工具的作用包括但不限于: 安装和卸载应用程序&…

2024年先进材料与清洁能源国际会议(ICAMCE 2024)

2024 International Conference on Advanced Materials and Clean Energy 【1】大会信息 大会时间:2024-07-15 大会地点:中国三亚 截稿时间:2024-07-01(以官网为准) 审稿通知:投稿后2-3日内通知 投稿邮箱&#xff1a…

揭秘相似矩阵:机器学习算法中的隐形“纽带”

在机器学习领域,数据的处理和分析至关重要。如何有效地从复杂的数据集中提取有价值的信息,是每一个机器学习研究者都在努力探索的问题。相似矩阵,作为衡量数据之间相似性的数学工具,在机器学习算法中扮演着不可或缺的角色。 相似矩…

Docker之路(三)docker安装nginx实现对springboot项目的负载均衡

Docker之路(三)dockernginxspringboot负载均衡 前言:一、安装docker二、安装nginx三、准备好我们的springboot项目四、将springboot项目分别build成docker镜像五、配置nginx并且启动六、nginx的负载均衡策略七、nginx的常用属性八、总结 前言…

SpringBoot高手之路02-全局异常处理器

RestControllerAdvice 可以将响应数据返回json格式然后响应 那么开始做全局异常处理器 首先先定义一个类 package com.healer.exception;import com.healer.common.Result; import org.springframework.web.bind.annotation.ExceptionHandler; import org.springframework.we…

高并发系统限流原理

短时间内巨大的访问流量,我们如何让系统在处理高并发的同时还能保证自身系统的稳定性?估计有人会说,增加机器就可以了,因为我的系统架构设计就是按照分布式思想进行架构设计的,所以可以只需要增加机器就可以解决问题了…

在Unity中配置Android项目以允许HTTP流量,解决AVPro在Android平台中无法播放http视频

解决方法快速通道:拉到底,看倒数第二张图 好记性不如烂笔头 最近在使用AVpro插件播放http视频,在Editor中一切正常,然而打包在Android平台下就播放不了 AVPro在Unity中的警告: 感觉只是个警告,没引起注意…

实验报告题目

1、Anaconda环境配置 主要包括环境介绍,安装,以及numpy,matplotlib,pandas包的使用 2、线性回归模型 线性回归模型的基本原理 线性模型代码编写 结果分析 3、线性判别模型 线性判别模型的基本原理 线性模型代码编写 结果分析 4、决策树模型 决策树模型…

为什么3D渲染有的模型是黑的?---模大狮模型网

在3D建模和渲染的过程中,遇到模型渲染出来是黑色的情况并不罕见。这可能是由于多种原因导致的,包括但不限于材质设置问题、灯光配置不当、渲染设置错误等。 一、材质设置问题 材质丢失或错误:如果模型在导入或创建时材质没有正确加载或设置&…

Qt QML学习资料链接

文章目录 1. Qt QML2. 学习和使用Qt QML目的3. Qt QML常用的一些模块举例4. 参考链接 1. Qt QML Qt QML模块为基于QML language开发应用程序和库提供了语言和基础建设。它定义并实现了QML的语法和基础框架,并提供API使得应用程序开发者可以使用定制的类型扩展QML&a…

TypeScript 中的 Mixins

1. Mixins 的概念 Mixins 是一种将多个类的功能组合到单个类中的技术。它允许一个类通过混合(Mix)其他类的功能来扩展自己的功能,从而实现代码的复用和组合。Mixins 可以用于解决传统继承带来的多重继承问题,并且更加灵活和可扩展…

webservice、WCF、webAPI、MVC权限认证

webservice 权限认证 》》soapHeader SOAPHeader案例 服务引用下生成的服务方法参数中会自动加入一个soapHeader的参数, WEB服务引用则没有,我感觉采用WEB服务引用基于这种验证比较方便, 因为只需将soapHeader实例赋值一次就可以多次调用不…

IDC数据中心内容合规检测

摘要:电信数据中心通过实施AI数据质检方案,自动检测反馈预警,提升服务质量,保障网络合规,降低风险,增强企业竞争力。 1. 客户背景 该客户是国内领先的IDC云计算综合服务提供商,专注于为客户提供…

Vue3中的常见组件通信之v-model

Vue3中的常见组件通信之v-model 概述 ​ 在vue3中常见的组件通信有props、mitt、v-model、 r e f s 、 refs、 refs、parent、provide、inject、pinia、slot等。不同的组件关系用不同的传递方式。常见的撘配形式如下表所示。 组件关系传递方式父传子1. props2. v-model3. $r…

jmap 和 jstack 的线上使用及操作过程示例

jmap 和 jstack 的线上使用及操作过程示例 一、jmap 1. 工具简介 jmap是JDK提供的一个命令行工具&#xff0c;主要用于生成Java堆的转储快照&#xff08;dump文件&#xff09;以及查看Java进程中的内存使用情况。 2. 命令语法 jmap [option] <pid>其中&#xff0c;o…

.NET周刊【6月第1期 2024-06-02】

国内文章 一文带你了解.NET能做什么&#xff1f; https://www.cnblogs.com/Can-daydayup/p/18214473 .NET是一个免费、开源、跨平台的开发平台框架&#xff0c;广泛应用于桌面、Web、移动、云服务、游戏、物联网、大数据和人工智能等领域开发。它支持C#、Visual Basic、F#等…

Redis 单线程问题 BigKey问题

前言 简单的redis基础类型以及常用操作我们都也已经介绍过了 现在今天我们来谈谈redis对于单线程是需要怎么理解的 以及redis假设遇见大key我们需要怎么去查询和删除呢??? redis单线程 假设有个人现在问你一个问题:redis是单线程的还是多线程的 这个问题本身就不严谨 就像问…

I.MX6ULL 串口格式化函数移植实验

系列文章目录 I.MX6ULL高精度延时实验 I.MX6ULL高精度延时实验 系列文章目录一、前言二、串口格式化函数简介三、硬件原理分析四、实验程序编写五、编译下载验证 一、前言 上一节实验实现了 UART1 基本的数据收发功能&#xff0c;虽然可以用来调试程序&#xff0c;但是功能太单…

vxe-table表格父级元素根据需求设置宽高

1. 问题 vxe-table表格父级元素未设置宽高&#xff0c;导致初始进入时表格内数据一列一列向下呈现动态递增&#xff0c;或者无表格数据时表格显示异常。 2.vxe-table 表格 vxe-table 是一个基于 Vue.js 的表格组件&#xff0c;它提供了丰富的功能和灵活的配置选项。在使用 vxe…

蒙层(css)

如何在 Vue 中实现一个包含图像和蒙层效果的组件&#xff1f;这个组件根据某个条件显示或隐藏蒙层&#xff0c;用于表示图像是否已读。 1. 创建基础模板 首先&#xff0c;我们在模板中使用 div 包裹我们的图像组件 GraphImage&#xff0c;并为最外层 div 设置 position: relat…