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

Input data from file "591-tomodel-phylip/chr3_6129.phylip" (7 taxa, 674 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	-2082.4069	4206.1001	0.6768	0.6768
GTRIG	21	-2082.4074	4208.2320	0.2331	0.9099
GTRI	20	-2084.4948	4210.2759	0.0839	0.9938
HKYG	16	-2091.7317	4216.2915	0.0041	0.9979
HKYIG	17	-2091.7413	4218.4155	0.0014	0.9993
HKYI	16	-2093.5744	4219.9769	0.0007	1.0000
GTR	19	-2094.5661	4228.2944	0.0000	1.0000
SYMG	17	-2098.7748	4232.4825	0.0000	1.0000
SYMIG	18	-2098.7747	4234.5936	0.0000	1.0000
SYMI	17	-2100.7726	4236.4781	0.0000	1.0000
K2PG	13	-2105.0349	4236.6213	0.0000	1.0000
HKY	15	-2103.0317	4236.7929	0.0000	1.0000
K2PIG	14	-2105.0336	4238.7046	0.0000	1.0000
K2PI	13	-2106.9042	4240.3599	0.0000	1.0000
SYM	16	-2110.5651	4253.9582	0.0000	1.0000
K2P	12	-2116.6078	4257.6876	0.0000	1.0000
F81G	15	-2132.1619	4295.0533	0.0000	1.0000
F81IG	16	-2132.1623	4297.1526	0.0000	1.0000
F81I	15	-2133.8025	4298.3346	0.0000	1.0000
JC69G	12	-2143.0160	4310.5040	0.0000	1.0000
JC69IG	13	-2143.0152	4312.5820	0.0000	1.0000
JC69I	12	-2144.6835	4313.8390	0.0000	1.0000
F81	14	-2143.3308	4315.2989	0.0000	1.0000
JC69	11	-2154.3799	4331.1586	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRG	20	-2082.4069	4204.8138	0.6671	0.6671
GTRIG	21	-2082.4074	4206.8148	0.2453	0.9124
GTRI	20	-2084.4948	4208.9895	0.0827	0.9950
HKYG	16	-2091.7317	4215.4635	0.0032	0.9983
HKYIG	17	-2091.7413	4217.4826	0.0012	0.9995
HKYI	16	-2093.5744	4219.1489	0.0005	1.0000
GTR	19	-2094.5661	4227.1323	0.0000	1.0000
SYMG	17	-2098.7748	4231.5496	0.0000	1.0000
SYMIG	18	-2098.7747	4233.5493	0.0000	1.0000
SYMI	17	-2100.7726	4235.5452	0.0000	1.0000
HKY	15	-2103.0317	4236.0634	0.0000	1.0000
K2PG	13	-2105.0349	4236.0697	0.0000	1.0000
K2PIG	14	-2105.0336	4238.0673	0.0000	1.0000
K2PI	13	-2106.9042	4239.8083	0.0000	1.0000
SYM	16	-2110.5651	4253.1302	0.0000	1.0000
K2P	12	-2116.6078	4257.2156	0.0000	1.0000
F81G	15	-2132.1619	4294.3238	0.0000	1.0000
F81IG	16	-2132.1623	4296.3246	0.0000	1.0000
F81I	15	-2133.8025	4297.6051	0.0000	1.0000
JC69G	12	-2143.0160	4310.0320	0.0000	1.0000
JC69IG	13	-2143.0152	4312.0304	0.0000	1.0000
JC69I	12	-2144.6835	4313.3670	0.0000	1.0000
F81	14	-2143.3308	4314.6616	0.0000	1.0000
JC69	11	-2154.3799	4330.7598	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	16	-2091.7317	4287.6752	0.7934	0.7934
HKYI	16	-2093.5744	4291.3606	0.1257	0.9190
HKYIG	17	-2091.7413	4294.2075	0.0303	0.9493
K2PG	13	-2105.0349	4294.7417	0.0232	0.9725
GTRG	20	-2082.4069	4295.0784	0.0196	0.9921
K2PI	13	-2106.9042	4298.4803	0.0036	0.9956
GTRI	20	-2084.4948	4299.2541	0.0024	0.9981
K2PIG	14	-2105.0336	4301.2525	0.0009	0.9990
GTRIG	21	-2082.4074	4301.5927	0.0008	0.9997
HKY	15	-2103.0317	4303.7619	0.0003	1.0000
SYMG	17	-2098.7748	4308.2745	0.0000	1.0000
K2P	12	-2116.6078	4311.3743	0.0000	1.0000
SYMI	17	-2100.7726	4312.2701	0.0000	1.0000
GTR	19	-2094.5661	4312.8837	0.0000	1.0000
SYMIG	18	-2098.7747	4314.7875	0.0000	1.0000
SYM	16	-2110.5651	4325.3418	0.0000	1.0000
F81G	15	-2132.1619	4362.0223	0.0000	1.0000
JC69G	12	-2143.0160	4364.1908	0.0000	1.0000
F81I	15	-2133.8025	4365.3036	0.0000	1.0000
JC69I	12	-2144.6835	4367.5258	0.0000	1.0000
F81IG	16	-2132.1623	4368.5362	0.0000	1.0000
JC69IG	13	-2143.0152	4370.7024	0.0000	1.0000
F81	14	-2143.3308	4377.8468	0.0000	1.0000
JC69	11	-2154.3799	4380.4053	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/chr3_6129.phylip = 1 - 674;
 Partition Dummy = 1:591-tomodel-phylip/chr3_6129.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/chr3_6129.phylip = 1 - 674;
 Partition Dummy = 1:591-tomodel-phylip/chr3_6129.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/chr3_6129.phylip = 1 - 674;
 Partition Dummy = 1:591-tomodel-phylip/chr3_6129.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

