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

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

Minimum AIC  model: HKYG
Minimum AICc model: HKYG
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
HKYG	16	-1406.4582	2846.0665	0.3743	0.3743
GTRG	20	-1402.7548	2847.3006	0.2019	0.5762
HKYI	16	-1407.4118	2847.9738	0.1442	0.7204
HKYIG	17	-1406.3477	2847.9920	0.1429	0.8633
GTRIG	21	-1402.5946	2849.1635	0.0796	0.9428
GTRI	20	-1404.0167	2849.8244	0.0572	1.0000
HKY	15	-1429.2333	2889.4793	0.0000	1.0000
GTR	19	-1427.0339	2893.6848	0.0000	1.0000
F81G	15	-1451.3949	2933.8025	0.0000	1.0000
F81I	15	-1452.3467	2935.7062	0.0000	1.0000
F81IG	16	-1451.2795	2935.7091	0.0000	1.0000
K2PG	13	-1454.9851	2936.7348	0.0000	1.0000
K2PI	13	-1455.7668	2938.2982	0.0000	1.0000
K2PIG	14	-1454.8341	2938.5524	0.0000	1.0000
SYMG	17	-1453.5424	2942.3815	0.0000	1.0000
SYMI	17	-1454.3679	2944.0325	0.0000	1.0000
SYMIG	18	-1453.3518	2944.1559	0.0000	1.0000
F81	14	-1471.6220	2972.1282	0.0000	1.0000
K2P	12	-1476.5748	2977.8037	0.0000	1.0000
SYM	16	-1475.4644	2984.0789	0.0000	1.0000
JC69G	12	-1490.2534	3005.1609	0.0000	1.0000
JC69IG	13	-1490.1053	3006.9753	0.0000	1.0000
JC69I	12	-1491.2352	3007.1245	0.0000	1.0000
JC69	11	-1510.7587	3044.0698	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
HKYG	16	-1406.4582	2844.9164	0.3256	0.3256
GTRG	20	-1402.7548	2845.5096	0.2421	0.5677
HKYIG	17	-1406.3477	2846.6953	0.1338	0.7015
HKYI	16	-1407.4118	2846.8237	0.1255	0.8270
GTRIG	21	-1402.5946	2847.1891	0.1045	0.9315
GTRI	20	-1404.0167	2848.0334	0.0685	1.0000
HKY	15	-1429.2333	2888.4666	0.0000	1.0000
GTR	19	-1427.0339	2892.0678	0.0000	1.0000
F81G	15	-1451.3949	2932.7898	0.0000	1.0000
F81IG	16	-1451.2795	2934.5590	0.0000	1.0000
F81I	15	-1452.3467	2934.6935	0.0000	1.0000
K2PG	13	-1454.9851	2935.9701	0.0000	1.0000
K2PI	13	-1455.7668	2937.5335	0.0000	1.0000
K2PIG	14	-1454.8341	2937.6682	0.0000	1.0000
SYMG	17	-1453.5424	2941.0849	0.0000	1.0000
SYMIG	18	-1453.3518	2942.7036	0.0000	1.0000
SYMI	17	-1454.3679	2942.7359	0.0000	1.0000
F81	14	-1471.6220	2971.2440	0.0000	1.0000
K2P	12	-1476.5748	2977.1496	0.0000	1.0000
SYM	16	-1475.4644	2982.9288	0.0000	1.0000
JC69G	12	-1490.2534	3004.5068	0.0000	1.0000
JC69IG	13	-1490.1053	3006.2106	0.0000	1.0000
JC69I	12	-1491.2352	3006.4704	0.0000	1.0000
JC69	11	-1510.7587	3043.5175	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	16	-1406.4582	2912.0268	0.6964	0.6964
HKYI	16	-1407.4118	2913.9341	0.2683	0.9647
HKYIG	17	-1406.3477	2918.0002	0.0351	0.9998
GTRG	20	-1402.7548	2929.3977	0.0001	1.0000
GTRI	20	-1404.0167	2931.9215	0.0000	1.0000
GTRIG	21	-1402.5946	2935.2717	0.0000	1.0000
HKY	15	-1429.2333	2951.3827	0.0000	1.0000
GTR	19	-1427.0339	2971.7615	0.0000	1.0000
K2PG	13	-1454.9851	2990.4974	0.0000	1.0000
K2PI	13	-1455.7668	2992.0608	0.0000	1.0000
F81G	15	-1451.3949	2995.7059	0.0000	1.0000
K2PIG	14	-1454.8341	2996.3899	0.0000	1.0000
F81I	15	-1452.3467	2997.6096	0.0000	1.0000
F81IG	16	-1451.2795	3001.6695	0.0000	1.0000
SYMG	17	-1453.5424	3012.3898	0.0000	1.0000
SYMI	17	-1454.3679	3014.0408	0.0000	1.0000
SYMIG	18	-1453.3518	3018.2029	0.0000	1.0000
K2P	12	-1476.5748	3027.4824	0.0000	1.0000
F81	14	-1471.6220	3029.9657	0.0000	1.0000
SYM	16	-1475.4644	3050.0393	0.0000	1.0000
JC69G	12	-1490.2534	3054.8397	0.0000	1.0000
JC69I	12	-1491.2352	3056.8033	0.0000	1.0000
JC69IG	13	-1490.1053	3060.7379	0.0000	1.0000
JC69	11	-1510.7587	3089.6559	0.0000	1.0000



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

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

BEGIN MRBAYES;
 Charset 591-tomodel-phylip/chr3_5845.phylip = 1 - 490;
 Partition Dummy = 1:591-tomodel-phylip/chr3_5845.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;


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

BEGIN MRBAYES;
 Charset 591-tomodel-phylip/chr3_5845.phylip = 1 - 490;
 Partition Dummy = 1:591-tomodel-phylip/chr3_5845.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;


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

BEGIN MRBAYES;
 Charset 591-tomodel-phylip/chr3_5845.phylip = 1 - 490;
 Partition Dummy = 1:591-tomodel-phylip/chr3_5845.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

