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

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

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


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

Model	df	lnL		AIC		wAIC	cwAIC
GTRIG	21	-3243.6145	6529.2290	0.7329	0.7329
GTRG	20	-3246.2204	6532.4408	0.1471	0.8799
GTRI	20	-3246.4273	6532.8546	0.1196	0.9995
SYMIG	18	-3254.0855	6544.1710	0.0004	1.0000
SYMG	17	-3257.8978	6549.7955	0.0000	1.0000
SYMI	17	-3258.3593	6550.7185	0.0000	1.0000
HKYIG	17	-3267.7972	6569.5944	0.0000	1.0000
HKYG	16	-3271.8948	6575.7896	0.0000	1.0000
HKYI	16	-3272.3362	6576.6723	0.0000	1.0000
K2PIG	14	-3288.0015	6604.0030	0.0000	1.0000
K2PG	13	-3291.5647	6609.1293	0.0000	1.0000
K2PI	13	-3291.8665	6609.7331	0.0000	1.0000
F81IG	16	-3341.4477	6714.8955	0.0000	1.0000
F81G	15	-3343.2608	6716.5216	0.0000	1.0000
F81I	15	-3343.7647	6717.5295	0.0000	1.0000
JC69IG	13	-3359.8194	6745.6388	0.0000	1.0000
JC69G	12	-3361.6648	6747.3295	0.0000	1.0000
JC69I	12	-3362.0604	6748.1208	0.0000	1.0000
GTR	19	-3382.2615	6802.5230	0.0000	1.0000
SYM	16	-3395.6456	6823.2913	0.0000	1.0000
HKY	15	-3404.6659	6839.3317	0.0000	1.0000
K2P	12	-3423.8556	6871.7113	0.0000	1.0000
F81	14	-3469.0359	6966.0717	0.0000	1.0000
JC69	11	-3487.1037	6996.2074	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
GTRG	20	-3246.2204	6633.1813	0.3142	0.3142
GTRI	20	-3246.4273	6633.5952	0.2555	0.5697
SYMIG	18	-3254.0855	6634.8375	0.1373	0.7070
GTRIG	21	-3243.6145	6635.0065	0.1262	0.8331
SYMG	17	-3257.8978	6635.4250	0.1023	0.9355
SYMI	17	-3258.3593	6636.3480	0.0645	1.0000
HKYIG	17	-3267.7972	6655.2239	0.0000	1.0000
HKYG	16	-3271.8948	6656.3820	0.0000	1.0000
HKYI	16	-3272.3362	6657.2648	0.0000	1.0000
K2PIG	14	-3288.0015	6674.5214	0.0000	1.0000
K2PG	13	-3291.5647	6674.6107	0.0000	1.0000
K2PI	13	-3291.8665	6675.2144	0.0000	1.0000
F81G	15	-3343.2608	6792.0771	0.0000	1.0000
F81I	15	-3343.7647	6793.0849	0.0000	1.0000
F81IG	16	-3341.4477	6795.4879	0.0000	1.0000
JC69G	12	-3361.6648	6807.7739	0.0000	1.0000
JC69I	12	-3362.0604	6808.5651	0.0000	1.0000
JC69IG	13	-3359.8194	6811.1202	0.0000	1.0000
GTR	19	-3382.2615	6898.2265	0.0000	1.0000
SYM	16	-3395.6456	6903.8837	0.0000	1.0000
HKY	15	-3404.6659	6914.8872	0.0000	1.0000
K2P	12	-3423.8556	6932.1556	0.0000	1.0000
F81	14	-3469.0359	7036.5901	0.0000	1.0000
JC69	11	-3487.1037	7051.6147	0.0000	1.0000



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

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

BEGIN MRBAYES;
 Charset 683-tomodel-phylip/chr2_24841.phylip = 1 - 1138;
 Partition Dummy = 1:683-tomodel-phylip/chr2_24841.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=6 rates=invgamma;
 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) pinvarpr=Uniform(0.0,1.0);
END;


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

BEGIN MRBAYES;
 Charset 683-tomodel-phylip/chr2_24841.phylip = 1 - 1138;
 Partition Dummy = 1:683-tomodel-phylip/chr2_24841.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=6 rates=invgamma;
 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) pinvarpr=Uniform(0.0,1.0);
END;


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

BEGIN MRBAYES;
 Charset 683-tomodel-phylip/chr2_24841.phylip = 1 - 1138;
 Partition Dummy = 1:683-tomodel-phylip/chr2_24841.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

