Simulation data and code

- PF_simulation_data.zipcontains Simulation data to create figure 2 of de Vries, Erten and Kokko- Code_PF.zip contains C++ code to create the data used to create figure 2 (see PF_simulation_data.zip for the datafiles produced), and it also contains the R script to create figure 2 from the data (Figure2_cloud_25.R). All code files were created by Pen, I., & Flatt, T. (2021). Asymmetry, division of labour and the evolution of ageing in multicellular organisms. Philosophical Transactions of the Royal Society B, 376(1823), 20190729. C++ code is slightly adjusted to change output. Note that the R script takes a long time to run (multiple days on our laptops), and uses a lot of swap memory, we advice running it on a server. Alternatively, you can edit the code to use less than the last 25 days bychanging this line: ddead% filter(t>4975)to for example ddead% filter(t>4998) to use the last 2 time steps only. However, note that there will be insufficient data at high ages to estimate mortality rates.

PF_simulation_data.zip 包含用于复现de Vries、Erten与Kokko论文中图2的仿真数据集；Code_PF.zip 内含生成该图2所用数据的C++代码（生成的数据文件可参见PF_simulation_data.zip），同时附带了用于从该数据绘制图2的R脚本Figure2_cloud_25.R。所有代码由Pen, I. 与 Flatt, T. (2021) 完成，其对应研究论文为《Asymmetry, division of labour and the evolution of ageing in multicellular organisms》，发表于《Philosophical Transactions of the Royal Society B》第376卷第1823期，文章编号20190729。本次提供的C++代码已做小幅调整以修改输出格式。
需注意，该R脚本运行耗时极长（在我们的笔记本电脑上需耗时数日），且会占用大量交换内存，建议在服务器上运行该脚本。此外，也可通过修改代码调整运行规模：将代码行`ddead %>% filter(t>4975)`修改为例如`ddead %>% filter(t>4998)`，仅使用最后2个时间步进行计算。但需注意，此时高龄阶段的样本数据将不足以准确估算死亡率。