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

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

Minimum AIC  model: GTRG
Minimum AICc model: GTRG
Minimum BIC  model: K2PG


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	-2558.7627	5158.6729	0.5908	0.5908
GTRIG	21	-2558.1372	5159.5383	0.3833	0.9741
GTRI	20	-2561.9620	5165.0715	0.0241	0.9982
SYMG	17	-2568.2407	5171.3140	0.0011	0.9992
SYMIG	18	-2567.6501	5172.2320	0.0007	0.9999
SYMI	17	-2570.7789	5176.3905	0.0001	1.0000
HKYG	16	-2574.1697	5181.0785	0.0000	1.0000
K2PG	13	-2577.7245	5181.9417	0.0000	1.0000
HKYIG	17	-2573.7928	5182.4184	0.0000	1.0000
K2PIG	14	-2577.3571	5183.2832	0.0000	1.0000
HKYI	16	-2575.8494	5184.4379	0.0000	1.0000
K2PI	13	-2579.3684	5185.2294	0.0000	1.0000
GTR	19	-2608.6236	5256.2841	0.0000	1.0000
SYM	16	-2618.4611	5269.6614	0.0000	1.0000
HKY	15	-2622.0643	5274.7799	0.0000	1.0000
K2P	12	-2626.7087	5277.8390	0.0000	1.0000
JC69G	12	-2627.4771	5279.3759	0.0000	1.0000
F81G	15	-2624.3659	5279.3830	0.0000	1.0000
JC69IG	13	-2627.2778	5281.0483	0.0000	1.0000
F81IG	16	-2624.1638	5281.0667	0.0000	1.0000
JC69I	12	-2629.1022	5282.6260	0.0000	1.0000
F81I	15	-2626.0136	5282.6785	0.0000	1.0000
F81	14	-2674.0318	5376.6326	0.0000	1.0000
JC69	11	-2677.6484	5377.6530	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRG	20	-2558.7627	5157.5254	0.5776	0.5776
GTRIG	21	-2558.1372	5158.2743	0.3972	0.9749
GTRI	20	-2561.9620	5163.9239	0.0236	0.9984
SYMG	17	-2568.2407	5170.4813	0.0009	0.9993
SYMIG	18	-2567.6501	5171.3002	0.0006	0.9999
SYMI	17	-2570.7789	5175.5578	0.0001	1.0000
HKYG	16	-2574.1697	5180.3393	0.0000	1.0000
K2PG	13	-2577.7245	5181.4491	0.0000	1.0000
HKYIG	17	-2573.7928	5181.5857	0.0000	1.0000
K2PIG	14	-2577.3571	5182.7141	0.0000	1.0000
HKYI	16	-2575.8494	5183.6987	0.0000	1.0000
K2PI	13	-2579.3684	5184.7368	0.0000	1.0000
GTR	19	-2608.6236	5255.2473	0.0000	1.0000
SYM	16	-2618.4611	5268.9222	0.0000	1.0000
HKY	15	-2622.0643	5274.1287	0.0000	1.0000
K2P	12	-2626.7087	5277.4173	0.0000	1.0000
F81G	15	-2624.3659	5278.7318	0.0000	1.0000
JC69G	12	-2627.4771	5278.9543	0.0000	1.0000
F81IG	16	-2624.1638	5280.3276	0.0000	1.0000
JC69IG	13	-2627.2778	5280.5557	0.0000	1.0000
F81I	15	-2626.0136	5282.0272	0.0000	1.0000
JC69I	12	-2629.1022	5282.2043	0.0000	1.0000
F81	14	-2674.0318	5376.0635	0.0000	1.0000
JC69	11	-2677.6484	5377.2968	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
K2PG	13	-2577.7245	5241.5619	0.7747	0.7747
K2PI	13	-2579.3684	5244.8497	0.1497	0.9244
K2PIG	14	-2577.3571	5247.4510	0.0408	0.9652
SYMG	17	-2568.2407	5249.0904	0.0180	0.9832
GTRG	20	-2558.7627	5250.0067	0.0114	0.9945
SYMI	17	-2570.7789	5254.1669	0.0014	0.9960
HKYG	16	-2574.1697	5254.3244	0.0013	0.9973
SYMIG	18	-2567.6501	5254.5333	0.0012	0.9984
GTRIG	21	-2558.1372	5255.3797	0.0008	0.9992
GTRI	20	-2561.9620	5256.4052	0.0005	0.9997
HKYI	16	-2575.8494	5257.6838	0.0002	0.9999
HKYIG	17	-2573.7928	5260.1948	0.0001	1.0000
K2P	12	-2626.7087	5332.9061	0.0000	1.0000
JC69G	12	-2627.4771	5334.4431	0.0000	1.0000
JC69I	12	-2629.1022	5337.6931	0.0000	1.0000
JC69IG	13	-2627.2778	5340.6685	0.0000	1.0000
SYM	16	-2618.4611	5342.9073	0.0000	1.0000
GTR	19	-2608.6236	5343.1045	0.0000	1.0000
HKY	15	-2622.0643	5343.4896	0.0000	1.0000
F81G	15	-2624.3659	5348.0927	0.0000	1.0000
F81I	15	-2626.0136	5351.3882	0.0000	1.0000
F81IG	16	-2624.1638	5354.3126	0.0000	1.0000
JC69	11	-2677.6484	5428.1615	0.0000	1.0000
F81	14	-2674.0318	5440.8004	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/chrUn_random_7370.phylip = 1 - 753;
 Partition Dummy = 1:591-tomodel-phylip/chrUn_random_7370.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/chrUn_random_7370.phylip = 1 - 753;
 Partition Dummy = 1:591-tomodel-phylip/chrUn_random_7370.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 (K2PG)]

BEGIN MRBAYES;
 Charset 591-tomodel-phylip/chrUn_random_7370.phylip = 1 - 753;
 Partition Dummy = 1:591-tomodel-phylip/chrUn_random_7370.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=2 rates=gamma;
 Prset applyto=(1) statefreqpr=Fixed(Equal) tratiopr=Beta(1.0,1.0) shapepr=Uniform(0.1,50.0);
END;


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

