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

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

Minimum AIC  model: GTRG
Minimum AICc model: GTRG
Minimum BIC  model: GTRG


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

Model	df	lnL		AIC		wAIC	cwAIC
GTRG	20	-2781.7048	5603.4096	0.5637	0.5637
GTRIG	21	-2781.3111	5604.6222	0.3074	0.8711
GTRI	20	-2783.1816	5606.3632	0.1287	0.9999
SYMG	17	-2793.7972	5621.5943	0.0001	0.9999
SYMIG	18	-2793.4479	5622.8958	0.0000	1.0000
SYMI	17	-2795.2374	5624.4748	0.0000	1.0000
HKYG	16	-2796.5380	5625.0759	0.0000	1.0000
HKYIG	17	-2796.3446	5626.6893	0.0000	1.0000
HKYI	16	-2797.6420	5627.2840	0.0000	1.0000
K2PG	13	-2808.5733	5643.1466	0.0000	1.0000
K2PIG	14	-2808.3847	5644.7693	0.0000	1.0000
K2PI	13	-2809.6607	5645.3215	0.0000	1.0000
GTR	19	-2814.4002	5666.8004	0.0000	1.0000
SYM	16	-2824.5748	5681.1496	0.0000	1.0000
HKY	15	-2826.1947	5682.3895	0.0000	1.0000
K2P	12	-2837.9916	5699.9831	0.0000	1.0000
F81G	15	-2840.0470	5710.0941	0.0000	1.0000
F81I	15	-2840.4414	5710.8829	0.0000	1.0000
F81IG	16	-2839.4985	5710.9969	0.0000	1.0000
JC69G	12	-2852.8578	5729.7157	0.0000	1.0000
JC69I	12	-2853.2557	5730.5114	0.0000	1.0000
JC69IG	13	-2852.3171	5730.6343	0.0000	1.0000
F81	14	-2870.2108	5768.4215	0.0000	1.0000
JC69	11	-2882.3689	5786.7378	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
GTRG	20	-2781.7048	5701.8822	0.4724	0.4724
HKYG	16	-2796.5380	5703.8540	0.1763	0.6486
GTRI	20	-2783.1816	5704.8357	0.1079	0.7565
SYMG	17	-2793.7972	5705.2960	0.0857	0.8422
HKYI	16	-2797.6420	5706.0621	0.0584	0.9007
K2PG	13	-2808.5733	5707.1537	0.0339	0.9345
GTRIG	21	-2781.3111	5708.0184	0.0220	0.9565
SYMI	17	-2795.2374	5708.1764	0.0203	0.9768
K2PI	13	-2809.6607	5709.3287	0.0114	0.9882
HKYIG	17	-2796.3446	5710.3910	0.0067	0.9949
SYMIG	18	-2793.4479	5711.5211	0.0038	0.9987
K2PIG	14	-2808.3847	5713.7001	0.0013	1.0000
HKY	15	-2826.1947	5756.2439	0.0000	1.0000
K2P	12	-2837.9916	5759.0667	0.0000	1.0000
SYM	16	-2824.5748	5759.9276	0.0000	1.0000
GTR	19	-2814.4002	5760.3493	0.0000	1.0000
F81G	15	-2840.0470	5783.9485	0.0000	1.0000
F81I	15	-2840.4414	5784.7373	0.0000	1.0000
JC69G	12	-2852.8578	5788.7992	0.0000	1.0000
JC69I	12	-2853.2557	5789.5949	0.0000	1.0000
F81IG	16	-2839.4985	5789.7750	0.0000	1.0000
JC69IG	13	-2852.3171	5794.6415	0.0000	1.0000
F81	14	-2870.2108	5837.3523	0.0000	1.0000
JC69	11	-2882.3689	5840.8977	0.0000	1.0000



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

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

BEGIN MRBAYES;
 Charset 591-tomodel-phylip/chr4_9746.phylip = 1 - 1016;
 Partition Dummy = 1:591-tomodel-phylip/chr4_9746.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=6 rates=gamma;
 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) shapepr=Uniform(0.1,50.0);
END;


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

BEGIN MRBAYES;
 Charset 591-tomodel-phylip/chr4_9746.phylip = 1 - 1016;
 Partition Dummy = 1:591-tomodel-phylip/chr4_9746.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=6 rates=gamma;
 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) shapepr=Uniform(0.1,50.0);
END;


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

BEGIN MRBAYES;
 Charset 591-tomodel-phylip/chr4_9746.phylip = 1 - 1016;
 Partition Dummy = 1:591-tomodel-phylip/chr4_9746.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=6 rates=gamma;
 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) shapepr=Uniform(0.1,50.0);
END;


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

