因果推论新思维反事实分析架构.ppt

《因果推论新思维反事实分析架构.ppt》由会员分享，可在线阅读，更多相关《因果推论新思维反事实分析架构.ppt（103页珍藏版）》请在三一办公上搜索。

1、因果推論新思維:反事實分析架構,關秉寅政治大學社會學系2010.01.16,A New Paradigm for Causal Inference:The Counterfactual Framework,06.09.06,1,因果推論：從古到今,參考Judea Pearl（professor of Computer Science and Statistics and director of the Cognitive Systems Laboratory,UCLA）http:/singapore.cs.ucla.edu/LECTURE/lecture_sec1.htm,06.09.06,1

2、,因果推論：從古到今,從亞當與夏娃說起：亞當與夏娃吃過智慧樹的果子後，他們已經是因果解釋的專家。當上帝問亞當你是否吃了智慧樹的果子時（詢問事實），亞當的回答是：祢給我的那個女人，從樹上拿了果子讓我吃了（不只是說明事實，還做了解釋）；上帝問夏娃，夏娃回答是：我是吃了，可是是蛇誘惑和欺騙了我。,06.09.06,1,因果推論：從古到今,聖經故事的意涵：因果解釋是一個man-made concept。因果解釋是拿來歸咎責任的。只有神（為了某些目的）、人與動物（有自有意志）可以讓事情發生，而不是東西（objects）或物理的過程（physical processes）。,06.09.06,1,因果推論

3、：從古到今,但當工程師開始建造有許多滑輪及繩纜的 system來幫人做事後，physical objects開始有了因果的特性。當系統無法運作時，歸咎於神或人並沒有用，比較有效的解釋是滑輪或繩纜壞了，只要加以替換後，系統就可再運作。至此，causes之概念有雙重意義：the targets of credit and blame the carriers of physical flow of control on the other,06.09.06,1,因果推論：從古到今,至文藝復興時，當上帝做為final cause逐漸被人的科學知識所取代時，這雙重意義遇到困難與挑戰。Galileo在1

4、638年出版 Discorsi（兩門新科學的對話）這本書後，有了革命性的轉變。此書有兩項重要的格言（maxims）：,06.09.06,1,因果推論：從古到今,1.先描述，後解釋（Description first,explanation second）:The how precedes the why.Ask not,said Galileo,whether an object falls because it is pulled from below or pushed from above.Ask how well you can predict the time it takes fo

5、r the object to travel a certain distance,and how that time will vary from object to object,and as the angle of the track changes.2.以數學（方程式）來描述，而不是語言：如d=t2。,06.09.06,1,因果推論：從古到今,從此，物理學充滿了有用但未被解釋的經驗法則，如Snell law,Hookes law,Ohms law,Joules law。另一項比預測實驗結果更重要的是代數方程式的運用，因為工程師除了可以問“how to”外，還可以問“what if”。

6、,06.09.06,1,因果推論：從古到今,至啟蒙時代，David Hume 將 Galileo的第一項格言發揮至極致，他認為 the WHY is not merely second to the HOW,but that the WHY is totally superfluous as it is subsumed by the HOW。On page 156 of Treatise of Human Nature:Thus we remember to have seen that species of object we call FLAME,and to have felt tha

7、t species of sensation we call HEAT.We likewise call to mind their constant conjunction in all past instances.Without any farther ceremony,we call the one CAUSE and the other EFFECT,and infer the existence of the one from that of the other.,06.09.06,1,因果推論：從古到今,從實證主義的角度來看，Hume是在說“A caused B”與“Whenev

8、er A occurs,then B does”相同。但是我們都知道雞啼是伴隨著朝陽，但雞啼並不是讓太陽出來的原因。大難題1：如果如Hume所說，我們的知識是來自於經驗，而經驗是以相關的形式存於我們的心靈中，那我們如何得到因果方面的知識？If regularity of succession is not sufficient;what WOULD be sufficient?,06.09.06,1,因果推論：從古到今,大難題2：我們知道某些關連有或沒有因果間的關係，有何差別嗎？知道因果又有何用？當然知道因果關係，就可以做某些事。如果雞啼是造成太陽升起的原因，那我們要縮短夜晚的話，就可以讓雞早

