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

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

Minimum AIC  model: HKYG
Minimum AICc model: HKYG
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
HKYG	16	-1390.1900	2813.4903	0.4510	0.4510
HKYI	16	-1390.8190	2814.7483	0.2404	0.6914
HKYIG	17	-1390.2095	2815.6706	0.1516	0.8430
GTRG	20	-1387.5425	2816.8133	0.0856	0.9286
GTRIG	21	-1387.2411	2818.3874	0.0390	0.9676
GTRI	20	-1388.5146	2818.7576	0.0324	1.0000
K2PG	13	-1411.1916	2849.1216	0.0000	1.0000
K2PI	13	-1411.6758	2850.0900	0.0000	1.0000
F81G	15	-1409.9564	2850.8904	0.0000	1.0000
K2PIG	14	-1411.3124	2851.4785	0.0000	1.0000
F81I	15	-1410.3593	2851.6962	0.0000	1.0000
F81IG	16	-1409.4531	2852.0164	0.0000	1.0000
HKY	15	-1411.2329	2853.4435	0.0000	1.0000
SYMG	17	-1409.3932	2854.0380	0.0000	1.0000
SYMI	17	-1409.9798	2855.2112	0.0000	1.0000
SYMIG	18	-1409.1138	2855.6293	0.0000	1.0000
GTR	19	-1408.8040	2857.1687	0.0000	1.0000
JC69G	12	-1428.3196	2881.2708	0.0000	1.0000
JC69I	12	-1428.7471	2882.1258	0.0000	1.0000
JC69IG	13	-1427.9858	2882.7099	0.0000	1.0000
F81	14	-1429.2699	2887.3935	0.0000	1.0000
K2P	12	-1431.8848	2888.4012	0.0000	1.0000
SYM	16	-1430.1924	2893.4950	0.0000	1.0000
JC69	11	-1447.5709	2917.6752	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
HKYG	16	-1390.1900	2812.3801	0.4204	0.4204
HKYI	16	-1390.8190	2813.6381	0.2241	0.6445
HKYIG	17	-1390.2095	2814.4190	0.1517	0.7961
GTRG	20	-1387.5425	2815.0849	0.1087	0.9048
GTRIG	21	-1387.2411	2816.4822	0.0541	0.9589
GTRI	20	-1388.5146	2817.0292	0.0411	1.0000
K2PG	13	-1411.1916	2848.3833	0.0000	1.0000
K2PI	13	-1411.6758	2849.3517	0.0000	1.0000
F81G	15	-1409.9564	2849.9128	0.0000	1.0000
K2PIG	14	-1411.3124	2850.6249	0.0000	1.0000
F81I	15	-1410.3593	2850.7186	0.0000	1.0000
F81IG	16	-1409.4531	2850.9062	0.0000	1.0000
HKY	15	-1411.2329	2852.4659	0.0000	1.0000
SYMG	17	-1409.3932	2852.7864	0.0000	1.0000
SYMI	17	-1409.9798	2853.9597	0.0000	1.0000
SYMIG	18	-1409.1138	2854.2277	0.0000	1.0000
GTR	19	-1408.8040	2855.6081	0.0000	1.0000
JC69G	12	-1428.3196	2880.6393	0.0000	1.0000
JC69I	12	-1428.7471	2881.4942	0.0000	1.0000
JC69IG	13	-1427.9858	2881.9716	0.0000	1.0000
F81	14	-1429.2699	2886.5399	0.0000	1.0000
K2P	12	-1431.8848	2887.7696	0.0000	1.0000
SYM	16	-1430.1924	2892.3848	0.0000	1.0000
JC69	11	-1447.5709	2917.1419	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	16	-1390.1900	2880.0363	0.6342	0.6342
HKYI	16	-1390.8190	2881.2942	0.3381	0.9723
HKYIG	17	-1390.2095	2886.3037	0.0276	0.9999
GTRG	20	-1387.5425	2899.6552	0.0000	1.0000
GTRI	20	-1388.5146	2901.5994	0.0000	1.0000
K2PG	13	-1411.1916	2903.3539	0.0000	1.0000
K2PI	13	-1411.6758	2904.3223	0.0000	1.0000
GTRIG	21	-1387.2411	2905.2809	0.0000	1.0000
K2PIG	14	-1411.3124	2909.8240	0.0000	1.0000
F81G	15	-1409.9564	2913.3405	0.0000	1.0000
F81I	15	-1410.3593	2914.1462	0.0000	1.0000
HKY	15	-1411.2329	2915.8935	0.0000	1.0000
F81IG	16	-1409.4531	2918.5623	0.0000	1.0000
SYMG	17	-1409.3932	2924.6711	0.0000	1.0000
SYMI	17	-1409.9798	2925.8443	0.0000	1.0000
SYMIG	18	-1409.1138	2930.3409	0.0000	1.0000
JC69G	12	-1428.3196	2931.3814	0.0000	1.0000
JC69I	12	-1428.7471	2932.2364	0.0000	1.0000
GTR	19	-1408.8040	2935.9498	0.0000	1.0000
JC69IG	13	-1427.9858	2936.9422	0.0000	1.0000
K2P	12	-1431.8848	2938.5118	0.0000	1.0000
F81	14	-1429.2699	2945.7390	0.0000	1.0000
SYM	16	-1430.1924	2960.0410	0.0000	1.0000
JC69	11	-1447.5709	2963.6555	0.0000	1.0000



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

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

BEGIN MRBAYES;
 Charset 591-tomodel-phylip/chr1_15365.phylip = 1 - 507;
 Partition Dummy = 1:591-tomodel-phylip/chr1_15365.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;


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

BEGIN MRBAYES;
 Charset 591-tomodel-phylip/chr1_15365.phylip = 1 - 507;
 Partition Dummy = 1:591-tomodel-phylip/chr1_15365.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;


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

BEGIN MRBAYES;
 Charset 591-tomodel-phylip/chr1_15365.phylip = 1 - 507;
 Partition Dummy = 1:591-tomodel-phylip/chr1_15365.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

