Identifying Dynamic Sequential Plans

Size: px
Start display at page:

Download "Identifying Dynamic Sequential Plans"

Transcription

1 Identfyng Dynamc Sequental Plans Jn Tan Department of Computer Scence Iowa State Unversty Ames, IA Abstract We address the problem of dentfyng dynamc sequental plans n the framework of causal Bayesan networks, and show that the problem s reduced to dentfyng causal effects, for whch there are complete dentfcaton algorthms avalable n the lterature. 1 Introducton Ths paper deals wth the problem of evaluatng the effects of sequental plans from a combnaton of nonexpermental data and qualtatve causal assumptons. The causal assumptons wll be represented n the form of an acyclc causal dagram [Sprtes et al., 1993, Heckerman and Shachter, 1995, Laurtzen, 2000, Pearl, 2000] n whch arrows represent the potental exstence of drect causal relatonshps between the correspondng varables. The causal dagram may contan unmeasured varables, and our task s to decde whether we can estmate the effects of a sequence of actons from the observed data. We motvate the study by consderng a medcal treatment problem dscussed n [Pearl and Robns, 1995, Dawd and Ddelez, 2005]. There are a sequence of medcal treatments (X 1,..., X k ) appled to a patent over tme. We have observatons Z before and between the treatments. Doctors may prescrbe a treatment based on prevous treatments and observatons. There s a outcome varable Y (say survval) of specal nterest and we want to estmate the effects of the sequental treatments on Y. In general there may be unobserved confounders that have nfluence on the observed varables. There are dfferent possble strateges for choosng the treatment acton (X ). The smplest acton nvolves fxng the value of X to a partcular value x, called atomc nterventon and denoted by do(x ) n [Pearl, 2000], e.g. fxng the dosage of a treatment rrespectve of any observatons on the patent. A uncondtonal plan conssts of a sequence of pre-defned atomc actons. The problem of dentfyng uncondtonal plans s to compute the dstrbuton of Y under atomc nterventons on a set X of acton varables, denoted by P x (y) = P (y do(x)), a quantty known as the causal effects of X on Y. Suffcent graphcal crtera for the dentfablty of uncondtonal plans are derved n [Pearl and Robns, 1995]. Recently the general problem of dentfyng causal effects P x (y) n a causal dagram contanng unobserved varables has been solved and complete algorthms for dentfcaton are gven n [Tan and Pearl, 2003, Shptser and Pearl, 2006b, Huang and Valtorta, 2006]. In general we may want to use dynamc treatment strateges n whch the values of acton varables (X ) are determned based on the prevously observed varables and prevously taken actons. Suffcent graphcal crtera for the dentfablty of dynamc plans are derved n [Dawd and Ddelez, 2005]. However the dentfablty problem s far from beng solved. In ths paper, we show that the problem of dentfyng dynamc sequental plans can be reduced to the well-studed problem of dentfyng causal effects and therefore essentally solved the sequental plan dentfcaton problem. Although Pearl (2000, Secton 4.2) has suggested that dynamc condtonal plans may be dentfed by dentfyng condtonal causal effects of the form P x (y c), for whch complete dentfcaton algorthms have been developed [Tan, 2004, Shptser and Pearl, 2006a], n ths paper, we wll show that ths gves a suffcent condton for dentfyng dynamc sequental plans but t s not necessary. The rest of the paper s organzed as follows. In Secton 2, we revew the work n [Dawd and Ddelez, 2005] and defne useful notaton. In Secton 3, we formulate the sequental plan problem n the framework of causal Bayesan networks. We

2 show how to reduce the problem of dentfyng dynamc sequental plans nto a problem of dentfyng causal effects n Secton 4, and dscuss n Secton 5 the problem versus that of dentfyng uncondtonal plans and condtonal causal effects. Secton 6 concludes the paper. 2 Prevous Work and Notaton Dawd and Ddelez (2005) formulated the problem of dentfyng dynamc sequental plans n the framework of regme ndcators and nfluence dagrams. An nfluence dagram (ID) s a DAG over a set V = {V 1,..., V n } of varables that also ncludes regme ndcators as specal nodes of ther own called decson nodes [Dawd, 2002]. We assume that all varables are dscrete. The DAG s assumed to represent condtonal ndependence assertons that each varable s ndependent of all ts non-descendants gven ts drect parents n the graph. 1 These assertons mply that the jont probablty functon P (v) = P (v 1,..., v n ) factorzes accordng to the product [Pearl, 1988] P (v) = P (v pa ) (1) where pa are (values of) the parents of varable V n the graph. 2 The queston of causal nference s consdered as a problem of nference across dfferent regmes, n whch we may ntervene n certan varables n certan ways and observe the behavor of other varables. Regme ndcators are used to represent dfferent types of nterventons. Here we wll roughly follow the notaton used n [Ddelez et al., 2006]. The regme ndcator for an nterventon n a varable V s denoted by σ V and can take values n a set of strateges. Under strategy σ V, the condtonal probablty P (v pa ) s changed to P (v pa ; σ V ). We wll consder the followng types of nterventons. Idle regme σ V = : No nterventon takes place, therefore P (v pa ; σ V = ) = P (v pa ). The dle regme s also called the observatonal regme under whch we wll assume that observatonal data has been collected. Therefore P (v pa ) can be estmated from data f V and P a are observed. 1 We use famly relatonshps such as parents, chldren, and ancestors to descrbe the obvous graphcal relatonshps. 2 We use uppercase letters to represent varables or sets of varables, and use correspondng lowercase letters to represent ther values (nstantatons). Atomc nterventon σ V = do(v ): The strategy of settng V to a fxed value v, denoted by do(v = v ) or smply do(v ) n Pearl (2000), such that P (v pa ; σ V = do(v )) = δ(v, v ), where δ(v, v ) s one f v = v and zero otherwse. Condtonal nterventon σ V = do(g(c)): In general, V may be made to respond n a specfed way to some set C of prevously observed varables, denoted by do(v = g(c)) n Pearl (2000), such that P (v pa ; σ V = do(g(c))) = δ(v, g(c)), where g(.) s a pre-specfed determnstc functon and the varables n C can not be descendants of V. Random nterventon σ V = d C : More generally, we may let V take on a random value accordng to some dstrbuton possbly dependng on some set C of prevously observed varables such that P (v pa ; σ V = d C ) = P (v c), where P (v c) s a pre-specfed probablty dstrbuton and the varables n C can not be descendants of V. In a sequental decson problem, we may ntervene, at least n prncple, n a set of varables X = {X } V, called control varables or acton varables, and are nterested n the response of a varable Y, called response varable or outcome varable. Let Z be the rest of observed varables whch are often called covarates. The varables are assumed to be ordered n a sequence (L 1, X 1,..., L K, X K, Y ) where L Z are the set of observed covarates after X 1 and before X. We denote L = (L 1,..., L ) and X = (X 1,..., X ). Gven an nterventon strategy σ X = {σ X }, under a condton called smple stablty whch says that the observed covarates L and the outcome Y are ndependent of how acton varables are generated once all earler observables ( L 1, X 1 for L ; X, Z for Y ) are gven, Dawd and Ddelez (2005) show that the postnterventon dstrbuton of Y s dentfed as P (y; σ X ) P (y x, z) P (l l 1, x 1 ) x,z P (x x 1, l ; σ X ), (2) where P (x x 1, l ; σ X ) are determned by the chosen regme and the other quanttes can be estmated from

