
Output from mraic.pl version 1.4.4 by Johan Nylander
-------------------------------------------------------------
16:12 04/04/2011

Input data from file "683-tomodel-phylip/chr2_23133.phylip" (7 taxa, 620 characters)

Minimum AIC  model: GTR
Minimum AICc model: GTR
Minimum BIC  model: HKY


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
GTR	19	-1345.5874	2730.4414	0.3688	0.3688
GTRI	20	-1344.8651	2731.1325	0.2611	0.6299
GTRG	20	-1344.9876	2731.3775	0.2310	0.8608
GTRIG	21	-1344.8679	2733.2809	0.0892	0.9500
HKY	15	-1352.7649	2736.3245	0.0195	0.9695
HKYI	16	-1352.0653	2737.0327	0.0137	0.9831
HKYG	16	-1352.1835	2737.2692	0.0121	0.9953
HKYIG	17	-1352.0695	2739.1557	0.0047	1.0000
F81	14	-1377.0580	2782.8102	0.0000	1.0000
F81I	15	-1376.3941	2783.5828	0.0000	1.0000
F81G	15	-1376.4941	2783.7829	0.0000	1.0000
F81IG	16	-1376.3961	2785.6944	0.0000	1.0000
K2P	12	-1386.5834	2797.6808	0.0000	1.0000
K2PI	13	-1385.5508	2797.7022	0.0000	1.0000
K2PG	13	-1385.6785	2797.9576	0.0000	1.0000
SYM	16	-1382.9920	2798.8861	0.0000	1.0000
SYMI	17	-1382.0461	2799.1088	0.0000	1.0000
SYMG	17	-1382.1895	2799.3956	0.0000	1.0000
K2PIG	14	-1385.5550	2799.8042	0.0000	1.0000
SYMIG	18	-1382.0503	2801.2386	0.0000	1.0000
JC69	11	-1408.2574	2838.9489	0.0000	1.0000
JC69I	12	-1407.4345	2839.3830	0.0000	1.0000
JC69G	12	-1407.5306	2839.5753	0.0000	1.0000
JC69IG	13	-1407.4379	2841.4765	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTR	19	-1345.5874	2729.1747	0.3550	0.3550
GTRI	20	-1344.8651	2729.7302	0.2689	0.6239
GTRG	20	-1344.9876	2729.9752	0.2379	0.8618
GTRIG	21	-1344.8679	2731.7358	0.0987	0.9605
HKY	15	-1352.7649	2735.5298	0.0148	0.9753
HKYI	16	-1352.0653	2736.1306	0.0110	0.9862
HKYG	16	-1352.1835	2736.3670	0.0097	0.9960
HKYIG	17	-1352.0695	2738.1390	0.0040	1.0000
F81	14	-1377.0580	2782.1160	0.0000	1.0000
F81I	15	-1376.3941	2782.7881	0.0000	1.0000
F81G	15	-1376.4941	2782.9882	0.0000	1.0000
F81IG	16	-1376.3961	2784.7922	0.0000	1.0000
K2PI	13	-1385.5508	2797.1015	0.0000	1.0000
K2P	12	-1386.5834	2797.1668	0.0000	1.0000
K2PG	13	-1385.6785	2797.3570	0.0000	1.0000
SYM	16	-1382.9920	2797.9840	0.0000	1.0000
SYMI	17	-1382.0461	2798.0921	0.0000	1.0000
SYMG	17	-1382.1895	2798.3790	0.0000	1.0000
K2PIG	14	-1385.5550	2799.1100	0.0000	1.0000
SYMIG	18	-1382.0503	2800.1005	0.0000	1.0000
JC69	11	-1408.2574	2838.5147	0.0000	1.0000
JC69I	12	-1407.4345	2838.8690	0.0000	1.0000
JC69G	12	-1407.5306	2839.0613	0.0000	1.0000
JC69IG	13	-1407.4379	2840.8758	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKY	15	-1352.7649	2801.9756	0.8622	0.8622
HKYI	16	-1352.0653	2807.0061	0.0697	0.9319
HKYG	16	-1352.1835	2807.2425	0.0619	0.9938
GTR	19	-1345.5874	2813.3394	0.0029	0.9967
HKYIG	17	-1352.0695	2813.4443	0.0028	0.9995
GTRI	20	-1344.8651	2818.3246	0.0002	0.9998
GTRG	20	-1344.9876	2818.5696	0.0002	1.0000
GTRIG	21	-1344.8679	2824.7599	0.0000	1.0000
F81	14	-1377.0580	2844.1320	0.0000	1.0000
F81I	15	-1376.3941	2849.2339	0.0000	1.0000
F81G	15	-1376.4941	2849.4340	0.0000	1.0000
K2P	12	-1386.5834	2850.3234	0.0000	1.0000
K2PI	13	-1385.5508	2854.6879	0.0000	1.0000
K2PG	13	-1385.6785	2854.9433	0.0000	1.0000
F81IG	16	-1376.3961	2855.6677	0.0000	1.0000
K2PIG	14	-1385.5550	2861.1260	0.0000	1.0000
SYM	16	-1382.9920	2868.8595	0.0000	1.0000
SYMI	17	-1382.0461	2873.3974	0.0000	1.0000
SYMG	17	-1382.1895	2873.6843	0.0000	1.0000
SYMIG	18	-1382.0503	2879.8355	0.0000	1.0000
JC69	11	-1408.2574	2887.2416	0.0000	1.0000
JC69I	12	-1407.4345	2892.0256	0.0000	1.0000
JC69G	12	-1407.5306	2892.2179	0.0000	1.0000
JC69IG	13	-1407.4379	2898.4622	0.0000	1.0000



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

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

BEGIN MRBAYES;
 Charset 683-tomodel-phylip/chr2_23133.phylip = 1 - 620;
 Partition Dummy = 1:683-tomodel-phylip/chr2_23133.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=6 rates=equal;
 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);
END;


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

BEGIN MRBAYES;
 Charset 683-tomodel-phylip/chr2_23133.phylip = 1 - 620;
 Partition Dummy = 1:683-tomodel-phylip/chr2_23133.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=6 rates=equal;
 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);
END;


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

BEGIN MRBAYES;
 Charset 683-tomodel-phylip/chr2_23133.phylip = 1 - 620;
 Partition Dummy = 1:683-tomodel-phylip/chr2_23133.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=2 rates=equal;
 Prset applyto=(1) statefreqpr=Dirichlet(1.0,1.0,1.0,1.0) tratiopr=Beta(1.0,1.0);
END;


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

