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

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

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


Note: number of parameters compared to sample size
is low (Nchar/Nparams < 40). Consider the use of AICc.


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

Model	df	lnL		AICc		wAICc	cwAICc
GTRG	20	-2668.7169	5378.4985	0.4131	0.4131
GTRI	20	-2668.9655	5378.9956	0.3222	0.7352
GTRIG	21	-2668.1982	5379.5690	0.2419	0.9771
HKYG	16	-2676.4172	5385.5203	0.0123	0.9894
HKYIG	17	-2676.2157	5387.2041	0.0053	0.9947
HKYI	16	-2677.2641	5387.2143	0.0053	1.0000
K2PG	13	-2705.1193	5436.6960	0.0000	1.0000
K2PIG	14	-2705.0160	5438.5604	0.0000	1.0000
K2PI	13	-2706.8651	5440.1875	0.0000	1.0000
SYMG	17	-2703.2987	5441.3702	0.0000	1.0000
SYMIG	18	-2703.1920	5443.2487	0.0000	1.0000
SYMI	17	-2704.8597	5444.4921	0.0000	1.0000
GTR	19	-2705.6140	5450.1901	0.0000	1.0000
HKY	15	-2713.3423	5457.2892	0.0000	1.0000
K2P	12	-2744.1713	5512.7341	0.0000	1.0000
SYM	16	-2741.9394	5516.5647	0.0000	1.0000
F81G	15	-2747.4420	5525.4885	0.0000	1.0000
F81I	15	-2748.1012	5526.8069	0.0000	1.0000
F81IG	16	-2747.2913	5527.2686	0.0000	1.0000
JC69G	12	-2770.4666	5565.3247	0.0000	1.0000
JC69IG	13	-2770.3998	5567.2569	0.0000	1.0000
JC69I	12	-2771.5815	5567.5545	0.0000	1.0000
F81	14	-2781.0514	5590.6312	0.0000	1.0000
JC69	11	-2805.4117	5633.1542	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRG	20	-2668.7169	5377.4339	0.4091	0.4091
GTRI	20	-2668.9655	5377.9310	0.3191	0.7282
GTRIG	21	-2668.1982	5378.3964	0.2528	0.9810
HKYG	16	-2676.4172	5384.8343	0.0101	0.9911
HKYIG	17	-2676.2157	5386.4314	0.0046	0.9957
HKYI	16	-2677.2641	5386.5283	0.0043	1.0000
K2PG	13	-2705.1193	5436.2387	0.0000	1.0000
K2PIG	14	-2705.0160	5438.0321	0.0000	1.0000
K2PI	13	-2706.8651	5439.7302	0.0000	1.0000
SYMG	17	-2703.2987	5440.5975	0.0000	1.0000
SYMIG	18	-2703.1920	5442.3840	0.0000	1.0000
SYMI	17	-2704.8597	5443.7194	0.0000	1.0000
GTR	19	-2705.6140	5449.2281	0.0000	1.0000
HKY	15	-2713.3423	5456.6847	0.0000	1.0000
K2P	12	-2744.1713	5512.3426	0.0000	1.0000
SYM	16	-2741.9394	5515.8787	0.0000	1.0000
F81G	15	-2747.4420	5524.8840	0.0000	1.0000
F81I	15	-2748.1012	5526.2024	0.0000	1.0000
F81IG	16	-2747.2913	5526.5826	0.0000	1.0000
JC69G	12	-2770.4666	5564.9332	0.0000	1.0000
JC69IG	13	-2770.3998	5566.7997	0.0000	1.0000
JC69I	12	-2771.5815	5567.1630	0.0000	1.0000
F81	14	-2781.0514	5590.1029	0.0000	1.0000
JC69	11	-2805.4117	5632.8233	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	16	-2676.4172	5459.9869	0.6766	0.6766
HKYI	16	-2677.2641	5461.6808	0.2901	0.9667
HKYIG	17	-2676.2157	5466.2810	0.0291	0.9958
GTRG	20	-2668.7169	5471.3745	0.0023	0.9981
GTRI	20	-2668.9655	5471.8716	0.0018	0.9999
GTRIG	21	-2668.1982	5477.0341	0.0001	1.0000
K2PG	13	-2705.1193	5497.3001	0.0000	1.0000
K2PI	13	-2706.8651	5500.7917	0.0000	1.0000
K2PIG	14	-2705.0160	5503.7905	0.0000	1.0000
SYMG	17	-2703.2987	5520.4471	0.0000	1.0000
SYMI	17	-2704.8597	5523.5689	0.0000	1.0000
SYMIG	18	-2703.1920	5526.9306	0.0000	1.0000
HKY	15	-2713.3423	5527.1402	0.0000	1.0000
GTR	19	-2705.6140	5538.4717	0.0000	1.0000
K2P	12	-2744.1713	5568.7070	0.0000	1.0000
SYM	16	-2741.9394	5591.0313	0.0000	1.0000
F81G	15	-2747.4420	5595.3395	0.0000	1.0000
F81I	15	-2748.1012	5596.6579	0.0000	1.0000
F81IG	16	-2747.2913	5601.7351	0.0000	1.0000
JC69G	12	-2770.4666	5621.2976	0.0000	1.0000
JC69I	12	-2771.5815	5623.5274	0.0000	1.0000
JC69IG	13	-2770.3998	5627.8611	0.0000	1.0000
F81	14	-2781.0514	5655.8614	0.0000	1.0000
JC69	11	-2805.4117	5684.4907	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/chr2_22616.phylip = 1 - 810;
 Partition Dummy = 1:591-tomodel-phylip/chr2_22616.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/chr2_22616.phylip = 1 - 810;
 Partition Dummy = 1:591-tomodel-phylip/chr2_22616.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/chr2_22616.phylip = 1 - 810;
 Partition Dummy = 1:591-tomodel-phylip/chr2_22616.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