3 observatonal data. Eq. (2) s known as the G-formula, and has been obtaned n [Robns, 1986, Robns, 1987] n the framework of potental response models. When there are unobserved confounders, the smple stablty may not hold. Dawd and Ddelez (2005) makes extended stablty assumpton whch essentally s (smple) stablty wth respect to the extended doman that ncludes unobserved U varables gnorng the dstncton between Z and U. The G-formula (2) no longer holds unless we nclude unobserved U varables, but then the condtonal probabltes nvolvng U varables can no longer be estmated from the data. Suffcent graphcal crtera for dentfyng P (y; σ X ) are derved. The crtera were obtaned by dentfyng graphcal condtons under whch the smple stablty can be reganed such that the G-formula can be used, and by extendng the work n [Pearl and Robns, 1995] to dynamc plans. 3 Problem Formulaton In ths paper, we wll formulate the sequental plan problem n the framework of causal Bayesan networks. A causal Bayesan network (CBN) conssts of a DAG G over a set V = {V 1,..., V n } of varables, called a causal dagram. The nterpretaton of such a graph has two components, probablstc and causal. The probablstc nterpretaton vews G as representng condtonal ndependence assertons such that the jont probablty functon P (v) = P (v 1,..., v n ) factorzes accordng to Eq. (1). The causal nterpretaton vews the drected edges n G as representng causal nfluences between the correspondng varables. In ths nterpretaton, the factorzaton of (1) stll holds, but the factors are further assumed to represent autonomous data-generaton processes, that s, each condtonal probablty P (v pa ) represents a stochastc process by whch the values of V are chosen n response to the values pa (prevously chosen for V s parents), and the stochastc varaton of ths assgnment s assumed ndependent of the varatons n all other assgnments. Moreover, each assgnment process remans nvarant to possble changes n the assgnment processes that govern other varables n the system. Ths modularty assumpton enables us to predct the effects of nterventons, whenever nterventons are descrbed as specfc modfcatons of some factors n the product of (1). We typcally assume that every varable V can potentally be manpulated by external nterventon. So we mght thnk of a CBN as an ID such that each node s (mplctly) ponted to by a correspondng regme/nterventon ndcator. In a sequental decson problem, we may ntervene n a set of acton varables X = {X } V, and are nterested n the response of a set of outcome varables Y. Assume that all the varables V are observed and let the rest of covarate varables be Z = {Z } = V \ (X Y ). Gven an nterventon strategy σ X = {σ X }, by modularty assumpton, we can predct the effects of σ X as P (v; σ X ) = P (y pa y ) P (z pa z ) P (x pa x ; σ X ), (3) where, by modularty assumpton, those condtonal probabltes correspondng to unmanpulated varables reman unaltered. We note that Dawd and Ddelez s (2005) smple stablty assumpton leads to Eq. (3) n the framework of CBNs. We see that, gven a CBN, whenever all varables n V are observed, the consequence of an nterventon strategy on the outcome varables Y s computed as P (y; σ X ) P (y pa y ) P (z pa z ) x,z P (x pa x ; σ X ), (4) where P (x pa x ; σ X ) are determned by the chosen regme and the other quanttes can be estmated from observatonal data. We note that the G-formula (2) can be reduced to Eq. (4) by usng the condtonal ndependence relatonshps mpled by the CBN that each varable s ndependent of all ts non-descendants gven ts parents. In general we may be concerned wth confoundng effects due to unobserved nfluental varables. In the presence of unobserved confounders, the dstrbuton over observed varables can no longer factorze accordng to (1). Lettng V = Y Z X and U = {U 1,..., U n } stand for the sets of observed and unobserved varables, respectvely, the observed probablty dstrbuton, P (v), becomes a mxture of products: P (v) P (v pa v ) P (u pa u ). u { V V } { U U} (5) We stll make modularty assumpton n the CBN wth unobserved varables, and the effects of an nterventon strategy σ X on the outcome varables Y can be expressed as P (y; σ X ) P (y pa y ) P (z pa z ) x,z,u P (x pa x ; σ X ) P (u pa u ). (6)

4 We note that Dawd and Ddelez s (2005) extended stablty assumpton leads to Eq. (6) n the framework of CBNs. In (6), the quanttes P (y pa y ) and P (z pa z ) (and P (u pa u )) may nvolve elements of U and may not be estmable from data. Then the queston of dentfablty arses,.e., whether t s possble to express P (y; σ X ) as a functon of the observed dstrbuton P (v). Defnton 1 [Plan Identfablty] A sequental plan s sad to be dentfable f P (y; σ X ) s unquely computable from the observed dstrbuton P (v). 4 Identfcaton of Sequental Plans Frst we make the followng assumpton about the type of nterventons we wll consder. Assumpton 1 P (x pa x ; σ X ) does not depend on the unobserved varable. That s, for condtonal nterventon σ X = do(g(c)) or random nterventon σ X = d C, we requre C X Z. Ths assumpton corresponds to Condton 6.6 or 7.2 n [Dawd and Ddelez, 2005]. Under Assumpton 1, Eq. (6) becomes P (y; σ X ) P (x pa x ; σ X ) P (y pa y ) x,z u P (z pa z ) P (u pa u ) (7) P (x pa x ; σ X )P x (y, z) (8) x,z Obvously a suffcent condton for P (y; σ X ) beng dentfable s that the causal effect P x (y, z) s dentfable. In partcular, a smple suffcent condton for P x (y, z) beng dentfable s f all the parents of acton (X) varables are observables, whch s Condton 6.3 n [Dawd and Ddelez, 2005]. Proposton 1 If all the parents of acton (X) varables are observables, then P (y; σ X ) s dentfable [Dawd and Ddelez, 2005]. Proof: If all the parents of acton (X) varables are observables, then P (x pa x ) contans no unobserved (U) varables, and Eq. (5) can be wrtten as P (v) = P (x pa x )P x (y, z), (9) from whch we obtan that P x (y, z) s dentfed as p(x, y, z) P x (y, z) = P (x (10) pa x ) = P (v v ), (11) { V Y Z} where we have used the chan rule assumng an order of V varables that s consstent wth the DAG and v denotes the V varables ordered ahead of V. Hence the sequental plan s dentfed as P (y; σ X ) x,z P (x pa x ; σ X ) whch s essentally the G-formula (2). { V Y Z} P (v v ), (12) In general P x (y, z) beng dentfable s not a necessary condton for P (y; σ X ) beng dentfable. Eq. (7) may be smplfed n that a factor P (z pa z ) may be summed out (usng z P (z pa z ) = 1) f Z does not appear n any other factors (graphcally, f Z does not have any chldren). We can derve stronger dentfcaton crteron by summng out as many factors as possble from Eq. (7). Before presentng our result, we frst ntroduce some notaton. Followng [Tan and Pearl, 2003], for any observed set S V of varables, we defne the quantty Q[S] to denote the post-nterventon dstrbuton of S under atomc nterventons to all other varables: Q[S](v) = P v\s (s) u { V S} P (v pa v ) { U U} P (u pa u ). (13) For convenence, we wll often wrte Q[S](v) as Q[S]. Eq. (7) can be wrtten as P (y; σ X ) x,z P (x pa x ; σ X )Q[Y Z]. (14) Let G σx denote the manpulated graph under the nterventon strategy σ X, whch can be constructed from the orgnal causal graph G as follows: For an atomc nterventon σ X = do(x ), cut off all the arrows enterng X ; For a condtonal nterventon σ X = do(g (c )) or a random nterventon σ X = d C, cut off all the arrows enterng X and then add an arrow enterng X from each varable n C.

5 Based on Eq. (14), we obtan the followng suffcent crteron for dentfyng P (y; σ X ). Theorem 1 Let Z D be the set of varables n Z that are ancestors of Y n G σx. P (y; σ X ) s dentfable f the causal effects Q[Y Z D ] = P x,z\zd (y, z D ) s dentfable. Proof: Let X D be the set of varables n X that are ancestors of Y n G σx. Then all the non-ancestors of Y can be summed out from Eq. (14) leadng to P (y; σ X ) P (x pa x ; σ X )Q[Y Z D ] x D,z D { X X D} x D,z D { X X D} P (x pa x ; σ X )P x,z\zd (y, z D ). (15) (16) We conjecture that the condton n Theorem 1 s also necessary. It mght appear that Eq. (15) can be further smplfed as follows. Let Z σxd be the set of Z varables that appear n the term { X P (x X D} pa x ; σ X ) (the set of condtonng varables n the strategy σ XD ). Then the term { X P (x X D} pa x ; σ X ) s a functon of X D and Z σxd. Eq. (15) becomes P (y; σ X ) x D,z σxd { X X D} x D,z σxd g(x D, z σxd ) P (x pa x ; σ X ) z D\z σxd Q[Y Z D ] (17) z D\z σxd Q[Y Z D ] (18) From Eq. (18), P (y; σ X ) s dentfable f z D\z σxd Q[Y Z D ] s dentfable, and ntutvely, the condton appears to be necessary too, snce the term g(x D, z σxd ) s specfed externally and no more factors can be summed out (as far as the functon g(.) s not ndependent of any varables n z σxd ). A strct proof of ths necessty s stll under study. On the other hand, due to the fact that none of the factors correspondng to the varables n Z D \ Z σxd can be summed out from Q[Y Z D ], t has been shown that z D\z σxd Q[Y Z D ] can be dentfed only va dentfyng Q[Y Z D ] and t s a f and only f condton (Lemma 11 n [Huang and Valtorta, 2006]). So from the pont of vew of dentfyng P (y; σ X ) the reducton from Eq. (15) to Eq. (17) s not necessary. U Z X 1 X 2 Y Fgure 1: An example causal graph We therefore have reduced the problem of dentfyng dynamc sequental plans P (y; σ X ) nto that of dentfyng causal effects Q[Y Z D ] whle the latter problem has been solved and complete algorthms are gven n [Tan and Pearl, 2003, Shptser and Pearl, 2006b, Huang and Valtorta, 2006]. We demonstrate the applcaton of Theorem 1 and the dentfcaton process wth an example. Consder the problem of dentfyng P (y; σ X1, σ X2 ) n Fgure 1, whch was studed n [Dawd and Ddelez, 2005] and s troublng to the methods presented theren. Theorem 1 calls for dentfyng Q[{Y, Z}] whch can be shown to be dentfable. We are gven the observatonal dstrbuton where P (v) = P (y x 1, x 2, z)p (x 2 )Q[{Z, X 1 }], (19) Q[{Z, X 1 }] u We want to compute P (z x 2, u)p (x 1 u)p (u). (20) P (y; σ X1, σ X2 ) P (x 1 ; σ X1 )P (x 2 ; σ X2 )P (y x 1, x 2, z)q[{z}], x 1,x 2,z (21) whch calls for computng Q[{Z}]. From Eq. (20), t s clear that From Eq. (19), we obtan Q[{Z}] x 1 Q[{Z, X 1 }]. (22) Q[{Z, X 1 }] = P (z, x 1 x 2 ) (23) Therefore Q[{Z}] s dentfed and we fnally obtan P (y; σ X1, σ X2 ) P (x 1 ; σ X1 )P (x 2 ; σ X2 )P (y x 1, x 2, z)p (z x 2 ). x 1,x 2,z (24)

