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.