define plan {
m = 0
for i=1 to 4 do {
for j=1 to ncircles[i] do {
exec sprintf("C[i,j]: circle thick 1.6 at Here \"\large\sf $%g\"",m+j)
if i>1 then {
if "$1"=="" then {
for k=1 to ncircles[i-1] do { arrow from C[i,j] to C[i-1,k] chop } }
}
if j!=ncircles[i] then { move right_ circlesep from C[i,j] }
}
m += ncircles[i]
case(i,
move to C[1,1]+(-circlesep/2,-rowsep),
move to C[2,1]+(0,-rowsep),
move to C[3,1]+(circlesep/2,-rowsep), )
}
}
gridsiz = 0.78/4
define grid { [ SW: 0,0
for i=1 to $1 do {
move to (0,(0.5-i)*gridsiz)
exec sprintf("fv = $%g",2+i)
for j=1 to $2 do { fl = 1-(int(fv / 2^(j-1)) % 2)
box wid gridsiz ht gridsiz outlined "gray" shaded rgbstring(fl,fl,fl) }
} ]}
define cxt {
m = 1
for i=2 to 4 do {
for j=1 to ncircles[i] do {
exec sprintf("narrows = $%g",m)
for k=1 to narrows do {
exec sprintf("cnum = $%g",m+k)
arrow from C[i,j] to C[i-1,cnum] chop
}
m += narrows+1
} } }
array(ncircles,3,4,4,3)
plan()
for j=1 to 3 do {
sprintf("$\hat{x}_%g$",j) at C[1,j].n above
sprintf("$x_%g$",j) at C[4,j].s+(0,-2bp__) below
}
"V" at 0.5 between C[2,1].w and C[1,1].w
"W$^2$" at 0.5 between C[3,1].w and C[2,1].w
"W$^1$" at 0.5 between C[3,1].w and C[4,1].w
GV: grid(3,4,0,15,6) with .w at C[2,4]+(circlerad*2,rowsep/2)
"$=\hbox{M}^{\hbox{V}}$" at GV.e ljust
G2: grid(4,4,13,0,0,13) with .w at (GV.w,C[3,4]+(circlerad*2,rowsep/2))
"$=\hbox{M}^{\hbox{V}^2}$" at G2.e ljust
G1: grid(4,3,1,5,1,1) with .w at (GV.w,C[4,3]+(circlerad*2,rowsep/2))
"$=\hbox{M}^{\hbox{V}^1}$" at G1.e ljust
move to (GV.x+3*circlesep,C[1,1].y)
plan(3,1,2,1,2,2,1,2,1,2,2,3,1,2)
for i=2 to 4 do {
rgbdraw(0.5,0.5,0.5,
cxt(2,1,3,0,0,2,1,3, 0,4,1,2,3,4,0,0, 0,4,1,2,3,4,0) )
cxt(0,1,1,1,1,0, 2,2,3,0,2,2,3,2,2,3, 0,0,3,1,3,4)
}
for j=1 to 3 do { sprintf("$x_%g$",j) at C[4,j].s+(0,-2bp__) below }
arrow dotted from C[1,2].n up circlesep/2; "$p(x_2)$" above
"$p(x_3|x_2)$" at (C[1,3]+(circlesep/2,0),last "") above
arrow dotted from C[1,3] to last "".s chop circlerad chop textoffset
"$p(x_1|x_2,x_3)$" at (C[1,1]+(-circlesep/2,0),last "") above
arrow dotted from C[1,1] to last "".s chop circlerad chop textoffset
