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

Input data from file "591-tomodel-phylip/chr11_3549.phylip" (7 taxa, 933 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	-3110.9550	6261.9100	0.7516	0.7516
GTRIG	21	-3111.6164	6265.2329	0.1427	0.8944
GTRI	20	-3112.9253	6265.8505	0.1048	0.9991
HKYG	16	-3122.3262	6276.6524	0.0005	0.9996
HKYI	16	-3123.1554	6278.3109	0.0002	0.9998
HKYIG	17	-3122.3384	6278.6769	0.0002	1.0000
SYMG	17	-3128.1979	6290.3957	0.0000	1.0000
SYMIG	18	-3128.8527	6293.7054	0.0000	1.0000
SYMI	17	-3129.9668	6293.9335	0.0000	1.0000
K2PG	13	-3138.8386	6303.6773	0.0000	1.0000
K2PI	13	-3139.5961	6305.1922	0.0000	1.0000
K2PIG	14	-3138.8572	6305.7144	0.0000	1.0000
GTR	19	-3157.1388	6352.2777	0.0000	1.0000
HKY	15	-3166.7734	6363.5469	0.0000	1.0000
SYM	16	-3171.6246	6375.2491	0.0000	1.0000
K2P	12	-3183.0140	6390.0279	0.0000	1.0000
F81G	15	-3208.9007	6447.8013	0.0000	1.0000
F81I	15	-3209.6197	6449.2394	0.0000	1.0000
F81IG	16	-3208.9957	6449.9914	0.0000	1.0000
JC69G	12	-3224.1675	6472.3350	0.0000	1.0000
JC69I	12	-3224.9255	6473.8511	0.0000	1.0000
JC69IG	13	-3224.2445	6474.4890	0.0000	1.0000
F81	14	-3249.8185	6527.6371	0.0000	1.0000
JC69	11	-3263.7314	6549.4628	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	16	-3122.3262	6354.0669	0.6301	0.6301
HKYI	16	-3123.1554	6355.7253	0.2750	0.9051
GTRG	20	-3110.9550	6358.6781	0.0628	0.9679
HKYIG	17	-3122.3384	6360.9297	0.0204	0.9883
GTRI	20	-3112.9253	6362.6186	0.0088	0.9971
K2PG	13	-3138.8386	6366.5765	0.0012	0.9983
GTRIG	21	-3111.6164	6366.8394	0.0011	0.9993
K2PI	13	-3139.5961	6368.0915	0.0006	0.9999
SYMG	17	-3128.1979	6372.6486	0.0001	1.0000
K2PIG	14	-3138.8572	6373.4521	0.0000	1.0000
SYMI	17	-3129.9668	6376.1864	0.0000	1.0000
SYMIG	18	-3128.8527	6380.7967	0.0000	1.0000
HKY	15	-3166.7734	6436.1230	0.0000	1.0000
GTR	19	-3157.1388	6444.2074	0.0000	1.0000
K2P	12	-3183.0140	6448.0888	0.0000	1.0000
SYM	16	-3171.6246	6452.6636	0.0000	1.0000
F81G	15	-3208.9007	6520.3774	0.0000	1.0000
F81I	15	-3209.6197	6521.8154	0.0000	1.0000
F81IG	16	-3208.9957	6527.4059	0.0000	1.0000
JC69G	12	-3224.1675	6530.3959	0.0000	1.0000
JC69I	12	-3224.9255	6531.9119	0.0000	1.0000
JC69IG	13	-3224.2445	6537.3882	0.0000	1.0000
F81	14	-3249.8185	6595.3747	0.0000	1.0000
JC69	11	-3263.7314	6602.6853	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/chr11_3549.phylip = 1 - 933;
 Partition Dummy = 1:591-tomodel-phylip/chr11_3549.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/chr11_3549.phylip = 1 - 933;
 Partition Dummy = 1:591-tomodel-phylip/chr11_3549.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/chr11_3549.phylip = 1 - 933;
 Partition Dummy = 1:591-tomodel-phylip/chr11_3549.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

