SlideShare a Scribd company logo
Learning LWF Chain Graphs:
A Markov Blanket Discovery Approach
Mohammad Ali Javidian
@ali javidian
Marco Valtorta
@MarcoGV2
Pooyan Jamshidi
@PooyanJamshidi
Department of Computer Science and Engineering
University of South Carolina
UAI 2020
1 / 9
LWF Chain Graphs (CGs)
Chain graphs:
admit both directed and undirected edges,
2 / 9
LWF Chain Graphs (CGs)
Chain graphs:
admit both directed and undirected edges,
there are no partially directed cycles.
2 / 9
LWF Chain Graphs (CGs)
Chain graphs:
admit both directed and undirected edges,
there are no partially directed cycles.
A Partially directed cycle: is a sequence of n distinct vertices
v1, v2, . . . , vn(n ≥ 3), and vn+1 ≡ v1, s.t.
for all i (1 ≤ i ≤ n) either vi − vi+1 or vi → vi+1, and
there exists a j (1 ≤ j ≤ n) such that vj ← vj+1.
2 / 9
LWF Chain Graphs (CGs)
Chain graphs:
admit both directed and undirected edges,
there are no partially directed cycles.
A Partially directed cycle: is a sequence of n distinct vertices
v1, v2, . . . , vn(n ≥ 3), and vn+1 ≡ v1, s.t.
for all i (1 ≤ i ≤ n) either vi − vi+1 or vi → vi+1, and
there exists a j (1 ≤ j ≤ n) such that vj ← vj+1.
Example:
2 / 9
LWF Chain Graphs (CGs)
Chain graphs:
admit both directed and undirected edges,
there are no partially directed cycles.
A Partially directed cycle: is a sequence of n distinct vertices
v1, v2, . . . , vn(n ≥ 3), and vn+1 ≡ v1, s.t.
for all i (1 ≤ i ≤ n) either vi − vi+1 or vi → vi+1, and
there exists a j (1 ≤ j ≤ n) such that vj ← vj+1.
Example:
Chain graphs under different interpretations: LWF, AMP, MVR, ...
Here, we focus on chain graphs (CGs) under the
Lauritzen-Wermuth-Frydenberg (LWF) interpretation.
2 / 9
Markov Blankets Enable Locality in Causal Structure Recovery
H
B
J
ADE
FC
I
OM N
K L
GT 𝑴𝒃(𝑇) H
B
J
ADE
FC
I
OM N
K L
GT
3 / 9
Markov Blankets Enable Locality in Causal Structure Recovery
H
B
J
ADE
FC
I
OM N
K L
GT 𝑴𝒃(𝑇) H
B
J
ADE
FC
I
OM N
K L
GT
3 / 9
Markov Blankets Enable Locality in Causal Structure Recovery
H
B
J
ADE
FC
I
OM N
K L
GT 𝑴𝒃(𝑇) H
B
J
ADE
FC
I
OM N
K L
GT
𝒄𝒉(𝑇)
3 / 9
Markov Blankets Enable Locality in Causal Structure Recovery
H
B
J
ADE
FC
I
OM N
K L
GT 𝑴𝒃(𝑇) H
B
J
ADE
FC
I
OM N
K L
GT
3 / 9
Markov Blankets Enable Locality in Causal Structure Recovery
H
B
J
ADE
FC
I
OM N
K L
GT 𝑴𝒃(𝑇) H
B
J
ADE
FC
I
OM N
K L
GT
Markov blankets can be used as a powerful tool in:
classification,
local causal discovery
3 / 9
Markov Blankets: a Missing Concept in the Context of Chain Graph Models
H
B
J
ADE
FC
I
OM N
K L
GT H
B
J
ADE
FC
I
OM N
K L
GT𝑴𝒃(𝑇)
4 / 9
Markov Blankets: a Missing Concept in the Context of Chain Graph Models
H
B
J
ADE
FC
I
OM N
K L
GT H
B
J
ADE
FC
I
OM N
K L
GT𝑴𝒃(𝑇)
𝒑𝒂(𝑇)
4 / 9
Markov Blankets: a Missing Concept in the Context of Chain Graph Models
H
B
J
ADE
FC
I
OM N
K L
GT H
B
J
ADE
FC
I
OM N
K L
GT𝑴𝒃(𝑇)
𝒄𝒉(𝑇)
4 / 9
Markov Blankets: a Missing Concept in the Context of Chain Graph Models
H
B
J
ADE
FC
I
OM N
K L
GT H
B
J
ADE
FC
I
OM N
K L
GT𝑴𝒃(𝑇)
𝒏𝒆(𝑇)
4 / 9
Markov Blankets: a Missing Concept in the Context of Chain Graph Models
H
B
J
ADE
FC
I
OM N
K L
GT H
B
J
ADE
FC
I
OM N
K L
GT𝑴𝒃(𝑇)
𝒄𝒔𝒑(𝑇)
4 / 9
Markov Blankets: a Missing Concept in the Context of Chain Graph Models
H
B
J
ADE
FC
I
OM N
K L
GT H
B
J
ADE
FC
I
OM N
K L
GT𝑴𝒃(𝑇)
4 / 9
Markov Blankets in LWF Chain Graphs: Main Results
Theorem (Characterization of Markov Blankets in LWF CGs)
Let G = (V, E, P) be an LWF chain graph model. Then, G entails conditional
independency T⊥⊥pV  {T, Mb(T)}|Mb(T), i.e., the Markov blanket of the target
variable T in an LWF CG probabilistically shields T from the rest of the variables.
5 / 9
Markov Blankets in LWF Chain Graphs: Main Results
Theorem (Characterization of Markov Blankets in LWF CGs)
Let G = (V, E, P) be an LWF chain graph model. Then, G entails conditional
independency T⊥⊥pV  {T, Mb(T)}|Mb(T), i.e., the Markov blanket of the target
variable T in an LWF CG probabilistically shields T from the rest of the variables.
Theorem (Standard algorithms for Markov blanket recovery in LWF CGs)
Given the Markov assumption, the faithfulness assumption, a graphical model
represented by an LWF CG, and i.i.d. sampling, in the large sample limit, the
Markov blanket recovery algorithms Grow-Shrink, Incremental Association
Markov blanket recovery, and its variants identify all Markov blankets for each
variable.
5 / 9
Markov Blankets in LWF Chain Graphs: Main Results
Theorem (Characterization of Markov Blankets in LWF CGs)
Let G = (V, E, P) be an LWF chain graph model. Then, G entails conditional
independency T⊥⊥pV  {T, Mb(T)}|Mb(T), i.e., the Markov blanket of the target
variable T in an LWF CG probabilistically shields T from the rest of the variables.
Theorem (Standard algorithms for Markov blanket recovery in LWF CGs)
Given the Markov assumption, the faithfulness assumption, a graphical model
represented by an LWF CG, and i.i.d. sampling, in the large sample limit, the
Markov blanket recovery algorithms Grow-Shrink, Incremental Association
Markov blanket recovery, and its variants identify all Markov blankets for each
variable.
The characterization of Markov blankets in chain graphs
enables us to develop new algorithms that are specifically
designed for learning Markov blankets in chain graphs.
5 / 9
MBC-CSP Algorithm: Markov Blanket Discovery in LWF CGs
H
B
J
A
E
FC
I
M
K L
GT
D
Grow
Phase: 
Step 1-1 
𝒂𝒅𝒋(𝑇)
H
B
J
A
E
FC
I
M
K L
GT
D
Grow
Phase: 
Step 1-2 
𝒂𝒅𝒋(𝑇)
H
B
J
A
E
FC
I
M
K L
GT
D
Grow
Phase: 
Step 2 
𝐜𝐬𝐩(𝑇)
6 / 9
MBC-CSP Algorithm: Markov Blanket Discovery in LWF CGs
H
B
J
A
E
FC
I
M
K L
GT
D
Grow
Phase: 
Step 1-1 
𝒂𝒅𝒋(𝑇)
H
B
J
A
E
FC
I
M
K L
GT
D
Grow
Phase: 
Step 1-2 
𝒂𝒅𝒋(𝑇)
H
B
J
A
E
FC
I
M
K L
GT
D
Grow
Phase: 
Step 2 
𝐜𝐬𝐩(𝑇)
H
B
J
A
E
FC
I
M
K L
GT
D
Grow
Phase: 
Step 1-1 
𝒂𝒅𝒋(𝑇)
H
B
J
A
E
FC
I
M
K L
GT
D
Grow
Phase: 
Step 1-2 
𝒂𝒅𝒋(𝑇)
H
B
J
A
E
FC
I
M
K L
GT
D
Grow
Phase: 
Step 2 
𝐜𝐬𝐩(𝑇)
6 / 9
MBC-CSP Algorithm: Markov Blanket Discovery in LWF CGs
H
B
J
A
E
FC
I
M
K L
GT
D
Grow
Phase: 
Step 1-1 
𝒂𝒅𝒋(𝑇)
H
B
J
A
E
FC
I
M
K L
GT
D
Grow
Phase: 
Step 1-2 
𝒂𝒅𝒋(𝑇)
H
B
J
A
E
FC
I
M
K L
GT
D
Grow
Phase: 
Step 2 
𝐜𝐬𝐩(𝑇)
H
B
J
A
E
FC
I
M
K L
GT
D
Grow
Phase: 
Step 1-1 
𝒂𝒅𝒋(𝑇)
H
B
J
A
E
FC
I
M
K L
GT
D
Grow
Phase: 
Step 1-2 
𝒂𝒅𝒋(𝑇)
H
B
J
A
E
FC
I
M
K L
GT
D
Grow
Phase: 
Step 2 
𝐜𝐬𝐩(𝑇)
Shrink
Phase:
𝐌𝐛(𝑇) H
B
J
A
E
FC
I
M
K L
GT
D
H
B
J
A
E
FC
I
M
K L
GT
D
6 / 9
MBC-CSP Algorithm: Markov Blanket Discovery in LWF CGs
H
B
J
A
E
FC
I
M
K L
GT
D
Grow
Phase: 
Step 1-1 
𝒂𝒅𝒋(𝑇)
H
B
J
A
E
FC
I
M
K L
GT
D
Grow
Phase: 
Step 1-2 
𝒂𝒅𝒋(𝑇)
H
B
J
A
E
FC
I
M
K L
GT
D
Grow
Phase: 
Step 2 
𝐜𝐬𝐩(𝑇)
H
B
J
A
E
FC
I
M
K L
GT
D
Grow
Phase: 
Step 1-1 
𝒂𝒅𝒋(𝑇)
H
B
J
A
E
FC
I
M
K L
GT
D
Grow
Phase: 
Step 1-2 
𝒂𝒅𝒋(𝑇)
H
B
J
A
E
FC
I
M
K L
GT
D
Grow
Phase: 
Step 2 
𝐜𝐬𝐩(𝑇)
Shrink
Phase:
𝐌𝐛(𝑇) H
B
J
A
E
FC
I
M
K L
GT
D
H
B
J
A
E
FC
I
M
K L
GT
D
Such Markov blanket discovery algorithms help
us to design new scalable algorithms for learning
chain graphs based on local structure discovery.
6 / 9
MbLWF Algorithm: Learning LWF CGs via Markov Blanket Discovery
Observational Data
Markov Blanket
Discovery Algorithm
A
DC E
B
𝐌𝐛 𝐴 = 𝐶, 𝐵 , 𝐒𝐞𝐩𝐬𝐞𝐭 𝐴, 𝐷 = 𝐶, 𝐵 , 𝐒𝐞𝐩𝐬𝐞𝐭 𝐴, 𝐸 = {𝐶, 𝐵}
𝐌𝐛 𝐵 = 𝐸, 𝐴 , 𝐒𝐞𝐩𝐬𝐞𝐭 𝐵, 𝐶 = 𝐸, 𝐴 , 𝐒𝐞𝐩𝐬𝐞𝐭 𝐵, 𝐷 = {𝐸, 𝐴}
𝐌𝐛 𝐶 = 𝐴, 𝐷 , 𝐒𝐞𝐩𝐬𝐞𝐭 𝐶, 𝐸 = 𝐴, 𝐷 , 𝐒𝐞𝐩𝐬𝐞𝐭 𝐶, 𝐵 = {𝐴, 𝐷}
𝐌𝐛 𝐷 = 𝐶, 𝐸 , 𝐒𝐞𝐩𝐬𝐞𝐭 𝐷, 𝐴 = 𝐶, 𝐵 , 𝐒𝐞𝐩𝐬𝐞𝐭 𝐷, 𝐵 = {𝐸, 𝐴}
𝐌𝐛 𝐸 = 𝐷, 𝐵 , 𝐒𝐞𝐩𝐬𝐞𝐭 𝐸, 𝐶 = 𝐴, 𝐷 , 𝐒𝐞𝐩𝐬𝐞𝐭 𝐸, 𝐴 = {𝐶, 𝐵}
7 / 9
MbLWF Algorithm: Learning LWF CGs via Markov Blanket Discovery
Observational Data
Markov Blanket
Discovery Algorithm
A
DC E
B
𝐌𝐛 𝐴 = 𝐶, 𝐵 , 𝐒𝐞𝐩𝐬𝐞𝐭 𝐴, 𝐷 = 𝐶, 𝐵 , 𝐒𝐞𝐩𝐬𝐞𝐭 𝐴, 𝐸 = {𝐶, 𝐵}
𝐌𝐛 𝐵 = 𝐸, 𝐴 , 𝐒𝐞𝐩𝐬𝐞𝐭 𝐵, 𝐶 = 𝐸, 𝐴 , 𝐒𝐞𝐩𝐬𝐞𝐭 𝐵, 𝐷 = {𝐸, 𝐴}
𝐌𝐛 𝐶 = 𝐴, 𝐷 , 𝐒𝐞𝐩𝐬𝐞𝐭 𝐶, 𝐸 = 𝐴, 𝐷 , 𝐒𝐞𝐩𝐬𝐞𝐭 𝐶, 𝐵 = {𝐴, 𝐷}
𝐌𝐛 𝐷 = 𝐶, 𝐸 , 𝐒𝐞𝐩𝐬𝐞𝐭 𝐷, 𝐴 = 𝐶, 𝐵 , 𝐒𝐞𝐩𝐬𝐞𝐭 𝐷, 𝐵 = {𝐸, 𝐴}
𝐌𝐛 𝐸 = 𝐷, 𝐵 , 𝐒𝐞𝐩𝐬𝐞𝐭 𝐸, 𝐶 = 𝐴, 𝐷 , 𝐒𝐞𝐩𝐬𝐞𝐭 𝐸, 𝐴 = {𝐶, 𝐵}
A
DC E
B
Super Skeleton
Recovery
7 / 9
MbLWF Algorithm: Learning LWF CGs via Markov Blanket Discovery
Observational Data
Markov Blanket
Discovery Algorithm
A
DC E
B
𝐌𝐛 𝐴 = 𝐶, 𝐵 , 𝐒𝐞𝐩𝐬𝐞𝐭 𝐴, 𝐷 = 𝐶, 𝐵 , 𝐒𝐞𝐩𝐬𝐞𝐭 𝐴, 𝐸 = {𝐶, 𝐵}
𝐌𝐛 𝐵 = 𝐸, 𝐴 , 𝐒𝐞𝐩𝐬𝐞𝐭 𝐵, 𝐶 = 𝐸, 𝐴 , 𝐒𝐞𝐩𝐬𝐞𝐭 𝐵, 𝐷 = {𝐸, 𝐴}
𝐌𝐛 𝐶 = 𝐴, 𝐷 , 𝐒𝐞𝐩𝐬𝐞𝐭 𝐶, 𝐸 = 𝐴, 𝐷 , 𝐒𝐞𝐩𝐬𝐞𝐭 𝐶, 𝐵 = {𝐴, 𝐷}
𝐌𝐛 𝐷 = 𝐶, 𝐸 , 𝐒𝐞𝐩𝐬𝐞𝐭 𝐷, 𝐴 = 𝐶, 𝐵 , 𝐒𝐞𝐩𝐬𝐞𝐭 𝐷, 𝐵 = {𝐸, 𝐴}
𝐌𝐛 𝐸 = 𝐷, 𝐵 , 𝐒𝐞𝐩𝐬𝐞𝐭 𝐸, 𝐶 = 𝐴, 𝐷 , 𝐒𝐞𝐩𝐬𝐞𝐭 𝐸, 𝐴 = {𝐶, 𝐵}
A
DC E
B
Super Skeleton
Recovery
𝐒𝐞𝐩𝐬𝐞𝐭 𝐴, 𝐷 = 𝐶, 𝐵
Skeleton
A
DC E
B
7 / 9
MbLWF Algorithm: Learning LWF CGs via Markov Blanket Discovery
Observational Data
Markov Blanket
Discovery Algorithm
A
DC E
B
𝐌𝐛 𝐴 = 𝐶, 𝐵 , 𝐒𝐞𝐩𝐬𝐞𝐭 𝐴, 𝐷 = 𝐶, 𝐵 , 𝐒𝐞𝐩𝐬𝐞𝐭 𝐴, 𝐸 = {𝐶, 𝐵}
𝐌𝐛 𝐵 = 𝐸, 𝐴 , 𝐒𝐞𝐩𝐬𝐞𝐭 𝐵, 𝐶 = 𝐸, 𝐴 , 𝐒𝐞𝐩𝐬𝐞𝐭 𝐵, 𝐷 = {𝐸, 𝐴}
𝐌𝐛 𝐶 = 𝐴, 𝐷 , 𝐒𝐞𝐩𝐬𝐞𝐭 𝐶, 𝐸 = 𝐴, 𝐷 , 𝐒𝐞𝐩𝐬𝐞𝐭 𝐶, 𝐵 = {𝐴, 𝐷}
𝐌𝐛 𝐷 = 𝐶, 𝐸 , 𝐒𝐞𝐩𝐬𝐞𝐭 𝐷, 𝐴 = 𝐶, 𝐵 , 𝐒𝐞𝐩𝐬𝐞𝐭 𝐷, 𝐵 = {𝐸, 𝐴}
𝐌𝐛 𝐸 = 𝐷, 𝐵 , 𝐒𝐞𝐩𝐬𝐞𝐭 𝐸, 𝐶 = 𝐴, 𝐷 , 𝐒𝐞𝐩𝐬𝐞𝐭 𝐸, 𝐴 = {𝐶, 𝐵}
A
DC E
B
Super Skeleton
Recovery
𝐒𝐞𝐩𝐬𝐞𝐭 𝐴, 𝐷 = 𝐶, 𝐵
Skeleton
A
DC E
B
A
DC E
BComplex
Recovery
7 / 9
Experimental Evaluation: Markov Blankets Make a Broad Range of
Inference/Learning Problems Computationally Tractable and More Precise.
●●
●●
●
●
●
●
●
●
●
●
● ●
●
●
●
●
●
●
●
●0.90
0.95
1.00 GSLWF
fastIAMBLWF
fdrIAMBLWF
interIAMBLWF
IAMBLWF
MBCCSPLWF
LCD
Precision
sample size
size = 200
size = 2000
alpha = 0.05
8 / 9
Experimental Evaluation: Markov Blankets Make a Broad Range of
Inference/Learning Problems Computationally Tractable and More Precise.
●●
●●
●
●
●
●
●
●
●
●
● ●
●
●
●
●
●
●
●
●0.90
0.95
1.00 GSLWF
fastIAMBLWF
fdrIAMBLWF
interIAMBLWF
IAMBLWF
MBCCSPLWF
LCD
Precision
sample size
size = 200
size = 2000
alpha = 0.05
● ● ● ●
●
●
0.5
0.6
0.7
0.8
0.9
1.0
GSLWF
fastIAMBLWF
fdrIAMBLWF
interIAMBLWF
IAMBLWF
MBCCSPLWF
LCD
Recall
alpha = 0.05
8 / 9
Markov Blankets: a Missing Concept in the Context of Chain Graph Models
H
B
J
ADE
FC
I
OM N
K L
GT H
B
J
ADE
FC
I
OM N
K L
GT𝑴𝒃(𝑇)
4 / 9
Markov Blankets: a Missing Concept in the Context of Chain Graph Models
H
B
J
ADE
FC
I
OM N
K L
GT H
B
J
ADE
FC
I
OM N
K L
GT𝑴𝒃(𝑇)
4 / 9
Markov Blankets in LWF Chain Graphs: Main Results
Theorem
Let G = (V, E, P) be an LWF chain graph model. Then, G entails conditional
independency T⊥⊥pV  {T, Mb(T)}|Mb(T).
Theorem
Given the Markov assumption, the faithfulness assumption, a graphical model
represented by an LWF CG, and i.i.d. sampling, in the large sample limit, the
Markov blanket recovery algorithms Grow-Shrink, Incremental Association
Markov blanket recovery, and its variants identify all Markov blankets for each
variable.
5 / 9
Markov Blankets: a Missing Concept in the Context of Chain Graph Models
H
B
J
ADE
FC
I
OM N
K L
GT H
B
J
ADE
FC
I
OM N
K L
GT𝑴𝒃(𝑇)
4 / 9
Markov Blankets in LWF Chain Graphs: Main Results
Theorem
Let G = (V, E, P) be an LWF chain graph model. Then, G entails conditional
independency T⊥⊥pV  {T, Mb(T)}|Mb(T).
Theorem
Given the Markov assumption, the faithfulness assumption, a graphical model
represented by an LWF CG, and i.i.d. sampling, in the large sample limit, the
Markov blanket recovery algorithms Grow-Shrink, Incremental Association
Markov blanket recovery, and its variants identify all Markov blankets for each
variable.
5 / 9
MBC-CSP Algorithm: Markov Blanket Discovery in LWF CGs
H
B
J
A
E
FC
I
M
K L
GT
D
Grow
Phase: 
Step 1-1 
𝒂𝒅𝒋(𝑇)
H
B
J
A
E
FC
I
M
K L
GT
D
Grow
Phase: 
Step 1-2 
𝒂𝒅𝒋(𝑇)
H
B
J
A
E
FC
I
M
K L
GT
D
Grow
Phase: 
Step 2 
𝐜𝐬𝐩(𝑇)
H
B
J
A
E
FC
I
M
K L
GT
D
Grow
Phase: 
Step 1-1 
𝒂𝒅𝒋(𝑇)
H
B
J
A
E
FC
I
M
K L
GT
D
Grow
Phase: 
Step 1-2 
𝒂𝒅𝒋(𝑇)
H
B
J
A
E
FC
I
M
K L
GT
D
Grow
Phase: 
Step 2 
𝐜𝐬𝐩(𝑇)
Shrink
Phase:
𝐌𝐛(𝑇) H
B
J
A
E
FC
I
M
K L
GT
D
H
B
J
A
E
FC
I
M
K L
GT
D
6 / 9
Markov Blankets: a Missing Concept in the Context of Chain Graph Models
H
B
J
ADE
FC
I
OM N
K L
GT H
B
J
ADE
FC
I
OM N
K L
GT𝑴𝒃(𝑇)
4 / 9
Markov Blankets in LWF Chain Graphs: Main Results
Theorem
Let G = (V, E, P) be an LWF chain graph model. Then, G entails conditional
independency T⊥⊥pV  {T, Mb(T)}|Mb(T).
Theorem
Given the Markov assumption, the faithfulness assumption, a graphical model
represented by an LWF CG, and i.i.d. sampling, in the large sample limit, the
Markov blanket recovery algorithms Grow-Shrink, Incremental Association
Markov blanket recovery, and its variants identify all Markov blankets for each
variable.
5 / 9
MBC-CSP Algorithm: Markov Blanket Discovery in LWF CGs
H
B
J
A
E
FC
I
M
K L
GT
D
Grow
Phase: 
Step 1-1 
𝒂𝒅𝒋(𝑇)
H
B
J
A
E
FC
I
M
K L
GT
D
Grow
Phase: 
Step 1-2 
𝒂𝒅𝒋(𝑇)
H
B
J
A
E
FC
I
M
K L
GT
D
Grow
Phase: 
Step 2 
𝐜𝐬𝐩(𝑇)
H
B
J
A
E
FC
I
M
K L
GT
D
Grow
Phase: 
Step 1-1 
𝒂𝒅𝒋(𝑇)
H
B
J
A
E
FC
I
M
K L
GT
D
Grow
Phase: 
Step 1-2 
𝒂𝒅𝒋(𝑇)
H
B
J
A
E
FC
I
M
K L
GT
D
Grow
Phase: 
Step 2 
𝐜𝐬𝐩(𝑇)
Shrink
Phase:
𝐌𝐛(𝑇) H
B
J
A
E
FC
I
M
K L
GT
D
H
B
J
A
E
FC
I
M
K L
GT
D
6 / 9
MbLWF Algorithm: Learning LWF CGs via Markov Blanket Discovery
Observational Data
Markov Blanket
Discovery Algorithm
A
DC E
B
𝐌𝐛 𝐴 = 𝐶, 𝐵 , 𝐒𝐞𝐩𝐬𝐞𝐭 𝐴, 𝐷 = 𝐶, 𝐵 , 𝐒𝐞𝐩𝐬𝐞𝐭 𝐴, 𝐸 = {𝐶, 𝐵}
𝐌𝐛 𝐵 = 𝐸, 𝐴 , 𝐒𝐞𝐩𝐬𝐞𝐭 𝐵, 𝐶 = 𝐸, 𝐴 , 𝐒𝐞𝐩𝐬𝐞𝐭 𝐵, 𝐷 = {𝐸, 𝐴}
𝐌𝐛 𝐶 = 𝐴, 𝐷 , 𝐒𝐞𝐩𝐬𝐞𝐭 𝐶, 𝐸 = 𝐴, 𝐷 , 𝐒𝐞𝐩𝐬𝐞𝐭 𝐶, 𝐵 = {𝐴, 𝐷}
𝐌𝐛 𝐷 = 𝐶, 𝐸 , 𝐒𝐞𝐩𝐬𝐞𝐭 𝐷, 𝐴 = 𝐶, 𝐵 , 𝐒𝐞𝐩𝐬𝐞𝐭 𝐷, 𝐵 = {𝐸, 𝐴}
𝐌𝐛 𝐸 = 𝐷, 𝐵 , 𝐒𝐞𝐩𝐬𝐞𝐭 𝐸, 𝐶 = 𝐴, 𝐷 , 𝐒𝐞𝐩𝐬𝐞𝐭 𝐸, 𝐴 = {𝐶, 𝐵}
A
DC E
B
Super Skeleton
Recovery
𝐒𝐞𝐩𝐬𝐞𝐭 𝐴, 𝐷 = 𝐶, 𝐵
Skeleton
A
DC E
B
A
DC E
BComplex
Recovery
7 / 9
Markov Blankets: a Missing Concept in the Context of Chain Graph Models
H
B
J
ADE
FC
I
OM N
K L
GT H
B
J
ADE
FC
I
OM N
K L
GT𝑴𝒃(𝑇)
4 / 9
Markov Blankets in LWF Chain Graphs: Main Results
Theorem
Let G = (V, E, P) be an LWF chain graph model. Then, G entails conditional
independency T⊥⊥pV  {T, Mb(T)}|Mb(T).
Theorem
Given the Markov assumption, the faithfulness assumption, a graphical model
represented by an LWF CG, and i.i.d. sampling, in the large sample limit, the
Markov blanket recovery algorithms Grow-Shrink, Incremental Association
Markov blanket recovery, and its variants identify all Markov blankets for each
variable.
5 / 9
MBC-CSP Algorithm: Markov Blanket Discovery in LWF CGs
H
B
J
A
E
FC
I
M
K L
GT
D
Grow
Phase: 
Step 1-1 
𝒂𝒅𝒋(𝑇)
H
B
J
A
E
FC
I
M
K L
GT
D
Grow
Phase: 
Step 1-2 
𝒂𝒅𝒋(𝑇)
H
B
J
A
E
FC
I
M
K L
GT
D
Grow
Phase: 
Step 2 
𝐜𝐬𝐩(𝑇)
H
B
J
A
E
FC
I
M
K L
GT
D
Grow
Phase: 
Step 1-1 
𝒂𝒅𝒋(𝑇)
H
B
J
A
E
FC
I
M
K L
GT
D
Grow
Phase: 
Step 1-2 
𝒂𝒅𝒋(𝑇)
H
B
J
A
E
FC
I
M
K L
GT
D
Grow
Phase: 
Step 2 
𝐜𝐬𝐩(𝑇)
Shrink
Phase:
𝐌𝐛(𝑇) H
B
J
A
E
FC
I
M
K L
GT
D
H
B
J
A
E
FC
I
M
K L
GT
D
6 / 9
MbLWF Algorithm: Learning LWF CGs via Markov Blanket Discovery
Observational Data
Markov Blanket
Discovery Algorithm
A
DC E
B
𝐌𝐛 𝐴 = 𝐶, 𝐵 , 𝐒𝐞𝐩𝐬𝐞𝐭 𝐴, 𝐷 = 𝐶, 𝐵 , 𝐒𝐞𝐩𝐬𝐞𝐭 𝐴, 𝐸 = {𝐶, 𝐵}
𝐌𝐛 𝐵 = 𝐸, 𝐴 , 𝐒𝐞𝐩𝐬𝐞𝐭 𝐵, 𝐶 = 𝐸, 𝐴 , 𝐒𝐞𝐩𝐬𝐞𝐭 𝐵, 𝐷 = {𝐸, 𝐴}
𝐌𝐛 𝐶 = 𝐴, 𝐷 , 𝐒𝐞𝐩𝐬𝐞𝐭 𝐶, 𝐸 = 𝐴, 𝐷 , 𝐒𝐞𝐩𝐬𝐞𝐭 𝐶, 𝐵 = {𝐴, 𝐷}
𝐌𝐛 𝐷 = 𝐶, 𝐸 , 𝐒𝐞𝐩𝐬𝐞𝐭 𝐷, 𝐴 = 𝐶, 𝐵 , 𝐒𝐞𝐩𝐬𝐞𝐭 𝐷, 𝐵 = {𝐸, 𝐴}
𝐌𝐛 𝐸 = 𝐷, 𝐵 , 𝐒𝐞𝐩𝐬𝐞𝐭 𝐸, 𝐶 = 𝐴, 𝐷 , 𝐒𝐞𝐩𝐬𝐞𝐭 𝐸, 𝐴 = {𝐶, 𝐵}
A
DC E
B
Super Skeleton
Recovery
𝐒𝐞𝐩𝐬𝐞𝐭 𝐴, 𝐷 = 𝐶, 𝐵
Skeleton
A
DC E
B
A
DC E
BComplex
Recovery
7 / 9
All code, data, and supplementary materials are available at:
https://majavid.github.io/structurelearning/blog/2020/uai/

