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

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

Minimum AIC  model: HKYI
Minimum AICc model: HKYI
Minimum BIC  model: HKYI


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
HKYI	16	-1587.8608	3208.4667	0.5058	0.5058
HKYG	16	-1588.3792	3209.5035	0.3012	0.8070
HKYIG	17	-1587.8607	3210.5609	0.1775	0.9845
GTRI	20	-1587.8031	3216.7633	0.0080	0.9925
GTRG	20	-1588.3177	3217.7924	0.0048	0.9972
GTRIG	21	-1587.8030	3218.8805	0.0028	1.0000
HKY	15	-1604.6965	3240.0496	0.0000	1.0000
GTR	19	-1604.6311	3248.3077	0.0000	1.0000
K2PI	13	-1617.1709	3260.8383	0.0000	1.0000
K2PG	13	-1617.2137	3260.9240	0.0000	1.0000
K2PIG	14	-1617.1709	3262.9156	0.0000	1.0000
SYMI	17	-1616.0287	3266.8970	0.0000	1.0000
SYMG	17	-1616.0769	3266.9933	0.0000	1.0000
SYMIG	18	-1616.0287	3268.9969	0.0000	1.0000
F81I	15	-1621.2022	3273.0611	0.0000	1.0000
F81G	15	-1621.2033	3273.0633	0.0000	1.0000
F81IG	16	-1621.2033	3275.1518	0.0000	1.0000
K2P	12	-1634.6750	3293.7751	0.0000	1.0000
SYM	16	-1633.7647	3300.2745	0.0000	1.0000
F81	14	-1636.2889	3301.1516	0.0000	1.0000
JC69G	12	-1648.3964	3321.2179	0.0000	1.0000
JC69I	12	-1648.3964	3321.2180	0.0000	1.0000
JC69IG	13	-1647.5362	3321.5691	0.0000	1.0000
JC69	11	-1663.7502	3349.8595	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
HKYI	16	-1587.8608	3207.7215	0.4996	0.4996
HKYG	16	-1588.3792	3208.7583	0.2975	0.7971
HKYIG	17	-1587.8607	3209.7214	0.1838	0.9809
GTRI	20	-1587.8031	3215.6063	0.0097	0.9906
GTRG	20	-1588.3177	3216.6353	0.0058	0.9964
GTRIG	21	-1587.8030	3217.6060	0.0036	1.0000
HKY	15	-1604.6965	3239.3930	0.0000	1.0000
GTR	19	-1604.6311	3247.2623	0.0000	1.0000
K2PI	13	-1617.1709	3260.3417	0.0000	1.0000
K2PG	13	-1617.2137	3260.4274	0.0000	1.0000
K2PIG	14	-1617.1709	3262.3418	0.0000	1.0000
SYMI	17	-1616.0287	3266.0575	0.0000	1.0000
SYMG	17	-1616.0769	3266.1538	0.0000	1.0000
SYMIG	18	-1616.0287	3268.0573	0.0000	1.0000
F81I	15	-1621.2022	3272.4044	0.0000	1.0000
F81G	15	-1621.2033	3272.4066	0.0000	1.0000
F81IG	16	-1621.2033	3274.4066	0.0000	1.0000
K2P	12	-1634.6750	3293.3500	0.0000	1.0000
SYM	16	-1633.7647	3299.5293	0.0000	1.0000
F81	14	-1636.2889	3300.5779	0.0000	1.0000
JC69G	12	-1648.3964	3320.7928	0.0000	1.0000
JC69I	12	-1648.3964	3320.7929	0.0000	1.0000
JC69IG	13	-1647.5362	3321.0725	0.0000	1.0000
JC69	11	-1663.7502	3349.5003	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYI	16	-1587.8608	3281.5786	0.6127	0.6127
HKYG	16	-1588.3792	3282.6154	0.3649	0.9776
HKYIG	17	-1587.8607	3288.1945	0.0224	1.0000
GTRI	20	-1587.8031	3307.9276	0.0000	1.0000
HKY	15	-1604.6965	3308.6340	0.0000	1.0000
GTRG	20	-1588.3177	3308.9566	0.0000	1.0000
GTRIG	21	-1587.8030	3314.5433	0.0000	1.0000
K2PI	13	-1617.1709	3320.3506	0.0000	1.0000
K2PG	13	-1617.2137	3320.4363	0.0000	1.0000
K2PIG	14	-1617.1709	3326.9667	0.0000	1.0000
GTR	19	-1604.6311	3334.9675	0.0000	1.0000
F81I	15	-1621.2022	3341.6454	0.0000	1.0000
F81G	15	-1621.2033	3341.6476	0.0000	1.0000
SYMI	17	-1616.0287	3344.5306	0.0000	1.0000
SYMG	17	-1616.0769	3344.6269	0.0000	1.0000
F81IG	16	-1621.2033	3348.2637	0.0000	1.0000
K2P	12	-1634.6750	3348.7428	0.0000	1.0000
SYMIG	18	-1616.0287	3351.1465	0.0000	1.0000
F81	14	-1636.2889	3365.2028	0.0000	1.0000
SYM	16	-1633.7647	3373.3864	0.0000	1.0000
JC69G	12	-1648.3964	3376.1856	0.0000	1.0000
JC69I	12	-1648.3964	3376.1857	0.0000	1.0000
JC69IG	13	-1647.5362	3381.0813	0.0000	1.0000
JC69	11	-1663.7502	3400.2770	0.0000	1.0000



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

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

BEGIN MRBAYES;
 Charset 591-tomodel-phylip/chr7_4885.phylip = 1 - 747;
 Partition Dummy = 1:591-tomodel-phylip/chr7_4885.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=2 rates=propinv;
 Prset applyto=(1) statefreqpr=Dirichlet(1.0,1.0,1.0,1.0) tratiopr=Beta(1.0,1.0) pinvarpr=Uniform(0.0,1.0);
END;


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

BEGIN MRBAYES;
 Charset 591-tomodel-phylip/chr7_4885.phylip = 1 - 747;
 Partition Dummy = 1:591-tomodel-phylip/chr7_4885.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=2 rates=propinv;
 Prset applyto=(1) statefreqpr=Dirichlet(1.0,1.0,1.0,1.0) tratiopr=Beta(1.0,1.0) pinvarpr=Uniform(0.0,1.0);
END;


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

BEGIN MRBAYES;
 Charset 591-tomodel-phylip/chr7_4885.phylip = 1 - 747;
 Partition Dummy = 1:591-tomodel-phylip/chr7_4885.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=2 rates=propinv;
 Prset applyto=(1) statefreqpr=Dirichlet(1.0,1.0,1.0,1.0) tratiopr=Beta(1.0,1.0) pinvarpr=Uniform(0.0,1.0);
END;


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

