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

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

Minimum AIC  model: HKYG
Minimum AICc model: HKYG
Minimum BIC  model: HKYG


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
HKYG	16	-1181.2854	2396.0214	0.3120	0.3120
HKYI	16	-1181.5519	2396.5544	0.2390	0.5511
GTRG	20	-1177.4876	2397.2393	0.1697	0.7208
GTRI	20	-1177.8527	2397.9696	0.1178	0.8386
HKYIG	17	-1181.3116	2398.2595	0.1019	0.9405
GTRIG	21	-1177.4197	2399.3366	0.0595	1.0000
F81G	15	-1198.4725	2428.2217	0.0000	1.0000
F81I	15	-1198.5359	2428.3483	0.0000	1.0000
GTR	19	-1194.6203	2429.2836	0.0000	1.0000
F81IG	16	-1198.4725	2430.3958	0.0000	1.0000
HKY	15	-1200.4382	2432.1531	0.0000	1.0000
K2PG	13	-1210.3986	2447.7602	0.0000	1.0000
K2PI	13	-1210.5056	2447.9742	0.0000	1.0000
K2PIG	14	-1210.3986	2449.9113	0.0000	1.0000
SYMG	17	-1208.1114	2451.8592	0.0000	1.0000
SYMI	17	-1208.1157	2451.8679	0.0000	1.0000
SYMIG	18	-1208.1618	2454.1575	0.0000	1.0000
F81	14	-1216.7466	2462.6072	0.0000	1.0000
JC69G	12	-1224.7628	2474.3488	0.0000	1.0000
JC69I	12	-1224.7681	2474.3594	0.0000	1.0000
JC69IG	13	-1224.7628	2476.4886	0.0000	1.0000
K2P	12	-1229.7711	2484.3654	0.0000	1.0000
SYM	16	-1227.1371	2487.7249	0.0000	1.0000
JC69	11	-1243.5159	2509.7266	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
HKYG	16	-1181.2854	2394.5707	0.2611	0.2611
GTRG	20	-1177.4876	2394.9752	0.2133	0.4744
HKYI	16	-1181.5519	2395.1037	0.2000	0.6744
GTRI	20	-1177.8527	2395.7054	0.1480	0.8225
HKYIG	17	-1181.3116	2396.6232	0.0936	0.9160
GTRIG	21	-1177.4197	2396.8393	0.0840	1.0000
F81G	15	-1198.4725	2426.9451	0.0000	1.0000
F81I	15	-1198.5359	2427.0717	0.0000	1.0000
GTR	19	-1194.6203	2427.2406	0.0000	1.0000
F81IG	16	-1198.4725	2428.9451	0.0000	1.0000
HKY	15	-1200.4382	2430.8765	0.0000	1.0000
K2PG	13	-1210.3986	2446.7973	0.0000	1.0000
K2PI	13	-1210.5056	2447.0113	0.0000	1.0000
K2PIG	14	-1210.3986	2448.7972	0.0000	1.0000
SYMG	17	-1208.1114	2450.2228	0.0000	1.0000
SYMI	17	-1208.1157	2450.2315	0.0000	1.0000
SYMIG	18	-1208.1618	2452.3237	0.0000	1.0000
F81	14	-1216.7466	2461.4931	0.0000	1.0000
JC69G	12	-1224.7628	2473.5256	0.0000	1.0000
JC69I	12	-1224.7681	2473.5362	0.0000	1.0000
JC69IG	13	-1224.7628	2475.5256	0.0000	1.0000
K2P	12	-1229.7711	2483.5422	0.0000	1.0000
SYM	16	-1227.1371	2486.2742	0.0000	1.0000
JC69	11	-1243.5159	2509.0318	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	16	-1181.2854	2458.1109	0.5507	0.5507
HKYI	16	-1181.5519	2458.6439	0.4219	0.9726
HKYIG	17	-1181.3116	2464.1346	0.0271	0.9997
GTRG	20	-1177.4876	2474.4004	0.0002	0.9999
GTRI	20	-1177.8527	2475.1307	0.0001	1.0000
GTRIG	21	-1177.4197	2480.2358	0.0000	1.0000
F81G	15	-1198.4725	2486.5140	0.0000	1.0000
F81I	15	-1198.5359	2486.6407	0.0000	1.0000
HKY	15	-1200.4382	2490.4454	0.0000	1.0000
F81IG	16	-1198.4725	2492.4853	0.0000	1.0000
K2PG	13	-1210.3986	2498.4237	0.0000	1.0000
K2PI	13	-1210.5056	2498.6377	0.0000	1.0000
GTR	19	-1194.6203	2502.6946	0.0000	1.0000
K2PIG	14	-1210.3986	2504.3949	0.0000	1.0000
F81	14	-1216.7466	2517.0908	0.0000	1.0000
SYMG	17	-1208.1114	2517.7343	0.0000	1.0000
SYMI	17	-1208.1157	2517.7430	0.0000	1.0000
JC69G	12	-1224.7628	2521.1807	0.0000	1.0000
JC69I	12	-1224.7681	2521.1913	0.0000	1.0000
SYMIG	18	-1208.1618	2523.8064	0.0000	1.0000
JC69IG	13	-1224.7628	2527.1520	0.0000	1.0000
K2P	12	-1229.7711	2531.1974	0.0000	1.0000
SYM	16	-1227.1371	2549.8144	0.0000	1.0000
JC69	11	-1243.5159	2552.7157	0.0000	1.0000



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

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

BEGIN MRBAYES;
 Charset 591-tomodel-phylip/chr4_7551.phylip = 1 - 392;
 Partition Dummy = 1:591-tomodel-phylip/chr4_7551.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=2 rates=gamma;
 Prset applyto=(1) statefreqpr=Dirichlet(1.0,1.0,1.0,1.0) tratiopr=Beta(1.0,1.0) shapepr=Uniform(0.1,50.0);
END;


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

BEGIN MRBAYES;
 Charset 591-tomodel-phylip/chr4_7551.phylip = 1 - 392;
 Partition Dummy = 1:591-tomodel-phylip/chr4_7551.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=2 rates=gamma;
 Prset applyto=(1) statefreqpr=Dirichlet(1.0,1.0,1.0,1.0) tratiopr=Beta(1.0,1.0) shapepr=Uniform(0.1,50.0);
END;


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

BEGIN MRBAYES;
 Charset 591-tomodel-phylip/chr4_7551.phylip = 1 - 392;
 Partition Dummy = 1:591-tomodel-phylip/chr4_7551.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=2 rates=gamma;
 Prset applyto=(1) statefreqpr=Dirichlet(1.0,1.0,1.0,1.0) tratiopr=Beta(1.0,1.0) shapepr=Uniform(0.1,50.0);
END;


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

