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

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

Minimum AIC  model: GTRI
Minimum AICc model: GTRI
Minimum BIC  model: HKYI


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

Model	df	lnL		AIC		wAIC	cwAIC
GTRI	20	-3269.6280	6579.2561	0.4834	0.4834
GTRG	20	-3270.0371	6580.0742	0.3211	0.8045
GTRIG	21	-3269.6288	6581.2576	0.1777	0.9822
HKYI	16	-3277.6621	6587.3243	0.0086	0.9908
HKYG	16	-3278.0118	6588.0236	0.0060	0.9968
HKYIG	17	-3277.6650	6589.3300	0.0031	0.9999
GTR	19	-3279.5401	6597.0803	0.0001	1.0000
HKY	15	-3287.0885	6604.1769	0.0000	1.0000
K2PI	13	-3327.9984	6681.9968	0.0000	1.0000
K2PG	13	-3328.1892	6682.3784	0.0000	1.0000
K2PIG	14	-3327.9807	6683.9614	0.0000	1.0000
SYMI	17	-3326.0438	6686.0876	0.0000	1.0000
SYMG	17	-3326.3915	6686.7829	0.0000	1.0000
SYMIG	18	-3326.0469	6688.0938	0.0000	1.0000
K2P	12	-3338.3128	6700.6256	0.0000	1.0000
SYM	16	-3336.2620	6704.5239	0.0000	1.0000
F81I	15	-3367.3070	6764.6139	0.0000	1.0000
F81G	15	-3367.7349	6765.4698	0.0000	1.0000
F81IG	16	-3367.3113	6766.6225	0.0000	1.0000
F81	14	-3376.0619	6780.1239	0.0000	1.0000
JC69I	12	-3407.6148	6839.2295	0.0000	1.0000
JC69G	12	-3407.9778	6839.9556	0.0000	1.0000
JC69IG	13	-3407.6189	6841.2377	0.0000	1.0000
JC69	11	-3416.9607	6855.9214	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYI	16	-3277.6621	6665.5740	0.5732	0.5732
HKYG	16	-3278.0118	6666.2733	0.4040	0.9772
HKYIG	17	-3277.6650	6672.4703	0.0182	0.9954
GTRI	20	-3269.6280	6677.0683	0.0018	0.9973
HKY	15	-3287.0885	6677.5360	0.0014	0.9987
GTRG	20	-3270.0371	6677.8863	0.0012	0.9999
GTRIG	21	-3269.6288	6683.9604	0.0001	1.0000
GTR	19	-3279.5401	6690.0019	0.0000	1.0000
K2PI	13	-3327.9984	6745.5748	0.0000	1.0000
K2PG	13	-3328.1892	6745.9563	0.0000	1.0000
K2PIG	14	-3327.9807	6752.4299	0.0000	1.0000
K2P	12	-3338.3128	6759.3129	0.0000	1.0000
SYMI	17	-3326.0438	6769.2280	0.0000	1.0000
SYMG	17	-3326.3915	6769.9233	0.0000	1.0000
SYMIG	18	-3326.0469	6776.1248	0.0000	1.0000
SYM	16	-3336.2620	6782.7736	0.0000	1.0000
F81I	15	-3367.3070	6837.9730	0.0000	1.0000
F81G	15	-3367.7349	6838.8290	0.0000	1.0000
F81IG	16	-3367.3113	6844.8723	0.0000	1.0000
F81	14	-3376.0619	6848.5924	0.0000	1.0000
JC69I	12	-3407.6148	6897.9168	0.0000	1.0000
JC69G	12	-3407.9778	6898.6429	0.0000	1.0000
JC69IG	13	-3407.6189	6904.8157	0.0000	1.0000
JC69	11	-3416.9607	6909.7181	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 683-tomodel-phylip/chr1_32309.phylip = 1 - 983;
 Partition Dummy = 1:683-tomodel-phylip/chr1_32309.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 683-tomodel-phylip/chr1_32309.phylip = 1 - 983;
 Partition Dummy = 1:683-tomodel-phylip/chr1_32309.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 683-tomodel-phylip/chr1_32309.phylip = 1 - 983;
 Partition Dummy = 1:683-tomodel-phylip/chr1_32309.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

