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

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

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


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

Model	df	lnL		AIC		wAIC	cwAIC
GTRG	20	-3208.2061	6456.4121	0.6236	0.6236
GTRIG	21	-3208.1975	6458.3950	0.2314	0.8550
GTRI	20	-3209.6856	6459.3712	0.1420	0.9971
HKYG	16	-3218.1245	6468.2489	0.0017	0.9987
HKYIG	17	-3218.0903	6470.1806	0.0006	0.9994
HKYI	16	-3219.1202	6470.2404	0.0006	1.0000
SYMG	17	-3228.8786	6491.7572	0.0000	1.0000
K2PG	13	-3233.0752	6492.1504	0.0000	1.0000
SYMIG	18	-3228.8375	6493.6751	0.0000	1.0000
SYMI	17	-3229.8539	6493.7078	0.0000	1.0000
K2PI	13	-3233.9294	6493.8589	0.0000	1.0000
K2PIG	14	-3232.9982	6493.9964	0.0000	1.0000
GTR	19	-3243.5183	6525.0366	0.0000	1.0000
HKY	15	-3254.8167	6539.6334	0.0000	1.0000
SYM	16	-3265.9383	6563.8767	0.0000	1.0000
K2P	12	-3271.7418	6567.4836	0.0000	1.0000
F81G	15	-3280.4386	6590.8772	0.0000	1.0000
F81I	15	-3280.6932	6591.3863	0.0000	1.0000
F81IG	16	-3280.4692	6592.9383	0.0000	1.0000
JC69G	12	-3292.8098	6609.6196	0.0000	1.0000
JC69I	12	-3293.0730	6610.1459	0.0000	1.0000
JC69IG	13	-3292.8292	6611.6584	0.0000	1.0000
F81	14	-3314.2239	6656.4478	0.0000	1.0000
JC69	11	-3327.6964	6677.3929	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	16	-3218.1245	6547.3852	0.6940	0.6940
HKYI	16	-3219.1202	6549.3767	0.2564	0.9504
HKYIG	17	-3218.0903	6554.2629	0.0223	0.9727
GTRG	20	-3208.2061	6555.3324	0.0131	0.9857
K2PG	13	-3233.0752	6556.4486	0.0075	0.9932
K2PI	13	-3233.9294	6558.1571	0.0032	0.9964
GTRI	20	-3209.6856	6558.2915	0.0030	0.9993
GTRIG	21	-3208.1975	6562.2613	0.0004	0.9997
K2PIG	14	-3232.9982	6563.2406	0.0003	1.0000
SYMG	17	-3228.8786	6575.8395	0.0000	1.0000
SYMI	17	-3229.8539	6577.7900	0.0000	1.0000
SYMIG	18	-3228.8375	6582.7033	0.0000	1.0000
HKY	15	-3254.8167	6613.8236	0.0000	1.0000
GTR	19	-3243.5183	6619.0108	0.0000	1.0000
K2P	12	-3271.7418	6626.8358	0.0000	1.0000
SYM	16	-3265.9383	6643.0129	0.0000	1.0000
F81G	15	-3280.4386	6665.0674	0.0000	1.0000
F81I	15	-3280.6932	6665.5765	0.0000	1.0000
JC69G	12	-3292.8098	6668.9718	0.0000	1.0000
JC69I	12	-3293.0730	6669.4981	0.0000	1.0000
F81IG	16	-3280.4692	6672.0746	0.0000	1.0000
JC69IG	13	-3292.8292	6675.9566	0.0000	1.0000
F81	14	-3314.2239	6725.6920	0.0000	1.0000
JC69	11	-3327.6964	6731.7991	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/chr5_3386.phylip = 1 - 1039;
 Partition Dummy = 1:591-tomodel-phylip/chr5_3386.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/chr5_3386.phylip = 1 - 1039;
 Partition Dummy = 1:591-tomodel-phylip/chr5_3386.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 (HKYG)]

BEGIN MRBAYES;
 Charset 591-tomodel-phylip/chr5_3386.phylip = 1 - 1039;
 Partition Dummy = 1:591-tomodel-phylip/chr5_3386.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