9、些啼。If causal information has an empirical meaning beyond regularity of succession,then that information should show up in the laws of physics.But it does not!,06.09.06,1,因果推論：從古到今,Russell（1913）認為“All philosophers imagine that causation is one of the fundamental axioms of science,yet oddly enough,in

10、advanced sciences,the word cause never occurs.The law of causality,I believe,is a relic of bygone age,surviving,like the monarchy,only because it is erroneously supposed to do no harm.”“It could not possibly be an abbreviation,because the laws of physics are all symmetrical,going both ways,while cau

11、sal relations are uni-directional,going from cause to effect.”,06.09.06,1,因果推論：從古到今,但另一位科學哲學家Patrick Suppes 則指出“There is scarcely an issue of Physical Review that does not contain at least one article using either cause or causality in its title.”物理學家一方面寫無因果意涵的方程式，但另一方面卻大談因果關係。統計學一百多年前發現相關（correlati

12、on）的概念後，卻無法忽略相關與因果間的區辨。,06.09.06,1,因果推論：從古到今,Francis Galton於1888年進行個人的前臂與其頭大小關係的測量，企圖瞭解一個數值預測另一個數值的程度時，發現到：If you plot one quantity against the other and scale the two axes properly,then the slope of the best-fit line has some nice mathematical properties:The slope is 1 only when one quantity can pr

13、edict the other precisely;it is zero whenever the prediction is no better than a random guess and,most remarkably,the slope is the same no matter if you plot X against Y or Y against X.我們開始可以根據資料客觀的測量兩個變項間的關係，而不是根據我們的意見或判斷。,06.09.06,1,因果推論：從古到今,Galton的發現震撼了其學生Karl Pearson（公認是現代統計學之父），使他終其一生認為我們只需要相關

14、這個更廣泛的概念，而不需要另一個獨立的因果關係的概念。他也從不在其論文內提到因果的概念。這一直要到Sir Ronald Fisher 建立 randomized experiment 的研究設計後，才成為唯一被主流統計學所接受，並認為這是唯一可從資料驗證因果間關係的科學方法。但這種謹慎的看法使得無法做實驗，而需靠統計分析為引導的社會科學處於近乎癱瘓的狀態。,06.09.06,1,因果推論的新典範,Pearl 認為這樣的困境是源自統計學的官方語言：機率的語言。因為 cause 並不是機率的字彙。我們無法以機率的語言表達：Mud does not cause rain。我們只能說兩者相關。Natu

15、rally,if we lack a language to express a certain concept explicitly,we cant expect to develop scientific activity around that concept.Scientific development requires that knowledge be transferred reliably from one study to another and,as Galileo has shown 350 years ago,such transference requires the

16、 precision and computational benefits of a formal language.,06.09.06,1,因果推論的新典範,當研究者開始企圖使用電腦來建立因果關係時，對於先前兩大難題有了新的想法。從概念層次來說，機器人所遇到的問題是與經濟學家企圖建立稅收及預算模式，或流行病學家企圖建立流行病的模式是一樣的。不論是機器人、經濟學家或流行病學家都需要在行動受限制的條件下，依據充滿雜音的資訊，從環境中找出因果的關係。這就是先前的第一個大難題:HOW?,06.09.06,1,因果推論的新典範,機器人的世界也與第二大難題有關。如果我們教導機器人關於因果的知識，機器人要

17、如何組織及運用這些知識呢？對機器人而言，這兩大難題是具體而實際的問題，就是如何從與環境的互動中得到因果的知識，以及如何運用從創造者/程式設計師得到的因果知識？,06.09.06,1,新典範與新語言,J.Pearl 的答案：第二個難題可以結合graphs與equations的方式解決，如此則第一個難題也比較容易解決。解題的主要關鍵概念是：(1)treating causation as a summary of behavior under interventions.(2)using equations and graphs as a mathematical language within

