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

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

Minimum AIC  model: GTRI
Minimum AICc model: GTRI
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
GTRI	20	-1181.1355	2403.6458	0.2849	0.2849
GTRG	20	-1181.4376	2404.2500	0.2106	0.4955
HKYI	16	-1185.7682	2404.4209	0.1933	0.6888
HKYG	16	-1186.0517	2404.9880	0.1456	0.8344
GTRIG	21	-1181.1353	2405.7854	0.0977	0.9321
HKYIG	17	-1185.7680	2406.5327	0.0673	0.9994
GTR	19	-1188.9645	2417.1709	0.0003	0.9997
HKY	15	-1193.3998	2417.5787	0.0003	1.0000
SYMI	17	-1198.2473	2431.4914	0.0000	1.0000
K2PI	13	-1202.4783	2431.5457	0.0000	1.0000
SYMG	17	-1198.5492	2432.0952	0.0000	1.0000
K2PG	13	-1202.7586	2432.1061	0.0000	1.0000
SYMIG	18	-1198.1587	2433.4333	0.0000	1.0000
K2PIG	14	-1202.4477	2433.5760	0.0000	1.0000
F81I	15	-1202.8075	2436.3943	0.0000	1.0000
F81G	15	-1202.9735	2436.7262	0.0000	1.0000
F81IG	16	-1202.7711	2438.4268	0.0000	1.0000
K2P	12	-1210.0507	2444.6054	0.0000	1.0000
SYM	16	-1205.9644	2444.8134	0.0000	1.0000
F81	14	-1209.9186	2448.5179	0.0000	1.0000
JC69I	12	-1218.4398	2461.3836	0.0000	1.0000
JC69G	12	-1218.6120	2461.7280	0.0000	1.0000
JC69IG	13	-1218.3695	2463.3280	0.0000	1.0000
JC69	11	-1225.5712	2473.5683	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRI	20	-1181.1355	2402.2710	0.3091	0.3091
GTRG	20	-1181.4376	2402.8752	0.2285	0.5375
HKYI	16	-1185.7682	2403.5363	0.1642	0.7017
HKYG	16	-1186.0517	2404.1034	0.1236	0.8253
GTRIG	21	-1181.1353	2404.2707	0.1137	0.9390
HKYIG	17	-1185.7680	2405.5359	0.0604	0.9994
GTR	19	-1188.9645	2415.9291	0.0003	0.9998
HKY	15	-1193.3998	2416.7995	0.0002	1.0000
SYMI	17	-1198.2473	2430.4946	0.0000	1.0000
K2PI	13	-1202.4783	2430.9567	0.0000	1.0000
SYMG	17	-1198.5492	2431.0985	0.0000	1.0000
K2PG	13	-1202.7586	2431.5171	0.0000	1.0000
SYMIG	18	-1198.1587	2432.3175	0.0000	1.0000
K2PIG	14	-1202.4477	2432.8953	0.0000	1.0000
F81I	15	-1202.8075	2435.6151	0.0000	1.0000
F81G	15	-1202.9735	2435.9470	0.0000	1.0000
F81IG	16	-1202.7711	2437.5423	0.0000	1.0000
SYM	16	-1205.9644	2443.9288	0.0000	1.0000
K2P	12	-1210.0507	2444.1013	0.0000	1.0000
F81	14	-1209.9186	2447.8372	0.0000	1.0000
JC69I	12	-1218.4398	2460.8796	0.0000	1.0000
JC69G	12	-1218.6120	2461.2240	0.0000	1.0000
JC69IG	13	-1218.3695	2462.7390	0.0000	1.0000
JC69	11	-1225.5712	2473.1425	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYI	16	-1185.7682	2474.7185	0.5534	0.5534
HKYG	16	-1186.0517	2475.2857	0.4167	0.9701
HKYIG	17	-1185.7680	2481.1670	0.0220	0.9921
HKY	15	-1193.3998	2483.5328	0.0067	0.9989
K2PI	13	-1202.4783	2488.7922	0.0005	0.9994
K2PG	13	-1202.7586	2489.3527	0.0004	0.9997
GTRI	20	-1181.1355	2491.2488	0.0001	0.9999
GTRG	20	-1181.4376	2491.8530	0.0001	1.0000
K2PIG	14	-1202.4477	2495.1798	0.0000	1.0000
K2P	12	-1210.0507	2497.4880	0.0000	1.0000
GTRIG	21	-1181.1353	2497.6974	0.0000	1.0000
GTR	19	-1188.9645	2500.4580	0.0000	1.0000
F81I	15	-1202.8075	2502.3484	0.0000	1.0000
F81G	15	-1202.9735	2502.6804	0.0000	1.0000
SYMI	17	-1198.2473	2506.1257	0.0000	1.0000
SYMG	17	-1198.5492	2506.7296	0.0000	1.0000
F81IG	16	-1202.7711	2508.7245	0.0000	1.0000
F81	14	-1209.9186	2510.1217	0.0000	1.0000
SYMIG	18	-1198.1587	2512.3975	0.0000	1.0000
JC69I	12	-1218.4398	2514.2663	0.0000	1.0000
JC69G	12	-1218.6120	2514.6106	0.0000	1.0000
SYM	16	-1205.9644	2515.1111	0.0000	1.0000
JC69IG	13	-1218.3695	2520.5745	0.0000	1.0000
JC69	11	-1225.5712	2522.0803	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/chr5_11541.phylip = 1 - 632;
 Partition Dummy = 1:591-tomodel-phylip/chr5_11541.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 (GTRI)]

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

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

