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

Input data from file "683-tomodel-phylip/chr3_209.phylip" (7 taxa, 1145 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	-3313.4135	6658.8270	0.5795	0.5795
HKYIG	17	-3313.3745	6660.7491	0.2216	0.8011
HKYI	16	-3314.7803	6661.5605	0.1477	0.9489
GTRG	20	-3312.3273	6664.6547	0.0314	0.9803
GTRIG	21	-3312.2867	6666.5735	0.0120	0.9924
GTRI	20	-3313.7426	6667.4852	0.0076	1.0000
K2PG	13	-3332.7162	6691.4325	0.0000	1.0000
K2PIG	14	-3332.7010	6693.4020	0.0000	1.0000
K2PI	13	-3333.8745	6693.7490	0.0000	1.0000
SYMG	17	-3331.8182	6697.6364	0.0000	1.0000
SYMIG	18	-3331.8117	6699.6234	0.0000	1.0000
SYMI	17	-3332.8873	6699.7746	0.0000	1.0000
HKY	15	-3344.8659	6719.7317	0.0000	1.0000
GTR	19	-3343.6966	6725.3932	0.0000	1.0000
K2P	12	-3362.9346	6749.8692	0.0000	1.0000
SYM	16	-3361.8330	6755.6660	0.0000	1.0000
F81G	15	-3406.2964	6842.5929	0.0000	1.0000
F81IG	16	-3406.2921	6844.5842	0.0000	1.0000
F81I	15	-3407.2991	6844.5982	0.0000	1.0000
JC69G	12	-3421.2108	6866.4216	0.0000	1.0000
JC69I	12	-3422.1415	6868.2829	0.0000	1.0000
JC69IG	13	-3421.2123	6868.4247	0.0000	1.0000
F81	14	-3434.9650	6897.9299	0.0000	1.0000
JC69	11	-3449.5111	6921.0221	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	16	-3313.4135	6739.5175	0.7777	0.7777
HKYI	16	-3314.7803	6742.2511	0.1983	0.9759
HKYIG	17	-3313.3745	6746.4828	0.0239	0.9998
K2PG	13	-3332.7162	6756.9936	0.0001	1.0000
K2PI	13	-3333.8745	6759.3101	0.0000	1.0000
K2PIG	14	-3332.7010	6764.0062	0.0000	1.0000
GTRG	20	-3312.3273	6765.5179	0.0000	1.0000
GTRI	20	-3313.7426	6768.3484	0.0000	1.0000
GTRIG	21	-3312.2867	6772.4798	0.0000	1.0000
SYMG	17	-3331.8182	6783.3701	0.0000	1.0000
SYMI	17	-3332.8873	6785.5083	0.0000	1.0000
SYMIG	18	-3331.8117	6790.4003	0.0000	1.0000
HKY	15	-3344.8659	6795.3791	0.0000	1.0000
K2P	12	-3362.9346	6810.3871	0.0000	1.0000
GTR	19	-3343.6966	6821.2132	0.0000	1.0000
SYM	16	-3361.8330	6836.3566	0.0000	1.0000
F81G	15	-3406.2964	6918.2403	0.0000	1.0000
F81I	15	-3407.2991	6920.2456	0.0000	1.0000
F81IG	16	-3406.2921	6925.2748	0.0000	1.0000
JC69G	12	-3421.2108	6926.9395	0.0000	1.0000
JC69I	12	-3422.1415	6928.8008	0.0000	1.0000
JC69IG	13	-3421.2123	6933.9858	0.0000	1.0000
F81	14	-3434.9650	6968.5342	0.0000	1.0000
JC69	11	-3449.5111	6976.4969	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 683-tomodel-phylip/chr3_209.phylip = 1 - 1145;
 Partition Dummy = 1:683-tomodel-phylip/chr3_209.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 683-tomodel-phylip/chr3_209.phylip = 1 - 1145;
 Partition Dummy = 1:683-tomodel-phylip/chr3_209.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 683-tomodel-phylip/chr3_209.phylip = 1 - 1145;
 Partition Dummy = 1:683-tomodel-phylip/chr3_209.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

