Matlab codes for running the maximum-likelihood estimation, data files for San Francisco, and simulated sighting records for the simulation exercise.
Ecological Archives E085-049-S1.
Authors
File list (downloadable)
Description
Andrew R. Solow
Marine Policy Center
Woods Hole Oceanographic Institution
Woods Hole, MA 02543
e-mail: asolow@whoi.eduChristopher J. Costello
Donald Bren School of Environmental Science and Management
University of California, Santa Barbara
Santa Barbara, CA 93106
e-mail: costello@bren.ucsb.edu
The main Matlab code (Count_Ecology_main.m) calls the remaining Matlab function files to run the maximum likelihood estimation of the parameters of the introduction and discovery process. The text file NumDis.txt contains the data to run the original analysis. The text files P1data.txt...P4data.txt contain the simulated sighting records over the 145 year sighting history (100 simulated sighting records for each of 4 parameter vectors, for a total of 400 independent 145-year sighting records).Count_Ecology_main.m - Main Matlab code for running the MLE.
count_lambda.m - Matlab function file for calculating lambda.
count_log_like.m - Matlab function file for calculating the log-likelihood value.
count_m.m - Matlab function file for calculating m.
count_p.m - Matlab function file for calculating p.
count_pi.m - Matlab function file for calculating pi.NumDis.txt - Text file containing the number of discoveries by year in San Francisco from 1851–1995.
P1data.txt - 100 simulated sighting records for first set of parameters (rows are years, columns are sighting records).
P2data.txt - 100 simulated sighting records for second set of parameters (rows are years, columns are sighting records).
P3data.txt - 100 simulated sighting records for third set of parameters (rows are years, columns are sighting records).
P4data.txt - 100 simulated sighting records for fourth set of parameters (rows are years, columns are sighting records).