6 X 1 Y U 1 X 2 Z (a) G U 2 X 1 Y U 1 X 2 Z (b) G σx Fgure 2: An example causal graph 5 Dscusson 5.1 Uncondtonal plans are easer In general dentfyng dynamc plans s more dffcult than dentfyng uncondtonal plans that nvolve only atomc nterventons. Let an nterventon strategy σ X consst of all atomc nterventons. Then the manpulated graph G σ X s a subgraph of G σx. Let Z D be the set of varables n Z that are ancestors of Y n G σ X. Then Z D s a subset of Z D. We have P (y; σ X) = P x (y) z D U 2 Q[Y Z D]. (25) Therefore dentfyng P x (y) calls for dentfyng Q[Y Z D ] whle P (y; σ X) calls for dentfyng Q[Y Z D ]. Now we have Q[Y Z D] = z D\z D Q[Y Z D ] (26) The factors of the varables n Z D \Z D are summed out from Q[Y Z D ] snce the varables n Z D \ Z D can be ancestors of Y only through X s. In general, whenever Q[Y Z D ] (and therefore P (y; σ X )) s dentfable, then P x (y) s dentfable. However, t s possble that P x (y) s dentfable but Q[Y Z D ] (and therefore P (y; σ X )) s not. We demonstrate ths pont wth an example. Consder the problem of dentfyng P (y; σ X1, σ X2 ) n Fgure 2(a), whch was studed n [Pearl and Robns, 1995]. If P (x 2 x 1, z; σ X2 ) depends on Z, say σ X2 = do(g(x 1, z)), then Z D = {Z}, and by Theorem 1, to dentfy P (y; σ X1, σ X2 ) we need to dentfy Q[{Y, Z}], whch can be shown to be not dentfable (by theorems n [Huang and Valtorta, 2006]). More specfcally, gven the observatonal dstrbuton P (v) = P (x 2 x 1, z)q[{x 1, Z, Y }], (27) we want to dentfy P (y; σ X1, σ X2 ) P (x 1 ; σ X1 )P (x 2 x 1, z; σ X2 )Q[{Y, Z}] (28) x 1,x 2,z From Eq. (28), we see that f P (x 2 x 1, z; σ X2 ) depends on Z, then the dentfablty of P (y; σ X1, σ X2 ) depends on the dentfablty of Q[{Y, Z}]. We therefore conclude that P (y; σ X1, σ X2 ) s not dentfable. On the other hand, f P (x 2 x 1, z; σ X2 ) s ndependent of Z, say P (x 2 x 1, z; σ X2 ) = P (x 2 x 1 ) (or σ X2 = do(x 2 )), then the set Z D of varables n Z that are ancestors of Y n G σx becomes empty (see Fgure 2(b)), and, by Theorem 1, the dentfablty of P (y; σ X1, σ X2 ) depends on the dentfablty of Q[{Y }]. In fact, n ths case, Eq. (28) becomes P (y; σ X1, σ X2 = d X1 ) P (x 1 ; σ X1 )P (x 2 x 1 ) Q[{Y, Z}] x 1,x 2 z P (x 1 ; σ X1 )P (x 2 x 1 )Q[{Y }] (29) x 1,x 2 From Eq. (27) we obtan Q[{X 1, Z, Y }] = P (v)/p (x 2 x 1, z) = P (y x 1, x 2, z)p (x 1, z). (30) It can be shown (or confrmed) that Q[{Y }] z Q[{X 1, Z, Y }] y,z Q[{X 1, Z, Y }] We obtan z P (y x 1, x 2, z)p (z x 1 ). (31) P (y; σ X1, σ X2 = d X1 ) P (x 1 ; σ X1 )P (x 2 x 1 ) P (y x 1, x 2, z)p (z x 1 ). x 1,x 2 z (32) And n partcular, the uncondtonal plan s dentfed as P x1,x 2 (y) = Q[{Y }] z P (y x 1, x 2, z)p (z x 1 ). 5.2 Identfcaton va condtonal causal effects? (33) Pearl (2000) has suggested that dynamc sequental plans nvolvng condtonal and random nterventons

7 may be dentfed by dentfyng condtonal causal effects of the form P x (y z). For nterventons on a sngle varable X, we can show [Pearl, 2000, Secton 4.2] that and P (y; σ X = do(g(z))) z P (y; σ X = d Z ) x,z P x (y z) x=g(z) P (z), P x (y z)p (x z)p (z). Therefore t appears that P (y; σ X ) s dentfable f and only f P x (y z) s dentfable. Ths dea was generalzed to dynamc sequental plans. Consder a plan nvolvng a sequence of condtonal nterventons σ X = do(g (C )). Let Z σx = Z ( C ) be the set of condtonng varables n the strategy σ X. Pearl (2006) shows that P (y; σ X ) z σx P xz (y z σx )P xz (z σx ), (34) where x z are the values attaned by X when the condtonng set Z σx takes the values z σx. Pearl then suggests that sequental condtonal plans can be dentfed by dentfyng condtonal causal effects P x (y z) and P x (z). Ths motvated the study of the dentfablty of condtonal causal effects and complete algorthms have been developed n [Tan, 2004, Shptser and Pearl, 2006a]. Next we show that although the dentfcaton of P xz (y z σx ) and P xz (z σx ) s suffcent for dentfyng P (y; σ X ), t s nonetheless not necessary. Rewrte Eq. (34) n the followng P (y; σ X ) δ(x, g (C ))P x (y, z σx ) (35) x,z σx δ(x, g (C )) Q[Y, Z] (36) x,z σx z\z σx Comparng the reducton from Eq. (14) nto (17) wth the reducton from (14) to (36), we obtan that f X D = X then (36) s equvalent to (17), otherwse Eq. (36) may be further reduced n that more factors could be summed out from Q[Y, Z]. We obtan the followng concluson If all the varables n X are ancestors of Y n G σx, then the sequental plan P (y; σ X ) can be dentfed by dentfyng the causal effects P x (y, z σx ), otherwse t s possble that P (y; σ X ) s dentfable even f P x (y, z σx ) s not. We demonstrate ths pont wth an example. Consder the problem of dentfyng P (y; σ X ) where Z 1 U 1 U 2 Z 3 X 1 Z 2 X 2 X 3 Y (a) G Z 1 U 1 U 2 Z 3 X 1 Z 2 X 2 X 3 Y (b) G σx Fgure 3: An example causal graph σ X = {σ X1 = do(g 1 (Z 1 )), σ X2 = do(g 2 (Z 2 )), σ X3 = do(g 3 (Z 3 ))} n Fgure 3(a). The graph G σx s shown n Fgure 3(b). We have Z D = {Z 1, Z 3 }, and Theorem 1 calls for dentfyng Q[{Y, Z 1, Z 3 }] whch can be shown to be dentfable. On the other hand, P x1x 2x 3 (y, z 1, z 2, z 3 ) = Q[{Y, Z 1, Z 2, Z 3 }] s not dentfable. More specfcally, gven the observatonal dstrbuton P (v) =P (y x 1, x 3, z 3 )P (x 3 x 2, z 3 )P (z 1 )Q[{X 2, Z 3 }] Q[{X 1, Z 2 }], (37) we want to dentfy P (y; σ X ) δ(x, g (z ))P (y x 1, x 3, z 3 )P (z 1 ) x,z Q[{Z 3 }]Q[{Z 2 }] (38) x 1,x 3,z 1,z 3 δ(x 1, g 1 (z 1 ))δ(x 3, g 3 (z 3 ))P (y x 1, x 3, z 3 ) P (z 1 )Q[{Z 3 }], (39) where Q[{Z 3 }] can be dentfed as We obtan Q[{Z 3 }] u 2 P (z 3 u 2 )P (u 2 ) = P (z 3 ). (40) P (y; σ X ) z 1,z 3 P (y g 1 (z 1 ), g 3 (z 3 ), z 3 )P (z 1 )P (z 3 ). On the other hand, P x1x 2x 3 (y, z 1, z 2, z 3 ) = P (y x 1, x 3, z 3 )P (z 1 )Q[{Z 3 }]Q[{Z 2 }] (41) = P (y x 1, x 3, z 3 )P (z 1 )P (z 3 )Q[{Z 2 }] (42)

