-
Notifications
You must be signed in to change notification settings - Fork 0
Expand file tree
/
Copy pathplot_results.R
More file actions
105 lines (83 loc) · 18.1 KB
/
Copy pathplot_results.R
File metadata and controls
105 lines (83 loc) · 18.1 KB
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
# plotting code for SLiM 6 performance results
#
# Benjamin C. Haller, 30 September 2026
# the data here were generated by `_performance_tests.txt` on a 2022 Mac Studio (Mac13,2; Apple M1 Ultra; 128 GB RAM).
# the SLiM 6 data were generated by SLiM 6.0b1 (but nothing changed for SLiM 6.0 final that would affect runtimes).
# note that the run lengths of the models are tailored to produce runtimes of ~30 s on that machine, so these plots use the same y axis.
# note also that timing results are stochastic and irreproducible, and will depend on the hardware and software environment.
plot_perf <- function(title, old, old_name, new, new_name, y_text=10.0)
{
df <- data.frame(old, new)
speedup <- mean(old) / mean(new)
roundedSpeedup <- round(speedup * 100) / 100
if (roundedSpeedup >= 1.0) {
speedupText <- paste(format(roundedSpeedup, nsmall=2), "× faster", sep="")
speedupColor <- "red"
} else {
speedupText <- paste(format(roundedSpeedup, nsmall=2), "× as fast", sep="")
speedupColor <- "royalblue"
}
range(old, new)
quartz(title, 3, 3)
par(mai=c(0.75, 0.75, 0.1, 0.1), mgp=c(2.5, 0.5, 0))
boxplot(df, ylim=c(0, 40), boxwex=0.5, pch=19, cex=0.5, cex.axis = 0.9, names=c(old_name, new_name), xlab="", ylab="runtime (seconds)")
text(1.5, y_text, speedupText, adj=c(0.5, 0.5), cex=1.5, col=speedupColor)
mtext("100 replicates", side=1, line=2.5, cex=0.75, font=3)
box()
}
# test 1A
title <- "1A: neutral WF"
d_slim52 <- c(29.7843,30.5407,29.071,30.2636,29.8749,31.7629,30.1756,29.8783,32.1261,30.5665,29.9706,30.3427,30.6645,30.5949,32.4576,34.3851,29.6036,30.9779,33.6778,30.8726,32.7185,29.9824,30.1753,30.4782,30.6812,30.3754,32.5635,29.5317,29.9737,31.8403,30.0222,30.6247,33.3275,29.8428,29.8174,30.9252,30.3561,28.9833,30.7283,29.4952,30.5553,35.1756,32.8974,30.9783,34.1368,30.1113,31.7746,32.1162,32.7996,30.0003,29.7677,29.495,30.7392,30.5468,30.3525,29.3986,32.9512,29.214,30.6868,32.54,30.4151,32.6366,32.5511,30.363,30.5928,28.7838,29.3711,31.258,29.9271,34.865,30.082,33.2546,33.5238,32.28,30.1285,28.7689,30.6782,29.8643,31.3079,32.4442,32.3466,30.1384,30.2768,30.2385,30.1351,30.964,30.5698,30.6748,31.7497,32.7833,30.518,32.046,32.4045,30.6283,31.2541,30.0068,30.5946,29.9239,29.7393,32.8781)
d_slim6 <- c(32.1493,33.1116,29.3838,28.4715,31.2654,29.2563,32.3199,31.373,31.3443,32.5122,32.7147,30.7866,32.7065,31.4148,29.1986,30.5709,29.2067,31.9844,31.1155,32.8675,31.1836,33.5308,28.6994,30.77,30.6971,29.5578,28.6731,30.1578,31.8769,29.5834,32.3055,31.9744,30.7366,31.5699,29.6743,31.4276,31.4321,33.2664,28.7529,30.9312,30.8503,29.3458,29.7034,31.0481,29.4589,29.2696,31.2016,31.0034,29.492,29.4821,33.4071,30.563,29.9611,32.5434,31.5056,31.6792,31.1865,29.4767,31.6671,29.2879,29.2137,31.7471,34.146,29.6614,31.244,33.2967,30.4446,29.5363,32.7121,29.0899,31.8523,31.981,29.1723,31.7639,31.8799,32.6311,28.3142,30.5419,32.3295,32.0356,31.9488,29.5575,29.395,28.8985,29.5596,28.2353,28.9844,29.0989,28.6256,31.6643,28.9833,30.7073,28.7438,28.9372,30.549,28.8288,28.4762,28.4828,28.8197,28.4518)
plot_perf(title, d_slim52, "SLiM 5.2", d_slim6, "SLiM 6")
# test 1B
title <- "1B: neutral nonWF"
d_slim52 <- c(29.5579,31.0637,30.4872,31.6075,31.8403,30.3561,30.4607,30.176,32.0803,28.367,30.7172,32.2587,32.3346,32.146,29.8827,31.3887,33.8508,31.3789,29.8044,31.1312,31.3504,30.3924,30.9655,31.9361,30.1852,30.3645,28.5644,30.8167,28.9153,31.2513,30.9313,30.8019,29.2257,30.4162,30.8729,31.2537,31.8647,29.4198,29.7183,30.3753,30.1114,29.3081,33.0164,29.6163,29.7477,28.9625,31.6622,29.6478,30.8001,32.515,28.2373,31.9463,30.0534,30.3382,29.1569,30.6976,31.037,29.4217,29.5382,30.9488,29.9687,29.545,28.5061,30.2224,29.4504,30.752,31.2475,29.5589,28.7394,31.3424,30.8563,30.7708,30.3058,32.1291,29.5597,31.6384,29.6289,28.5404,29.9398,29.8032,30.6761,30.4217,29.6483,29.1975,32.3885,31.7487,30.5942,31.2687,31.6461,30.4858,32.616,31.9198,29.0379,31.2086,30.1019,30.2941,29.1182,29.5563,31.3382,32.4427)
d_slim6 <- c(31.5573,29.4562,30.1721,29.4684,28.33,28.6412,30.8707,29.9008,29.29,30.3619,28.8323,28.7624,30.189,30.1096,28.7145,30.8238,29.9815,30.2642,30.3498,29.4181,30.2358,29.6243,28.3685,28.3482,29.012,30.1452,28.4033,28.1137,31.3756,30.4138,30.7218,30.6429,28.0598,31.3003,32.2767,31.0938,28.9476,31.1109,29.9246,30.9746,28.393,30.7704,28.5797,30.629,29.4444,29.8029,29.9915,28.7616,31.6417,30.5602,31.5102,28.5582,31.1012,30.0358,30.8657,28.79,29.8893,29.2253,27.7042,30.6763,27.411,28.9095,30.0091,30.3973,29.3858,29.018,27.7518,29.4229,30.2313,31.235,31.5852,29.6337,32.2927,30.2011,30.9678,29.0955,30.344,29.3424,28.6644,29.7338,29.5482,30.9758,30.8228,28.7872,30.8049,28.5608,30.2374,28.2357,29.43,31.7734,30.0926,27.704,32.5994,31.6661,28.5987,30.6921,30.2337,27.7775,31.1764,28.8927)
plot_perf(title, d_slim52, "SLiM 5.2", d_slim6, "SLiM 6")
# test 1C
title <- "1C: non-neutral WF"
d_slim52 <- c(29.5422,33.5925,31.5651,33.5189,33.0997,34.1133,33.9837,32.5114,34.7346,37.7572,31.803,30.5734,33.6575,31.2878,33.8586,33.9287,31.2307,29.8,34.0934,31.0648,32.242,34.6699,32.4524,33.0189,32.6346,34.815,28.9066,34.6633,33.5221,31.0624,33.385,33.4361,34.1557,31.1843,33.2506,34.4816,33.2535,33.6875,32.4132,32.364,33.0777,34.4972,31.1549,33.4763,32.915,33.9022,32.1707,34.5732,32.9301,30.6513,32.4231,32.3248,31.8556,33.6236,33.9785,30.934,32.803,33.3467,33.4587,33.178,32.7529,33.7967,34.5359,35.2095,33.3717,31.9859,34.5213,34.7866,34.4875,33.4126,31.2254,32.7533,33.9583,33.8044,32.2437,30.2963,33.5166,32.7695,34.7585,30.3444,34.5663,33.0596,33.9081,32.8251,29.6267,30.7109,33.4057,30.4835,30.8878,33.2564,35.1825,31.7843,32.5203,34.1076,31.4458,34.2135,32.4749,33.1985,31.8156,34.1232)
d_slim6 <- c(28.4484,29.7142,29.9173,29.5195,30.2613,28.8164,28.2483,30.2676,29.1324,28.7391,28.9862,28.5175,28.38,28.0904,29.2682,29.3491,27.7419,29.4322,30.0566,29.9683,28.9821,30.0025,29.3522,28.2637,28.4797,29.524,27.2163,28.3444,29.1627,29.3517,29.9917,28.9576,29.2234,28.8053,29.8522,29.3794,27.8617,29.6714,28.487,28.5644,28.4624,30.5796,29.6577,30.1287,29.3969,29.6541,29.4906,29.3534,29.5328,28.803,30.3172,29.1156,29.9272,30.2998,29.9529,29.418,32.2415,32.2055,28.2458,31.1973,29.3201,30.5889,29.7752,28.7191,29.3105,28.7268,30.1709,29.9291,29.0994,30.5288,29.0428,29.2047,30.4725,30.5556,29.0555,29.3469,29.8292,29.5503,30.1808,29.2664,29.5099,30.1305,29.7128,29.9259,29.3512,28.3299,31.9286,30.3482,29.8519,29.8647,29.6271,30.3277,30.5227,28.8286,29.3386,30.9476,28.0173,30.1051,29.1143,29.4745)
plot_perf(title, d_slim52, "SLiM 5.2", d_slim6, "SLiM 6")
# test 1D
title <- "1D: non-neutral nonWF"
d_slim52 <- c(30.8866,32.971,34.4197,30.2932,34.6693,31.5377,30.877,34.4187,32.376,36.1689,30.7312,31.1956,33.141,33.5939,34.4577,30.4911,33.9181,32.6292,31.6086,33.071,31.5393,32.2511,29.6564,30.0178,31.3387,30.8666,32.612,32.2594,30.6054,31.9316,34.796,33.6886,32.9838,30.5096,31.2621,32.746,31.9021,31.683,30.5157,31.7835,33.4214,29.6857,31.278,33.0108,33.9972,33.4525,35.6491,35.9101,33.5811,30.2924,30.425,31.8682,30.5731,31.4639,31.7611,29.9713,30.2481,31.8928,33.3341,30.9102,32.7929,31.8356,33.0815,33.4886,31.53,32.1285,31.8081,32.756,33.1937,32.0444,31.116,33.5255,30.9469,32.6189,32.254,30.8088,31.9657,31.6674,31.8534,34.7836,30.5327,34.1314,31.524,30.4131,34.586,31.5669,31.9076,32.3678,30.5654,31.021,32.5594,30.6565,34.5822,32.4742,34.0015,33.0956,32.7224,32.495,31.5813,31.922)
d_slim6 <- c(17.9087,19.4156,20.4257,18.4271,20.508,18.9306,18.768,20.3361,19.2114,21.3876,18.3037,18.3713,19.665,19.8389,20.4215,18.0736,19.8847,19.2484,18.6134,20.092,19.0189,19.0812,17.9445,17.9867,18.3812,18.1607,19.2516,19.0175,18.0431,18.4072,20.4779,19.8695,19.349,17.4485,18.3551,19.0777,18.8601,18.6254,18.1196,18.8651,20.2677,17.5886,18.2093,19.2084,20.1122,19.7779,19.6737,19.2896,19.6803,17.9101,17.9764,18.9658,18.2062,18.7629,19.0348,17.9582,18.1259,18.8013,19.7837,18.1855,19.1189,18.7888,19.2718,19.744,18.7238,18.937,18.8004,19.2812,19.4551,19.2073,18.7837,19.7486,17.9529,19.0684,18.8127,18.333,19.4322,19.0162,19.0052,20.6048,18.4289,20.2882,18.5862,18.4236,20.1539,18.5892,18.5589,19.3063,18.056,18.335,18.871,18.0852,20.1489,19.0465,19.7914,19.3964,19.338,19.3707,18.6652,18.8371)
plot_perf(title, d_slim52, "SLiM 5.2", d_slim6, "SLiM 6")
# test 2
title <- "2: one quantitative trait"
d_slim52 <- c(30.1059,30.619,32.3251,30.304,29.6394,32.0794,30.9515,30.5428,32.3475,30.511,33.0333,30.2843,32.6404,31.3811,31.9103,31.2664,31.9405,31.4758,29.4495,32.4737,30.5145,28.8982,29.91,32.5369,33.0834,31.1411,30.7486,29.7624,30.6599,30.844,29.7804,31.5459,30.4904,30.5153,30.3721,31.714,32.1251,30.7356,30.6408,28.8653,29.8476,31.4457,31.7779,30.8493,32.3238,30.6987,31.1592,32.3557,30.6644,31.6665,31.5751,30.445,31.1484,32.3584,32.2811,29.1574,30.8752,30.3381,30.2985,31.1848,30.2116,30.4814,32.6629,33.0036,31.4838,31.4824,29.9524,31.7541,30.6577,31.7034,32.0947,29.5622,31.7632,31.9633,30.5004,30.0097,29.372,29.3214,30.3327,30.475,31.5192,30.9829,33.7084,31.9331,32.5706,30.1375,31.5643,32.232,29.4734,32.3733,31.273,30.0791,33.2788,32.3032,33.2352,30.6521,31.291,31.2346,32.451,29.2628)
d_slim6 <- c(2.70822,2.71071,2.70565,2.71186,2.71269,2.70676,2.70792,2.7086,2.71501,2.7084,2.73572,2.70302,2.70599,2.71503,2.71248,2.71493,2.72117,2.72231,2.7168,2.7179,2.71706,2.71042,2.70342,2.71324,2.70689,2.72956,2.70145,2.71343,2.71694,2.70408,2.72377,2.70723,2.72569,2.71215,2.71037,2.70852,2.71752,2.70159,2.71743,2.72375,2.72051,2.72153,2.70839,2.72496,2.72934,2.72922,2.71755,2.71574,2.73031,2.79091,2.81083,2.74082,2.72198,2.72451,2.72819,2.78271,2.81873,2.73335,2.73087,2.72776,2.72313,2.72765,2.7298,2.73203,2.73383,2.7207,2.73685,2.73392,2.72291,2.73078,2.7314,2.73024,2.72906,2.73178,2.72533,2.72707,2.73479,2.73349,2.73168,2.7197,2.72803,2.72774,2.72514,2.72085,2.72782,2.7264,2.73242,2.73754,2.72304,2.73897,2.7548,2.74363,2.7282,2.72969,2.72783,2.72424,2.72924,2.73718,2.71857,2.73407)
plot_perf(title, d_slim52, "SLiM 5.2", d_slim6, "SLiM 6")
# test 3
title <- "3: two quantitative traits"
d_slim52 <- c(37.0258,36.6948,29.8178,33.6511,35.1247,29.3844,32.3656,33.3357,33.9673,32.9387,37.4043,37.7035,30.3256,35.562,35.544,31.3566,35.8093,37.2946,33.2609,35.7463,36.293,33.7692,33.6676,35.2264,32.2234,30.858,37.1562,33.6402,37.3571,34.8356,34.5319,41.6703,34.2033,30.1924,35.0184,35.9591,33.5006,36.3299,36.6132,36.0138,30.9334,36.8524,39.0747,31.1706,35.7253,33.8848,38.0398,34.422,33.4737,36.0479,35.251,35.3788,37.9555,34.1363,34.2876,33.4927,32.923,36.1787,34.8382,29.9904,35.0849,32.3007,33.7516,29.9042,41.603,35.6173,34.7265,32.6911,33.4913,39.7704,36.1424,33.1774,34.9319,35.4965,35.525,39.2577,37.0979,39.1197,35.8009,36.3435,37.3288,30.2871,35.7801,32.1636,33.1997,29.0881,33.4397,34.7982,33.7176,30.4469,37.1609,37.9122,36.0826,35.2765,34.8106,36.9119,32.7658,38.1318,32.5747,34.7493)
d_slim6 <- c(1.01834,1.01538,1.01989,1.02054,1.01136,1.01917,1.02426,1.00899,1.01748,1.0213,1.01782,1.01178,1.02026,1.01873,1.01226,1.01799,1.02004,1.01605,1.01377,1.02029,1.02017,1.01135,1.02031,1.01798,1.0106,1.02032,1.02428,1.00804,1.0176,1.02043,1.01486,1.01617,1.01149,1.01058,1.01707,1.01273,1.01154,1.01268,1.01234,1.01285,1.00846,1.01036,1.01138,1.01354,1.01064,1.00992,1.01805,1.00829,1.0176,1.01659,1.01919,1.02447,1.01734,1.01492,1.01984,1.02076,1.01223,1.01927,1.01886,1.01067,1.02176,1.02116,1.01692,1.0116,1.02091,1.01826,1.01093,1.0217,1.0192,1.01183,1.01631,1.01786,1.01696,1.01448,1.02214,1.0226,1.01272,1.02027,1.02095,1.01784,1.0183,1.0218,1.0213,1.0116,1.01915,1.01672,1.01114,1.02117,1.01871,1.01828,1.01119,1.01954,1.01993,1.01106,1.04597,1.04997,1.05416,1.05383,1.05013,1.03574)
plot_perf(title, d_slim52, "SLiM 5.2", d_slim6, "SLiM 6")
# test 4
title <- "4: per-mutation dominance"
d_slim52 <- c(32.2939,25.3197,22.8152,23.2937,20.8524,19.1245,21.8223,32.3391,26.8895,34.5446,40.9653,22.8506,23.0378,25.1731,30.2291,23.7791,27.832,21.9099,19.6343,27.1798,30.0543,23.2934,26.0966,21.1344,23.2445,27.7886,24.9254,20.213,30.6084,22.9272,21.3234,28.0688,26.589,23.1304,28.2689,24.9578,21.4875,27.1234,22.7099,22.7723,30.9642,22.1579,22.5505,18.35,23.7056,22.5925,20.7491,20.8485,19.8363,25.304,21.0302,25.2881,23.0973,23.1007,23.8579,26.225,22.5129,22.7898,26.3547,23.178,30.7722,24.2164,23.5798,22.9353,21.7585,20.4048,22.4678,21.4712,44.9814,21.2915,29.1027,25.9859,29.9075,29.5028,24.5162,21.2802,22.9435,21.8098,22.8677,28.1437,18.7497,23.7063,18.6813,23.0742,22.8742,19.35,24.8259,22.7449,22.0106,31.8595,23.9189,21.2064,23.6814,23.1276,22.4387,26.1209,20.2369,24.7857,27.1354,29.6132)
d_slim6 <- c(3.11036,2.85937,2.83693,3.02276,2.82598,2.69732,2.79441,2.92622,2.90608,2.81454,2.85593,2.83602,2.82043,2.90333,2.76678,2.78747,2.7256,2.78694,2.84923,2.81058,2.86647,2.91391,2.75641,2.92895,2.7199,3.41248,2.77714,2.92812,2.93305,2.72553,2.95386,2.90018,3.03192,2.80183,2.80676,3.0221,2.93146,2.86805,2.71795,2.73556,2.83317,3.21387,3.24548,2.77582,3.13537,2.89495,2.76256,2.88613,2.86882,2.96597,2.79331,2.75692,2.76729,2.92092,2.77773,2.82843,2.89782,2.68262,2.93727,3.08257,2.82628,2.81043,2.73012,3.02861,2.97824,2.7851,2.72073,3.04182,2.68107,3.11084,3.96743,2.79598,2.8866,2.87379,2.72261,2.90509,2.75521,2.74122,3.65407,2.84204,2.78031,2.72998,2.868,2.89179,2.81343,2.85041,2.8833,2.84177,2.86481,2.84261,2.74336,2.86955,2.80203,2.82296,2.74133,2.91681,2.84643,3.06657,2.87235,2.95706)
plot_perf(title, d_slim52, "SLiM 5.2", d_slim6, "SLiM 6")
# test 5
title <- "5: independent dominance"
d_intermediate <-
c(31.5227,31.2336,31.0492,31.1269,31.2052,31.2736,31.1945,31.558,31.4461,31.1124,31.4651,31.357,31.2479,31.2121,31.341,31.609,31.1276,30.8803,31.3749,31.1305,30.9711,31.0128,31.5328,31.9604,31.0795,31.1947,31.0326,31.234,30.8077,31.2784,31.1044,31.4614,31.3735,31.3608,31.1754,31.6137,31.4675,31.2404,31.1746,31.4045,31.5295,31.1071,31.0486,31.202,30.9929,31.4524,31.0288,31.3505,31.3496,30.93,31.4036,31.6074,31.2835,31.0053,30.6949,31.079,31.3337,31.4021,31.6197,31.3635,31.5599,31.4975,31.7946,31.0457,31.5998,31.5598,31.2413,31.2763,31.0164,31.3999,31.6109,31.3001,31.1475,31.0388,31.179,31.2922,31.8731,31.3501,31.4165,31.378,31.0824,31.1458,31.2849,31.4981,31.4329,31.146,31.5733,31.3354,31.5517,31.3762,31.5999,31.5116,31.2429,31.0908,31.2612,30.8931,31.4176,31.6591,31.2468,31.2586)
d_independent <-
c(13.7587,13.7046,13.6332,13.6633,13.7623,13.5134,13.7241,13.5658,13.793,13.9706,14.0229,13.4333,13.6032,13.7376,13.6661,13.9284,13.7072,13.8063,13.7692,13.5931,13.7737,13.6819,13.6776,13.715,13.6345,13.6645,13.4524,13.7673,13.7938,13.565,13.658,13.6755,13.5336,13.5418,13.5536,13.6407,13.647,13.7672,13.7686,13.6179,13.7885,13.8218,13.6229,13.8386,13.6142,13.5209,13.6323,13.589,13.6587,13.6802,13.6362,13.6685,13.7442,13.6536,13.6718,13.7329,13.9149,13.7142,13.6003,13.6414,13.5367,13.7213,13.6494,13.526,13.7669,13.6106,13.5314,13.6436,13.5571,13.716,13.5889,13.6133,13.7212,13.7213,13.8185,13.771,13.7886,13.8228,13.6179,13.606,13.7353,13.9159,13.7272,13.605,13.6513,13.6328,13.7814,13.5905,13.7877,13.7201,13.834,13.9071,13.5573,13.5123,13.832,13.683,13.6608,13.8829,13.6933,13.5559)
plot_perf(title, d_intermediate, "intermediate", d_independent, "independent", y_text=8.0)
# test 6
title <- "6: substitution accumulation"
d_without <-
c(32.8924,34.983,33.1125,35.0127,34.5132,35.0983,34.8623,35.0428,33.6825,34.0884,35.0261,32.1158,34.6895,32.707,34.4929,36.8461,33.5463,35.1163,35.4066,33.3596,34.5157,34.3507,34.2629,34.3959,33.1476,31.4371,33.0486,35.4053,32.7745,35.8137,35.5723,34.1355,34.6091,33.39,33.346,34.3644,35.0916,32.051,35.1855,34.4962,34.4089,33.8094,33.3005,33.4259,34.8716,36.0845,33.3175,32.8225,37.2751,35.0732,35.7545,35.4348,34.1948,34.4582,35.1564,33.8989,33.3253,34.1358,36.65,31.9519,33.8004,36.1909,33.8342,34.6404,36.2999,33.7832,37.2155,36.7351,34.9805,36.5729,36.0557,32.4001,33.9455,35.0137,36.73,32.262,33.0253,33.4449,32.9028,34.4026,34.2655,35.501,33.8612,36.5041,37.3979,36.9617,35.8771,35.3915,33.1406,35.2521,35.9076,32.6116,37.5301,32.7174,33.5999,36.78,32.8476,33.8845,36.1825,33.1793)
d_with <-
c(5.36533,5.36009,5.01326,5.45619,5.42684,5.23493,5.28235,5.2664,5.25902,5.14212,5.44959,5.22785,5.32007,5.47138,5.41213,5.4841,5.25674,5.40874,5.21404,5.3932,5.48653,5.37517,5.12869,5.24407,5.20933,5.19343,5.19176,5.44021,5.4285,5.19889,5.57975,5.24773,5.71959,5.36737,5.36754,5.26934,5.57734,5.10334,5.50042,5.46263,5.50098,5.63512,5.2435,5.25861,5.31473,5.38091,5.21502,5.21582,5.41185,5.43247,5.23701,5.27865,5.31382,5.22599,5.43975,5.28216,5.31727,5.18497,5.43872,5.56291,5.49747,5.472,5.61547,5.66761,5.54903,5.34811,5.58148,5.67878,5.3935,5.51593,5.29879,5.27304,5.4321,5.44394,5.56942,5.36311,5.55799,5.49202,5.3963,5.2432,5.08824,5.27418,5.26971,5.3765,5.41008,5.34223,5.23431,5.51839,5.1588,5.40417,5.43116,5.25808,5.43543,5.35008,5.21981,5.33175,5.36932,5.36105,5.30937,5.43429)
plot_perf(title, d_without, "without", d_with, "with", y_text=15.0)
# test 7
title <- "7: mutation effect accumulation"
d_without <-
c(31.6215,32.9004,33.0603,33.7637,32.7159,32.155,31.6858,33.7141,33.239,32.9221,32.7364,31.9916,32.7332,32.1732,31.9633,32.3895,32.994,33.1138,33.2161,32.8082,33.3671,33.1271,33.3538,33.5308,33.2711,32.0978,32.7468,32.9912,33.2723,33.5391,32.3651,34.0971,33.5884,33.8409,33.6548,33.7762,32.4828,33.4161,33.1612,34.3357,35.7392,35.0524,32.0486,33.0793,32.9051,33.0939,33.0531,32.8028,32.214,33.1656,33.0018,32.5518,33.2895,33.1297,32.5971,33.4813,33.8837,33.0554,32.6723,32.7431,33.4951,32.8753,31.9493,34.2303,33.8153,32.7046,32.7856,31.9744,33.8893,33.2069,32.0732,34.0653,33.6379,31.5844,32.6449,32.7103,33.2933,32.7741,34.5781,35.1723,33.7547,33.2643,32.2045,33.2073,32.3229,32.7945,34.2226,32.3511,32.4234,33.1157,33.1886,32.1748,34.0987,32.0875,33.0042,32.2953,33.4912,32.3546,33.6642,33.7959)
d_with <-
c(8.37121,8.3632,8.36356,8.37135,8.38949,8.68661,8.67912,8.75519,8.70175,8.77573,8.6815,8.69717,8.541,8.38061,8.46785,8.38448,8.38176,8.41062,8.42988,8.52306,8.39713,8.35994,8.38632,8.38437,8.36701,8.37411,8.38583,8.36787,8.36458,8.37527,8.39244,8.57228,8.41035,8.55396,8.36308,8.37644,8.37159,8.36625,8.36977,8.32101,8.32516,8.32468,8.33053,8.3132,8.33142,8.32286,8.31268,8.32308,8.32625,8.33378,8.33176,8.31513,8.32681,8.32989,8.3275,8.32068,8.33589,8.3584,8.37812,8.36501,8.62263,8.35952,8.36733,8.36801,8.3557,8.36613,8.37244,8.37854,8.6542,8.37125,8.3701,8.36029,8.3768,8.3636,8.36177,8.47192,8.45054,8.3576,8.3905,8.36239,8.35107,8.4665,8.36613,8.48425,8.36186,8.36249,8.3596,8.3812,8.3747,8.49057,8.37397,8.36016,8.35785,8.37063,8.37384,8.37058,8.37875,8.44408,8.37262,8.47262)
plot_perf(title, d_without, "without", d_with, "with", y_text=15.0)