18、which causal thoughts can be represented and manipulated.(3)Treating interventions as a surgery over equations.,06.09.06,1,新典範與新語言,社會科學家過去75年來不時使用graphs，如structural equations modeling 及 path diagrams。但是由於代數方程式的便利性，因而壓抑了圖型的使用，也進而失去了圖型帶來的好處。這些diagrams事實上捕捉了因果的本質預測不正常情況或是新的操弄狀態下，會產生什麼結果。這種預測是代數或相關分析無法做

19、到的。,06.09.06,1,新典範與新語言,06.09.06,1,新典範與新語言,從這樣的角度來看因果關係，可以理解為何科學家如此熱衷於因果解釋，因為建立因果模式會得到一種“deep understanding”及“being in control”的感覺。Deep understanding的意思是“knowing,not merely how things behaved yesterday,but also how things will behave under new hypothetical circumstances,control being one such circums

20、tance”.,06.09.06,1,新典範與新語言,即使我們無法實際上控制事情，我們也因理解因果關係而得到“in control”的感覺。例如，我們無法控制星體的運轉，但萬有引力的理論，讓我們能夠理解並獲得假設性控制（hypothetical control）的感覺。我們也可以預測當萬有引力改變時，對潮汐會產生什麼影響。因果模式也是做為區辨有意識的論證（deliberate reasoning）及被動或本能的反應（reactive or instinctive response）的試金石。前者可在即使不真正去從事新的操弄（manipulation）下，也可以預期新的操弄可以得到什麼結果。,0

21、6.09.06,1,新典範與新語言,Equations vs.Diagrams(J.Pearl),06.09.06,1,新典範與新語言,Definition of Causation:Y is a cause of Z if we can change Z by manipulating Y,namely,if after surgically removing the equation for Y,the solution for Z will depend on the new value we substitute for Y.THE DIAGRAM TELLS US WHICH EQU

22、ATION IS TO BE DELETED WHEN WE MANIPULATE Y.INTERVENTION AMOUNTS TO A SURGERY ON EQUATIONS,GUIDED BY A DIAGRAM,AND CAUSATION MEANS PREDICTING THE CONSEQUENCES OF SUCH A SURGERY.,06.09.06,1,新典範與新語言,INTERVENTION AS SURGERY-CONTROLLED EXPERIMENTS 假定我們要研究某種藥物是否能幫助病人從某些疾病復原。但影響復原的因素不只是藥物治療，還有其他的因素，如社經地位、

23、年齡、生活方式等。這對我們判斷藥物的效果來說是個問題，因為我們不知道影響復原的程度是那個因素造成的。因此，我們希望能夠比較背景相同的病人，而這就是Fisher的Randomized Experiment能夠做到的。但隨機分派的實驗設計如何做到這種比較呢？,06.09.06,1,新典範與新語言,這種實驗設計實際上有兩個部份：Randomization and Intervention.Intervention 就是我們將藥物給一些在正常情況下不會需要這種治療的病人，另一方面我們給一些尋求治療的病人安慰劑（placebo）。這也是先前提及的SURGERY 的概念，因為我們將一個functional

24、 link（如SES）切斷，而用另一種替代。Fisher的偉大貢獻是 connecting the new link to a random coin flip，以此保證我們想要切斷的 link 是確實被切斷了，因為我們可以假定這個 random coin 是不受任何我們可以測量到的因素所影響的。,06.09.06,1,新典範與新語言,06.09.06,1,新典範與新思維,在新典範下，我們可以如何在無法從事實驗的情況下，利用observational studies（如調查或病歷）的資料思考因果關係呢？,06.09.06,1,發問方式,兩類探索因果關係的研究問題第一類研究問題：X 對 Y 的影

25、響為何？研究目的是要知道：X 對 Y 有影響嗎？有的話，影響有多大？這是實驗設計所問的研究問題第二類研究問題：影響 Y 的因素有哪些？研究目的是要知道：所有可能解釋 Y 的因素有哪些。,06.09.06,1,如何認定因果關係?,參考 Morgan,Stephen L.&Christopher Winship(2007).Counterfactuals and Causal Inference:Methods and Principles for Social Research.New York,NY:Cambridge University Press.,06.09.06,1,如何認定因果關係