8 s not dentfable snce Q[{Z 2 }] s not dentfable. In fact the condtonal causal effect P x1x 2x 3 (y z 1, z 2, z 3 ) = P (y x 1, x 3, z 3 ) (43) s dentfable but the causal effect P x1x 2x 3 (z 1, z 2, z 3 ) = P (z 1 )P (z 3 )Q[{Z 2 }] (44) s not dentfable. 6 Concluson We present a method for dentfyng dynamc sequental plans. A closed-form expresson for the probablty of the outcome varables under a dynamc plan can be obtaned n terms of the observed dstrbuton, by usng the algorthms for dentfyng causal effects avalable n the lterature. Acknowledgments Ths research was partly supported by NSF grant IIS References [Dawd and Ddelez, 2005] A.P. Dawd and V. Ddelez. Identfyng the consequences of dynamc treatment strateges. Techncal report, Department of Statstcal Scence, Unversty College London, UK, [Dawd, 2002] A.P. Dawd. Influence dagrams for causal modellng and nference. Internatonal Statstcal Revew, 70(2), [Ddelez et al., 2006] V. Ddelez, A.P. Dawd, and S. Genelett. Drect and ndrect effects of sequental treatments. In Proceedngs of the 22nd Annual Conference on Uncertanty n Artfcal Intellgence, pages , [Heckerman and Shachter, 1995] D. Heckerman and R. Shachter. Decson-theoretc foundatons for causal reasonng. Journal of Artfcal Intellgence Research, 3: , [Huang and Valtorta, 2006] Y. Huang and M. Valtorta. Identfablty n causal bayesan networks: A sound and complete algorthm. In Proceedngs of the Twenty-Frst Natonal Conference on Artfcal Intellgence, pages , Menlo Park, CA, July AAAI Press. [Laurtzen, 2000] S. Laurtzen. Graphcal models for causal nference. In O.E. Barndorff-Nelsen, D. Cox, and C. Kluppelberg, edtors, Complex Stochastc Systems, chapter 2, pages Chapman and Hall/CRC Press, London/Boca Raton, [Pearl and Robns, 1995] J. Pearl and J.M. Robns. Probablstc evaluaton of sequental plans from causal models wth hdden varables. In P. Besnard and S. Hanks, edtors, Uncertanty n Artfcal Intellgence 11, pages Morgan Kaufmann, San Francsco, [Pearl, 1988] J. Pearl. Probablstc Reasonng n Intellgent Systems. Morgan Kaufmann, San Mateo, CA, [Pearl, 2000] J. Pearl. Causalty: Models, Reasonng, and Inference. Cambrdge Unversty Press, NY, [Pearl, 2006] J. Pearl. Comment on dentfyng condtonal plans [Robns, 1986] J.M. Robns. A new approach to causal nference n mortalty studes wth a sustaned exposure perod applcatons to control of the healthy workers survvor effect. Mathematcal Modelng, 7: , [Robns, 1987] J.M. Robns. A graphcal approach to the dentfcaton and estmaton of causal parameters n mortalty studes wth sustaned exposure perods. Journal of Chronc Dseases, 40(Suppl 2):139S 161S, [Shptser and Pearl, 2006a] I. Shptser and J. Pearl. Identfcaton of condtonal nterventonal dstrbutons. In R. Dechter and T.S. Rchardson, edtors, Proceedngs of the Twenty-Second Conference on Uncertanty n Artfcal Intellgence, pages AUAI Press, July [Shptser and Pearl, 2006b] I. Shptser and J. Pearl. Identfcaton of jont nterventonal dstrbutons n recursve sem-markovan causal models. In Proceedngs of the Twenty-Frst Natonal Conference on Artfcal Intellgence, pages , Menlo Park, CA, July AAAI Press. [Sprtes et al., 1993] P. Sprtes, C. Glymour, and R. Schenes. Causaton, Predcton, and Search. Sprnger-Verlag, New York, [Tan and Pearl, 2003] J. Tan and J. Pearl. On the dentfcaton of causal effects. Techncal Report R-290-L, Department of Computer Scence, Unversty of Calforna, Los Angeles, jtan/r290-l.pdf. [Tan, 2004] J. Tan. Identfyng condtonal causal effects. In Proceedngs of the Conference on Uncertanty n Artfcal Intellgence (UAI), 2004.

A General Identification Condition for Causal Effects

A General Identification Condition for Causal Effects From: AAAI-02 Proceedngs. Copyrght 2002, AAAI (www.aaa.org). All rghts reserved. A General Identfcaton Condton for Causal Effects Jn Tan and Judea Pearl Cogntve Systems Laboratory Computer Scence Department

More information

The Multiple Classical Linear Regression Model (CLRM): Specification and Assumptions. 1. Introduction

The Multiple Classical Linear Regression Model (CLRM): Specification and Assumptions. 1. Introduction ECONOMICS 5* -- NOTE (Summary) ECON 5* -- NOTE The Multple Classcal Lnear Regresson Model (CLRM): Specfcaton and Assumptons. Introducton CLRM stands for the Classcal Lnear Regresson Model. The CLRM s also

More information

Outline. Bayesian Networks: Maximum Likelihood Estimation and Tree Structure Learning. Our Model and Data. Outline

Outline. Bayesian Networks: Maximum Likelihood Estimation and Tree Structure Learning. Our Model and Data. Outline Outlne Bayesan Networks: Maxmum Lkelhood Estmaton and Tree Structure Learnng Huzhen Yu janey.yu@cs.helsnk.f Dept. Computer Scence, Unv. of Helsnk Probablstc Models, Sprng, 200 Notces: I corrected a number

More information

Hidden Markov Models & The Multivariate Gaussian (10/26/04)

Hidden Markov Models & The Multivariate Gaussian (10/26/04) CS281A/Stat241A: Statstcal Learnng Theory Hdden Markov Models & The Multvarate Gaussan (10/26/04) Lecturer: Mchael I. Jordan Scrbes: Jonathan W. Hu 1 Hdden Markov Models As a bref revew, hdden Markov models

More information

Hidden Markov Models

Hidden Markov Models CM229S: Machne Learnng for Bonformatcs Lecture 12-05/05/2016 Hdden Markov Models Lecturer: Srram Sankararaman Scrbe: Akshay Dattatray Shnde Edted by: TBD 1 Introducton For a drected graph G we can wrte

More information

2E Pattern Recognition Solutions to Introduction to Pattern Recognition, Chapter 2: Bayesian pattern classification

2E Pattern Recognition Solutions to Introduction to Pattern Recognition, Chapter 2: Bayesian pattern classification E395 - Pattern Recognton Solutons to Introducton to Pattern Recognton, Chapter : Bayesan pattern classfcaton Preface Ths document s a soluton manual for selected exercses from Introducton to Pattern Recognton

More information

Time-Varying Systems and Computations Lecture 6

Time-Varying Systems and Computations Lecture 6 Tme-Varyng Systems and Computatons Lecture 6 Klaus Depold 14. Januar 2014 The Kalman Flter The Kalman estmaton flter attempts to estmate the actual state of an unknown dscrete dynamcal system, gven nosy

More information

Using T.O.M to Estimate Parameter of distributions that have not Single Exponential Family

Using T.O.M to Estimate Parameter of distributions that have not Single Exponential Family IOSR Journal of Mathematcs IOSR-JM) ISSN: 2278-5728. Volume 3, Issue 3 Sep-Oct. 202), PP 44-48 www.osrjournals.org Usng T.O.M to Estmate Parameter of dstrbutons that have not Sngle Exponental Famly Jubran

