source("http://192.38.117.59/~linearpredictors/datafiles/readFever.R")

## Table of fever against fetal deaths.
(t <- table(fever$episodes,dnn = c("Number of fever episodes")))

## total no of episodes is the sum of "cases * no of episodes for that case".
no.episodes <-  sum(t*as.numeric(dimnames(t)[[1]]))
no.women <- nrow(fever[!is.na(fever$episodes),])

## estimated probability 
p.hat <- no.episodes/(no.women*14)

## expected numbers of episodes
round(dbinom(0:10, 14, 0.0143)*no.women, digits=1)