These files accompany the paper "Improving estimations of life history parameters of small animals in mesocosm experiments: A case study on mosquitoes".
The code provided enables replication of our results.
If you wish to apply this methodology to your own research project then this code should provide a good starting point.

To run our model and the optimisation, use one of the following:
- normal_model_plus_optimisation.py
- mixed_model_plus_optimisation.py
The normal version should run on a regular laptop. The mixed version will require a supercomputer or computing cluster.
The mixed version requires Python3.7 or later.

The mixed version uses initial conditions defined in a separate text tile (initial_contidions_for_mixed.txt). 
We originally used very basic initial conditions(x0 = [0.5 for i in range(0,96)]+[1,1,1,2]), but found that this gave a non-optimal result for treatment 1.
To fix this we ran the optimisation again, with initial conditions for treatments 2-12 as the optimised values we had just calculated, and the initial 
conditions for treatment 1 as the optimised values from the all-normal approach.

The optimised probabilities for each approach are stored in:
- normal_probabilities.txt
- mixed_probabilities.txt

To run the validation use:
- validation.py

To make graphs of adult emergence and survival and development probabilities use:
- make_plots.py

Data for the main and validation experiments is stored in:
- main_data.csv
- validation_data.csv

To save the predicted populations and optimised probabilities to a csv file use:
- write_outputs_normal.py
- write_outputs_mixed.py