
Output from mraic.pl version 1.4.4 by Johan Nylander
-------------------------------------------------------------
00:45 02/17/2011

Input data from file "183-phylip/chr2_23281.phylip" (29 taxa, 834 characters)

Minimum AIC  model: GTRIG
Minimum AICc model: GTRIG
Minimum BIC  model: GTRG


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
GTRIG	65	-6592.9624	13327.0967	0.8707	0.8707
GTRG	64	-6596.0460	13330.9112	0.1293	1.0000
SYMIG	62	-6613.5752	13361.2826	0.0000	1.0000
SYMG	61	-6615.9481	13363.6941	0.0000	1.0000
HKYG	60	-6622.6335	13374.7365	0.0000	1.0000
HKYIG	61	-6621.5815	13374.9610	0.0000	1.0000
K2PG	57	-6626.7819	13376.0844	0.0000	1.0000
K2PIG	58	-6626.2073	13377.2456	0.0000	1.0000
GTRI	64	-6720.3468	13579.5129	0.0000	1.0000
SYMI	61	-6734.0040	13599.8059	0.0000	1.0000
HKYI	60	-6737.1617	13603.7930	0.0000	1.0000
K2PI	57	-6742.2204	13606.9613	0.0000	1.0000
JC69IG	57	-6784.4631	13691.4467	0.0000	1.0000
JC69G	56	-6786.9026	13694.0213	0.0000	1.0000
F81IG	60	-6783.4727	13696.4150	0.0000	1.0000
F81G	59	-6785.9542	13699.0556	0.0000	1.0000
JC69I	56	-6895.9706	13912.1574	0.0000	1.0000
F81I	59	-6894.6842	13916.5156	0.0000	1.0000
GTR	63	-7021.5829	14179.6386	0.0000	1.0000
SYM	60	-7040.2094	14209.8884	0.0000	1.0000
HKY	59	-7050.1794	14227.5061	0.0000	1.0000
K2P	56	-7056.5520	14233.3202	0.0000	1.0000
JC69	55	-7200.0947	14518.1071	0.0000	1.0000
F81	58	-7198.5189	14521.8689	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRIG	65	-6592.9624	13315.9248	0.8893	0.8893
GTRG	64	-6596.0460	13320.0920	0.1107	1.0000
SYMIG	62	-6613.5752	13351.1503	0.0000	1.0000
SYMG	61	-6615.9481	13353.8962	0.0000	1.0000
HKYIG	61	-6621.5815	13365.1630	0.0000	1.0000
HKYG	60	-6622.6335	13365.2669	0.0000	1.0000
K2PG	57	-6626.7819	13367.5638	0.0000	1.0000
K2PIG	58	-6626.2073	13368.4146	0.0000	1.0000
GTRI	64	-6720.3468	13568.6937	0.0000	1.0000
SYMI	61	-6734.0040	13590.0080	0.0000	1.0000
HKYI	60	-6737.1617	13594.3234	0.0000	1.0000
K2PI	57	-6742.2204	13598.4407	0.0000	1.0000
JC69IG	57	-6784.4631	13682.9261	0.0000	1.0000
JC69G	56	-6786.9026	13685.8051	0.0000	1.0000
F81IG	60	-6783.4727	13686.9454	0.0000	1.0000
F81G	59	-6785.9542	13689.9083	0.0000	1.0000
JC69I	56	-6895.9706	13903.9411	0.0000	1.0000
F81I	59	-6894.6842	13907.3683	0.0000	1.0000
GTR	63	-7021.5829	14169.1659	0.0000	1.0000
SYM	60	-7040.2094	14200.4188	0.0000	1.0000
HKY	59	-7050.1794	14218.3589	0.0000	1.0000
K2P	56	-7056.5520	14225.1040	0.0000	1.0000
JC69	55	-7200.0947	14510.1894	0.0000	1.0000
F81	58	-7198.5189	14513.0379	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
GTRG	64	-6596.0460	13622.5709	0.5691	0.5691
GTRIG	65	-6592.9624	13623.1300	0.4304	0.9995
K2PG	57	-6626.7819	13636.9591	0.0004	0.9999
SYMG	61	-6615.9481	13642.1965	0.0000	1.0000
K2PIG	58	-6626.2073	13642.5362	0.0000	1.0000
SYMIG	62	-6613.5752	13644.1768	0.0000	1.0000
HKYG	60	-6622.6335	13648.8409	0.0000	1.0000
HKYIG	61	-6621.5815	13653.4633	0.0000	1.0000
K2PI	57	-6742.2204	13867.8360	0.0000	1.0000
GTRI	64	-6720.3468	13871.1726	0.0000	1.0000
HKYI	60	-6737.1617	13877.8974	0.0000	1.0000
SYMI	61	-6734.0040	13878.3082	0.0000	1.0000
JC69G	56	-6786.9026	13950.4742	0.0000	1.0000
JC69IG	57	-6784.4631	13952.3214	0.0000	1.0000
F81G	59	-6785.9542	13968.7561	0.0000	1.0000
F81IG	60	-6783.4727	13970.5194	0.0000	1.0000
JC69I	56	-6895.9706	14168.6102	0.0000	1.0000
F81I	59	-6894.6842	14186.2161	0.0000	1.0000
GTR	63	-7021.5829	14466.9186	0.0000	1.0000
SYM	60	-7040.2094	14483.9928	0.0000	1.0000
K2P	56	-7056.5520	14489.7731	0.0000	1.0000
HKY	59	-7050.1794	14497.2066	0.0000	1.0000
JC69	55	-7200.0947	14770.1322	0.0000	1.0000
F81	58	-7198.5189	14787.1594	0.0000	1.0000



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

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

BEGIN MRBAYES;
 Charset 183-phylip/chr2_23281.phylip = 1 - 834;
 Partition Dummy = 1:183-phylip/chr2_23281.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=6 rates=invgamma;
 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) pinvarpr=Uniform(0.0,1.0);
END;


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

BEGIN MRBAYES;
 Charset 183-phylip/chr2_23281.phylip = 1 - 834;
 Partition Dummy = 1:183-phylip/chr2_23281.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=6 rates=invgamma;
 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) pinvarpr=Uniform(0.0,1.0);
END;


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

BEGIN MRBAYES;
 Charset 183-phylip/chr2_23281.phylip = 1 - 834;
 Partition Dummy = 1:183-phylip/chr2_23281.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;


-------------------------------------------------------------
End of Output