26、?,X Y(Had X taken a different value,then Y would have taken a different value)The causal relationship between X and Y is confounded if:Z X Y,06.09.06,1,如何認定因果關係?,Z A X YZ A(unobserved)X Y,06.09.06,1,如何認定因果關係?,X Y C(collider)Z A(unobserved)X Y,06.09.06,1,Statistical Relations vs.Causal Relations,Stat

27、istical dependence may reflectRandom fluctuation(c.i.collider bias),06.09.06,1,如何推估 the effect of X on Y?,當X與Y有common causes時，如能認定X與Y間的因果路徑，可以因有效的控制一組共變項 Z，且 Z 內沒有受X影響時，則我們可以清楚的推估 the effect of X on Y。Pearls Back-door Criterion,06.09.06,1,如何推估the effect of D(treatment)on Y?,V(unobserved)A F GU(unobs

28、erved)B D Y C,06.09.06,1,Pearls Back-door Criterion,If one or more back-door paths connects the causal variable to the outcome variable,Pearl shows that the causal effect is identified by conditioning on a set of variables Z if and only if all back-door paths between the causal variable and the outc

29、ome variable are blocked after conditioning on Z.,06.09.06,1,Pearls Back-door Criterion,A back-door path of D and Y is blocked by Z if and only if the back-door path satisfies any one of the following:contains a chain of mediation A Z B,or contains a fork of mutual dependence A Z B;contains an inver

