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

Input data from file "591-tomodel-phylip/chr2_7410.phylip" (7 taxa, 701 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	-1248.1423	2529.0798	0.3477	0.3477
HKYG	16	-1248.2859	2529.3671	0.3012	0.6489
HKYIG	17	-1248.1416	2531.1792	0.1217	0.7706
GTRI	20	-1245.1742	2531.5838	0.0994	0.8700
GTRG	20	-1245.3513	2531.9380	0.0833	0.9533
GTRIG	21	-1245.1740	2533.7088	0.0344	0.9876
HKY	15	-1252.7692	2536.2391	0.0097	0.9973
GTR	19	-1250.0240	2539.1641	0.0022	0.9996
F81I	15	-1256.7011	2544.1029	0.0002	0.9998
F81G	15	-1256.7891	2544.2789	0.0002	0.9999
F81IG	16	-1256.7005	2546.1964	0.0001	1.0000
F81	14	-1261.0685	2550.7492	0.0000	1.0000
K2PI	13	-1271.5635	2569.6569	0.0000	1.0000
K2PG	13	-1271.6349	2569.7997	0.0000	1.0000
K2PIG	14	-1271.5631	2571.7384	0.0000	1.0000
SYMI	17	-1269.1157	2573.1274	0.0000	1.0000
SYMG	17	-1269.1766	2573.2492	0.0000	1.0000
SYMIG	18	-1269.1160	2575.2349	0.0000	1.0000
K2P	12	-1275.8488	2576.1511	0.0000	1.0000
SYM	16	-1273.3862	2579.5677	0.0000	1.0000
JC69I	12	-1279.2803	2583.0141	0.0000	1.0000
JC69G	12	-1279.3226	2583.0986	0.0000	1.0000
JC69IG	13	-1279.2796	2585.0891	0.0000	1.0000
JC69	11	-1283.3898	2589.1627	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
HKYI	16	-1248.1423	2528.2845	0.3273	0.3273
HKYG	16	-1248.2859	2528.5718	0.2835	0.6107
HKYIG	17	-1248.1416	2530.2831	0.1205	0.7312
GTRI	20	-1245.1742	2530.3485	0.1166	0.8478
GTRG	20	-1245.3513	2530.7027	0.0977	0.9455
GTRIG	21	-1245.1740	2532.3479	0.0429	0.9884
HKY	15	-1252.7692	2535.5384	0.0087	0.9971
GTR	19	-1250.0240	2538.0481	0.0025	0.9996
F81I	15	-1256.7011	2543.4021	0.0002	0.9998
F81G	15	-1256.7891	2543.5782	0.0002	0.9999
F81IG	16	-1256.7005	2545.4011	0.0001	1.0000
F81	14	-1261.0685	2550.1370	0.0000	1.0000
K2PI	13	-1271.5635	2569.1271	0.0000	1.0000
K2PG	13	-1271.6349	2569.2699	0.0000	1.0000
K2PIG	14	-1271.5631	2571.1261	0.0000	1.0000
SYMI	17	-1269.1157	2572.2314	0.0000	1.0000
SYMG	17	-1269.1766	2572.3531	0.0000	1.0000
SYMIG	18	-1269.1160	2574.2319	0.0000	1.0000
K2P	12	-1275.8488	2575.6976	0.0000	1.0000
SYM	16	-1273.3862	2578.7724	0.0000	1.0000
JC69I	12	-1279.2803	2582.5606	0.0000	1.0000
JC69G	12	-1279.3226	2582.6452	0.0000	1.0000
JC69IG	13	-1279.2796	2584.5593	0.0000	1.0000
JC69	11	-1283.3898	2588.7795	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYI	16	-1248.1423	2601.1246	0.4598	0.4598
HKYG	16	-1248.2859	2601.4119	0.3983	0.8581
HKY	15	-1252.7692	2603.8260	0.1191	0.9772
HKYIG	17	-1248.1416	2607.6758	0.0174	0.9946
F81I	15	-1256.7011	2611.6898	0.0023	0.9969
F81G	15	-1256.7891	2611.8658	0.0021	0.9991
F81	14	-1261.0685	2613.8721	0.0008	0.9999
F81IG	16	-1256.7005	2618.2412	0.0001	1.0000
GTRI	20	-1245.1742	2621.3986	0.0000	1.0000
GTRG	20	-1245.3513	2621.7528	0.0000	1.0000
GTR	19	-1250.0240	2624.5457	0.0000	1.0000
GTRIG	21	-1245.1740	2627.9506	0.0000	1.0000
K2PI	13	-1271.5635	2628.3097	0.0000	1.0000
K2PG	13	-1271.6349	2628.4525	0.0000	1.0000
K2P	12	-1275.8488	2630.3277	0.0000	1.0000
K2PIG	14	-1271.5631	2634.8612	0.0000	1.0000
JC69I	12	-1279.2803	2637.1907	0.0000	1.0000
JC69G	12	-1279.3226	2637.2753	0.0000	1.0000
JC69	11	-1283.3898	2638.8571	0.0000	1.0000
JC69IG	13	-1279.2796	2643.7419	0.0000	1.0000
SYMI	17	-1269.1157	2649.6240	0.0000	1.0000
SYMG	17	-1269.1766	2649.7458	0.0000	1.0000
SYM	16	-1273.3862	2651.6125	0.0000	1.0000
SYMIG	18	-1269.1160	2656.1771	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/chr2_7410.phylip = 1 - 701;
 Partition Dummy = 1:591-tomodel-phylip/chr2_7410.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/chr2_7410.phylip = 1 - 701;
 Partition Dummy = 1:591-tomodel-phylip/chr2_7410.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/chr2_7410.phylip = 1 - 701;
 Partition Dummy = 1:591-tomodel-phylip/chr2_7410.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

