
Output from mraic.pl version 1.4.4 by Johan Nylander
-------------------------------------------------------------
21:42 04/04/2011

Input data from file "591-tomodel-phylip/chr7_10269.phylip" (7 taxa, 831 characters)

Minimum AIC  model: GTRIG
Minimum AICc model: GTRIG
Minimum BIC  model: HKYG


Note: number of parameters compared to sample size
is low (Nchar/Nparams < 40). Consider the use of AICc.


(Output sorted by AICc. w: Akaike weight, cw: cumulative w.)

Model	df	lnL		AICc		wAICc	cwAICc
GTRIG	21	-2415.9917	4875.1256	0.4934	0.4934
GTRG	20	-2417.1367	4875.3105	0.4498	0.9431
GTRI	20	-2419.2385	4879.5140	0.0550	0.9981
HKYIG	17	-2425.9558	4886.6644	0.0015	0.9997
HKYG	16	-2428.6750	4890.0184	0.0003	1.0000
HKYI	16	-2430.5878	4893.8440	0.0000	1.0000
K2PIG	14	-2444.2836	4917.0819	0.0000	1.0000
SYMG	17	-2441.1947	4917.1422	0.0000	1.0000
SYMIG	18	-2440.2974	4917.4372	0.0000	1.0000
SYMI	17	-2442.7019	4920.1565	0.0000	1.0000
K2PG	13	-2446.9639	4920.3733	0.0000	1.0000
K2PI	13	-2448.6272	4923.7000	0.0000	1.0000
F81G	15	-2456.0838	4942.7566	0.0000	1.0000
F81IG	16	-2455.1983	4943.0650	0.0000	1.0000
F81I	15	-2457.9264	4946.4418	0.0000	1.0000
JC69G	12	-2471.7281	4967.8375	0.0000	1.0000
JC69IG	13	-2470.9117	4968.2690	0.0000	1.0000
JC69I	12	-2473.4274	4971.2361	0.0000	1.0000
GTR	19	-2471.4414	4981.8199	0.0000	1.0000
HKY	15	-2477.9352	4986.4593	0.0000	1.0000
K2P	12	-2496.4279	5017.2373	0.0000	1.0000
SYM	16	-2493.3928	5019.4538	0.0000	1.0000
F81	14	-2503.4780	5035.4708	0.0000	1.0000
JC69	11	-2519.4721	5061.2665	0.0000	1.0000



(Output sorted by AIC. w: Akaike weight, cw: cumulative w.)

Model	df	lnL		AIC		wAIC	cwAIC
GTRIG	21	-2415.9917	4873.9834	0.5066	0.5066
GTRG	20	-2417.1367	4874.2734	0.4382	0.9449
GTRI	20	-2419.2385	4878.4770	0.0536	0.9984
HKYIG	17	-2425.9558	4885.9117	0.0013	0.9997
HKYG	16	-2428.6750	4889.3501	0.0002	1.0000
HKYI	16	-2430.5878	4893.1757	0.0000	1.0000
SYMG	17	-2441.1947	4916.3894	0.0000	1.0000
K2PIG	14	-2444.2836	4916.5672	0.0000	1.0000
SYMIG	18	-2440.2974	4916.5948	0.0000	1.0000
SYMI	17	-2442.7019	4919.4038	0.0000	1.0000
K2PG	13	-2446.9639	4919.9278	0.0000	1.0000
K2PI	13	-2448.6272	4923.2544	0.0000	1.0000
F81G	15	-2456.0838	4942.1676	0.0000	1.0000
F81IG	16	-2455.1983	4942.3967	0.0000	1.0000
F81I	15	-2457.9264	4945.8528	0.0000	1.0000
JC69G	12	-2471.7281	4967.4561	0.0000	1.0000
JC69IG	13	-2470.9117	4967.8234	0.0000	1.0000
JC69I	12	-2473.4274	4970.8547	0.0000	1.0000
GTR	19	-2471.4414	4980.8828	0.0000	1.0000
HKY	15	-2477.9352	4985.8704	0.0000	1.0000
K2P	12	-2496.4279	5016.8559	0.0000	1.0000
SYM	16	-2493.3928	5018.7855	0.0000	1.0000
F81	14	-2503.4780	5034.9561	0.0000	1.0000
JC69	11	-2519.4721	5060.9442	0.0000	1.0000