30、ted fork of mutual causation A C*B,where C*and all its descendants are not in Z.,06.09.06,1,如何推估the effect of D on Y?(控制 B&F 即可，Why?）,V(unobserved)A F GU(unobserved)B D Y C,06.09.06,1,Pearls Back-door Criterion(continued),從 Pearl 的 Back-door Criterion 來看，並不是控制越多變項就好，因為要是控制了 colliders 反而會有問題，因為本來X與沒有

31、相關或因果關係的，控制這類變項後，反而會產生相關。Example如果collider是申請入學時是否被一所菁英學校接受（adm:1-接受；0-拒絕）是否被接受是根據兩個獨立變項：SAT及面試時對動機的評估(Motivation)因此：adm 是的兩個 causes 是 SAT 及 Motivation，而 SAT 及 Motivation 間是獨立的。,06.09.06,1,Example of controlling a collider,06.09.06,1,Example of controlling a collider,06.09.06,1,Example of controllin

32、g a collider,06.09.06,1,如何推估the effect of D on Y?,A F G B D Y C H,06.09.06,1,如何有效的控制 confounding variables(Z)呢？,如果D只有兩個值：0及1，我們可以用實驗設計的方式，將兩組相同特性的人，隨機分配到控制組（D0）及實驗組（D1），即可達到前述的狀態。Why?如果無法用實驗設計，而是用調查方法蒐集資料時，如何推估呢？條件性控制（conditioning）或是配對（matching）：by holding constant or by balancing/homogenizing the t

33、reatment&control groups.,06.09.06,1,The Counterfactual Framework,反事實因果推論的想像,06.09.06,1,The Counterfactual Framework,反事實分析架構的想像可看成是一種thought experiment。要想像的是同一個個體或群體在不同的狀態下，會有什麼可能的結果（potential outcomes）。這些可能結果間的差異，即為不同狀態（因）的效果。Counterfactuals should be reasonable!,06.09.06,1,The Counterfactual Framew

34、ork,Q：什麼是unreasonable 的 counterfactuals 呢？有什麼狀態不適合看成為 causes 的嗎？有什麼樣的結果不適合想像counterfactual情況的嗎？,06.09.06,1,The Counterfactual Framework,個人層次的真正因果效應：i=Yi1 Yi0The Fundamental Problem of Causal Inference：無法觀察同一個人同時在實驗組及控制組。加上一些假定，如 SUTVA，則可推估群層次的因果效應。,06.09.06,1,The Counterfactual Framework,SUTVA：The S

35、table Unit Treatment Value Assumption a priori assumption that the value of Y for unit u when exposed to treatment t will be the same no matter what mechanism is used to assign treatment t to unit u and no matter what treatments the other units receive.,06.09.06,1,The Counterfactual Framework,實驗設計是假

36、設我們能夠將觀察到的替代無法觀察到的。如隨機分派到實驗組與控制組的個體的特性相同，則我們可以假定：如果實驗組的個人沒有接受treatment的話，其結果與控制組觀察到的相同；如果控制組的個人接受treatment的話，其結果與實驗組觀察到的相同。,06.09.06,1,The Counterfactual Framework,當使用調查方法得到資料時，即observational data，個人為何會接受或不接受treatment，往往不是一個隨機的現象。Observational data通常有兩個問題：接受treatment者與不接受者有baseline differences，以及het

37、erogeneity of treatment effect.可能有些影響接受treatment與否的變項，並未觀察到，亦即omitted variables的問題。,06.09.06,1,The Counterfactual Framework,06.09.06,1,The Counterfactual Framework,以反事實架構的觀點來看，母群體層次的真正因果效應（ATE；Average Treatment Effect）為：E=EY1 Y0=EY1 EY0=EY1|D=1+(1)EY1|D=0 EY0|D=1+(1)EY0|D=0=EY1|D=1 EY0|D=1+(1)EY1|D=

38、0 EY0|D=0=E|D=1+(1)E|D=0,06.09.06,1,The Counterfactual Framework,：母群體中接受 treatment 的比例不同組的因果效應：ATT（Average Treatment Effect on the Treated）:EY1|D=1 EY0|D=1，即 E|D=1 ATU（Average Treatment Effect on the Untreated）:EY1|D=0 EY0|D=0，即 E|D=0Q：我們可以假定不同組的人有同樣的因果效應嗎？如果我們能夠做此假定的話，則 EY1|D=0=EY1|D=1,06.09.06,1,T

39、he Counterfactual Framework,基準線的差異：EY0|D=1 EY0|D=0Q：我們可以假定不同組的人在未接受 treatment 前是一樣的嗎？如果我們能夠做此假定的話，則 EY0|D=1=EY0|D=0,06.09.06,1,The Counterfactual Framework,如果我們只以觀察到接受 treatment 的組與觀察到未接受 treatment 的組之間的差異做為 Causal Effect 的估計時，此估計是一種 Nave Estimate：Nave Estimate EY1|D=1 EY0|D=0=average causal effect=

40、E()+baseline bias+E(Y0|D=1)E(Y0|D=0)+differential effect bias+E(|D=1)E(|D=0)(1),06.09.06,1,The Counterfactual Framework:A Review,反事實分析架構的五個關鍵概念：Potential/Hypothetical States&Outcomes:因果效應（causal effect）是利用“potential”或“hypothetical”的概念，而不是只用到 actual observations。.The ceteris paribus condition其他條件相同的條

41、件下，也就是將其他因素控制成等同（equal）、固定不變（fixed）或是constant。,06.09.06,1,The Counterfactual Framework:A Review,Heterogeneity:個人對於treatment的反應是因人而異的。亦即因果效應在個人層次即被認定是有差異的。每個人的因果效應是：potential outcome under the potential treatment state potential outcome under the potential control stateFundamental Problem of Causal I

42、nference:由於 the counterfactual definition of causal effect 意涵著評估個人層次的因果效應會有 missing data 的問題。但是如果我們願意做一些假定的話，我們可以評估幾種 Average Causal Effects。,06.09.06,1,The Counterfactual Framework:A Review,Basic Parameters of Interest:ATT:Average Treatment effect on the Treated ATU:Average Treatment effect on the

43、UntreatedATE:Average Treatment Effect the most basic one is ATT,and there are other meaningful causal parameters of interest than these three.,06.09.06,1,反事實分析架構下分析觀察資料的方法,假設控制影響treatment之共變項後，就能達成ignorability，這也稱為selection on observables）propensity score matching如果此假定不成立的話：instrumental variableHeck

44、man selection model利用長期追蹤資料的特性，使用如 fixed effect model，change score model不同的分析方法，要做不同的假定,06.09.06,1,傾向分數配對法介紹,關秉寅政治大學社會學系2010.01.16,Introduction to Propensity Score Matching,06.09.06,1,OLS迴歸分析的問題,一般多元OLS迴歸分析，是一種ATE的估計，其作法假定其控制足夠的共變項後，可以消除基準線差異，並假定接受 treatment 者的因果效應與未接受者相同。這些假定合理嗎？OLS迴歸分析通常無法克服自我選擇的問

45、題。OLS 迴歸分析可能將接受 treatment 及沒接受 treatment 兩組中無法比較的人納入分析。如果兩組人的特性（基準線）相當不同，則 OLS的推估會有大問題，因為其無法比較的部份是以imputation 的方式來推估。,06.09.06,1,反事實分析架構下分析觀察資料的方法,假設只有 omitted variables，且是可以用控制觀察到之變項來控制的話（即控制影響 treatment之共變項後，就能達成 ignorability，這也稱為selection on observables）propensity score matching（PSM）是這類假定下發展出來的分析

46、方法。,06.09.06,1,反事實分析架構下分析觀察資料的方法,selection on observables Z D Y U(unobservable),06.09.06,1,Propensity Score Matching（PSM）,假定：如果接受及不接受 treatment 兩種人的差異能夠被一組共變項（Z）完美解釋，那麼我們就可以用這些共變項進行分層配對，使得每層內有兩種人：接受者及不接受者，且這些人在各層中唯一的不同是他們是否接受 treatment。然後，我們觀察這兩種人在 outcome 的差異為何。我們再進一步將各分層的差異，以及分層所佔的比例做適當加權，則可得到好的AT

47、E，ATT，ATU 的估計。,06.09.06,1,Propensity Score Matching（PSM）,(Y1,Y0)D|Z實際上如果有許多共變項時，配對過程很麻煩，且在樣本有限的情況下，有些分層會沒有個案。Paul Rosenbaum及Donald Rubin在一系列的論文中，證明用所謂的傾向分數（propensity score）將是否會接受視為一種機率，然後以此分數來從事分層配對是可行的。,06.09.06,1,Propensity Score Matching（PSM）,如何得到傾向分數？找到有意義可解釋是否會接受 treatment 的共變項，然後做Logit 或 Prob

48、it 迴歸，應變項為是否接受treatment。根據傾向分數將接受者及不接受者進行配對，並找到 common support。在 common support 的樣本中，兩組人的分配是 balanced（檢視matching quality）。,06.09.06,1,Propensity Score Matching（PSM）,實際從事PSM的運算方法有四大類：Exact Matching Nearest Neighbor Matching Interval Matching Kernel Matching 不同運算方法的差異：With or without replacementHow ma

49、ny units to match,06.09.06,1,Propensity Score Matching（PSM）,選擇不同運算法的兩難要能 maximize 精準的配對（如用 strictly“nearest”or common-support region），可能會有比較多的樣本個案會被排除於分析之外。要能包括比較多的配對樣本個案（如擴大配對的範圍），就會有比較多不精準的配對。用bootstrapping方式求得PSM估計值的standard errors,06.09.06,1,Propensity Score Matching（PSM）,實際可從事PSM的程式：Stata:psmat

50、ch2 等SPSS:SPSS Macro for Propensity Score Matching(http:/ssw.unc.edu/VRC/Lectures/index.htm)SAS:“GREEDY”Macro(http:/“MatchIt”(http:/gking.harvard.edu/matchit/),06.09.06,1,運用PSM的實例：補習數學有用嗎？,研究問題：國三補習數學有用嗎？ATT：如果國三補數學的人，沒補的話，數學成就會有差異嗎？ATU：如果國三沒補數學的人，補習的話，數學成就會有差異嗎？ATE：如果所有的人國三都補數學的話，數學成就會有差異嗎？參考：關秉寅、李