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

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

Minimum AIC  model: GTRI
Minimum AICc model: GTRI
Minimum BIC  model: GTRI


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
GTRI	20	-2821.5011	5684.0368	0.6699	0.6699
GTRIG	21	-2821.5022	5686.1437	0.2336	0.9036
GTRG	20	-2823.4394	5687.9132	0.0964	1.0000
HKYI	16	-2840.2174	5713.1015	0.0000	1.0000
HKYIG	17	-2840.2247	5715.2004	0.0000	1.0000
HKYG	16	-2841.7796	5716.2259	0.0000	1.0000
GTR	19	-2847.5168	5733.9683	0.0000	1.0000
HKY	15	-2865.4566	5761.5008	0.0000	1.0000
SYMI	17	-2893.5701	5821.8912	0.0000	1.0000
SYMIG	18	-2893.5529	5823.9462	0.0000	1.0000
SYMG	17	-2894.6222	5823.9953	0.0000	1.0000
K2PI	13	-2901.1784	5828.8012	0.0000	1.0000
K2PG	13	-2901.9623	5830.3690	0.0000	1.0000
K2PIG	14	-2901.1297	5830.7729	0.0000	1.0000
F81I	15	-2914.3823	5859.3520	0.0000	1.0000
F81IG	16	-2914.3950	5861.4567	0.0000	1.0000
F81G	15	-2916.0048	5862.5971	0.0000	1.0000
SYM	16	-2920.3176	5873.3019	0.0000	1.0000
K2P	12	-2927.6322	5879.6448	0.0000	1.0000
F81	14	-2937.3518	5903.2171	0.0000	1.0000
JC69I	12	-2962.0655	5948.5115	0.0000	1.0000
JC69IG	13	-2962.0672	5950.5789	0.0000	1.0000
JC69G	12	-2963.2847	5950.9498	0.0000	1.0000
JC69	11	-2986.1629	5994.6473	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRI	20	-2821.5011	5683.0023	0.6616	0.6616
GTRIG	21	-2821.5022	5685.0044	0.2431	0.9048
GTRG	20	-2823.4394	5686.8787	0.0952	1.0000
HKYI	16	-2840.2174	5712.4348	0.0000	1.0000
HKYIG	17	-2840.2247	5714.4494	0.0000	1.0000
HKYG	16	-2841.7796	5715.5592	0.0000	1.0000
GTR	19	-2847.5168	5733.0335	0.0000	1.0000
HKY	15	-2865.4566	5760.9133	0.0000	1.0000
SYMI	17	-2893.5701	5821.1403	0.0000	1.0000
SYMIG	18	-2893.5529	5823.1059	0.0000	1.0000
SYMG	17	-2894.6222	5823.2443	0.0000	1.0000
K2PI	13	-2901.1784	5828.3568	0.0000	1.0000
K2PG	13	-2901.9623	5829.9245	0.0000	1.0000
K2PIG	14	-2901.1297	5830.2594	0.0000	1.0000
F81I	15	-2914.3823	5858.7645	0.0000	1.0000
F81IG	16	-2914.3950	5860.7901	0.0000	1.0000
F81G	15	-2916.0048	5862.0096	0.0000	1.0000
SYM	16	-2920.3176	5872.6352	0.0000	1.0000
K2P	12	-2927.6322	5879.2644	0.0000	1.0000
F81	14	-2937.3518	5902.7037	0.0000	1.0000
JC69I	12	-2962.0655	5948.1310	0.0000	1.0000
JC69IG	13	-2962.0672	5950.1345	0.0000	1.0000
JC69G	12	-2963.2847	5950.5693	0.0000	1.0000
JC69	11	-2986.1629	5994.3258	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
GTRI	20	-2821.5011	5777.5030	0.8439	0.8439
GTRG	20	-2823.4394	5781.3794	0.1215	0.9654
GTRIG	21	-2821.5022	5784.2301	0.0292	0.9946
HKYI	16	-2840.2174	5788.0353	0.0044	0.9989
HKYG	16	-2841.7796	5791.1597	0.0009	0.9999
HKYIG	17	-2840.2247	5794.7750	0.0001	1.0000
GTR	19	-2847.5168	5822.8092	0.0000	1.0000
HKY	15	-2865.4566	5831.7888	0.0000	1.0000
K2PI	13	-2901.1784	5889.7822	0.0000	1.0000
K2PG	13	-2901.9623	5891.3500	0.0000	1.0000
K2PIG	14	-2901.1297	5896.4099	0.0000	1.0000
SYMI	17	-2893.5701	5901.4659	0.0000	1.0000
SYMG	17	-2894.6222	5903.5699	0.0000	1.0000
SYMIG	18	-2893.5529	5908.1565	0.0000	1.0000
F81I	15	-2914.3823	5929.6400	0.0000	1.0000
F81G	15	-2916.0048	5932.8851	0.0000	1.0000
K2P	12	-2927.6322	5935.9648	0.0000	1.0000
F81IG	16	-2914.3950	5936.3906	0.0000	1.0000
SYM	16	-2920.3176	5948.2358	0.0000	1.0000
F81	14	-2937.3518	5968.8542	0.0000	1.0000
JC69I	12	-2962.0655	6004.8314	0.0000	1.0000
JC69G	12	-2963.2847	6007.2697	0.0000	1.0000
JC69IG	13	-2962.0672	6011.5599	0.0000	1.0000
JC69	11	-2986.1629	6046.3012	0.0000	1.0000



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

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

BEGIN MRBAYES;
 Charset 591-tomodel-phylip/chr1_31743.phylip = 1 - 833;
 Partition Dummy = 1:591-tomodel-phylip/chr1_31743.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=6 rates=propinv;
 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) pinvarpr=Uniform(0.0,1.0);
END;


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

BEGIN MRBAYES;
 Charset 591-tomodel-phylip/chr1_31743.phylip = 1 - 833;
 Partition Dummy = 1:591-tomodel-phylip/chr1_31743.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=6 rates=propinv;
 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) pinvarpr=Uniform(0.0,1.0);
END;


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

BEGIN MRBAYES;
 Charset 591-tomodel-phylip/chr1_31743.phylip = 1 - 833;
 Partition Dummy = 1:591-tomodel-phylip/chr1_31743.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=6 rates=propinv;
 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) pinvarpr=Uniform(0.0,1.0);
END;


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

