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

Input data from file "591-tomodel-phylip/chr1_5281.phylip" (7 taxa, 1148 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	-3368.2685	6768.5369	0.4016	0.4016
GTRG	20	-3364.9390	6769.8781	0.2054	0.6070
HKYIG	17	-3368.1173	6770.2346	0.1718	0.7788
HKYI	16	-3369.7368	6771.4735	0.0925	0.8713
GTRIG	21	-3364.8271	6771.6542	0.0845	0.9558
GTRI	20	-3366.4757	6772.9514	0.0442	1.0000
K2PG	13	-3406.0747	6838.1493	0.0000	1.0000
K2PIG	14	-3405.8605	6839.7210	0.0000	1.0000
K2PI	13	-3407.5311	6841.0622	0.0000	1.0000
SYMG	17	-3405.8766	6845.7533	0.0000	1.0000
SYMIG	18	-3405.6757	6847.3514	0.0000	1.0000
SYMI	17	-3407.2602	6848.5203	0.0000	1.0000
HKY	15	-3441.8718	6913.7436	0.0000	1.0000
GTR	19	-3438.6164	6915.2329	0.0000	1.0000
F81G	15	-3455.5111	6941.0222	0.0000	1.0000
F81IG	16	-3455.4504	6942.9008	0.0000	1.0000
F81I	15	-3456.5992	6943.1985	0.0000	1.0000
K2P	12	-3478.1603	6980.3206	0.0000	1.0000
SYM	16	-3477.4755	6986.9509	0.0000	1.0000
JC69G	12	-3485.3126	6994.6252	0.0000	1.0000
JC69IG	13	-3485.2511	6996.5021	0.0000	1.0000
JC69I	12	-3486.4040	6996.8080	0.0000	1.0000
F81	14	-3521.6565	7071.3131	0.0000	1.0000
JC69	11	-3551.0089	7124.0179	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	16	-3368.2685	6849.2694	0.7907	0.7907
HKYI	16	-3369.7368	6852.2060	0.1821	0.9728
HKYIG	17	-3368.1173	6856.0128	0.0271	1.0000
GTRG	20	-3364.9390	6870.7936	0.0000	1.0000
GTRI	20	-3366.4757	6873.8669	0.0000	1.0000
GTRIG	21	-3364.8271	6877.6155	0.0000	1.0000
K2PG	13	-3406.0747	6903.7444	0.0000	1.0000
K2PI	13	-3407.5311	6906.6573	0.0000	1.0000
K2PIG	14	-3405.8605	6910.3619	0.0000	1.0000
SYMG	17	-3405.8766	6931.5315	0.0000	1.0000
SYMI	17	-3407.2602	6934.2985	0.0000	1.0000
SYMIG	18	-3405.6757	6938.1753	0.0000	1.0000
HKY	15	-3441.8718	6989.4302	0.0000	1.0000
GTR	19	-3438.6164	7011.1027	0.0000	1.0000
F81G	15	-3455.5111	7016.7088	0.0000	1.0000
F81I	15	-3456.5992	7018.8851	0.0000	1.0000
F81IG	16	-3455.4504	7023.6332	0.0000	1.0000
K2P	12	-3478.1603	7040.8699	0.0000	1.0000
JC69G	12	-3485.3126	7055.1745	0.0000	1.0000
JC69I	12	-3486.4040	7057.3573	0.0000	1.0000
JC69IG	13	-3485.2511	7062.0972	0.0000	1.0000
SYM	16	-3477.4755	7067.6833	0.0000	1.0000
F81	14	-3521.6565	7141.9539	0.0000	1.0000
JC69	11	-3551.0089	7179.5214	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/chr1_5281.phylip = 1 - 1148;
 Partition Dummy = 1:591-tomodel-phylip/chr1_5281.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/chr1_5281.phylip = 1 - 1148;
 Partition Dummy = 1:591-tomodel-phylip/chr1_5281.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/chr1_5281.phylip = 1 - 1148;
 Partition Dummy = 1:591-tomodel-phylip/chr1_5281.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