(Output sorted by BIC. w: Akaike weight, cw: cumulative w.)

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	16	-2428.6750	4964.9121	0.5384	0.5384
HKYIG	17	-2425.9558	4966.1964	0.2833	0.8218
GTRG	20	-2417.1367	4968.7260	0.0800	0.9017
HKYI	16	-2430.5878	4968.7377	0.0795	0.9813
GTRI	20	-2419.2385	4972.9296	0.0098	0.9910
GTRIG	21	-2415.9917	4973.1586	0.0087	0.9998
K2PG	13	-2446.9639	4981.3219	0.0001	0.9999
K2PIG	14	-2444.2836	4982.6840	0.0001	1.0000
K2PI	13	-2448.6272	4984.6486	0.0000	1.0000
SYMG	17	-2441.1947	4996.6741	0.0000	1.0000
SYMI	17	-2442.7019	4999.6885	0.0000	1.0000
SYMIG	18	-2440.2974	5001.6022	0.0000	1.0000
F81G	15	-2456.0838	5013.0070	0.0000	1.0000
F81I	15	-2457.9264	5016.6922	0.0000	1.0000
F81IG	16	-2455.1983	5017.9587	0.0000	1.0000
JC69G	12	-2471.7281	5024.1277	0.0000	1.0000
JC69I	12	-2473.4274	5027.5263	0.0000	1.0000
JC69IG	13	-2470.9117	5029.2176	0.0000	1.0000
HKY	15	-2477.9352	5056.7098	0.0000	1.0000
GTR	19	-2471.4414	5070.6127	0.0000	1.0000
K2P	12	-2496.4279	5073.5274	0.0000	1.0000
SYM	16	-2493.3928	5094.3476	0.0000	1.0000
F81	14	-2503.4780	5101.0729	0.0000	1.0000
JC69	11	-2519.4721	5112.8931	0.0000	1.0000



To specify default values for best AIC, AICc, or BIC models in MrBayes:

[Mrbayes block for the best AIC model (GTRIG)]

BEGIN MRBAYES;
 Charset 591-tomodel-phylip/chr7_10269.phylip = 1 - 831;
 Partition Dummy = 1:591-tomodel-phylip/chr7_10269.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=6 rates=invgamma;
 Prset applyto=(1) revmatpr=Dirichlet(1.0,1.0,1.0,1.0,1.0,1.0) statefreqpr=Dirichlet(1.0,1.0,1.0,1.0) shapepr=Uniform(0.1,50.0) pinvarpr=Uniform(0.0,1.0);
END;


[Mrbayes block for the best AICc model (GTRIG)]

BEGIN MRBAYES;
 Charset 591-tomodel-phylip/chr7_10269.phylip = 1 - 831;
 Partition Dummy = 1:591-tomodel-phylip/chr7_10269.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=6 rates=invgamma;
 Prset applyto=(1) revmatpr=Dirichlet(1.0,1.0,1.0,1.0,1.0,1.0) statefreqpr=Dirichlet(1.0,1.0,1.0,1.0) shapepr=Uniform(0.1,50.0) pinvarpr=Uniform(0.0,1.0);
END;


[Mrbayes block for the best BIC model (HKYG)]

BEGIN MRBAYES;
 Charset 591-tomodel-phylip/chr7_10269.phylip = 1 - 831;
 Partition Dummy = 1:591-tomodel-phylip/chr7_10269.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=2 rates=gamma;
 Prset applyto=(1) statefreqpr=Dirichlet(1.0,1.0,1.0,1.0) tratiopr=Beta(1.0,1.0) shapepr=Uniform(0.1,50.0);
END;


-------------------------------------------------------------
End of Output

