Hern\(\'{a}\)n M. and Robins J. Causal Inference: What If.

这一章主要讨论的是, 观测得到的数据(而非随机实验)在什么条件下可以视为是随机试验.

outcome predictors: 一些会导致\(Y\)发生的诱因

3.1

  1. 我们所考虑的\(A\)和实验中实际的采取的手段\(A\)是相一致的.

  2. 采取何种手段\(A\)仅仅与\(L\)有关(这里考虑, \(L, A, Y\)三个元素).

  3. \(\mathrm{Pr}(A|L) > 0\), 即正定性.

下面是一点一点的分析这三个点的重要性.

3.2 Exchangeability

这个对应的是第二点, 即我们要探究是否\(A\)仅仅与\(L\)有关, 从而有可交换性:

\[Y^a \amalg A |L.
\]

一旦遇到上面的情况, 往往就没有上述可交换性的保证了.

3.3 Positivity

设想\(L\)代表的是一个人是否吸烟, 倘若一个医生仅仅给不吸烟的人进行心脏迁移手术, 即

\[\mathrm{Pr}[A=1|L=1] = 0,
\]

则我们就完全丢失了这部分信息, 自然也没办法计算casual effect, 因为

\[\mathrm{Pr}[Y|A=1, L=1]
\]

压根没有定义.

3.4 Consistency

一致性分类预期结果的一致性, 以及结果和观测数据的一致性

First

现在假设\(A \in \{0, 1\}\), 即代表是否进行心脏移植手术, 但是在实际中, \(A\)并非如此纯粹的0, 1.

实际上, 取决于器材, 外科医生的差别会衍生出不同版本的\(A\).

当然了, 这么讨论下去只会导致不可知论, 我们可以在某种程度上假设, 不过对\(A\)的描述越细致, 即越细分, 最后的结论也会更加精准.

Second

这个一致性, 用公式就是

\[Y^a = Y, A=a,
\]

这个很重要, 因为我们在计算causal effect的时候有这么一步

\[\mathrm{Pr}[Y|A=a, L] = \mathrm{Pr}[Y^a|A=a, L].
\]

这个一致性, 个人的理解是, 我们所观察的\(A=a\)有很多版本, 可能与我们所希望的\(Y^a\)并不一致, 导致\(Y^a \not = Y\).

这里有一个微妙的东西, 实在是不知道如何描述了.

Fine Point

3.1 Identifiability of causal effects

指, 倘若不是随机实验, 我们需要一些额外的假设来得以计算causal effect.

3.2 Crossover randomized experiments

p32

这个讨论的是在不同的时间点\(t=0, t=1\).

3.3 Possible worlds

p35

3.4 Attributable fraction

p38

Technical Point

3.1 Positivity for standardization and IP weighting

p32

上一章讲了利用standardization 和 IP weighting 在条件可交换的假定下, 我们可以计算causal effect.

但是, 实际上这同时是需要positivity的假定的.

standardization:

\[\sum_l \mathbb{E} [Y|A=a, L=l] \mathrm{Pr} [L=l],
\]

这个式子需要\(\mathbb{E}[Y|A=a, L=l]\), 但是这个在某些\(P(A=a|L=l)=0\)的情况下是没有定义的.

另一方面, IP weighting

\[\mathbb{E} [\frac{I(A=a)Y}{f(A|L)}] = \mathrm{Pr}[L \in Q(a)]\sum_{l} \mathbb{E} [Y|A=a, L=l, L\in Q(a)] \mathrm{Pr} [L=l|L \in Q(a)],
\]

其中\(Q(a) = \{l; \mathrm{Pr} (A=a|L=l)>0\}\).

相当于, 认为地目标的集合缩小了.

里头还说, 上述的与

\[\mathbb{E} [\frac{I(A=a)Y}{f(a|L)}]
\]

不同, 而且说后者是undefined的, 可是后决定后者才是等价于上面所说的啊.

不过我倒是觉得无所谓的, 毕竟我们应该关心我们所关心的, 限定在\(f(a|L)\not = 0\)才是合适的区域.

3.2 Cheating consistency

p40

Chapter 3 Observational Studies的更多相关文章

  1. descriptive statistics|inferential statistics|Observational Studies| Designed Experiments

    descriptive statistics:组织和总结信息,为自身(可以是population也可以是sample)审视和探索, inferential statistics.从sample中推论p ...

  2. Weighted Effect Coding: Dummy coding when size matters

    If your regression model contains a categorical predictor variable, you commonly test the significan ...

  3. hbase官方文档(转)

    FROM:http://www.just4e.com/hbase.html Apache HBase™ 参考指南  HBase 官方文档中文版 Copyright © 2012 Apache Soft ...

  4. HBase官方文档

    HBase官方文档 目录 序 1. 入门 1.1. 介绍 1.2. 快速开始 2. Apache HBase (TM)配置 2.1. 基础条件 2.2. HBase 运行模式: 独立和分布式 2.3. ...

  5. 【统计】Causal Inference

    [统计]Causal Inference 原文传送门 http://www.stat.cmu.edu/~larry/=sml/Causation.pdf 过程 一.Prediction 和 causa ...

  6. Propensity Scores

    目录 基本的概念 重要的结果 应用 Propensity Score Matching Stratification on the Propensity Score Inverse Probabili ...

  7. R数据分析:样本量计算的底层逻辑与实操,pwr包

    样本量问题真的是好多人的老大难,是很多同学科研入门第一个拦路虎,今天给本科同学改大创标书又遇到这个问题,我想想不止是本科生对这个问题不会,很多同学从上研究生到最后脱离科研估计也没能把这个问题弄得很明白 ...

  8. Modern C++ CHAPTER 2(读书笔记)

    CHAPTER 2 Recipe 2-1. Initializing Variables Recipe 2-2. Initializing Objects with Initializer Lists ...

  9. Android Programming: Pushing the Limits -- Chapter 7:Android IPC -- ApiWrapper

    前面两片文章讲解了通过AIDL和Messenger两种方式实现Android IPC.而本文所讲的并不是第三种IPC方式,而是对前面两种方式进行封装,这样我们就不用直接把Aidl文件,java文件拷贝 ...

随机推荐

  1. C++构造函数和析构函数初步认识(2)

    构造函数的三个作用1.构造对象2.对象初始化3.类型转换  //Test1.h #include<iostream> using namespace std; //构造对象 //初始化对象 ...

  2. 【STM32】基于正点原子『探索者』开发板的烧录

    项目需要一个功能,开发板范例正好有,就买了一块,不过还是有点贵 我手边没有J-Link 用的都是串口烧录 烧录时,先打开右上的开关 如果是仿真器烧录,它无法供电,需要接12V适配器或是杜邦线供电 然后 ...

  3. 【Linux】【Basis】块存储,文件存储,对象存储

    1. 块存储: 定义:这种接口通常以QEMU Driver或者Kernel Module的方式存在,这种接口需要实现Linux的Block Device的接口或者QEMU提供的Block Driver ...

  4. Mave 下载与安装

    一,Maven 介绍 我们在开发中经常需要依赖第三方的包,包与包之间存在依赖关系,版本间还有兼容性问题,有时还需要将旧的包升级或降级,当项目复杂到一定程度时包管理变得非常重要.Maven是当前最受欢迎 ...

  5. 【C/C++】子数组的最大累加和问题

    #include <bits/stdc++.h> using namespace std; class Solution { public: /** * max sum of the su ...

  6. kubernetes list/watch设计原理

    overview kubernetes的设计里面大致上分为3部分: API驱动型的特点 (API-driven) 控制循环(control loops)与 条件触发 (Level Trigger) A ...

  7. Spring Cloud Eureka源码分析之服务注册的流程与数据存储设计!

    Spring Cloud是一个生态,它提供了一套标准,这套标准可以通过不同的组件来实现,其中就包含服务注册/发现.熔断.负载均衡等,在spring-cloud-common这个包中,org.sprin ...

  8. Springboot MVC 自动配置

    Springboot MVC 自动配置 官方文档阅读 https://docs.spring.io/spring-boot/docs/current/reference/html/web.html#w ...

  9. AT266 迷子のCDケース 题解

    Content 有 \(n+1\) 个碟,编号为 \(0\sim n\),一开始 \(0\) 号碟在播放机上,其他的碟依次放进了 \(n\) 个盒子里面.现在有 \(m\) 次操作,每次操作找到当前在 ...

  10. 【RTOS】FreeRTOS中的任务堆栈溢出检测机制

    目录 前言 任务堆栈 堆栈溢出 任务堆栈溢出检测机制 API 两种堆栈溢出检测方式 堆栈溢出钩子函数 内核何时检测任务堆栈溢出 任务堆栈溢出检测存在的局限性 前言 注意:本笔记发布时可能忘记补充查看d ...