More Related Content

PDF
Libxc a library of exchange and correlation functionals
PDF
The all-electron GW method based on WIEN2k: Implementation and applications.
PDF
A Closed-Form Expression for Queuing Delay in Rayleigh Fading Channels Using ...
PDF
Presentation on Probability Genrating Function
PDF
Guarding Polygons via CSP
PDF
Graph Kernelpdf
PDF
Metodo Monte Carlo -Wang Landau
PDF
Polynomial Kernel for Interval Vertex Deletion
Libxc a library of exchange and correlation functionals
The all-electron GW method based on WIEN2k: Implementation and applications.
A Closed-Form Expression for Queuing Delay in Rayleigh Fading Channels Using ...
Presentation on Probability Genrating Function
Guarding Polygons via CSP
Graph Kernelpdf
Metodo Monte Carlo -Wang Landau
Polynomial Kernel for Interval Vertex Deletion

What's hot (19)

PDF
Kernel for Chordal Vertex Deletion
PDF
Wien2k getting started
PDF
Guarding Terrains though the Lens of Parameterized Complexity
PDF
Mathematics Colloquium, UCSC
PDF
Graph Modification: Beyond the known Boundaries
PDF
Fine Grained Complexity
DOCX
Question bank of module iv packet switching networks
PDF
Message Passing
PDF
vasp-gpu on Balena: Usage and Some Benchmarks
PDF
Fine Grained Complexity of Rainbow Coloring and its Variants
PDF
Split Contraction: The Untold Story
PDF
1 s2.0-s2090447917300291-main
PDF
Phonons & Phonopy: Pro Tips (2015)
PDF
Ultrasound Modular Architecture
PDF
Research Inventy : International Journal of Engineering and Science is publis...
PPTX
Blaha krakow 2004
PDF
ECE 467 Mini project 1
PDF
Density functional theory (DFT) and the concepts of the augmented-plane-wave ...
PDF
Bayesian Inference and Uncertainty Quantification for Inverse Problems
Kernel for Chordal Vertex Deletion
Wien2k getting started
Guarding Terrains though the Lens of Parameterized Complexity
Mathematics Colloquium, UCSC
Graph Modification: Beyond the known Boundaries
Fine Grained Complexity
Question bank of module iv packet switching networks
Message Passing
vasp-gpu on Balena: Usage and Some Benchmarks
Fine Grained Complexity of Rainbow Coloring and its Variants
Split Contraction: The Untold Story
1 s2.0-s2090447917300291-main
Phonons & Phonopy: Pro Tips (2015)
Ultrasound Modular Architecture
Research Inventy : International Journal of Engineering and Science is publis...
Blaha krakow 2004
ECE 467 Mini project 1
Density functional theory (DFT) and the concepts of the augmented-plane-wave ...
Bayesian Inference and Uncertainty Quantification for Inverse Problems
Ad

