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

Input data from file "591-tomodel-phylip/chr2_17003.phylip" (7 taxa, 614 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	-1593.5000	3228.4166	0.4561	0.4561
GTRI	20	-1593.6614	3228.7393	0.3882	0.8443
GTRIG	21	-1593.5112	3230.5832	0.1544	0.9987
HKYG	16	-1604.3661	3241.6434	0.0006	0.9993
HKYI	16	-1604.6097	3242.1306	0.0005	0.9998
HKYIG	17	-1604.3668	3243.7605	0.0002	1.0000
GTR	19	-1604.3878	3248.0550	0.0000	1.0000
HKY	15	-1614.2014	3259.2054	0.0000	1.0000
SYMG	17	-1625.5286	3286.0841	0.0000	1.0000
SYMI	17	-1625.6439	3286.3147	0.0000	1.0000
SYMIG	18	-1625.5348	3288.2193	0.0000	1.0000
K2PG	13	-1631.0018	3288.6102	0.0000	1.0000
K2PI	13	-1631.2680	3289.1426	0.0000	1.0000
K2PIG	14	-1630.9978	3290.6969	0.0000	1.0000
SYM	16	-1635.3325	3303.5762	0.0000	1.0000
F81G	15	-1637.5913	3305.9854	0.0000	1.0000
K2P	12	-1640.9499	3306.4190	0.0000	1.0000
F81I	15	-1637.9284	3306.6594	0.0000	1.0000
F81IG	16	-1637.5836	3308.0784	0.0000	1.0000
F81	14	-1646.5246	3321.7503	0.0000	1.0000
JC69G	12	-1660.6933	3345.9057	0.0000	1.0000
JC69I	12	-1661.0467	3346.6126	0.0000	1.0000
JC69IG	13	-1660.6852	3347.9770	0.0000	1.0000
JC69	11	-1669.9076	3362.2537	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRG	20	-1593.5000	3227.0001	0.4510	0.4510
GTRI	20	-1593.6614	3227.3228	0.3838	0.8349
GTRIG	21	-1593.5112	3229.0224	0.1641	0.9990
HKYG	16	-1604.3661	3240.7321	0.0005	0.9994
HKYI	16	-1604.6097	3241.2193	0.0004	0.9998
HKYIG	17	-1604.3668	3242.7337	0.0002	1.0000
GTR	19	-1604.3878	3246.7755	0.0000	1.0000
HKY	15	-1614.2014	3258.4027	0.0000	1.0000
SYMG	17	-1625.5286	3285.0572	0.0000	1.0000
SYMI	17	-1625.6439	3285.2878	0.0000	1.0000
SYMIG	18	-1625.5348	3287.0697	0.0000	1.0000
K2PG	13	-1631.0018	3288.0035	0.0000	1.0000
K2PI	13	-1631.2680	3288.5360	0.0000	1.0000
K2PIG	14	-1630.9978	3289.9957	0.0000	1.0000
SYM	16	-1635.3325	3302.6650	0.0000	1.0000
F81G	15	-1637.5913	3305.1827	0.0000	1.0000
F81I	15	-1637.9284	3305.8567	0.0000	1.0000
K2P	12	-1640.9499	3305.8999	0.0000	1.0000
F81IG	16	-1637.5836	3307.1672	0.0000	1.0000
F81	14	-1646.5246	3321.0492	0.0000	1.0000
JC69G	12	-1660.6933	3345.3865	0.0000	1.0000
JC69I	12	-1661.0467	3346.0935	0.0000	1.0000
JC69IG	13	-1660.6852	3347.3704	0.0000	1.0000
JC69	11	-1669.9076	3361.8152	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	16	-1604.3661	3311.4521	0.4789	0.4789
HKYI	16	-1604.6097	3311.9393	0.3754	0.8542
GTRG	20	-1593.5000	3315.4000	0.0665	0.9208
GTRI	20	-1593.6614	3315.7227	0.0566	0.9774
HKYIG	17	-1604.3668	3317.8736	0.0193	0.9967
GTRIG	21	-1593.5112	3321.8423	0.0027	0.9993
HKY	15	-1614.2014	3324.7026	0.0006	1.0000
GTR	19	-1604.3878	3330.7554	0.0000	1.0000
K2PG	13	-1631.0018	3345.4635	0.0000	1.0000
K2PI	13	-1631.2680	3345.9959	0.0000	1.0000
K2PIG	14	-1630.9978	3351.8756	0.0000	1.0000
K2P	12	-1640.9499	3358.9398	0.0000	1.0000
SYMG	17	-1625.5286	3360.1971	0.0000	1.0000
SYMI	17	-1625.6439	3360.4278	0.0000	1.0000
SYMIG	18	-1625.5348	3366.6296	0.0000	1.0000
F81G	15	-1637.5913	3371.4826	0.0000	1.0000
F81I	15	-1637.9284	3372.1566	0.0000	1.0000
SYM	16	-1635.3325	3373.3849	0.0000	1.0000
F81IG	16	-1637.5836	3377.8871	0.0000	1.0000
F81	14	-1646.5246	3382.9291	0.0000	1.0000
JC69G	12	-1660.6933	3398.4265	0.0000	1.0000
JC69I	12	-1661.0467	3399.1334	0.0000	1.0000
JC69IG	13	-1660.6852	3404.8303	0.0000	1.0000
JC69	11	-1669.9076	3410.4351	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_17003.phylip = 1 - 614;
 Partition Dummy = 1:591-tomodel-phylip/chr2_17003.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_17003.phylip = 1 - 614;
 Partition Dummy = 1:591-tomodel-phylip/chr2_17003.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_17003.phylip = 1 - 614;
 Partition Dummy = 1:591-tomodel-phylip/chr2_17003.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

