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

Input data from file "591-tomodel-phylip/chr2_13006.phylip" (7 taxa, 737 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	-2359.5567	4760.2866	0.5294	0.5294
GTRIG	21	-2359.3957	4762.0837	0.2156	0.7450
GTRI	20	-2361.1834	4763.5400	0.1041	0.8491
HKYG	16	-2365.5355	4763.8266	0.0902	0.9393
HKYIG	17	-2365.4606	4765.7723	0.0341	0.9733
HKYI	16	-2366.7563	4766.2681	0.0266	0.9999
SYMG	17	-2373.0909	4781.0330	0.0000	1.0000
K2PG	13	-2377.3742	4781.2518	0.0000	1.0000
SYMIG	18	-2372.9561	4782.8649	0.0000	1.0000
K2PIG	14	-2377.3093	4783.2002	0.0000	1.0000
K2PI	13	-2378.5230	4783.5494	0.0000	1.0000
SYMI	17	-2374.5979	4784.0469	0.0000	1.0000
GTR	19	-2391.8100	4822.6800	0.0000	1.0000
HKY	15	-2396.8260	4824.3177	0.0000	1.0000
K2P	12	-2409.5040	4843.4389	0.0000	1.0000
SYM	16	-2405.6985	4844.1525	0.0000	1.0000
F81G	15	-2415.3176	4861.3010	0.0000	1.0000
F81I	15	-2416.0981	4862.8619	0.0000	1.0000
F81IG	16	-2415.3093	4863.3741	0.0000	1.0000
JC69G	12	-2425.7773	4875.9856	0.0000	1.0000
JC69I	12	-2426.6087	4877.6483	0.0000	1.0000
JC69IG	13	-2425.7663	4878.0361	0.0000	1.0000
F81	14	-2443.7009	4915.9835	0.0000	1.0000
JC69	11	-2454.5924	4931.5490	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRG	20	-2359.5567	4759.1134	0.5369	0.5369
GTRIG	21	-2359.3957	4760.7914	0.2320	0.7689
GTRI	20	-2361.1834	4762.3669	0.1055	0.8744
HKYG	16	-2365.5355	4763.0711	0.0742	0.9486
HKYIG	17	-2365.4606	4764.9211	0.0294	0.9781
HKYI	16	-2366.7563	4765.5126	0.0219	1.0000
SYMG	17	-2373.0909	4780.1818	0.0000	1.0000
K2PG	13	-2377.3742	4780.7483	0.0000	1.0000
SYMIG	18	-2372.9561	4781.9123	0.0000	1.0000
K2PIG	14	-2377.3093	4782.6185	0.0000	1.0000
K2PI	13	-2378.5230	4783.0460	0.0000	1.0000
SYMI	17	-2374.5979	4783.1957	0.0000	1.0000
GTR	19	-2391.8100	4821.6201	0.0000	1.0000
HKY	15	-2396.8260	4823.6520	0.0000	1.0000
K2P	12	-2409.5040	4843.0080	0.0000	1.0000
SYM	16	-2405.6985	4843.3969	0.0000	1.0000
F81G	15	-2415.3176	4860.6353	0.0000	1.0000
F81I	15	-2416.0981	4862.1961	0.0000	1.0000
F81IG	16	-2415.3093	4862.6185	0.0000	1.0000
JC69G	12	-2425.7773	4875.5547	0.0000	1.0000
JC69I	12	-2426.6087	4877.2174	0.0000	1.0000
JC69IG	13	-2425.7663	4877.5326	0.0000	1.0000
F81	14	-2443.7009	4915.4018	0.0000	1.0000
JC69	11	-2454.5924	4931.1849	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	16	-2365.5355	4836.7125	0.6529	0.6529
HKYI	16	-2366.7563	4839.1540	0.1926	0.8455
K2PG	13	-2377.3742	4840.5819	0.0943	0.9399
K2PI	13	-2378.5230	4842.8796	0.0299	0.9698
HKYIG	17	-2365.4606	4843.1651	0.0259	0.9957
K2PIG	14	-2377.3093	4847.0548	0.0037	0.9994
GTRG	20	-2359.5567	4851.1652	0.0005	0.9999
GTRI	20	-2361.1834	4854.4186	0.0001	1.0000
GTRIG	21	-2359.3957	4857.4458	0.0000	1.0000
SYMG	17	-2373.0909	4858.4258	0.0000	1.0000
SYMI	17	-2374.5979	4861.4397	0.0000	1.0000
SYMIG	18	-2372.9561	4864.7589	0.0000	1.0000
HKY	15	-2396.8260	4892.6908	0.0000	1.0000
K2P	12	-2409.5040	4898.2390	0.0000	1.0000
GTR	19	-2391.8100	4909.0692	0.0000	1.0000
SYM	16	-2405.6985	4917.0383	0.0000	1.0000
F81G	15	-2415.3176	4929.6741	0.0000	1.0000
JC69G	12	-2425.7773	4930.7857	0.0000	1.0000
F81I	15	-2416.0981	4931.2350	0.0000	1.0000
JC69I	12	-2426.6087	4932.4485	0.0000	1.0000
F81IG	16	-2415.3093	4936.2599	0.0000	1.0000
JC69IG	13	-2425.7663	4937.3663	0.0000	1.0000
F81	14	-2443.7009	4979.8380	0.0000	1.0000
JC69	11	-2454.5924	4981.8133	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_13006.phylip = 1 - 737;
 Partition Dummy = 1:591-tomodel-phylip/chr2_13006.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_13006.phylip = 1 - 737;
 Partition Dummy = 1:591-tomodel-phylip/chr2_13006.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_13006.phylip = 1 - 737;
 Partition Dummy = 1:591-tomodel-phylip/chr2_13006.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