Similar to Learning LWF Chain Graphs: A Markov Blanket Discovery Approach (20)

PPTX
Algorithm_NP-Completeness Proof
PDF
Learning Convolutional Neural Networks for Graphs
PDF
System overflow blocking-transients-for-queues-with-batch-arrivals-using-a-fa...
PDF
SLAM of Multi-Robot System Considering Its Network Topology
PPTX
Lec5(Polynomial Regression) & Support vector regression.pptx
PDF
PPTX
Krunal_lbm.pptx
PDF
TMPA-2015: Implementing the MetaVCG Approach in the C-light System
PPT
about power system operation and control13197214.ppt
PPTX
Generative models
PDF
lecture6.pdf
DOC
Report-Implementation of Quantum Gates using Verilog
PDF
Node Unique Label Cover
PDF
Extracting linear information from 
nonlinear large-scale structure observations
PDF
Section 4 3_the_scattering_matrix_package
PDF
Tuning of PID, SVFB and LQ Controllers Using Genetic Algorithms
PDF
Paper.pdf
PDF
Iclr2016 vaeまとめ
PPT
NIPS2007: structured prediction
Algorithm_NP-Completeness Proof
Learning Convolutional Neural Networks for Graphs
System overflow blocking-transients-for-queues-with-batch-arrivals-using-a-fa...
SLAM of Multi-Robot System Considering Its Network Topology
Lec5(Polynomial Regression) & Support vector regression.pptx
Krunal_lbm.pptx
TMPA-2015: Implementing the MetaVCG Approach in the C-light System
about power system operation and control13197214.ppt
Generative models
lecture6.pdf
Report-Implementation of Quantum Gates using Verilog
Node Unique Label Cover
Extracting linear information from 
nonlinear large-scale structure observations
Section 4 3_the_scattering_matrix_package
Tuning of PID, SVFB and LQ Controllers Using Genetic Algorithms
Paper.pdf
Iclr2016 vaeまとめ
NIPS2007: structured prediction
Ad