More information

A New Evolutionary Computation Based Approach for Learning Bayesian Network

A New Evolutionary Computation Based Approach for Learning Bayesian Network Avalable onlne at www.scencedrect.com Proceda Engneerng 15 (2011) 4026 4030 Advanced n Control Engneerng and Informaton Scence A New Evolutonary Computaton Based Approach for Learnng Bayesan Network Yungang

More information

Polynomial Constraints in Causal Bayesian Networks

Polynomial Constraints in Causal Bayesian Networks Polynomal Constrants n Causal Bayesan Networks Changsung Kang Department of Computer Scence Iowa State Unversty Ames, IA 50011 cskang@astate.edu Jn Tan Department of Computer Scence Iowa State Unversty

More information

3.1 Expectation of Functions of Several Random Variables. )' be a k-dimensional discrete or continuous random vector, with joint PMF p (, E X E X1 E X

3.1 Expectation of Functions of Several Random Variables. )' be a k-dimensional discrete or continuous random vector, with joint PMF p (, E X E X1 E X Statstcs 1: Probablty Theory II 37 3 EPECTATION OF SEVERAL RANDOM VARIABLES As n Probablty Theory I, the nterest n most stuatons les not on the actual dstrbuton of a random vector, but rather on a number

More information

Inner Product. Euclidean Space. Orthonormal Basis. Orthogonal

Inner Product. Euclidean Space. Orthonormal Basis. Orthogonal Inner Product Defnton 1 () A Eucldean space s a fnte-dmensonal vector space over the reals R, wth an nner product,. Defnton 2 (Inner Product) An nner product, on a real vector space X s a symmetrc, blnear,

More information

Econ107 Applied Econometrics Topic 3: Classical Model (Studenmund, Chapter 4)

Econ107 Applied Econometrics Topic 3: Classical Model (Studenmund, Chapter 4) I. Classcal Assumptons Econ7 Appled Econometrcs Topc 3: Classcal Model (Studenmund, Chapter 4) We have defned OLS and studed some algebrac propertes of OLS. In ths topc we wll study statstcal propertes

More information

Bayesian Networks. Course: CS40022 Instructor: Dr. Pallab Dasgupta

Bayesian Networks. Course: CS40022 Instructor: Dr. Pallab Dasgupta Bayesan Networks Course: CS40022 Instructor: Dr. Pallab Dasgupta Department of Computer Scence & Engneerng Indan Insttute of Technology Kharagpur Example Burglar alarm at home Farly relable at detectng

More information

4 Analysis of Variance (ANOVA) 5 ANOVA. 5.1 Introduction. 5.2 Fixed Effects ANOVA

4 Analysis of Variance (ANOVA) 5 ANOVA. 5.1 Introduction. 5.2 Fixed Effects ANOVA 4 Analyss of Varance (ANOVA) 5 ANOVA 51 Introducton ANOVA ANOVA s a way to estmate and test the means of multple populatons We wll start wth one-way ANOVA If the populatons ncluded n the study are selected

More information

Artificial Intelligence Bayesian Networks

Artificial Intelligence Bayesian Networks Artfcal Intellgence Bayesan Networks Adapted from sldes by Tm Fnn and Mare desjardns. Some materal borrowed from Lse Getoor. 1 Outlne Bayesan networks Network structure Condtonal probablty tables Condtonal

More information

NUMERICAL DIFFERENTIATION

NUMERICAL DIFFERENTIATION NUMERICAL DIFFERENTIATION 1 Introducton Dfferentaton s a method to compute the rate at whch a dependent output y changes wth respect to the change n the ndependent nput x. Ths rate of change s called the

More information

Homework Assignment 3 Due in class, Thursday October 15

Homework Assignment 3 Due in class, Thursday October 15 Homework Assgnment 3 Due n class, Thursday October 15 SDS 383C Statstcal Modelng I 1 Rdge regresson and Lasso 1. Get the Prostrate cancer data from http://statweb.stanford.edu/~tbs/elemstatlearn/ datasets/prostate.data.

More information

Conjugacy and the Exponential Family

Conjugacy and the Exponential Family CS281B/Stat241B: Advanced Topcs n Learnng & Decson Makng Conjugacy and the Exponental Famly Lecturer: Mchael I. Jordan Scrbes: Bran Mlch 1 Conjugacy In the prevous lecture, we saw conjugate prors for the

More information

The Expectation-Maximization Algorithm

The Expectation-Maximization Algorithm The Expectaton-Maxmaton Algorthm Charles Elan elan@cs.ucsd.edu November 16, 2007 Ths chapter explans the EM algorthm at multple levels of generalty. Secton 1 gves the standard hgh-level verson of the algorthm.

More information

Hidden Markov Models

Hidden Markov Models Hdden Markov Models Namrata Vaswan, Iowa State Unversty Aprl 24, 204 Hdden Markov Model Defntons and Examples Defntons:. A hdden Markov model (HMM) refers to a set of hdden states X 0, X,..., X t,...,

More information

ANSWERS. Problem 1. and the moment generating function (mgf) by. defined for any real t. Use this to show that E( U) var( U)

ANSWERS. Problem 1. and the moment generating function (mgf) by. defined for any real t. Use this to show that E( U) var( U) Econ 413 Exam 13 H ANSWERS Settet er nndelt 9 deloppgaver, A,B,C, som alle anbefales å telle lkt for å gøre det ltt lettere å stå. Svar er gtt . Unfortunately, there s a prntng error n the hnt of

More information

A Bayes Algorithm for the Multitask Pattern Recognition Problem Direct Approach

A Bayes Algorithm for the Multitask Pattern Recognition Problem Direct Approach A Bayes Algorthm for the Multtask Pattern Recognton Problem Drect Approach Edward Puchala Wroclaw Unversty of Technology, Char of Systems and Computer etworks, Wybrzeze Wyspanskego 7, 50-370 Wroclaw, Poland

More information

Bayesian predictive Configural Frequency Analysis

Bayesian predictive Configural Frequency Analysis Psychologcal Test and Assessment Modelng, Volume 54, 2012 (3), 285-292 Bayesan predctve Confgural Frequency Analyss Eduardo Gutérrez-Peña 1 Abstract Confgural Frequency Analyss s a method for cell-wse

More information

A new construction of 3-separable matrices via an improved decoding of Macula s construction

A new construction of 3-separable matrices via an improved decoding of Macula s construction Dscrete Optmzaton 5 008 700 704 Contents lsts avalable at ScenceDrect Dscrete Optmzaton journal homepage: wwwelsevercom/locate/dsopt A new constructon of 3-separable matrces va an mproved decodng of Macula

More information

THE SUMMATION NOTATION Ʃ

THE SUMMATION NOTATION Ʃ Sngle Subscrpt otaton THE SUMMATIO OTATIO Ʃ Most of the calculatons we perform n statstcs are repettve operatons on lsts of numbers. For example, we compute the sum of a set of numbers, or the sum of the

More information

Methods Lunch Talk: Causal Mediation Analysis

Methods Lunch Talk: Causal Mediation Analysis Methods Lunch Talk: Causal Medaton Analyss Taeyong Park Washngton Unversty n St. Lous Aprl 9, 2015 Park (Wash U.) Methods Lunch Aprl 9, 2015 1 / 1 References Baron and Kenny. 1986. The Moderator-Medator

More information

Global Sensitivity. Tuesday 20 th February, 2018

Global Sensitivity. Tuesday 20 th February, 2018 Global Senstvty Tuesday 2 th February, 28 ) Local Senstvty Most senstvty analyses [] are based on local estmates of senstvty, typcally by expandng the response n a Taylor seres about some specfc values

More information

COMPARISON OF SOME RELIABILITY CHARACTERISTICS BETWEEN REDUNDANT SYSTEMS REQUIRING SUPPORTING UNITS FOR THEIR OPERATIONS

COMPARISON OF SOME RELIABILITY CHARACTERISTICS BETWEEN REDUNDANT SYSTEMS REQUIRING SUPPORTING UNITS FOR THEIR OPERATIONS Avalable onlne at http://sck.org J. Math. Comput. Sc. 3 (3), No., 6-3 ISSN: 97-537 COMPARISON OF SOME RELIABILITY CHARACTERISTICS BETWEEN REDUNDANT SYSTEMS REQUIRING SUPPORTING UNITS FOR THEIR OPERATIONS

More information

CIS587 - Artificial Intellgence. Bayesian Networks CIS587 - AI. KB for medical diagnosis. Example.

