Supplement 1. The data and Mathematica codes used in the age-dependent survival analysis of Glanville fritillary butterflies.
File List
survival_modelI.nb – Mathematica code for model I of survival analysis survival_modelII.nb – Mathematica code for model II of survival analysis individual_analysis.nb – Mathematica code for the analysis of individual errors in death rates MCMCAlgorithms.m – Mathematica package needed to run the above programs MCMCAlgorithms.nb – Used to generate MCMCAlgorithms.m realdata.txt – Data y for each individual at each census covariates.txt – Individual covariates weather.txt – Temperature at each census jlast.txt – The last recapture census of each individual femalemobility.txt – Mobility of females measured by the number of moves between different cells in the closed cage readme.txt – Description of the data files all.zip – A zip of all files
Description
The files survival_modelI.nb and survival_modelII.nb are Mathematica programs that can be used to estimate the parameters of the main model, and the file individual_analysis.nb performs the regressions relating to individual error terms. All the data files are described in readme.txt.
Download all the files in the same directory. Change in the programs "workdirectory" to the current directory. The data files realdata.txt, covariates.txt, weather.txt, and jlast.txt are needed to run the program survival_xxx.nb. Change the random seed to obtain different starting values. The posterior distributions are saved in the results_xxx.txt files.
In the program individual_analysis.nb, the estimated distribution of individual errors is imported from the file individual_errors.txt. Set “infile” to femalemobility.txt if individual errors are regressed against the individual’s mobility. The posterior distribution is saved to the file "outfile".