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

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

Minimum AIC  model: GTRI
Minimum AICc model: HKYI
Minimum BIC  model: HKYI


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
HKYI	16	-1433.4725	2899.9734	0.3268	0.3268
GTRI	20	-1429.4423	2900.4845	0.2531	0.5799
HKYG	16	-1434.0821	2901.1926	0.1776	0.7576
GTRG	20	-1430.0682	2901.7364	0.1353	0.8929
HKYIG	17	-1434.0822	2903.3236	0.0612	0.9541
GTRIG	21	-1430.0682	2903.8998	0.0459	1.0000
GTR	19	-1456.2710	2951.9869	0.0000	1.0000
F81I	15	-1461.2502	2953.4060	0.0000	1.0000
F81G	15	-1461.6336	2954.1728	0.0000	1.0000
HKY	15	-1461.7144	2954.3344	0.0000	1.0000
F81IG	16	-1461.6337	2956.2957	0.0000	1.0000
F81	14	-1485.2565	2999.3040	0.0000	1.0000
SYMI	17	-1482.6449	3000.4489	0.0000	1.0000
SYMG	17	-1482.8578	3000.8748	0.0000	1.0000
SYMIG	18	-1482.8579	3003.0137	0.0000	1.0000
K2PI	13	-1488.9950	3004.6742	0.0000	1.0000
K2PG	13	-1489.3216	3005.3274	0.0000	1.0000
K2PIG	14	-1489.3217	3007.4343	0.0000	1.0000
JC69I	12	-1511.6404	3047.8662	0.0000	1.0000
JC69G	12	-1511.9047	3048.3949	0.0000	1.0000
SYM	16	-1508.4580	3049.9443	0.0000	1.0000
JC69IG	13	-1511.9048	3050.4939	0.0000	1.0000
K2P	12	-1515.3310	3055.2474	0.0000	1.0000
JC69	11	-1536.1643	3094.8229	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRI	20	-1429.4423	2898.8845	0.2921	0.2921
HKYI	16	-1433.4725	2898.9450	0.2834	0.5756
GTRG	20	-1430.0682	2900.1364	0.1562	0.7318
HKYG	16	-1434.0821	2900.1643	0.1541	0.8859
GTRIG	21	-1430.0682	2902.1365	0.0575	0.9433
HKYIG	17	-1434.0822	2902.1645	0.0567	1.0000
GTR	19	-1456.2710	2950.5420	0.0000	1.0000
F81I	15	-1461.2502	2952.5003	0.0000	1.0000
F81G	15	-1461.6336	2953.2672	0.0000	1.0000
HKY	15	-1461.7144	2953.4288	0.0000	1.0000
F81IG	16	-1461.6337	2955.2674	0.0000	1.0000
F81	14	-1485.2565	2998.5130	0.0000	1.0000
SYMI	17	-1482.6449	2999.2898	0.0000	1.0000
SYMG	17	-1482.8578	2999.7157	0.0000	1.0000
SYMIG	18	-1482.8579	3001.7158	0.0000	1.0000
K2PI	13	-1488.9950	3003.9900	0.0000	1.0000
K2PG	13	-1489.3216	3004.6432	0.0000	1.0000
K2PIG	14	-1489.3217	3006.6434	0.0000	1.0000
JC69I	12	-1511.6404	3047.2809	0.0000	1.0000
JC69G	12	-1511.9047	3047.8095	0.0000	1.0000
SYM	16	-1508.4580	3048.9159	0.0000	1.0000
JC69IG	13	-1511.9048	3049.8097	0.0000	1.0000
K2P	12	-1515.3310	3054.6621	0.0000	1.0000
JC69	11	-1536.1643	3094.3285	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYI	16	-1433.4725	2967.7869	0.6381	0.6381
HKYG	16	-1434.0821	2969.0062	0.3469	0.9850
HKYIG	17	-1434.0822	2975.3090	0.0148	0.9998
GTRI	20	-1429.4423	2984.9369	0.0001	0.9999
GTRG	20	-1430.0682	2986.1888	0.0001	1.0000
GTRIG	21	-1430.0682	2992.4915	0.0000	1.0000
F81I	15	-1461.2502	3017.0396	0.0000	1.0000
F81G	15	-1461.6336	3017.8065	0.0000	1.0000
HKY	15	-1461.7144	3017.9680	0.0000	1.0000
F81IG	16	-1461.6337	3024.1093	0.0000	1.0000
GTR	19	-1456.2710	3032.2918	0.0000	1.0000
F81	14	-1485.2565	3058.7497	0.0000	1.0000
K2PI	13	-1488.9950	3059.9240	0.0000	1.0000
K2PG	13	-1489.3216	3060.5773	0.0000	1.0000
K2PIG	14	-1489.3217	3066.8800	0.0000	1.0000
SYMI	17	-1482.6449	3072.4344	0.0000	1.0000
SYMG	17	-1482.8578	3072.8602	0.0000	1.0000
SYMIG	18	-1482.8579	3079.1630	0.0000	1.0000
JC69I	12	-1511.6404	3098.9123	0.0000	1.0000
JC69G	12	-1511.9047	3099.4409	0.0000	1.0000
JC69IG	13	-1511.9048	3105.7437	0.0000	1.0000
K2P	12	-1515.3310	3106.2935	0.0000	1.0000
SYM	16	-1508.4580	3117.7578	0.0000	1.0000
JC69	11	-1536.1643	3141.6573	0.0000	1.0000



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

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

BEGIN MRBAYES;
 Charset 591-tomodel-phylip/chr13_2779.phylip = 1 - 546;
 Partition Dummy = 1:591-tomodel-phylip/chr13_2779.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=6 rates=propinv;
 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) pinvarpr=Uniform(0.0,1.0);
END;


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

BEGIN MRBAYES;
 Charset 591-tomodel-phylip/chr13_2779.phylip = 1 - 546;
 Partition Dummy = 1:591-tomodel-phylip/chr13_2779.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=2 rates=propinv;
 Prset applyto=(1) statefreqpr=Dirichlet(1.0,1.0,1.0,1.0) tratiopr=Beta(1.0,1.0) pinvarpr=Uniform(0.0,1.0);
END;


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

BEGIN MRBAYES;
 Charset 591-tomodel-phylip/chr13_2779.phylip = 1 - 546;
 Partition Dummy = 1:591-tomodel-phylip/chr13_2779.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=2 rates=propinv;
 Prset applyto=(1) statefreqpr=Dirichlet(1.0,1.0,1.0,1.0) tratiopr=Beta(1.0,1.0) pinvarpr=Uniform(0.0,1.0);
END;


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

