
Output from mraic.pl version 1.4.4 by Johan Nylander
-------------------------------------------------------------
11:42 08/01/2011

Input data from file "261-aggressive-phylip/5719_hg18.phylip" (29 taxa, 1238 characters)

Minimum AIC  model: GTRIG
Minimum AICc model: GTRIG
Minimum BIC  model: GTRIG


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
GTRIG	65	-5859.8604	11857.0416	1.0000	1.0000
SYMIG	62	-5877.6822	11886.0128	0.0000	1.0000
HKYIG	61	-5888.8438	11906.1196	0.0000	1.0000
GTRG	64	-5890.7407	11916.5743	0.0000	1.0000
K2PIG	58	-5901.9957	11925.7963	0.0000	1.0000
HKYG	60	-5908.8687	11943.9565	0.0000	1.0000
SYMG	61	-5910.3333	11949.0985	0.0000	1.0000
K2PG	57	-5920.3849	11960.3732	0.0000	1.0000
HKYI	60	-5946.0177	12018.2545	0.0000	1.0000
K2PI	57	-5956.9157	12033.4348	0.0000	1.0000
SYMI	61	-5953.8412	12036.1144	0.0000	1.0000
GTRI	64	-5951.0823	12037.2576	0.0000	1.0000
F81IG	60	-5998.4244	12123.0679	0.0000	1.0000
F81G	59	-6015.7153	12155.4408	0.0000	1.0000
JC69IG	57	-6019.6259	12158.8553	0.0000	1.0000
JC69G	56	-6030.3836	12178.1728	0.0000	1.0000
F81I	59	-6038.3501	12200.7103	0.0000	1.0000
JC69I	56	-6052.2334	12221.8723	0.0000	1.0000
GTR	63	-6319.9434	12772.7556	0.0000	1.0000
HKY	59	-6334.6021	12793.2145	0.0000	1.0000
SYM	60	-6335.0340	12796.2872	0.0000	1.0000
K2P	56	-6346.8773	12811.1602	0.0000	1.0000
F81	58	-6438.0064	12997.8178	0.0000	1.0000
JC69	55	-6451.6921	13018.5956	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRIG	65	-5859.8604	11849.7208	1.0000	1.0000
SYMIG	62	-5877.6822	11879.3643	0.0000	1.0000
HKYIG	61	-5888.8438	11899.6876	0.0000	1.0000
GTRG	64	-5890.7407	11909.4814	0.0000	1.0000
K2PIG	58	-5901.9957	11919.9914	0.0000	1.0000
HKYG	60	-5908.8687	11937.7373	0.0000	1.0000
SYMG	61	-5910.3333	11942.6665	0.0000	1.0000
K2PG	57	-5920.3849	11954.7698	0.0000	1.0000
HKYI	60	-5946.0177	12012.0353	0.0000	1.0000
K2PI	57	-5956.9157	12027.8314	0.0000	1.0000
SYMI	61	-5953.8412	12029.6825	0.0000	1.0000
GTRI	64	-5951.0823	12030.1647	0.0000	1.0000
F81IG	60	-5998.4244	12116.8487	0.0000	1.0000
F81G	59	-6015.7153	12149.4306	0.0000	1.0000
JC69IG	57	-6019.6259	12153.2519	0.0000	1.0000
JC69G	56	-6030.3836	12172.7672	0.0000	1.0000
F81I	59	-6038.3501	12194.7002	0.0000	1.0000
JC69I	56	-6052.2334	12216.4667	0.0000	1.0000
GTR	63	-6319.9434	12765.8868	0.0000	1.0000
HKY	59	-6334.6021	12787.2043	0.0000	1.0000
SYM	60	-6335.0340	12790.0680	0.0000	1.0000
K2P	56	-6346.8773	12805.7547	0.0000	1.0000
F81	58	-6438.0064	12992.0129	0.0000	1.0000
JC69	55	-6451.6921	13013.3841	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
GTRIG	65	-5859.8604	12182.6022	0.9992	0.9992
SYMIG	62	-5877.6822	12196.8820	0.0008	1.0000
HKYIG	61	-5888.8438	12212.0840	0.0000	1.0000
K2PIG	58	-5901.9957	12217.0240	0.0000	1.0000
GTRG	64	-5890.7407	12237.2415	0.0000	1.0000
HKYG	60	-5908.8687	12245.0125	0.0000	1.0000
K2PG	57	-5920.3849	12246.6812	0.0000	1.0000
SYMG	61	-5910.3333	12255.0629	0.0000	1.0000
HKYI	60	-5946.0177	12319.3105	0.0000	1.0000
K2PI	57	-5956.9157	12319.7428	0.0000	1.0000
SYMI	61	-5953.8412	12342.0789	0.0000	1.0000
GTRI	64	-5951.0823	12357.9249	0.0000	1.0000
F81IG	60	-5998.4244	12424.1239	0.0000	1.0000
JC69IG	57	-6019.6259	12445.1633	0.0000	1.0000
F81G	59	-6015.7153	12451.5845	0.0000	1.0000
JC69G	56	-6030.3836	12459.5574	0.0000	1.0000
F81I	59	-6038.3501	12496.8541	0.0000	1.0000
JC69I	56	-6052.2334	12503.2569	0.0000	1.0000
GTR	63	-6319.9434	13088.5257	0.0000	1.0000
HKY	59	-6334.6021	13089.3582	0.0000	1.0000
K2P	56	-6346.8773	13092.5448	0.0000	1.0000
SYM	60	-6335.0340	13097.3432	0.0000	1.0000
F81	58	-6438.0064	13289.0455	0.0000	1.0000
JC69	55	-6451.6921	13295.0530	0.0000	1.0000



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

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

BEGIN MRBAYES;
 Charset 261-aggressive-phylip/5719_hg18.phylip = 1 - 1238;
 Partition Dummy = 1:261-aggressive-phylip/5719_hg18.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=6 rates=invgamma;
 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) pinvarpr=Uniform(0.0,1.0);
END;


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

BEGIN MRBAYES;
 Charset 261-aggressive-phylip/5719_hg18.phylip = 1 - 1238;
 Partition Dummy = 1:261-aggressive-phylip/5719_hg18.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=6 rates=invgamma;
 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) pinvarpr=Uniform(0.0,1.0);
END;


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

BEGIN MRBAYES;
 Charset 261-aggressive-phylip/5719_hg18.phylip = 1 - 1238;
 Partition Dummy = 1:261-aggressive-phylip/5719_hg18.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=6 rates=invgamma;
 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) pinvarpr=Uniform(0.0,1.0);
END;


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