More from Pooyan Jamshidi (20)

PDF
A Framework for Robust Control of Uncertainty in Self-Adaptive Software Conn...
PDF
Machine Learning Meets Quantitative Planning: Enabling Self-Adaptation in Aut...
PDF
Ensembles of Many Diverse Weak Defenses can be Strong: Defending Deep Neural ...
PDF
Transfer Learning for Performance Analysis of Machine Learning Systems
PDF
Transfer Learning for Performance Analysis of Configurable Systems: A Causal ...
PDF
Machine Learning meets DevOps
PDF
Learning to Sample
PDF
Integrated Model Discovery and Self-Adaptation of Robots
PDF
Transfer Learning for Performance Analysis of Highly-Configurable Software
PDF
Architectural Tradeoff in Learning-Based Software
PDF
Production-Ready Machine Learning for the Software Architect
PDF
Transfer Learning for Software Performance Analysis: An Exploratory Analysis
PDF
Architecting for Scale
PDF
Learning Software Performance Models for Dynamic and Uncertain Environments
PDF
Sensitivity Analysis for Building Adaptive Robotic Software
PDF
Transfer Learning for Improving Model Predictions in Highly Configurable Soft...
PDF
Transfer Learning for Improving Model Predictions in Robotic Systems
PDF
Machine Learning meets DevOps
PDF
An Uncertainty-Aware Approach to Optimal Configuration of Stream Processing S...
PDF
Configuration Optimization Tool
A Framework for Robust Control of Uncertainty in Self-Adaptive Software Conn...
Machine Learning Meets Quantitative Planning: Enabling Self-Adaptation in Aut...
Ensembles of Many Diverse Weak Defenses can be Strong: Defending Deep Neural ...
Transfer Learning for Performance Analysis of Machine Learning Systems
Transfer Learning for Performance Analysis of Configurable Systems: A Causal ...
Machine Learning meets DevOps
Learning to Sample
Integrated Model Discovery and Self-Adaptation of Robots
Transfer Learning for Performance Analysis of Highly-Configurable Software
Architectural Tradeoff in Learning-Based Software
Production-Ready Machine Learning for the Software Architect
Transfer Learning for Software Performance Analysis: An Exploratory Analysis
Architecting for Scale
Learning Software Performance Models for Dynamic and Uncertain Environments
Sensitivity Analysis for Building Adaptive Robotic Software
Transfer Learning for Improving Model Predictions in Highly Configurable Soft...
Transfer Learning for Improving Model Predictions in Robotic Systems
Machine Learning meets DevOps
An Uncertainty-Aware Approach to Optimal Configuration of Stream Processing S...
Configuration Optimization Tool