CIS587 - Artificial Intellgence. Bayesian Networks CIS587 - AI. KB for medical diagnosis. Example. CIS587 - Artfcal Intellgence Bayesan Networks KB for medcal dagnoss. Example. We want to buld a KB system for the dagnoss of pneumona. Problem descrpton: Dsease: pneumona Patent symptoms (fndngs, lab tests):

More information

Module 9. Lecture 6. Duality in Assignment Problems

Module 9. Lecture 6. Duality in Assignment Problems Module 9 1 Lecture 6 Dualty n Assgnment Problems In ths lecture we attempt to answer few other mportant questons posed n earler lecture for (AP) and see how some of them can be explaned through the concept

More information

Chapter 8 Indicator Variables

Chapter 8 Indicator Variables Chapter 8 Indcator Varables In general, e explanatory varables n any regresson analyss are assumed to be quanttatve n nature. For example, e varables lke temperature, dstance, age etc. are quanttatve n

More information

A Robust Method for Calculating the Correlation Coefficient

A Robust Method for Calculating the Correlation Coefficient A Robust Method for Calculatng the Correlaton Coeffcent E.B. Nven and C. V. Deutsch Relatonshps between prmary and secondary data are frequently quantfed usng the correlaton coeffcent; however, the tradtonal

More information

Constant Best-Response Functions: Interpreting Cournot

Constant Best-Response Functions: Interpreting Cournot Internatonal Journal of Busness and Economcs, 009, Vol. 8, No., -6 Constant Best-Response Functons: Interpretng Cournot Zvan Forshner Department of Economcs, Unversty of Hafa, Israel Oz Shy * Research

More information

Limited Dependent Variables and Panel Data. Tibor Hanappi

Limited Dependent Variables and Panel Data. Tibor Hanappi Lmted Dependent Varables and Panel Data Tbor Hanapp 30.06.2010 Lmted Dependent Varables Dscrete: Varables that can take onl a countable number of values Censored/Truncated: Data ponts n some specfc range

More information

Vapnik-Chervonenkis theory

Vapnik-Chervonenkis theory Vapnk-Chervonenks theory Rs Kondor June 13, 2008 For the purposes of ths lecture, we restrct ourselves to the bnary supervsed batch learnng settng. We assume that we have an nput space X, and an unknown

More information

Uncertainty and auto-correlation in. Measurement

Uncertainty and auto-correlation in. Measurement Uncertanty and auto-correlaton n arxv:1707.03276v2 [physcs.data-an] 30 Dec 2017 Measurement Markus Schebl Federal Offce of Metrology and Surveyng (BEV), 1160 Venna, Austra E-mal: markus.schebl@bev.gv.at

More information

x = , so that calculated

x = , so that calculated Stat 4, secton Sngle Factor ANOVA notes by Tm Plachowsk n chapter 8 we conducted hypothess tests n whch we compared a sngle sample s mean or proporton to some hypotheszed value Chapter 9 expanded ths to

More information

Difference Equations

Difference Equations Dfference Equatons c Jan Vrbk 1 Bascs Suppose a sequence of numbers, say a 0,a 1,a,a 3,... s defned by a certan general relatonshp between, say, three consecutve values of the sequence, e.g. a + +3a +1

More information

Linear Regression Analysis: Terminology and Notation

Linear Regression Analysis: Terminology and Notation ECON 35* -- Secton : Basc Concepts of Regresson Analyss (Page ) Lnear Regresson Analyss: Termnology and Notaton Consder the generc verson of the smple (two-varable) lnear regresson model. It s represented

More information

Open Systems: Chemical Potential and Partial Molar Quantities Chemical Potential

Open Systems: Chemical Potential and Partial Molar Quantities Chemical Potential Open Systems: Chemcal Potental and Partal Molar Quanttes Chemcal Potental For closed systems, we have derved the followng relatonshps: du = TdS pdv dh = TdS + Vdp da = SdT pdv dg = VdP SdT For open systems,

More information

Comparison of Regression Lines

Comparison of Regression Lines STATGRAPHICS Rev. 9/13/2013 Comparson of Regresson Lnes Summary... 1 Data Input... 3 Analyss Summary... 4 Plot of Ftted Model... 6 Condtonal Sums of Squares... 6 Analyss Optons... 7 Forecasts... 8 Confdence

More information

MATH 829: Introduction to Data Mining and Analysis The EM algorithm (part 2)

MATH 829: Introduction to Data Mining and Analysis The EM algorithm (part 2) 1/16 MATH 829: Introducton to Data Mnng and Analyss The EM algorthm (part 2) Domnque Gullot Departments of Mathematcal Scences Unversty of Delaware Aprl 20, 2016 Recall 2/16 We are gven ndependent observatons

More information

STAT 3008 Applied Regression Analysis

STAT 3008 Applied Regression Analysis STAT 3008 Appled Regresson Analyss Tutoral : Smple Lnear Regresson LAI Chun He Department of Statstcs, The Chnese Unversty of Hong Kong 1 Model Assumpton To quantfy the relatonshp between two factors,

More information

Psychology 282 Lecture #24 Outline Regression Diagnostics: Outliers

Psychology 282 Lecture #24 Outline Regression Diagnostics: Outliers Psychology 282 Lecture #24 Outlne Regresson Dagnostcs: Outlers In an earler lecture we studed the statstcal assumptons underlyng the regresson model, ncludng the followng ponts: Formal statement of assumptons.

More information

Reasoning under Uncertainty

