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

Input data from file "591-tomodel-phylip/chr3_5533.phylip" (7 taxa, 639 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	-1647.6911	3328.2569	0.2566	0.2566
HKYI	16	-1647.7619	3328.3985	0.2391	0.4957
GTRI	20	-1643.8131	3328.9854	0.1783	0.6740
GTRG	20	-1643.8156	3328.9904	0.1778	0.8518
HKYIG	17	-1647.7155	3330.4165	0.0872	0.9389
GTRIG	21	-1643.8157	3331.1289	0.0610	1.0000
HKY	15	-1659.6305	3350.0314	0.0000	1.0000
GTR	19	-1655.6547	3350.5373	0.0000	1.0000
K2PG	13	-1674.4171	3375.4166	0.0000	1.0000
K2PI	13	-1674.4739	3375.5301	0.0000	1.0000
K2PIG	14	-1674.4388	3377.5507	0.0000	1.0000
SYMG	17	-1672.0945	3379.1744	0.0000	1.0000
SYMI	17	-1672.1531	3379.2917	0.0000	1.0000
SYMIG	18	-1672.1122	3381.3276	0.0000	1.0000
F81G	15	-1676.3737	3383.5179	0.0000	1.0000
F81I	15	-1676.4846	3383.7396	0.0000	1.0000
F81IG	16	-1676.3865	3385.6477	0.0000	1.0000
K2P	12	-1686.2839	3397.0662	0.0000	1.0000
SYM	16	-1683.7527	3400.3801	0.0000	1.0000
F81	14	-1687.2725	3403.2181	0.0000	1.0000
JC69G	12	-1699.6333	3423.7650	0.0000	1.0000
JC69I	12	-1699.7372	3423.9728	0.0000	1.0000
JC69IG	13	-1699.6446	3425.8716	0.0000	1.0000
JC69	11	-1710.6866	3443.7943	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
HKYG	16	-1647.6911	3327.3823	0.2281	0.2281
HKYI	16	-1647.7619	3327.5239	0.2125	0.4406
GTRI	20	-1643.8131	3327.6262	0.2019	0.6426
GTRG	20	-1643.8156	3327.6312	0.2014	0.8440
HKYIG	17	-1647.7155	3329.4310	0.0819	0.9259
GTRIG	21	-1643.8157	3329.6313	0.0741	1.0000
HKY	15	-1659.6305	3349.2610	0.0000	1.0000
GTR	19	-1655.6547	3349.3095	0.0000	1.0000
K2PG	13	-1674.4171	3374.8342	0.0000	1.0000
K2PI	13	-1674.4739	3374.9477	0.0000	1.0000
K2PIG	14	-1674.4388	3376.8777	0.0000	1.0000
SYMG	17	-1672.0945	3378.1889	0.0000	1.0000
SYMI	17	-1672.1531	3378.3062	0.0000	1.0000
SYMIG	18	-1672.1122	3380.2244	0.0000	1.0000
F81G	15	-1676.3737	3382.7474	0.0000	1.0000
F81I	15	-1676.4846	3382.9692	0.0000	1.0000
F81IG	16	-1676.3865	3384.7731	0.0000	1.0000
K2P	12	-1686.2839	3396.5678	0.0000	1.0000
SYM	16	-1683.7527	3399.5055	0.0000	1.0000
F81	14	-1687.2725	3402.5451	0.0000	1.0000
JC69G	12	-1699.6333	3423.2666	0.0000	1.0000
JC69I	12	-1699.7372	3423.4744	0.0000	1.0000
JC69IG	13	-1699.6446	3425.2892	0.0000	1.0000
JC69	11	-1710.6866	3443.3732	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	16	-1647.6911	3398.7407	0.5074	0.5074
HKYI	16	-1647.7619	3398.8824	0.4728	0.9802
HKYIG	17	-1647.7155	3405.2494	0.0196	0.9998
HKY	15	-1659.6305	3416.1595	0.0001	0.9999
GTRI	20	-1643.8131	3416.8243	0.0001	0.9999
GTRG	20	-1643.8156	3416.8293	0.0001	1.0000
GTRIG	21	-1643.8157	3423.2893	0.0000	1.0000
K2PG	13	-1674.4171	3432.8129	0.0000	1.0000
K2PI	13	-1674.4739	3432.9265	0.0000	1.0000
GTR	19	-1655.6547	3434.0477	0.0000	1.0000
K2PIG	14	-1674.4388	3439.3163	0.0000	1.0000
F81G	15	-1676.3737	3449.6460	0.0000	1.0000
F81I	15	-1676.4846	3449.8677	0.0000	1.0000
K2P	12	-1686.2839	3450.0867	0.0000	1.0000
SYMG	17	-1672.0945	3454.0073	0.0000	1.0000
SYMI	17	-1672.1531	3454.1246	0.0000	1.0000
F81IG	16	-1676.3865	3456.1316	0.0000	1.0000
SYMIG	18	-1672.1122	3460.5027	0.0000	1.0000
F81	14	-1687.2725	3464.9837	0.0000	1.0000
SYM	16	-1683.7527	3470.8640	0.0000	1.0000
JC69G	12	-1699.6333	3476.7855	0.0000	1.0000
JC69I	12	-1699.7372	3476.9933	0.0000	1.0000
JC69IG	13	-1699.6446	3483.2679	0.0000	1.0000
JC69	11	-1710.6866	3492.4321	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/chr3_5533.phylip = 1 - 639;
 Partition Dummy = 1:591-tomodel-phylip/chr3_5533.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/chr3_5533.phylip = 1 - 639;
 Partition Dummy = 1:591-tomodel-phylip/chr3_5533.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/chr3_5533.phylip = 1 - 639;
 Partition Dummy = 1:591-tomodel-phylip/chr3_5533.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

