6.16 As an exercise in verification, simulate appropriate multinomial data using the simulation program, ecosim, below, and then obtain the maximum likelihood estimates of the parameters using the function ecolik.

(a)

% Program to simulate multinomial data corresponding to

% a simple bird survival model, for 4 years of recovery.

% There are three probabilites : the annual survival probabilities,

% ‘phi_1’ and ‘phi_2’, and the reporting probability, ‘lam’

% ‘n’ denotes the number of birds in the cohort.

% The ‘data’ vector is global, as it may be used by the function

% ECOLIK, which forms the corresponding negative log-likelihood,

% to be minimised by fminsearch, for example.

% Calls MULTSIM

%___________________________________________________________________