Reasoning under Uncertainty Reasonng under Uncertanty Course: CS40022 Instructor: Dr. Pallab Dasgupta Department of Computer Scence & Engneerng Indan Insttute of Technology Kharagpur Handlng uncertan knowledge p p Symptom(p, Toothache

More information

Introduction to Vapor/Liquid Equilibrium, part 2. Raoult s Law:

Introduction to Vapor/Liquid Equilibrium, part 2. Raoult s Law: CE304, Sprng 2004 Lecture 4 Introducton to Vapor/Lqud Equlbrum, part 2 Raoult s Law: The smplest model that allows us do VLE calculatons s obtaned when we assume that the vapor phase s an deal gas, and

More information

Convergence of random processes

Convergence of random processes DS-GA 12 Lecture notes 6 Fall 216 Convergence of random processes 1 Introducton In these notes we study convergence of dscrete random processes. Ths allows to characterze phenomena such as the law of large

More information

Winter 2008 CS567 Stochastic Linear/Integer Programming Guest Lecturer: Xu, Huan

Winter 2008 CS567 Stochastic Linear/Integer Programming Guest Lecturer: Xu, Huan Wnter 2008 CS567 Stochastc Lnear/Integer Programmng Guest Lecturer: Xu, Huan Class 2: More Modelng Examples 1 Capacty Expanson Capacty expanson models optmal choces of the tmng and levels of nvestments

More information

Lecture 12: Discrete Laplacian

Lecture 12: Discrete Laplacian Lecture 12: Dscrete Laplacan Scrbe: Tanye Lu Our goal s to come up wth a dscrete verson of Laplacan operator for trangulated surfaces, so that we can use t n practce to solve related problems We are mostly

More information

An identification algorithm of model kinetic parameters of the interfacial layer growth in fiber composites

An identification algorithm of model kinetic parameters of the interfacial layer growth in fiber composites IOP Conference Seres: Materals Scence and Engneerng PAPER OPE ACCESS An dentfcaton algorthm of model knetc parameters of the nterfacal layer growth n fber compostes o cte ths artcle: V Zubov et al 216

More information

Lecture Space-Bounded Derandomization

Lecture Space-Bounded Derandomization Notes on Complexty Theory Last updated: October, 2008 Jonathan Katz Lecture Space-Bounded Derandomzaton 1 Space-Bounded Derandomzaton We now dscuss derandomzaton of space-bounded algorthms. Here non-trval

More information

Lecture 3 Stat102, Spring 2007

Lecture 3 Stat102, Spring 2007 Lecture 3 Stat0, Sprng 007 Chapter 3. 3.: Introducton to regresson analyss Lnear regresson as a descrptve technque The least-squares equatons Chapter 3.3 Samplng dstrbuton of b 0, b. Contnued n net lecture

More information

Module 3 LOSSY IMAGE COMPRESSION SYSTEMS. Version 2 ECE IIT, Kharagpur

Module 3 LOSSY IMAGE COMPRESSION SYSTEMS. Version 2 ECE IIT, Kharagpur Module 3 LOSSY IMAGE COMPRESSION SYSTEMS Verson ECE IIT, Kharagpur Lesson 6 Theory of Quantzaton Verson ECE IIT, Kharagpur Instructonal Objectves At the end of ths lesson, the students should be able to:

More information

Week 5: Neural Networks

Week 5: Neural Networks Week 5: Neural Networks Instructor: Sergey Levne Neural Networks Summary In the prevous lecture, we saw how we can construct neural networks by extendng logstc regresson. Neural networks consst of multple

More information

The Feynman path integral

The Feynman path integral The Feynman path ntegral Aprl 3, 205 Hesenberg and Schrödnger pctures The Schrödnger wave functon places the tme dependence of a physcal system n the state, ψ, t, where the state s a vector n Hlbert space

More information

Power law and dimension of the maximum value for belief distribution with the max Deng entropy

Power law and dimension of the maximum value for belief distribution with the max Deng entropy Power law and dmenson of the maxmum value for belef dstrbuton wth the max Deng entropy Bngy Kang a, a College of Informaton Engneerng, Northwest A&F Unversty, Yanglng, Shaanx, 712100, Chna. Abstract Deng

More information

Marginal Models for categorical data.

Marginal Models for categorical data. Margnal Models for categorcal data Applcaton to condtonal ndependence and graphcal models Wcher Bergsma 1 Marcel Croon 2 Jacques Hagenaars 2 Tamas Rudas 3 1 London School of Economcs and Poltcal Scence

More information

Finding Dense Subgraphs in G(n, 1/2)

Finding Dense Subgraphs in G(n, 1/2) Fndng Dense Subgraphs n Gn, 1/ Atsh Das Sarma 1, Amt Deshpande, and Rav Kannan 1 Georga Insttute of Technology,atsh@cc.gatech.edu Mcrosoft Research-Bangalore,amtdesh,annan@mcrosoft.com Abstract. Fndng

More information

Assortment Optimization under MNL

Assortment Optimization under MNL Assortment Optmzaton under MNL Haotan Song Aprl 30, 2017 1 Introducton The assortment optmzaton problem ams to fnd the revenue-maxmzng assortment of products to offer when the prces of products are fxed.

More information

FACTORIZATION IN KRULL MONOIDS WITH INFINITE CLASS GROUP

FACTORIZATION IN KRULL MONOIDS WITH INFINITE CLASS GROUP C O L L O Q U I U M M A T H E M A T I C U M VOL. 80 1999 NO. 1 FACTORIZATION IN KRULL MONOIDS WITH INFINITE CLASS GROUP BY FLORIAN K A I N R A T H (GRAZ) Abstract. Let H be a Krull monod wth nfnte class

More information

Module 2. Random Processes. Version 2 ECE IIT, Kharagpur

Module 2. Random Processes. Version 2 ECE IIT, Kharagpur Module Random Processes Lesson 6 Functons of Random Varables After readng ths lesson, ou wll learn about cdf of functon of a random varable. Formula for determnng the pdf of a random varable. Let, X be

More information

NEW ASTERISKS IN VERSION 2.0 OF ACTIVEPI

NEW ASTERISKS IN VERSION 2.0 OF ACTIVEPI NEW ASTERISKS IN VERSION 2.0 OF ACTIVEPI ASTERISK ADDED ON LESSON PAGE 3-1 after the second sentence under Clncal Trals Effcacy versus Effectveness versus Effcency The apprasal of a new or exstng healthcare

More information

More metrics on cartesian products

More metrics on cartesian products More metrcs on cartesan products If (X, d ) are metrc spaces for 1 n, then n Secton II4 of the lecture notes we defned three metrcs on X whose underlyng topologes are the product topology The purpose of

More information

Primer on High-Order Moment Estimators

Primer on High-Order Moment Estimators Prmer on Hgh-Order Moment Estmators Ton M. Whted July 2007 The Errors-n-Varables Model We wll start wth the classcal EIV for one msmeasured regressor. The general case s n Erckson and Whted Econometrc

More information

6. Stochastic processes (2)

6. Stochastic processes (2) 6. Stochastc processes () Lect6.ppt S-38.45 - Introducton to Teletraffc Theory Sprng 5 6. Stochastc processes () Contents Markov processes Brth-death processes 6. Stochastc processes () Markov process

More information

Department of Quantitative Methods & Information Systems. Time Series and Their Components QMIS 320. Chapter 6

Department of Quantitative Methods & Information Systems. Time Series and Their Components QMIS 320. Chapter 6 Department of Quanttatve Methods & Informaton Systems Tme Seres and Ther Components QMIS 30 Chapter 6 Fall 00 Dr. Mohammad Zanal These sldes were modfed from ther orgnal source for educatonal purpose only.

More information

6. Stochastic processes (2)

6. Stochastic processes (2) Contents Markov processes Brth-death processes Lect6.ppt S-38.45 - Introducton to Teletraffc Theory Sprng 5 Markov process Consder a contnuous-tme and dscrete-state stochastc process X(t) wth state space

More information

A Refinement of the Common Cause Principle

A Refinement of the Common Cause Principle A Refnement of the Common Cause Prncple Nhat Ay SFI WORKING PAPER: 2008-01-001 SFI Workng Papers contan accounts of scentfc work of the author(s) and do not necessarly represent the vews of the Santa Fe

More information

Chapter 9: Statistical Inference and the Relationship between Two Variables

Chapter 9: Statistical Inference and the Relationship between Two Variables Chapter 9: Statstcal Inference and the Relatonshp between Two Varables Key Words The Regresson Model The Sample Regresson Equaton The Pearson Correlaton Coeffcent Learnng Outcomes After studyng ths chapter,

More information

8 : Learning in Fully Observed Markov Networks. 1 Why We Need to Learn Undirected Graphical Models. 2 Structural Learning for Completely Observed MRF

8 : Learning in Fully Observed Markov Networks. 1 Why We Need to Learn Undirected Graphical Models. 2 Structural Learning for Completely Observed MRF 10-708: Probablstc Graphcal Models 10-708, Sprng 2014 8 : Learnng n Fully Observed Markov Networks Lecturer: Erc P. Xng Scrbes: Meng Song, L Zhou 1 Why We Need to Learn Undrected Graphcal Models In the

More information

j) = 1 (note sigma notation) ii. Continuous random variable (e.g. Normal distribution) 1. density function: f ( x) 0 and f ( x) dx = 1

j) = 1 (note sigma notation) ii. Continuous random variable (e.g. Normal distribution) 1. density function: f ( x) 0 and f ( x) dx = 1 Random varables Measure of central tendences and varablty (means and varances) Jont densty functons and ndependence Measures of assocaton (covarance and correlaton) Interestng result Condtonal dstrbutons

More information

Probability Theory (revisited)

Probability Theory (revisited) Probablty Theory (revsted) Summary Probablty v.s. plausblty Random varables Smulaton of Random Experments Challenge The alarm of a shop rang. Soon afterwards, a man was seen runnng n the street, persecuted

More information

EM and Structure Learning

EM and Structure Learning EM and Structure Learnng Le Song Machne Learnng II: Advanced Topcs CSE 8803ML, Sprng 2012 Partally observed graphcal models Mxture Models N(μ 1, Σ 1 ) Z X N N(μ 2, Σ 2 ) 2 Gaussan mxture model Consder

More information

COS 521: Advanced Algorithms Game Theory and Linear Programming

COS 521: Advanced Algorithms Game Theory and Linear Programming COS 521: Advanced Algorthms Game Theory and Lnear Programmng Moses Charkar February 27, 2013 In these notes, we ntroduce some basc concepts n game theory and lnear programmng (LP). We show a connecton

More information

Maximizing the number of nonnegative subsets

Maximizing the number of nonnegative subsets Maxmzng the number of nonnegatve subsets Noga Alon Hao Huang December 1, 213 Abstract Gven a set of n real numbers, f the sum of elements of every subset of sze larger than k s negatve, what s the maxmum

More information

Linear Approximation with Regularization and Moving Least Squares

Linear Approximation with Regularization and Moving Least Squares Lnear Approxmaton wth Regularzaton and Movng Least Squares Igor Grešovn May 007 Revson 4.6 (Revson : March 004). 5 4 3 0.5 3 3.5 4 Contents: Lnear Fttng...4. Weghted Least Squares n Functon Approxmaton...

More information

