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

Input data from file "591-tomodel-phylip/chr5_17310.phylip" (7 taxa, 728 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	-2133.9876	4309.1633	0.5364	0.5364
GTRIG	21	-2133.1597	4309.6283	0.4251	0.9616
GTRI	20	-2136.6544	4314.4969	0.0373	0.9988
HKYG	16	-2144.6388	4322.0428	0.0009	0.9997
HKYIG	17	-2144.7902	4324.4424	0.0003	1.0000
HKYI	16	-2147.5467	4327.8585	0.0000	1.0000
SYMG	17	-2165.4447	4365.7514	0.0000	1.0000
SYMIG	18	-2164.6007	4366.1662	0.0000	1.0000
SYMI	17	-2167.9609	4370.7838	0.0000	1.0000
K2PG	13	-2173.7806	4374.0710	0.0000	1.0000
K2PIG	14	-2173.4410	4375.4710	0.0000	1.0000
K2PI	13	-2176.3928	4379.2955	0.0000	1.0000
GTR	19	-2175.6990	4390.4714	0.0000	1.0000
HKY	15	-2186.2113	4403.0968	0.0000	1.0000
F81IG	16	-2209.7486	4452.2623	0.0000	1.0000
F81G	15	-2210.8730	4452.4202	0.0000	1.0000
F81I	15	-2212.0027	4454.6796	0.0000	1.0000
SYM	16	-2211.2958	4455.3567	0.0000	1.0000
K2P	12	-2218.7655	4461.9674	0.0000	1.0000
JC69IG	13	-2234.3913	4495.2924	0.0000	1.0000
JC69G	12	-2235.7360	4495.9084	0.0000	1.0000
JC69I	12	-2236.6031	4497.6426	0.0000	1.0000
F81	14	-2250.0071	4528.6032	0.0000	1.0000
JC69	11	-2276.2923	4574.9534	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRG	20	-2133.9876	4307.9751	0.5227	0.5227
GTRIG	21	-2133.1597	4308.3195	0.4400	0.9628
GTRI	20	-2136.6544	4313.3087	0.0363	0.9991
HKYG	16	-2144.6388	4321.2776	0.0007	0.9997
HKYIG	17	-2144.7902	4323.5805	0.0002	1.0000
HKYI	16	-2147.5467	4327.0933	0.0000	1.0000
SYMG	17	-2165.4447	4364.8894	0.0000	1.0000
SYMIG	18	-2164.6007	4365.2015	0.0000	1.0000
SYMI	17	-2167.9609	4369.9218	0.0000	1.0000
K2PG	13	-2173.7806	4373.5612	0.0000	1.0000
K2PIG	14	-2173.4410	4374.8820	0.0000	1.0000
K2PI	13	-2176.3928	4378.7857	0.0000	1.0000
GTR	19	-2175.6990	4389.3980	0.0000	1.0000
HKY	15	-2186.2113	4402.4226	0.0000	1.0000
F81IG	16	-2209.7486	4451.4972	0.0000	1.0000
F81G	15	-2210.8730	4451.7460	0.0000	1.0000
F81I	15	-2212.0027	4454.0055	0.0000	1.0000
SYM	16	-2211.2958	4454.5916	0.0000	1.0000
K2P	12	-2218.7655	4461.5311	0.0000	1.0000
JC69IG	13	-2234.3913	4494.7826	0.0000	1.0000
JC69G	12	-2235.7360	4495.4720	0.0000	1.0000
JC69I	12	-2236.6031	4497.2063	0.0000	1.0000
F81	14	-2250.0071	4528.0141	0.0000	1.0000
JC69	11	-2276.2923	4574.5847	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	16	-2144.6388	4394.7225	0.8486	0.8486
GTRG	20	-2133.9876	4399.7812	0.0676	0.9162
HKYI	16	-2147.5467	4400.5382	0.0463	0.9625
HKYIG	17	-2144.7902	4401.6156	0.0270	0.9896
GTRIG	21	-2133.1597	4404.7158	0.0057	0.9953
GTRI	20	-2136.6544	4405.1148	0.0047	1.0000
K2PG	13	-2173.7806	4433.2351	0.0000	1.0000
K2PI	13	-2176.3928	4438.4596	0.0000	1.0000
K2PIG	14	-2173.4410	4439.1462	0.0000	1.0000
SYMG	17	-2165.4447	4442.9245	0.0000	1.0000
SYMIG	18	-2164.6007	4447.8269	0.0000	1.0000
SYMI	17	-2167.9609	4447.9569	0.0000	1.0000
HKY	15	-2186.2113	4471.2772	0.0000	1.0000
GTR	19	-2175.6990	4476.6137	0.0000	1.0000
K2P	12	-2218.7655	4516.6147	0.0000	1.0000
F81G	15	-2210.8730	4520.6005	0.0000	1.0000
F81I	15	-2212.0027	4522.8600	0.0000	1.0000
F81IG	16	-2209.7486	4524.9420	0.0000	1.0000
SYM	16	-2211.2958	4528.0364	0.0000	1.0000
JC69G	12	-2235.7360	4550.5557	0.0000	1.0000
JC69I	12	-2236.6031	4552.2899	0.0000	1.0000
JC69IG	13	-2234.3913	4554.4566	0.0000	1.0000
F81	14	-2250.0071	4592.2784	0.0000	1.0000
JC69	11	-2276.2923	4625.0780	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/chr5_17310.phylip = 1 - 728;
 Partition Dummy = 1:591-tomodel-phylip/chr5_17310.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/chr5_17310.phylip = 1 - 728;
 Partition Dummy = 1:591-tomodel-phylip/chr5_17310.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/chr5_17310.phylip = 1 - 728;
 Partition Dummy = 1:591-tomodel-phylip/chr5_17310.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

