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

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

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


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

Model	df	lnL		AIC		wAIC	cwAIC
HKYG	16	-3300.9911	6633.9822	0.4950	0.4950
GTRG	20	-3297.9810	6635.9620	0.1840	0.6790
HKYIG	17	-3301.0445	6636.0891	0.1726	0.8516
HKYI	16	-3302.9306	6637.8612	0.0712	0.9228
GTRIG	21	-3297.9811	6637.9622	0.0677	0.9904
GTRI	20	-3300.9389	6641.8778	0.0096	1.0000
K2PG	13	-3318.1068	6662.2136	0.0000	1.0000
K2PIG	14	-3318.1514	6664.3028	0.0000	1.0000
K2PI	13	-3319.7415	6665.4831	0.0000	1.0000
SYMG	17	-3316.2486	6666.4972	0.0000	1.0000
SYMIG	18	-3316.2626	6668.5253	0.0000	1.0000
SYMI	17	-3318.0813	6670.1625	0.0000	1.0000
HKY	15	-3358.9753	6747.9507	0.0000	1.0000
GTR	19	-3357.0488	6752.0977	0.0000	1.0000
F81G	15	-3364.0771	6758.1543	0.0000	1.0000
F81I	15	-3364.7497	6759.4994	0.0000	1.0000
F81IG	16	-3364.0006	6760.0011	0.0000	1.0000
K2P	12	-3375.3519	6774.7037	0.0000	1.0000
JC69G	12	-3377.2596	6778.5192	0.0000	1.0000
SYM	16	-3373.5979	6779.1959	0.0000	1.0000
JC69I	12	-3377.8396	6779.6792	0.0000	1.0000
JC69IG	13	-3377.1541	6780.3081	0.0000	1.0000
F81	14	-3419.2138	6866.4276	0.0000	1.0000
JC69	11	-3432.1494	6886.2987	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	16	-3300.9911	6711.1199	0.8501	0.8501
HKYI	16	-3302.9306	6714.9989	0.1222	0.9723
HKYIG	17	-3301.0445	6718.0479	0.0266	0.9989
K2PG	13	-3318.1068	6724.8880	0.0009	0.9998
K2PI	13	-3319.7415	6728.1575	0.0002	1.0000
K2PIG	14	-3318.1514	6731.7983	0.0000	1.0000
GTRG	20	-3297.9810	6732.3841	0.0000	1.0000
GTRI	20	-3300.9389	6738.2999	0.0000	1.0000
GTRIG	21	-3297.9811	6739.2054	0.0000	1.0000
SYMG	17	-3316.2486	6748.4561	0.0000	1.0000
SYMI	17	-3318.0813	6752.1214	0.0000	1.0000
SYMIG	18	-3316.2626	6755.3052	0.0000	1.0000
HKY	15	-3358.9753	6820.2673	0.0000	1.0000
F81G	15	-3364.0771	6830.4709	0.0000	1.0000
F81I	15	-3364.7497	6831.8160	0.0000	1.0000
K2P	12	-3375.3519	6832.5570	0.0000	1.0000
JC69G	12	-3377.2596	6836.3725	0.0000	1.0000
F81IG	16	-3364.0006	6837.1388	0.0000	1.0000
JC69I	12	-3377.8396	6837.5325	0.0000	1.0000
JC69IG	13	-3377.1541	6842.9825	0.0000	1.0000
GTR	19	-3357.0488	6843.6987	0.0000	1.0000
SYM	16	-3373.5979	6856.3336	0.0000	1.0000
F81	14	-3419.2138	6933.9231	0.0000	1.0000
JC69	11	-3432.1494	6939.3309	0.0000	1.0000



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

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

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


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

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


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

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

