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

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

Minimum AIC  model: HKYI
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	-1313.7353	2660.3341	0.3806	0.3806
HKYG	16	-1313.7988	2660.4611	0.3572	0.7378
HKYIG	17	-1313.8014	2662.5758	0.1241	0.8619
GTRI	20	-1311.4220	2664.1859	0.0555	0.9173
GTRG	20	-1311.5227	2664.3872	0.0502	0.9675
GTRIG	21	-1311.5226	2666.5235	0.0172	0.9847
HKY	15	-1318.3074	2667.3754	0.0113	0.9960
GTR	19	-1315.7271	2670.6663	0.0022	0.9982
F81I	15	-1321.2754	2673.3115	0.0006	0.9987
F81G	15	-1321.3557	2673.4722	0.0005	0.9993
K2PI	13	-1324.4066	2675.3883	0.0002	0.9995
K2PG	13	-1324.4717	2675.5184	0.0002	0.9997
F81IG	16	-1321.3557	2675.5749	0.0002	0.9999
K2PIG	14	-1324.5100	2677.6845	0.0001	0.9999
F81	14	-1325.0394	2678.7434	0.0000	1.0000
K2P	12	-1328.3837	2681.2596	0.0000	1.0000
SYMI	17	-1323.2947	2681.5625	0.0000	1.0000
SYMG	17	-1323.3789	2681.7309	0.0000	1.0000
SYMIG	18	-1323.4277	2683.9446	0.0000	1.0000
SYM	16	-1327.2648	2687.3931	0.0000	1.0000
JC69I	12	-1331.4806	2687.4533	0.0000	1.0000
JC69G	12	-1331.5618	2687.6157	0.0000	1.0000
JC69IG	13	-1331.6005	2689.7760	0.0000	1.0000
JC69	11	-1335.1692	2692.7542	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
HKYI	16	-1313.7353	2659.4706	0.3655	0.3655
HKYG	16	-1313.7988	2659.5976	0.3430	0.7084
HKYIG	17	-1313.8014	2661.6028	0.1258	0.8343
GTRI	20	-1311.4220	2662.8440	0.0677	0.9019
GTRG	20	-1311.5227	2663.0453	0.0612	0.9631
GTRIG	21	-1311.5226	2665.0451	0.0225	0.9856
HKY	15	-1318.3074	2666.6147	0.0103	0.9959
GTR	19	-1315.7271	2669.4542	0.0025	0.9984
F81I	15	-1321.2754	2672.5508	0.0005	0.9989
F81G	15	-1321.3557	2672.7115	0.0005	0.9994
F81IG	16	-1321.3557	2674.7114	0.0002	0.9995
K2PI	13	-1324.4066	2674.8132	0.0002	0.9997
K2PG	13	-1324.4717	2674.9433	0.0002	0.9999
K2PIG	14	-1324.5100	2677.0199	0.0001	0.9999
F81	14	-1325.0394	2678.0788	0.0000	1.0000
SYMI	17	-1323.2947	2680.5895	0.0000	1.0000
SYMG	17	-1323.3789	2680.7579	0.0000	1.0000
K2P	12	-1328.3837	2680.7675	0.0000	1.0000
SYMIG	18	-1323.4277	2682.8554	0.0000	1.0000
SYM	16	-1327.2648	2686.5296	0.0000	1.0000
JC69I	12	-1331.4806	2686.9612	0.0000	1.0000
JC69G	12	-1331.5618	2687.1236	0.0000	1.0000
JC69IG	13	-1331.6005	2689.2010	0.0000	1.0000
JC69	11	-1335.1692	2692.3384	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYI	16	-1313.7353	2731.0282	0.3097	0.3097
HKYG	16	-1313.7988	2731.1551	0.2906	0.6003
K2PI	13	-1324.4066	2732.9537	0.1182	0.7185
K2PG	13	-1324.4717	2733.0838	0.1108	0.8293
HKY	15	-1318.3074	2733.6999	0.0814	0.9108
K2P	12	-1328.3837	2734.4357	0.0564	0.9671
HKYIG	17	-1313.8014	2737.6327	0.0114	0.9785
K2PIG	14	-1324.5100	2739.6328	0.0042	0.9827
F81I	15	-1321.2754	2739.6360	0.0042	0.9869
F81G	15	-1321.3557	2739.7967	0.0039	0.9908
JC69I	12	-1331.4806	2740.6294	0.0025	0.9933
F81	14	-1325.0394	2740.6917	0.0025	0.9958
JC69G	12	-1331.5618	2740.7918	0.0023	0.9981
JC69	11	-1335.1692	2741.5342	0.0016	0.9997
F81IG	16	-1321.3557	2746.2690	0.0002	0.9999
JC69IG	13	-1331.6005	2747.3415	0.0001	1.0000
GTRI	20	-1311.4220	2752.2909	0.0000	1.0000
GTRG	20	-1311.5227	2752.4922	0.0000	1.0000
GTR	19	-1315.7271	2754.4288	0.0000	1.0000
SYMI	17	-1323.2947	2756.6194	0.0000	1.0000
SYMG	17	-1323.3789	2756.7878	0.0000	1.0000
SYM	16	-1327.2648	2758.0871	0.0000	1.0000
GTRIG	21	-1311.5226	2758.9644	0.0000	1.0000
SYMIG	18	-1323.4277	2763.3576	0.0000	1.0000



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

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

BEGIN MRBAYES;
 Charset 591-tomodel-phylip/chr1_32307.phylip = 1 - 647;
 Partition Dummy = 1:591-tomodel-phylip/chr1_32307.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 AICc model (HKYI)]

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