C/CS/Phy191 Problem Set 3 Solutions Out: Oct 1, 2008., where ( 00. ), so the overall state of the system is ) ( ( ( ( 00 ± 11 ), Φ ± = 1

C/CS/Phy191 Problem Set 3 Solutions Out: Oct 1, 2008., where ( 00. ), so the overall state of the system is ) ( ( ( ( 00 ± 11 ), Φ ± = 1 C/CS/Phy9 Problem Set 3 Solutons Out: Oct, 8 Suppose you have two qubts n some arbtrary entangled state ψ You apply the teleportaton protocol to each of the qubts separately What s the resultng state obtaned

More information

Bayesian belief networks

Bayesian belief networks CS 1571 Introducton to I Lecture 24 ayesan belef networks los Hauskrecht mlos@cs.ptt.edu 5329 Sennott Square CS 1571 Intro to I dmnstraton Homework assgnment 10 s out and due next week Fnal exam: December

More information

Bayesian decision theory. Nuno Vasconcelos ECE Department, UCSD

Bayesian decision theory. Nuno Vasconcelos ECE Department, UCSD Bayesan decson theory Nuno Vasconcelos ECE Department UCSD Notaton the notaton n DHS s qute sloppy e.. show that error error z z dz really not clear what ths means we wll use the follown notaton subscrpts

More information

Transfer Functions. Convenient representation of a linear, dynamic model. A transfer function (TF) relates one input and one output: ( ) system

Transfer Functions. Convenient representation of a linear, dynamic model. A transfer function (TF) relates one input and one output: ( ) system Transfer Functons Convenent representaton of a lnear, dynamc model. A transfer functon (TF) relates one nput and one output: x t X s y t system Y s The followng termnology s used: x y nput output forcng

More information

Chapter 13: Multiple Regression

Chapter 13: Multiple Regression Chapter 13: Multple Regresson 13.1 Developng the multple-regresson Model The general model can be descrbed as: It smplfes for two ndependent varables: The sample ft parameter b 0, b 1, and b are used to

More information

Physics 5153 Classical Mechanics. D Alembert s Principle and The Lagrangian-1

Physics 5153 Classical Mechanics. D Alembert s Principle and The Lagrangian-1 P. Guterrez Physcs 5153 Classcal Mechancs D Alembert s Prncple and The Lagrangan 1 Introducton The prncple of vrtual work provdes a method of solvng problems of statc equlbrum wthout havng to consder the

More information

Graph Reconstruction by Permutations

Graph Reconstruction by Permutations Graph Reconstructon by Permutatons Perre Ille and Wllam Kocay* Insttut de Mathémathques de Lumny CNRS UMR 6206 163 avenue de Lumny, Case 907 13288 Marselle Cedex 9, France e-mal: lle@ml.unv-mrs.fr Computer

More information

Lecture 3: Probability Distributions

Lecture 3: Probability Distributions Lecture 3: Probablty Dstrbutons Random Varables Let us begn by defnng a sample space as a set of outcomes from an experment. We denote ths by S. A random varable s a functon whch maps outcomes nto the

More information

PHYS 705: Classical Mechanics. Calculus of Variations II

PHYS 705: Classical Mechanics. Calculus of Variations II 1 PHYS 705: Classcal Mechancs Calculus of Varatons II 2 Calculus of Varatons: Generalzaton (no constrant yet) Suppose now that F depends on several dependent varables : We need to fnd such that has a statonary

More information

VARIATION OF CONSTANT SUM CONSTRAINT FOR INTEGER MODEL WITH NON UNIFORM VARIABLES

VARIATION OF CONSTANT SUM CONSTRAINT FOR INTEGER MODEL WITH NON UNIFORM VARIABLES VARIATION OF CONSTANT SUM CONSTRAINT FOR INTEGER MODEL WITH NON UNIFORM VARIABLES BÂRZĂ, Slvu Faculty of Mathematcs-Informatcs Spru Haret Unversty barza_slvu@yahoo.com Abstract Ths paper wants to contnue

More information

Complete subgraphs in multipartite graphs

Complete subgraphs in multipartite graphs Complete subgraphs n multpartte graphs FLORIAN PFENDER Unverstät Rostock, Insttut für Mathematk D-18057 Rostock, Germany Floran.Pfender@un-rostock.de Abstract Turán s Theorem states that every graph G

More information

Introduction to Hidden Markov Models

Introduction to Hidden Markov Models Introducton to Hdden Markov Models Alperen Degrmenc Ths document contans dervatons and algorthms for mplementng Hdden Markov Models. The content presented here s a collecton of my notes and personal nsghts

More information

MASSACHUSETTS INSTITUTE OF TECHNOLOGY 6.265/15.070J Fall 2013 Lecture 12 10/21/2013. Martingale Concentration Inequalities and Applications

MASSACHUSETTS INSTITUTE OF TECHNOLOGY 6.265/15.070J Fall 2013 Lecture 12 10/21/2013. Martingale Concentration Inequalities and Applications MASSACHUSETTS INSTITUTE OF TECHNOLOGY 6.65/15.070J Fall 013 Lecture 1 10/1/013 Martngale Concentraton Inequaltes and Applcatons Content. 1. Exponental concentraton for martngales wth bounded ncrements.

More information

A be a probability space. A random vector

A be a probability space. A random vector Statstcs 1: Probablty Theory II 8 1 JOINT AND MARGINAL DISTRIBUTIONS In Probablty Theory I we formulate the concept of a (real) random varable and descrbe the probablstc behavor of ths random varable by

More information

Maximum Likelihood Estimation of Binary Dependent Variables Models: Probit and Logit. 1. General Formulation of Binary Dependent Variables Models

Maximum Likelihood Estimation of Binary Dependent Variables Models: Probit and Logit. 1. General Formulation of Binary Dependent Variables Models ECO 452 -- OE 4: Probt and Logt Models ECO 452 -- OE 4 Maxmum Lkelhood Estmaton of Bnary Dependent Varables Models: Probt and Logt hs note demonstrates how to formulate bnary dependent varables models

More information

Generalized Linear Methods

Generalized Linear Methods Generalzed Lnear Methods 1 Introducton In the Ensemble Methods the general dea s that usng a combnaton of several weak learner one could make a better learner. More formally, assume that we have a set

More information

Lecture 12: Classification

Lecture 12: Classification Lecture : Classfcaton g Dscrmnant functons g The optmal Bayes classfer g Quadratc classfers g Eucldean and Mahalanobs metrcs g K Nearest Neghbor Classfers Intellgent Sensor Systems Rcardo Guterrez-Osuna

More information

1. Inference on Regression Parameters a. Finding Mean, s.d and covariance amongst estimates. 2. Confidence Intervals and Working Hotelling Bands

1. Inference on Regression Parameters a. Finding Mean, s.d and covariance amongst estimates. 2. Confidence Intervals and Working Hotelling Bands Content. Inference on Regresson Parameters a. Fndng Mean, s.d and covarance amongst estmates.. Confdence Intervals and Workng Hotellng Bands 3. Cochran s Theorem 4. General Lnear Testng 5. Measures of

More information

Double Layered Fuzzy Planar Graph

Double Layered Fuzzy Planar Graph Global Journal of Pure and Appled Mathematcs. ISSN 0973-768 Volume 3, Number 0 07), pp. 7365-7376 Research Inda Publcatons http://www.rpublcaton.com Double Layered Fuzzy Planar Graph J. Jon Arockaraj Assstant

More information

Department of Computer Science Artificial Intelligence Research Laboratory. Iowa State University MACHINE LEARNING

Department of Computer Science Artificial Intelligence Research Laboratory. Iowa State University MACHINE LEARNING MACHINE LEANING Vasant Honavar Bonformatcs and Computatonal Bology rogram Center for Computatonal Intellgence, Learnng, & Dscovery Iowa State Unversty honavar@cs.astate.edu www.cs.astate.edu/~honavar/

More information

Perfect Competition and the Nash Bargaining Solution

Perfect Competition and the Nash Bargaining Solution Perfect Competton and the Nash Barganng Soluton Renhard John Department of Economcs Unversty of Bonn Adenauerallee 24-42 53113 Bonn, Germany emal: rohn@un-bonn.de May 2005 Abstract For a lnear exchange

More information

Computation of Higher Order Moments from Two Multinomial Overdispersion Likelihood Models

Computation of Higher Order Moments from Two Multinomial Overdispersion Likelihood Models Computaton of Hgher Order Moments from Two Multnomal Overdsperson Lkelhood Models BY J. T. NEWCOMER, N. K. NEERCHAL Department of Mathematcs and Statstcs, Unversty of Maryland, Baltmore County, Baltmore,

More information

Anti-van der Waerden numbers of 3-term arithmetic progressions.

Anti-van der Waerden numbers of 3-term arithmetic progressions. Ant-van der Waerden numbers of 3-term arthmetc progressons. Zhanar Berkkyzy, Alex Schulte, and Mchael Young Aprl 24, 2016 Abstract The ant-van der Waerden number, denoted by aw([n], k), s the smallest

More information