prob earthquake: 1/(365 + 0.25 + 0.01 + 0.0025)
prob burglar:  1/(365 + 0.25 + 0.01 + 0.0025)

% P(i1 | earthquake):
i1 <- earthquake & i1WhenEartquake
prob i1WhenEartquake: 0.2
i1 <- ~earthquake & i1WhenNotEartquake
prob i1WhenNotEartquake: 0

% P(i2 | burglar):
i2 <- burglar & i2WhenBurglar 
prob i2WhenBurglar: 0.95
i2 <- ~burglar & i2WhenNotBurglar
prob i2WhenNotBurglar: 0

% P(Alarm | i1, i2):
Alarm <- i1 & i2 & i1i2both
prob i1i2both: 1
Alarm <- ~i1 & i2 & i1i2i2
prob i1i2i2: 1
Alarm <- i1 & ~i2 & i1i2i1
prob i1i2i1: 1
Alarm <- i1 & ~i2 & i1i2none
prob i1i2none: 0

% P(Gibbons calls | Alarm)
Gibbons <- Alarm & GibbonsWhenAlarm
prob GibbonsWhenAlarm: 0.99
Gibbons <- ~Alarm & GibbonsWhenNoAlarm
prob GibbonsWhenNoAlarm: 0.04

% P(Watson calls | Alarm)
Watson <- Alarm & WatsonWhenAlarm
prob WatsonWhenAlarm: 0.8
Watson <- ~Alarm & WatsonWhenNoAlarm
prob WatsonWhenNoAlarm: 0.4

% P(Radio | earthquake)
Radio <- earthquake & RadioWhenEarthquake
prob RadioWhenEarthquake: 0.9998
Radio <- ~earthquake & RadioWhenNoEarthquake
prob RadioWhenNoEarthquake: 0.0002