Recently uploaded (20)

PDF
GENETICS IN BIOLOGY IN SECONDARY LEVEL FORM 3
PPTX
Institutional Correction lecture only . . .
PDF
FourierSeries-QuestionsWithAnswers(Part-A).pdf
PDF
O7-L3 Supply Chain Operations - ICLT Program
PDF
STATICS OF THE RIGID BODIES Hibbelers.pdf
PDF
Abdominal Access Techniques with Prof. Dr. R K Mishra
PPTX
Introduction-to-Literarature-and-Literary-Studies-week-Prelim-coverage.pptx
PPTX
GDM (1) (1).pptx small presentation for students
PDF
O5-L3 Freight Transport Ops (International) V1.pdf
PPTX
Cell Types and Its function , kingdom of life
PDF
3rd Neelam Sanjeevareddy Memorial Lecture.pdf
PPTX
Final Presentation General Medicine 03-08-2024.pptx
PDF
Complications of Minimal Access Surgery at WLH
PDF
Saundersa Comprehensive Review for the NCLEX-RN Examination.pdf
PDF
Supply Chain Operations Speaking Notes -ICLT Program
PDF
102 student loan defaulters named and shamed – Is someone you know on the list?
PDF
The Lost Whites of Pakistan by Jahanzaib Mughal.pdf
PPTX
Microbial diseases, their pathogenesis and prophylaxis
PPTX
1st Inaugural Professorial Lecture held on 19th February 2020 (Governance and...
PDF
2.FourierTransform-ShortQuestionswithAnswers.pdf
GENETICS IN BIOLOGY IN SECONDARY LEVEL FORM 3
Institutional Correction lecture only . . .
FourierSeries-QuestionsWithAnswers(Part-A).pdf
O7-L3 Supply Chain Operations - ICLT Program
STATICS OF THE RIGID BODIES Hibbelers.pdf
Abdominal Access Techniques with Prof. Dr. R K Mishra
Introduction-to-Literarature-and-Literary-Studies-week-Prelim-coverage.pptx
GDM (1) (1).pptx small presentation for students
O5-L3 Freight Transport Ops (International) V1.pdf
Cell Types and Its function , kingdom of life
3rd Neelam Sanjeevareddy Memorial Lecture.pdf
Final Presentation General Medicine 03-08-2024.pptx
Complications of Minimal Access Surgery at WLH
Saundersa Comprehensive Review for the NCLEX-RN Examination.pdf
Supply Chain Operations Speaking Notes -ICLT Program
102 student loan defaulters named and shamed – Is someone you know on the list?
The Lost Whites of Pakistan by Jahanzaib Mughal.pdf
Microbial diseases, their pathogenesis and prophylaxis
1st Inaugural Professorial Lecture held on 19th February 2020 (Governance and...
2.FourierTransform-ShortQuestionswithAnswers.pdf

Learning LWF Chain Graphs: A Markov Blanket Discovery Approach

  • 1. Learning LWF Chain Graphs: A Markov Blanket Discovery Approach Mohammad Ali Javidian @ali javidian Marco Valtorta @MarcoGV2 Pooyan Jamshidi @PooyanJamshidi Department of Computer Science and Engineering University of South Carolina UAI 2020 1 / 9
  • 2. LWF Chain Graphs (CGs) Chain graphs: admit both directed and undirected edges, 2 / 9
  • 3. LWF Chain Graphs (CGs) Chain graphs: admit both directed and undirected edges, there are no partially directed cycles. 2 / 9
  • 4. LWF Chain Graphs (CGs) Chain graphs: admit both directed and undirected edges, there are no partially directed cycles. A Partially directed cycle: is a sequence of n distinct vertices v1, v2, . . . , vn(n ≥ 3), and vn+1 ≡ v1, s.t. for all i (1 ≤ i ≤ n) either vi − vi+1 or vi → vi+1, and there exists a j (1 ≤ j ≤ n) such that vj ← vj+1. 2 / 9
  • 5. LWF Chain Graphs (CGs) Chain graphs: admit both directed and undirected edges, there are no partially directed cycles. A Partially directed cycle: is a sequence of n distinct vertices v1, v2, . . . , vn(n ≥ 3), and vn+1 ≡ v1, s.t. for all i (1 ≤ i ≤ n) either vi − vi+1 or vi → vi+1, and there exists a j (1 ≤ j ≤ n) such that vj ← vj+1. Example: 2 / 9
  • 6. LWF Chain Graphs (CGs) Chain graphs: admit both directed and undirected edges, there are no partially directed cycles. A Partially directed cycle: is a sequence of n distinct vertices v1, v2, . . . , vn(n ≥ 3), and vn+1 ≡ v1, s.t. for all i (1 ≤ i ≤ n) either vi − vi+1 or vi → vi+1, and there exists a j (1 ≤ j ≤ n) such that vj ← vj+1. Example: Chain graphs under different interpretations: LWF, AMP, MVR, ... Here, we focus on chain graphs (CGs) under the Lauritzen-Wermuth-Frydenberg (LWF) interpretation. 2 / 9
  • 7. Markov Blankets Enable Locality in Causal Structure Recovery H B J ADE FC I OM N K L GT 𝑴𝒃(𝑇) H B J ADE FC I OM N K L GT 3 / 9
  • 8. Markov Blankets Enable Locality in Causal Structure Recovery H B J ADE FC I OM N K L GT 𝑴𝒃(𝑇) H B J ADE FC I OM N K L GT 3 / 9
  • 9. Markov Blankets Enable Locality in Causal Structure Recovery H B J ADE FC I OM N K L GT 𝑴𝒃(𝑇) H B J ADE FC I OM N K L GT 𝒄𝒉(𝑇) 3 / 9
  • 10. Markov Blankets Enable Locality in Causal Structure Recovery H B J ADE FC I OM N K L GT 𝑴𝒃(𝑇) H B J ADE FC I OM N K L GT 3 / 9
  • 11. Markov Blankets Enable Locality in Causal Structure Recovery H B J ADE FC I OM N K L GT 𝑴𝒃(𝑇) H B J ADE FC I OM N K L GT Markov blankets can be used as a powerful tool in: classification, local causal discovery 3 / 9
  • 12. Markov Blankets: a Missing Concept in the Context of Chain Graph Models H B J ADE FC I OM N K L GT H B J ADE FC I OM N K L GT𝑴𝒃(𝑇) 4 / 9
  • 13. Markov Blankets: a Missing Concept in the Context of Chain Graph Models H B J ADE FC I OM N K L GT H B J ADE FC I OM N K L GT𝑴𝒃(𝑇) 𝒑𝒂(𝑇) 4 / 9
  • 14. Markov Blankets: a Missing Concept in the Context of Chain Graph Models H B J ADE FC I OM N K L GT H B J ADE FC I OM N K L GT𝑴𝒃(𝑇) 𝒄𝒉(𝑇) 4 / 9
  • 15. Markov Blankets: a Missing Concept in the Context of Chain Graph Models H B J ADE FC I OM N K L GT H B J ADE FC I OM N K L GT𝑴𝒃(𝑇) 𝒏𝒆(𝑇) 4 / 9
  • 16. Markov Blankets: a Missing Concept in the Context of Chain Graph Models H B J ADE FC I OM N K L GT H B J ADE FC I OM N K L GT𝑴𝒃(𝑇) 𝒄𝒔𝒑(𝑇) 4 / 9
  • 17. Markov Blankets: a Missing Concept in the Context of Chain Graph Models H B J ADE FC I OM N K L GT H B J ADE FC I OM N K L GT𝑴𝒃(𝑇) 4 / 9
  • 18. Markov Blankets in LWF Chain Graphs: Main Results Theorem (Characterization of Markov Blankets in LWF CGs) Let G = (V, E, P) be an LWF chain graph model. Then, G entails conditional independency T⊥⊥pV {T, Mb(T)}|Mb(T), i.e., the Markov blanket of the target variable T in an LWF CG probabilistically shields T from the rest of the variables. 5 / 9
  • 19. Markov Blankets in LWF Chain Graphs: Main Results Theorem (Characterization of Markov Blankets in LWF CGs) Let G = (V, E, P) be an LWF chain graph model. Then, G entails conditional independency T⊥⊥pV {T, Mb(T)}|Mb(T), i.e., the Markov blanket of the target variable T in an LWF CG probabilistically shields T from the rest of the variables. Theorem (Standard algorithms for Markov blanket recovery in LWF CGs) Given the Markov assumption, the faithfulness assumption, a graphical model represented by an LWF CG, and i.i.d. sampling, in the large sample limit, the Markov blanket recovery algorithms Grow-Shrink, Incremental Association Markov blanket recovery, and its variants identify all Markov blankets for each variable. 5 / 9
  • 20. Markov Blankets in LWF Chain Graphs: Main Results Theorem (Characterization of Markov Blankets in LWF CGs) Let G = (V, E, P) be an LWF chain graph model. Then, G entails conditional independency T⊥⊥pV {T, Mb(T)}|Mb(T), i.e., the Markov blanket of the target variable T in an LWF CG probabilistically shields T from the rest of the variables. Theorem (Standard algorithms for Markov blanket recovery in LWF CGs) Given the Markov assumption, the faithfulness assumption, a graphical model represented by an LWF CG, and i.i.d. sampling, in the large sample limit, the Markov blanket recovery algorithms Grow-Shrink, Incremental Association Markov blanket recovery, and its variants identify all Markov blankets for each variable. The characterization of Markov blankets in chain graphs enables us to develop new algorithms that are specifically designed for learning Markov blankets in chain graphs. 5 / 9
  • 21. MBC-CSP Algorithm: Markov Blanket Discovery in LWF CGs H B J A E FC I M K L GT D Grow Phase:  Step 1-1  𝒂𝒅𝒋(𝑇) H B J A E FC I M K L GT D Grow Phase:  Step 1-2  𝒂𝒅𝒋(𝑇) H B J A E FC I M K L GT D Grow Phase:  Step 2  𝐜𝐬𝐩(𝑇) 6 / 9
  • 22. MBC-CSP Algorithm: Markov Blanket Discovery in LWF CGs H B J A E FC I M K L GT D Grow Phase:  Step 1-1  𝒂𝒅𝒋(𝑇) H B J A E FC I M K L GT D Grow Phase:  Step 1-2  𝒂𝒅𝒋(𝑇) H B J A E FC I M K L GT D Grow Phase:  Step 2  𝐜𝐬𝐩(𝑇) H B J A E FC I M K L GT D Grow Phase:  Step 1-1  𝒂𝒅𝒋(𝑇) H B J A E FC I M K L GT D Grow Phase:  Step 1-2  𝒂𝒅𝒋(𝑇) H B J A E FC I M K L GT D Grow Phase:  Step 2  𝐜𝐬𝐩(𝑇) 6 / 9
  • 23. MBC-CSP Algorithm: Markov Blanket Discovery in LWF CGs H B J A E FC I M K L GT D Grow Phase:  Step 1-1  𝒂𝒅𝒋(𝑇) H B J A E FC I M K L GT D Grow Phase:  Step 1-2  𝒂𝒅𝒋(𝑇) H B J A E FC I M K L GT D Grow Phase:  Step 2  𝐜𝐬𝐩(𝑇) H B J A E FC I M K L GT D Grow Phase:  Step 1-1  𝒂𝒅𝒋(𝑇) H B J A E FC I M K L GT D Grow Phase:  Step 1-2  𝒂𝒅𝒋(𝑇) H B J A E FC I M K L GT D Grow Phase:  Step 2  𝐜𝐬𝐩(𝑇) Shrink Phase: 𝐌𝐛(𝑇) H B J A E FC I M K L GT D H B J A E FC I M K L GT D 6 / 9
  • 24. MBC-CSP Algorithm: Markov Blanket Discovery in LWF CGs H B J A E FC I M K L GT D Grow Phase:  Step 1-1  𝒂𝒅𝒋(𝑇) H B J A E FC I M K L GT D Grow Phase:  Step 1-2  𝒂𝒅𝒋(𝑇) H B J A E FC I M K L GT D Grow Phase:  Step 2  𝐜𝐬𝐩(𝑇) H B J A E FC I M K L GT D Grow Phase:  Step 1-1  𝒂𝒅𝒋(𝑇) H B J A E FC I M K L GT D Grow Phase:  Step 1-2  𝒂𝒅𝒋(𝑇) H B J A E FC I M K L GT D Grow Phase:  Step 2  𝐜𝐬𝐩(𝑇) Shrink Phase: 𝐌𝐛(𝑇) H B J A E FC I M K L GT D H B J A E FC I M K L GT D Such Markov blanket discovery algorithms help us to design new scalable algorithms for learning chain graphs based on local structure discovery. 6 / 9
  • 25. MbLWF Algorithm: Learning LWF CGs via Markov Blanket Discovery Observational Data Markov Blanket Discovery Algorithm A DC E B 𝐌𝐛 𝐴 = 𝐶, 𝐵 , 𝐒𝐞𝐩𝐬𝐞𝐭 𝐴, 𝐷 = 𝐶, 𝐵 , 𝐒𝐞𝐩𝐬𝐞𝐭 𝐴, 𝐸 = {𝐶, 𝐵} 𝐌𝐛 𝐵 = 𝐸, 𝐴 , 𝐒𝐞𝐩𝐬𝐞𝐭 𝐵, 𝐶 = 𝐸, 𝐴 , 𝐒𝐞𝐩𝐬𝐞𝐭 𝐵, 𝐷 = {𝐸, 𝐴} 𝐌𝐛 𝐶 = 𝐴, 𝐷 , 𝐒𝐞𝐩𝐬𝐞𝐭 𝐶, 𝐸 = 𝐴, 𝐷 , 𝐒𝐞𝐩𝐬𝐞𝐭 𝐶, 𝐵 = {𝐴, 𝐷} 𝐌𝐛 𝐷 = 𝐶, 𝐸 , 𝐒𝐞𝐩𝐬𝐞𝐭 𝐷, 𝐴 = 𝐶, 𝐵 , 𝐒𝐞𝐩𝐬𝐞𝐭 𝐷, 𝐵 = {𝐸, 𝐴} 𝐌𝐛 𝐸 = 𝐷, 𝐵 , 𝐒𝐞𝐩𝐬𝐞𝐭 𝐸, 𝐶 = 𝐴, 𝐷 , 𝐒𝐞𝐩𝐬𝐞𝐭 𝐸, 𝐴 = {𝐶, 𝐵} 7 / 9
  • 26. MbLWF Algorithm: Learning LWF CGs via Markov Blanket Discovery Observational Data Markov Blanket Discovery Algorithm A DC E B 𝐌𝐛 𝐴 = 𝐶, 𝐵 , 𝐒𝐞𝐩𝐬𝐞𝐭 𝐴, 𝐷 = 𝐶, 𝐵 , 𝐒𝐞𝐩𝐬𝐞𝐭 𝐴, 𝐸 = {𝐶, 𝐵} 𝐌𝐛 𝐵 = 𝐸, 𝐴 , 𝐒𝐞𝐩𝐬𝐞𝐭 𝐵, 𝐶 = 𝐸, 𝐴 , 𝐒𝐞𝐩𝐬𝐞𝐭 𝐵, 𝐷 = {𝐸, 𝐴} 𝐌𝐛 𝐶 = 𝐴, 𝐷 , 𝐒𝐞𝐩𝐬𝐞𝐭 𝐶, 𝐸 = 𝐴, 𝐷 , 𝐒𝐞𝐩𝐬𝐞𝐭 𝐶, 𝐵 = {𝐴, 𝐷} 𝐌𝐛 𝐷 = 𝐶, 𝐸 , 𝐒𝐞𝐩𝐬𝐞𝐭 𝐷, 𝐴 = 𝐶, 𝐵 , 𝐒𝐞𝐩𝐬𝐞𝐭 𝐷, 𝐵 = {𝐸, 𝐴} 𝐌𝐛 𝐸 = 𝐷, 𝐵 , 𝐒𝐞𝐩𝐬𝐞𝐭 𝐸, 𝐶 = 𝐴, 𝐷 , 𝐒𝐞𝐩𝐬𝐞𝐭 𝐸, 𝐴 = {𝐶, 𝐵} A DC E B Super Skeleton Recovery 7 / 9
  • 27. MbLWF Algorithm: Learning LWF CGs via Markov Blanket Discovery Observational Data Markov Blanket Discovery Algorithm A DC E B 𝐌𝐛 𝐴 = 𝐶, 𝐵 , 𝐒𝐞𝐩𝐬𝐞𝐭 𝐴, 𝐷 = 𝐶, 𝐵 , 𝐒𝐞𝐩𝐬𝐞𝐭 𝐴, 𝐸 = {𝐶, 𝐵} 𝐌𝐛 𝐵 = 𝐸, 𝐴 , 𝐒𝐞𝐩𝐬𝐞𝐭 𝐵, 𝐶 = 𝐸, 𝐴 , 𝐒𝐞𝐩𝐬𝐞𝐭 𝐵, 𝐷 = {𝐸, 𝐴} 𝐌𝐛 𝐶 = 𝐴, 𝐷 , 𝐒𝐞𝐩𝐬𝐞𝐭 𝐶, 𝐸 = 𝐴, 𝐷 , 𝐒𝐞𝐩𝐬𝐞𝐭 𝐶, 𝐵 = {𝐴, 𝐷} 𝐌𝐛 𝐷 = 𝐶, 𝐸 , 𝐒𝐞𝐩𝐬𝐞𝐭 𝐷, 𝐴 = 𝐶, 𝐵 , 𝐒𝐞𝐩𝐬𝐞𝐭 𝐷, 𝐵 = {𝐸, 𝐴} 𝐌𝐛 𝐸 = 𝐷, 𝐵 , 𝐒𝐞𝐩𝐬𝐞𝐭 𝐸, 𝐶 = 𝐴, 𝐷 , 𝐒𝐞𝐩𝐬𝐞𝐭 𝐸, 𝐴 = {𝐶, 𝐵} A DC E B Super Skeleton Recovery 𝐒𝐞𝐩𝐬𝐞𝐭 𝐴, 𝐷 = 𝐶, 𝐵 Skeleton A DC E B 7 / 9
  • 28. MbLWF Algorithm: Learning LWF CGs via Markov Blanket Discovery Observational Data Markov Blanket Discovery Algorithm A DC E B 𝐌𝐛 𝐴 = 𝐶, 𝐵 , 𝐒𝐞𝐩𝐬𝐞𝐭 𝐴, 𝐷 = 𝐶, 𝐵 , 𝐒𝐞𝐩𝐬𝐞𝐭 𝐴, 𝐸 = {𝐶, 𝐵} 𝐌𝐛 𝐵 = 𝐸, 𝐴 , 𝐒𝐞𝐩𝐬𝐞𝐭 𝐵, 𝐶 = 𝐸, 𝐴 , 𝐒𝐞𝐩𝐬𝐞𝐭 𝐵, 𝐷 = {𝐸, 𝐴} 𝐌𝐛 𝐶 = 𝐴, 𝐷 , 𝐒𝐞𝐩𝐬𝐞𝐭 𝐶, 𝐸 = 𝐴, 𝐷 , 𝐒𝐞𝐩𝐬𝐞𝐭 𝐶, 𝐵 = {𝐴, 𝐷} 𝐌𝐛 𝐷 = 𝐶, 𝐸 , 𝐒𝐞𝐩𝐬𝐞𝐭 𝐷, 𝐴 = 𝐶, 𝐵 , 𝐒𝐞𝐩𝐬𝐞𝐭 𝐷, 𝐵 = {𝐸, 𝐴} 𝐌𝐛 𝐸 = 𝐷, 𝐵 , 𝐒𝐞𝐩𝐬𝐞𝐭 𝐸, 𝐶 = 𝐴, 𝐷 , 𝐒𝐞𝐩𝐬𝐞𝐭 𝐸, 𝐴 = {𝐶, 𝐵} A DC E B Super Skeleton Recovery 𝐒𝐞𝐩𝐬𝐞𝐭 𝐴, 𝐷 = 𝐶, 𝐵 Skeleton A DC E B A DC E BComplex Recovery 7 / 9
  • 29. Experimental Evaluation: Markov Blankets Make a Broad Range of Inference/Learning Problems Computationally Tractable and More Precise. ●● ●● ● ● ● ● ● ● ● ● ● ● ● ● ● ● ● ● ● ●0.90 0.95 1.00 GSLWF fastIAMBLWF fdrIAMBLWF interIAMBLWF IAMBLWF MBCCSPLWF LCD Precision sample size size = 200 size = 2000 alpha = 0.05 8 / 9
  • 30. Experimental Evaluation: Markov Blankets Make a Broad Range of Inference/Learning Problems Computationally Tractable and More Precise. ●● ●● ● ● ● ● ● ● ● ● ● ● ● ● ● ● ● ● ● ●0.90 0.95 1.00 GSLWF fastIAMBLWF fdrIAMBLWF interIAMBLWF IAMBLWF MBCCSPLWF LCD Precision sample size size = 200 size = 2000 alpha = 0.05 ● ● ● ● ● ● 0.5 0.6 0.7 0.8 0.9 1.0 GSLWF fastIAMBLWF fdrIAMBLWF interIAMBLWF IAMBLWF MBCCSPLWF LCD Recall alpha = 0.05 8 / 9
  • 31. Markov Blankets: a Missing Concept in the Context of Chain Graph Models H B J ADE FC I OM N K L GT H B J ADE FC I OM N K L GT𝑴𝒃(𝑇) 4 / 9
  • 32. Markov Blankets: a Missing Concept in the Context of Chain Graph Models H B J ADE FC I OM N K L GT H B J ADE FC I OM N K L GT𝑴𝒃(𝑇) 4 / 9 Markov Blankets in LWF Chain Graphs: Main Results Theorem Let G = (V, E, P) be an LWF chain graph model. Then, G entails conditional independency T⊥⊥pV {T, Mb(T)}|Mb(T). Theorem Given the Markov assumption, the faithfulness assumption, a graphical model represented by an LWF CG, and i.i.d. sampling, in the large sample limit, the Markov blanket recovery algorithms Grow-Shrink, Incremental Association Markov blanket recovery, and its variants identify all Markov blankets for each variable. 5 / 9
  • 33. Markov Blankets: a Missing Concept in the Context of Chain Graph Models H B J ADE FC I OM N K L GT H B J ADE FC I OM N K L GT𝑴𝒃(𝑇) 4 / 9 Markov Blankets in LWF Chain Graphs: Main Results Theorem Let G = (V, E, P) be an LWF chain graph model. Then, G entails conditional independency T⊥⊥pV {T, Mb(T)}|Mb(T). Theorem Given the Markov assumption, the faithfulness assumption, a graphical model represented by an LWF CG, and i.i.d. sampling, in the large sample limit, the Markov blanket recovery algorithms Grow-Shrink, Incremental Association Markov blanket recovery, and its variants identify all Markov blankets for each variable. 5 / 9 MBC-CSP Algorithm: Markov Blanket Discovery in LWF CGs H B J A E FC I M K L GT D Grow Phase:  Step 1-1  𝒂𝒅𝒋(𝑇) H B J A E FC I M K L GT D Grow Phase:  Step 1-2  𝒂𝒅𝒋(𝑇) H B J A E FC I M K L GT D Grow Phase:  Step 2  𝐜𝐬𝐩(𝑇) H B J A E FC I M K L GT D Grow Phase:  Step 1-1  𝒂𝒅𝒋(𝑇) H B J A E FC I M K L GT D Grow Phase:  Step 1-2  𝒂𝒅𝒋(𝑇) H B J A E FC I M K L GT D Grow Phase:  Step 2  𝐜𝐬𝐩(𝑇) Shrink Phase: 𝐌𝐛(𝑇) H B J A E FC I M K L GT D H B J A E FC I M K L GT D 6 / 9
  • 34. Markov Blankets: a Missing Concept in the Context of Chain Graph Models H B J ADE FC I OM N K L GT H B J ADE FC I OM N K L GT𝑴𝒃(𝑇) 4 / 9 Markov Blankets in LWF Chain Graphs: Main Results Theorem Let G = (V, E, P) be an LWF chain graph model. Then, G entails conditional independency T⊥⊥pV {T, Mb(T)}|Mb(T). Theorem Given the Markov assumption, the faithfulness assumption, a graphical model represented by an LWF CG, and i.i.d. sampling, in the large sample limit, the Markov blanket recovery algorithms Grow-Shrink, Incremental Association Markov blanket recovery, and its variants identify all Markov blankets for each variable. 5 / 9 MBC-CSP Algorithm: Markov Blanket Discovery in LWF CGs H B J A E FC I M K L GT D Grow Phase:  Step 1-1  𝒂𝒅𝒋(𝑇) H B J A E FC I M K L GT D Grow Phase:  Step 1-2  𝒂𝒅𝒋(𝑇) H B J A E FC I M K L GT D Grow Phase:  Step 2  𝐜𝐬𝐩(𝑇) H B J A E FC I M K L GT D Grow Phase:  Step 1-1  𝒂𝒅𝒋(𝑇) H B J A E FC I M K L GT D Grow Phase:  Step 1-2  𝒂𝒅𝒋(𝑇) H B J A E FC I M K L GT D Grow Phase:  Step 2  𝐜𝐬𝐩(𝑇) Shrink Phase: 𝐌𝐛(𝑇) H B J A E FC I M K L GT D H B J A E FC I M K L GT D 6 / 9 MbLWF Algorithm: Learning LWF CGs via Markov Blanket Discovery Observational Data Markov Blanket Discovery Algorithm A DC E B 𝐌𝐛 𝐴 = 𝐶, 𝐵 , 𝐒𝐞𝐩𝐬𝐞𝐭 𝐴, 𝐷 = 𝐶, 𝐵 , 𝐒𝐞𝐩𝐬𝐞𝐭 𝐴, 𝐸 = {𝐶, 𝐵} 𝐌𝐛 𝐵 = 𝐸, 𝐴 , 𝐒𝐞𝐩𝐬𝐞𝐭 𝐵, 𝐶 = 𝐸, 𝐴 , 𝐒𝐞𝐩𝐬𝐞𝐭 𝐵, 𝐷 = {𝐸, 𝐴} 𝐌𝐛 𝐶 = 𝐴, 𝐷 , 𝐒𝐞𝐩𝐬𝐞𝐭 𝐶, 𝐸 = 𝐴, 𝐷 , 𝐒𝐞𝐩𝐬𝐞𝐭 𝐶, 𝐵 = {𝐴, 𝐷} 𝐌𝐛 𝐷 = 𝐶, 𝐸 , 𝐒𝐞𝐩𝐬𝐞𝐭 𝐷, 𝐴 = 𝐶, 𝐵 , 𝐒𝐞𝐩𝐬𝐞𝐭 𝐷, 𝐵 = {𝐸, 𝐴} 𝐌𝐛 𝐸 = 𝐷, 𝐵 , 𝐒𝐞𝐩𝐬𝐞𝐭 𝐸, 𝐶 = 𝐴, 𝐷 , 𝐒𝐞𝐩𝐬𝐞𝐭 𝐸, 𝐴 = {𝐶, 𝐵} A DC E B Super Skeleton Recovery 𝐒𝐞𝐩𝐬𝐞𝐭 𝐴, 𝐷 = 𝐶, 𝐵 Skeleton A DC E B A DC E BComplex Recovery 7 / 9
  • 35. Markov Blankets: a Missing Concept in the Context of Chain Graph Models H B J ADE FC I OM N K L GT H B J ADE FC I OM N K L GT𝑴𝒃(𝑇) 4 / 9 Markov Blankets in LWF Chain Graphs: Main Results Theorem Let G = (V, E, P) be an LWF chain graph model. Then, G entails conditional independency T⊥⊥pV {T, Mb(T)}|Mb(T). Theorem Given the Markov assumption, the faithfulness assumption, a graphical model represented by an LWF CG, and i.i.d. sampling, in the large sample limit, the Markov blanket recovery algorithms Grow-Shrink, Incremental Association Markov blanket recovery, and its variants identify all Markov blankets for each variable. 5 / 9 MBC-CSP Algorithm: Markov Blanket Discovery in LWF CGs H B J A E FC I M K L GT D Grow Phase:  Step 1-1  𝒂𝒅𝒋(𝑇) H B J A E FC I M K L GT D Grow Phase:  Step 1-2  𝒂𝒅𝒋(𝑇) H B J A E FC I M K L GT D Grow Phase:  Step 2  𝐜𝐬𝐩(𝑇) H B J A E FC I M K L GT D Grow Phase:  Step 1-1  𝒂𝒅𝒋(𝑇) H B J A E FC I M K L GT D Grow Phase:  Step 1-2  𝒂𝒅𝒋(𝑇) H B J A E FC I M K L GT D Grow Phase:  Step 2  𝐜𝐬𝐩(𝑇) Shrink Phase: 𝐌𝐛(𝑇) H B J A E FC I M K L GT D H B J A E FC I M K L GT D 6 / 9 MbLWF Algorithm: Learning LWF CGs via Markov Blanket Discovery Observational Data Markov Blanket Discovery Algorithm A DC E B 𝐌𝐛 𝐴 = 𝐶, 𝐵 , 𝐒𝐞𝐩𝐬𝐞𝐭 𝐴, 𝐷 = 𝐶, 𝐵 , 𝐒𝐞𝐩𝐬𝐞𝐭 𝐴, 𝐸 = {𝐶, 𝐵} 𝐌𝐛 𝐵 = 𝐸, 𝐴 , 𝐒𝐞𝐩𝐬𝐞𝐭 𝐵, 𝐶 = 𝐸, 𝐴 , 𝐒𝐞𝐩𝐬𝐞𝐭 𝐵, 𝐷 = {𝐸, 𝐴} 𝐌𝐛 𝐶 = 𝐴, 𝐷 , 𝐒𝐞𝐩𝐬𝐞𝐭 𝐶, 𝐸 = 𝐴, 𝐷 , 𝐒𝐞𝐩𝐬𝐞𝐭 𝐶, 𝐵 = {𝐴, 𝐷} 𝐌𝐛 𝐷 = 𝐶, 𝐸 , 𝐒𝐞𝐩𝐬𝐞𝐭 𝐷, 𝐴 = 𝐶, 𝐵 , 𝐒𝐞𝐩𝐬𝐞𝐭 𝐷, 𝐵 = {𝐸, 𝐴} 𝐌𝐛 𝐸 = 𝐷, 𝐵 , 𝐒𝐞𝐩𝐬𝐞𝐭 𝐸, 𝐶 = 𝐴, 𝐷 , 𝐒𝐞𝐩𝐬𝐞𝐭 𝐸, 𝐴 = {𝐶, 𝐵} A DC E B Super Skeleton Recovery 𝐒𝐞𝐩𝐬𝐞𝐭 𝐴, 𝐷 = 𝐶, 𝐵 Skeleton A DC E B A DC E BComplex Recovery 7 / 9 All code, data, and supplementary materials are available at: https://majavid.github.io/structurelearning/blog/2020/uai/