3.14 Set up m-files to produce the geometric model likelihoods as explained in Section 2.2 for the two sets of fecundability data. Obtain the maximum likelihood estimates of p numerically, using the Matlab routine fminbnd as follows. Suppose one of the negative log-likelihoods is specified in the m-file: geo.m, then, in Matlab, use the command:

p = fminbnd (‘geo

, 0, 1)

Note the format for the m-file:

function y = geo(p)

data = ...

.

.

.

y = ...