{ "optimal_k": 6, "communities": { "0": { "community_id": 0, "taxa_count": 22, "timespan": [ 2014.8027397260273, 2022.1835616438357 ], "ols": { "method": "OLS", "status": "OK", "mu": 0.0014184720296959896, "d0": 0.009782068469976002, "t_ref": 2019.031256082042, "t_mrca": 2012.1350549126241, "se_mu": 8.448441371739223e-05, "se_d0": 0.0001822527019846743, "se_mrca": 0.4303659179698413, "ci_fieller": [ 2011.115764309522, 2012.9381130818283 ], "fieller_g": 0.015435653609799926, "ci_delta": [ 2011.23732733879, 2013.0327824864582 ], "ci_bootstrap": null, "ci_mrca": [ 2011.115764309522, 2012.9381130818283 ], "ci_method": "fieller", "r": 0.9663084095040785, "r2": 0.9337519422783018, "p_value": 2.9598545836506673e-13, "sigma2": 7.30753042375719e-07, "rmse": 0.0008150588614414851, "residuals": "[ 9.80286666e-04 9.80286666e-04 -9.00247392e-04 7.04357253e-04\n 1.47142715e-03 -8.82736155e-04 -1.11858639e-03 1.19669143e-04\n -3.72817457e-04 1.00154463e-03 2.82079228e-04 -7.36111915e-04\n -9.78636927e-04 -2.94955511e-04 -7.34074802e-05 -6.52146461e-04\n 2.60679232e-04 -4.19611463e-04 8.14479385e-04 -4.71357245e-04\n 1.31919854e-03 -1.03339349e-03]", "fitted": "[0.00996702 0.00996702 0.0124274 0.00965612 0.01231501 0.01302548\n 0.00378404 0.0045574 0.01425353 0.00560246 0.01007972 0.01414083\n 0.01018076 0.01006806 0.00499265 0.01134223 0.01009527 0.01077535\n 0.00838144 0.00560246 0.01375998 0.01023128]", "times": "[2019.16164384 2019.16164384 2020.89617486 2018.94246575 2020.81693989\n 2021.31780822 2014.80273973 2015.34794521 2022.18356164 2016.08469945\n 2019.24109589 2022.10410959 2019.31232877 2019.23287671 2015.65479452\n 2020.13114754 2019.25205479 2019.73150685 2018.04383562 2016.08469945\n 2021.83561644 2019.34794521]", "n": 22 }, "pgls": {}, "spline": { "method": "RESTRICTED_SPLINE", "t_mrca": 2012.5837824212535, "ci_mrca": [ 2012.5837824212535, 2012.5837824212535 ], "rate_ancestral": 0.001587314620409123, "ci_rate_ancestral": [ 0.001587314620409123, 0.001587314620409123 ], "rate_recent": 0.0012344117695646081, "ci_rate_recent": [ 0.0012344117695646081, 0.0012344117695646081 ], "rate_ratio": 0.7776730294630841, "beta_0": -3.194603662635549, "beta_1": 0.001587314620409123, "beta_2": -0.00018479896305730978, "beta": [ -3.194603662635549, 0.001587314620409123, -0.00018479896305730978 ], "ci_beta_2": [ -0.00018479896305730978, -0.00018479896305730978 ], "knots": [ 2014.8027397260273, 2019.2465753424658, 2021.7838356164382 ], "rss": 1.4049283818017545e-05, "aic": -307.8075949942697, "rss_linear": 1.4615060847514376e-05, "aic_linear": -308.93900988696504, "delta_aic": 1.131414892695318, "aic_reduction": -1.131414892695318, "f_stat": 0.7651467291630707, "p_f_test": 0.39264592492254147, "is_nonlinear_preferred": false, "r2": 0.9363165316220466, "rmse": 0.000799126904759118, "fitted": "[0.0101271 0.0101271 0.01238331 0.00982422 0.01228414 0.01290775\n 0.00352218 0.00438698 0.01397712 0.00554907 0.01023572 0.01387904\n 0.01033258 0.01022452 0.00487232 0.01141362 0.01025066 0.01089287\n 0.00853772 0.00554907 0.01354761 0.01038081]", "residuals": "[ 8.20210417e-04 8.20210417e-04 -8.56153842e-04 5.36257299e-04\n 1.50229490e-03 -7.65009496e-04 -8.56733469e-04 2.90082673e-04\n -9.64121680e-05 1.05493684e-03 1.26075495e-04 -4.74330592e-04\n -1.13045081e-03 -4.51409904e-04 4.69279059e-05 -7.23529988e-04\n 1.05286992e-04 -5.37122210e-04 6.58197790e-04 -4.17965039e-04\n 1.53156094e-03 -1.18292415e-03]", "n": 22 }, "calibrated_rate": 0.0014184720296959896, "calibrated_tmrca": 2012.1350549126241, "ci_mrca": [ 2011.115764309522, 2012.9381130818283 ], "r2": 0.9337519422783018, "pagel_lambda": null }, "1": { "community_id": 1, "taxa_count": 14, "timespan": [ 1987.0, 2018.2465753424658 ], "ols": { "method": "OLS", "status": "OK", "mu": 0.0010935447862130598, "d0": 0.016892183945858383, "t_ref": 1995.1757338551859, "t_mrca": 1979.7285531280427, "se_mu": 0.00011476133977873895, "se_d0": 0.0009519306525281072, "se_mrca": 1.8400317000484763, "ci_fieller": [ 1974.6709484807375, 1983.0818095658824 ], "fieller_g": 0.0522826893163999, "ci_delta": [ 1975.7194684529898, 1983.7376378030956 ], "ci_bootstrap": null, "ci_mrca": [ 1974.6709484807375, 1983.0818095658824 ], "ci_method": "fieller", "r": 0.9398231493088111, "r2": 0.8832675519767318, "p_value": 6.017210251174632e-07, "sigma2": 1.2686407541116233e-05, "rmse": 0.003297584510921753, "residuals": "[-1.10746778e-04 7.35836289e-03 -7.89738586e-04 6.01059159e-03\n -1.12205970e-03 -1.63279285e-03 -1.04912380e-03 2.74059798e-04\n 1.04596850e-03 -1.38164620e-03 -1.88855692e-03 4.07630658e-04\n -8.74645013e-05 -7.03448410e-03]", "fitted": "[0.01232583 0.02772834 0.01013874 0.02777627 0.01123229 0.01232583\n 0.01232583 0.01232583 0.01232583 0.0188871 0.00795165 0.01560647\n 0.01341938 0.04212118]", "times": "[1991. 2005.08493151 1989. 2005.12876712 1990.\n 1991. 1991. 1991. 1991. 1997.\n 1987. 1994. 1992. 2018.24657534]", "n": 14 }, "pgls": {}, "spline": { "method": "RESTRICTED_SPLINE", "t_mrca": 1985.2127712326528, "ci_mrca": [ 1985.2127712326528, 1985.2127712326528 ], "rate_ancestral": 0.0022410058067572937, "ci_rate_ancestral": [ 0.0022410058067572937, 0.0022410058067572937 ], "rate_recent": 0.0007858892535023718, "ci_rate_recent": [ 0.0007858892535023718, 0.0007858892535023718 ], "rate_ratio": 0.35068595142979275, "beta_0": -4.448873347981114, "beta_1": 0.0022410058067572937, "beta_2": -0.002196694445989066, "beta": [ -4.448873347981114, 0.0022410058067572937, -0.002196694445989066 ], "ci_beta_2": [ -0.002196694445989066, -0.002196694445989066 ], "knots": [ 1987.0, 1991.0, 2005.115616438356 ], "rss": 0.000109852647058465, "aic": -158.57599188816516, "rss_linear": 0.00015223689049339485, "aic_linear": -156.00782125762674, "delta_aic": -2.568170630538418, "aic_reduction": 2.568170630538418, "f_stat": 4.244109634755515, "p_f_test": 0.0638517595719145, "is_nonlinear_preferred": false, "r2": 0.91576701040457, "rmse": 0.002801181473420316, "fitted": "[0.01254082 0.02894488 0.00843365 0.02897933 0.01054748 0.01254082\n 0.01254082 0.01254082 0.01254082 0.02157713 0.00400519 0.01762825\n 0.0143821 0.03928847]", "residuals": "[-3.25734078e-04 6.14182177e-03 9.15351394e-04 4.80753668e-03\n -4.37251268e-04 -1.84778015e-03 -1.26411110e-03 5.90724983e-05\n 8.30981195e-04 -4.07167768e-03 2.05790585e-03 -1.61415212e-03\n -1.05019024e-03 -4.20177274e-03]", "n": 14 }, "calibrated_rate": 0.0010935447862130598, "calibrated_tmrca": 1979.7285531280427, "ci_mrca": [ 1974.6709484807375, 1983.0818095658824 ], "r2": 0.8832675519767318, "pagel_lambda": null }, "2": { "community_id": 2, "taxa_count": 12, "timespan": [ 1978.0, 1992.3387978142077 ], "ols": { "method": "OLS", "status": "OK", "mu": 0.002146671325446895, "d0": 0.012986609068794732, "t_ref": 1982.6948998178505, "t_mrca": 1976.64525033735, "se_mu": 0.00013649887570203593, "se_d0": 0.0005093562301261866, "se_mrca": 0.4519681215425009, "ci_fieller": [ 1975.4964988747404, 1977.5461580281751 ], "fieller_g": 0.02007296563655218, "ci_delta": [ 1975.6382026058916, 1977.6522980688085 ], "ci_bootstrap": null, "ci_mrca": [ 1975.4964988747404, 1977.5461580281751 ], "ci_method": "fieller", "r": 0.980377000504938, "r2": 0.9611390631190593, "p_value": 2.2172803371134364e-08, "sigma2": 3.1133252300203294e-06, "rmse": 0.001610725827595624, "residuals": "[ 3.99752340e-04 2.93441706e-03 7.99379028e-05 1.29983772e-03\n 1.65775379e-03 -7.53978698e-04 -7.42439815e-04 1.48688846e-03\n 2.12209473e-04 -2.16170130e-03 -1.94735116e-03 -2.46532576e-03]", "fitted": "[0.02008157 0.0179349 0.00720154 0.0179349 0.00720154 0.03368889\n 0.00720154 0.0029082 0.00934822 0.00934822 0.00934822 0.01364156]", "times": "[1986. 1985. 1980. 1985. 1980.\n 1992.33879781 1980. 1978. 1981. 1981.\n 1981. 1983. ]", "n": 12 }, "pgls": {}, "spline": { "method": "RESTRICTED_SPLINE", "t_mrca": 1975.876134808063, "ci_mrca": [ 1975.876134808063, 1975.876134808063 ], "rate_ancestral": 0.0017509589062900505, "ci_rate_ancestral": [ 0.0017509589062900505, 0.0017509589062900505 ], "rate_recent": 0.002256227346342468, "ci_rate_recent": [ 0.002256227346342468, 0.002256227346342468 ], "rate_ratio": 1.2885667037857476, "beta_0": -3.4596779159681383, "beta_1": 0.0017509589062900505, "beta_2": 0.0004435134084904606, "beta": [ -3.4596779159681383, 0.0017509589062900505, 0.0004435134084904606 ], "ci_beta_2": [ 0.0004435134084904606, 0.0004435134084904606 ], "knots": [ 1978.0, 1981.0, 1985.9 ], "rss": 2.976855546092809e-05, "aic": -148.88357468659632, "rss_linear": 3.1133252300203294e-05, "aic_linear": -150.34568904759405, "delta_aic": 1.462114360997731, "aic_reduction": -1.462114360997731, "f_stat": 0.41259212492180164, "p_f_test": 0.5366736404028668, "is_nonlinear_preferred": false, "r2": 0.962842495745414, "rmse": 0.0015750279219992707, "fitted": "[0.01993281 0.01767976 0.00727757 0.01767976 0.00727757 0.03423458\n 0.00727757 0.0037188 0.00916355 0.00916355 0.00916355 0.01327024]", "residuals": "[ 5.48514419e-04 3.18956336e-03 3.91273893e-06 1.55498402e-03\n 1.58172863e-03 -1.29967009e-03 -8.18464979e-04 6.76290042e-04\n 3.96874220e-04 -1.97703656e-03 -1.76268642e-03 -2.09400939e-03]", "n": 12 }, "calibrated_rate": 0.002146671325446895, "calibrated_tmrca": 1976.64525033735, "ci_mrca": [ 1975.4964988747404, 1977.5461580281751 ], "r2": 0.9611390631190593, "pagel_lambda": null }, "3": { "community_id": 3, "taxa_count": 18, "timespan": [ 2007.0, 2020.9453551912568 ], "ols": { "method": "OLS", "status": "OK", "mu": 0.00110249491183843, "d0": 0.007170431117603913, "t_ref": 2011.1414655953954, "t_mrca": 2004.6376431826034, "se_mu": 0.00015342119876596844, "se_d0": 0.0005958205370308166, "se_mrca": 1.054133635523376, "ci_fieller": [ 2001.5981672195874, 2006.4372071093746 ], "fieller_g": 0.08702629147364789, "ci_delta": [ 2002.4029797025703, 2006.8723066626364 ], "ci_bootstrap": null, "ci_mrca": [ 2001.5981672195874, 2006.4372071093746 ], "ci_method": "fieller", "r": 0.873757404514679, "r2": 0.7634520019442284, "p_value": 2.1669642934485367e-06, "sigma2": 6.390038022258433e-06, "rmse": 0.0023832821481232667, "residuals": "[ 0.00043672 -0.00101048 0.00029121 -0.00042234 0.00292058 -0.00101048\n -0.00027205 0.00078213 -0.00331916 0.00318617 0.00703423 -0.00139834\n 0.00059827 -0.00105712 -0.00263806 -0.00101048 -0.00290271 -0.00020806]", "fitted": "[0.01797917 0.00757326 0.00334452 0.00880677 0.00363449 0.00757326\n 0.00260449 0.00260449 0.00977523 0.00367376 0.01388618 0.01143653\n 0.00329619 0.00370698 0.00977523 0.00757326 0.00921946 0.00260449]", "times": "[2020.94535519 2011.50684932 2007.67123288 2012.62568306 2007.93424658\n 2011.50684932 2007. 2007. 2013.50410959 2007.96986301\n 2017.23287671 2015.0109589 2007.62739726 2008. 2013.50410959\n 2011.50684932 2013. 2007. ]", "n": 18 }, "pgls": {}, "spline": { "method": "RESTRICTED_SPLINE", "t_mrca": 1995.433692721854, "ci_mrca": [ 1995.433692721854, 1995.433692721854 ], "rate_ancestral": 0.0003199457495929273, "ci_rate_ancestral": [ 0.0003199457495929273, 0.0003199457495929273 ], "rate_recent": 0.0016485124360763016, "ci_rate_recent": [ 0.0016485124360763016, 0.0016485124360763016 ], "rate_ratio": 5.152474874799035, "beta_0": -0.6384305285808766, "beta_1": 0.0003199457495929273, "beta_2": 0.0008526786760078758, "beta": [ -0.6384305285808766, 0.0003199457495929273, 0.0008526786760078758 ], "ci_beta_2": [ 0.0008526786760078758, 0.0008526786760078758 ], "knots": [ 2007.0, 2011.5068493150684, 2015.6775342465753 ], "rss": 7.710725407946997e-05, "aic": -216.4923291554971, "rss_linear": 0.00010224060835613493, "aic_linear": -213.41396074952718, "delta_aic": -3.0783684059699112, "aic_reduction": 3.0783684059699112, "f_stat": 4.889297623819182, "p_f_test": 0.04296858865532803, "is_nonlinear_preferred": true, "r2": 0.8216015448133642, "rmse": 0.002069718903075568, "fitted": "[0.0208509 0.00617914 0.00391877 0.00748364 0.00400873 0.00617914\n 0.00370059 0.00370059 0.00870955 0.00402122 0.01473084 0.0110716\n 0.00390412 0.00403186 0.00870955 0.00617914 0.00798778 0.00370059]", "residuals": "[-2.43501422e-03 3.83641780e-04 -2.83050192e-04 9.00794722e-04\n 2.54633325e-03 3.83641780e-04 -1.36815544e-03 -3.13975715e-04\n -2.25347736e-03 2.83869564e-03 6.18956874e-03 -1.03341494e-03\n -9.66075418e-06 -1.38199482e-03 -1.57237715e-03 3.83641780e-04\n -1.67102785e-03 -1.30416924e-03]", "n": 18 }, "calibrated_rate": 0.00110249491183843, "calibrated_tmrca": 2004.6376431826034, "ci_mrca": [ 2001.5981672195874, 2006.4372071093746 ], "r2": 0.7634520019442284, "pagel_lambda": null }, "4": { "community_id": 4, "taxa_count": 18, "timespan": [ 1999.0, 2015.7205479452055 ], "ols": { "method": "OLS", "status": "OK", "mu": 0.0017368827016212994, "d0": 0.015606508928132309, "t_ref": 2007.2261251257996, "t_mrca": 1998.2407693999558, "se_mu": 9.653292120190284e-05, "se_d0": 0.00042018367105737273, "se_mrca": 0.554901140532882, "ci_fieller": [ 1996.9229588878763, 1999.0 ], "fieller_g": 0.013881705882045871, "ci_delta": [ 1997.0644315315963, 1999.0 ], "ci_bootstrap": null, "ci_mrca": [ 1996.9229588878763, 1999.0 ], "ci_method": "fieller", "r": 0.9761682815161132, "r2": 0.9529045138381217, "p_value": 4.855338353593197e-12, "sigma2": 3.1779777136185075e-06, "rmse": 0.001680734684229491, "residuals": "[ 4.32954710e-04 7.24670255e-04 -9.71675923e-04 -7.59968948e-05\n -1.39441849e-03 -1.45681419e-03 -1.62759629e-03 1.62201960e-03\n 4.92246704e-04 2.57333446e-03 2.25406899e-03 7.18838637e-04\n -7.59968948e-05 -4.50938562e-03 -7.59968948e-05 -1.39888290e-03\n 2.17198125e-03 5.96649476e-04]", "fitted": "[0.00947491 0.02433002 0.00131869 0.01173999 0.00826623 0.00652934\n 0.03036032 0.0223088 0.02404528 0.01868752 0.01521376 0.01686974\n 0.01173999 0.02302259 0.01173999 0.01306288 0.00826623 0.02394088]", "times": "[2003.69589041 2012.24863388 1999. 2005. 2003.\n 2002. 2015.72054795 2011.08493151 2012.08469945 2009.\n 2007. 2007.95342466 2005. 2011.49589041 2005.\n 2005.76164384 2003. 2012.02459016]", "n": 18 }, "pgls": {}, "spline": { "method": "RESTRICTED_SPLINE", "t_mrca": 1999.1342954242266, "ci_mrca": [ 1999.1342954242266, 1999.1342954242266 ], "rate_ancestral": 0.0021648344605581372, "ci_rate_ancestral": [ 0.0021648344605581372, 0.0021648344605581372 ], "rate_recent": 0.001494453645087632, "ci_rate_recent": [ 0.001494453645087632, 0.001494453645087632 ], "rate_ratio": 0.6903316037856918, "beta_0": -4.327794814017977, "beta_1": 0.0021648344605581372, "beta_2": -0.00039763868831088796, "beta": [ -4.327794814017977, 0.0021648344605581372, -0.00039763868831088796 ], "ci_beta_2": [ -0.00039763868831088796, -0.00039763868831088796 ], "knots": [ 1999.0, 2006.3808219178081, 2012.1338797814208 ], "rss": 4.485673106675283e-05, "aic": -226.2433558636232, "rss_linear": 5.084764341789612e-05, "aic_linear": -225.98687368827134, "delta_aic": -0.2564821753518629, "aic_reduction": 0.2564821753518629, "f_stat": 2.003348954105019, "p_f_test": 0.17737844658125063, "is_nonlinear_preferred": false, "r2": 0.958453343848028, "rmse": 0.0015786198449340492, "fitted": "[ 0.0096364 0.02409296 -0.00029073 0.01220036 0.00822108 0.00614154\n 0.02928158 0.02235045 0.02384797 0.01914701 0.01584895 0.01745791\n 0.01220036 0.02296725 0.01220036 0.01363449 0.00822108 0.02375814]", "residuals": "[ 0.00027146 0.00096173 0.00063775 -0.00053637 -0.00134927 -0.00106901\n -0.00054885 0.00158037 0.00068956 0.00211385 0.00161887 0.00013067\n -0.00053637 -0.00445405 -0.00053637 -0.0019705 0.00221713 0.00077939]", "n": 18 }, "calibrated_rate": 0.0017368827016212994, "calibrated_tmrca": 1998.2407693999558, "ci_mrca": [ 1996.9229588878763, 1999.0 ], "r2": 0.9529045138381217, "pagel_lambda": null }, "5": { "community_id": 5, "taxa_count": 4, "timespan": [ 1963.0, 1972.0 ], "ols": { "method": "OLS", "status": "MRCA_AFTER_EARLIEST_SAMPLE", "mu": 0.00722192079456151, "d0": 0.03215874636458593, "t_ref": 1968.75, "t_mrca": NaN, "se_mu": 0.0037741366595937815, "se_d0": 0.01317573756813896, "se_mrca": NaN, "ci_fieller": [ NaN, NaN ], "fieller_g": NaN, "ci_delta": [ NaN, NaN ], "ci_bootstrap": null, "ci_mrca": [ NaN, NaN ], "ci_method": "fieller", "r": 0.8042030847767735, "r2": 0.6467426015644784, "p_value": 0.1957969152236283, "sigma2": 0.0006944002418578733, "rmse": 0.018633306763130818, "residuals": "[ 0.0093673 -0.03222084 0.01235679 0.01049675]", "fitted": "[-0.0093673 0.03396423 0.04840807 0.05562999]", "times": "[1963. 1969. 1971. 1972.]", "n": 4 }, "pgls": {}, "spline": {}, "calibrated_rate": 0.00722192079456151, "calibrated_tmrca": NaN, "ci_mrca": [ NaN, NaN ], "r2": 0.6467426015644784, "pagel_lambda": null } }, "max_eigengap_value": 0.7277904748916626, "eigengaps": [ 0.7277904748916626, 0.14378488063812256, 0.06051468849182129, 0.04968225955963135, 0.01899111270904541, 0.002577066421508789, 0.0078507661819458, 0.0014787912368774414, 0.0008057355880737305, 0.000674128532409668 ], "sus_count": 2, "rooting_mode": "convex_decay", "contemporaneous_dyads_enabled": false, "n_contemporaneous_clusters": 0, "n_contemporaneous_taxa": 0, "summary_path": "/home/sweaver/programming/axomeme/reviews/moncla_h3nx_hoststrategy_REVIEW/REPRODUCTION/results/clock_communities_equine/autoclock_summary.json", "classified_metadata_path": "/home/sweaver/programming/axomeme/reviews/moncla_h3nx_hoststrategy_REVIEW/REPRODUCTION/results/clock_communities_equine/autoclock_classified_metadata.csv", "model_selection_path": "/home/sweaver/programming/axomeme/reviews/moncla_h3nx_hoststrategy_REVIEW/REPRODUCTION/results/clock_communities_equine/autoclock_model_selection.csv", "triage_path": "/home/sweaver/programming/axomeme/reviews/moncla_h3nx_hoststrategy_REVIEW/REPRODUCTION/results/clock_communities_equine/autoclock_sequence_triage.csv" }