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

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

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


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

Model	df	lnL		AIC		wAIC	cwAIC
GTRG	20	-3018.4178	6076.8357	0.6933	0.6933
GTRIG	21	-3018.3697	6078.7393	0.2677	0.9610
GTRI	20	-3021.7812	6083.5625	0.0240	0.9850
HKYG	16	-3027.0324	6086.0648	0.0069	0.9919
SYMG	17	-3026.7082	6087.4164	0.0035	0.9954
HKYIG	17	-3027.0098	6088.0195	0.0026	0.9979
SYMIG	18	-3026.6722	6089.3445	0.0013	0.9993
HKYI	16	-3030.0009	6092.0017	0.0004	0.9996
K2PG	13	-3033.7697	6093.5394	0.0002	0.9998
SYMI	17	-3029.9873	6093.9746	0.0001	0.9999
K2PIG	14	-3033.7454	6095.4908	0.0001	1.0000
K2PI	13	-3036.6846	6099.3691	0.0000	1.0000
GTR	19	-3057.9013	6153.8026	0.0000	1.0000
HKY	15	-3064.8762	6159.7525	0.0000	1.0000
SYM	16	-3065.0257	6162.0513	0.0000	1.0000
K2P	12	-3071.2974	6166.5948	0.0000	1.0000
F81G	15	-3081.9692	6193.9384	0.0000	1.0000
F81IG	16	-3081.9716	6195.9432	0.0000	1.0000
F81I	15	-3083.9497	6197.8994	0.0000	1.0000
JC69G	12	-3087.1523	6198.3045	0.0000	1.0000
JC69IG	13	-3087.1586	6200.3172	0.0000	1.0000
JC69I	12	-3089.1426	6202.2853	0.0000	1.0000
F81	14	-3116.3302	6260.6605	0.0000	1.0000
JC69	11	-3121.4838	6264.9676	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
K2PG	13	-3033.7697	6156.0999	0.8902	0.8902
K2PI	13	-3036.6846	6161.9296	0.0483	0.9385
K2PIG	14	-3033.7454	6162.8637	0.0303	0.9687
HKYG	16	-3027.0324	6163.0623	0.0274	0.9961
HKYI	16	-3030.0009	6168.9992	0.0014	0.9975
SYMG	17	-3026.7082	6169.2262	0.0013	0.9988
HKYIG	17	-3027.0098	6169.8294	0.0009	0.9997
GTRG	20	-3018.4178	6173.0826	0.0002	0.9999
SYMI	17	-3029.9873	6175.7844	0.0000	0.9999
SYMIG	18	-3026.6722	6175.9667	0.0000	1.0000
GTRIG	21	-3018.3697	6179.7986	0.0000	1.0000
GTRI	20	-3021.7812	6179.8094	0.0000	1.0000
K2P	12	-3071.2974	6224.3430	0.0000	1.0000
HKY	15	-3064.8762	6231.9376	0.0000	1.0000
SYM	16	-3065.0257	6239.0488	0.0000	1.0000
GTR	19	-3057.9013	6245.2371	0.0000	1.0000
JC69G	12	-3087.1523	6256.0527	0.0000	1.0000
JC69I	12	-3089.1426	6260.0334	0.0000	1.0000
JC69IG	13	-3087.1586	6262.8777	0.0000	1.0000
F81G	15	-3081.9692	6266.1235	0.0000	1.0000
F81I	15	-3083.9497	6270.0845	0.0000	1.0000
F81IG	16	-3081.9716	6272.9407	0.0000	1.0000
JC69	11	-3121.4838	6317.9034	0.0000	1.0000
F81	14	-3116.3302	6328.0333	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/chr3_3073.phylip = 1 - 909;
 Partition Dummy = 1:591-tomodel-phylip/chr3_3073.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/chr3_3073.phylip = 1 - 909;
 Partition Dummy = 1:591-tomodel-phylip/chr3_3073.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 (K2PG)]

BEGIN MRBAYES;
 Charset 591-tomodel-phylip/chr3_3073.phylip = 1 - 909;
 Partition Dummy = 1:591-tomodel-phylip/chr3_3073.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=2 rates=gamma;
 Prset applyto=(1) statefreqpr=Fixed(Equal) tratiopr=Beta(1.0,1.0) shapepr=Uniform(0.1,50.0);
END;


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

