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

Input data from file "591-tomodel-phylip/chr3_3116.phylip" (7 taxa, 591 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	-1490.8635	3023.2006	0.3053	0.3053
HKYG	16	-1495.5527	3024.0531	0.1994	0.5047
GTRI	20	-1491.4092	3024.2921	0.1769	0.6816
HKYI	16	-1496.0532	3025.0541	0.1209	0.8024
GTRIG	21	-1490.7281	3025.0800	0.1193	0.9217
HKYIG	17	-1495.4330	3025.9340	0.0778	0.9996
K2PG	13	-1505.7601	3038.1509	0.0002	0.9997
K2PI	13	-1506.2498	3039.1304	0.0001	0.9998
K2PIG	14	-1505.6450	3040.0192	0.0001	0.9999
SYMG	17	-1502.8620	3040.7921	0.0000	0.9999
SYMI	17	-1503.3940	3041.8560	0.0000	1.0000
SYMIG	18	-1502.7377	3042.6712	0.0000	1.0000
HKY	15	-1507.0682	3044.9713	0.0000	1.0000
GTR	19	-1502.8672	3045.0654	0.0000	1.0000
K2P	12	-1516.9498	3058.4395	0.0000	1.0000
SYM	16	-1514.2542	3061.4561	0.0000	1.0000
F81G	15	-1520.5860	3072.0067	0.0000	1.0000
F81I	15	-1520.9260	3072.6867	0.0000	1.0000
F81IG	16	-1520.5395	3074.0267	0.0000	1.0000
JC69G	12	-1529.6546	3083.8491	0.0000	1.0000
JC69I	12	-1529.9995	3084.5388	0.0000	1.0000
JC69IG	13	-1529.6066	3085.8441	0.0000	1.0000
F81	14	-1530.4044	3089.5381	0.0000	1.0000
JC69	11	-1539.5378	3101.5315	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRG	20	-1490.8635	3021.7270	0.3316	0.3316
GTRI	20	-1491.4092	3022.8184	0.1921	0.5237
HKYG	16	-1495.5527	3023.1053	0.1664	0.6901
GTRIG	21	-1490.7281	3023.4561	0.1397	0.8297
HKYI	16	-1496.0532	3024.1063	0.1009	0.9306
HKYIG	17	-1495.4330	3024.8660	0.0690	0.9997
K2PG	13	-1505.7601	3037.5201	0.0001	0.9998
K2PI	13	-1506.2498	3038.4995	0.0001	0.9999
K2PIG	14	-1505.6450	3039.2900	0.0001	0.9999
SYMG	17	-1502.8620	3039.7241	0.0000	0.9999
SYMI	17	-1503.3940	3040.7879	0.0000	1.0000
SYMIG	18	-1502.7377	3041.4754	0.0000	1.0000
GTR	19	-1502.8672	3043.7344	0.0000	1.0000
HKY	15	-1507.0682	3044.1365	0.0000	1.0000
K2P	12	-1516.9498	3057.8997	0.0000	1.0000
SYM	16	-1514.2542	3060.5083	0.0000	1.0000
F81G	15	-1520.5860	3071.1720	0.0000	1.0000
F81I	15	-1520.9260	3071.8519	0.0000	1.0000
F81IG	16	-1520.5395	3073.0790	0.0000	1.0000
JC69G	12	-1529.6546	3083.3093	0.0000	1.0000
JC69I	12	-1529.9995	3083.9990	0.0000	1.0000
JC69IG	13	-1529.6066	3085.2132	0.0000	1.0000
F81	14	-1530.4044	3088.8089	0.0000	1.0000
JC69	11	-1539.5378	3101.0756	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	16	-1495.5527	3093.2144	0.3948	0.3948
HKYI	16	-1496.0532	3094.2154	0.2393	0.6341
K2PG	13	-1505.7601	3094.4837	0.2093	0.8434
K2PI	13	-1506.2498	3095.4631	0.1283	0.9717
HKYIG	17	-1495.4330	3099.3569	0.0183	0.9900
K2PIG	14	-1505.6450	3100.6354	0.0097	0.9996
GTRG	20	-1490.8635	3109.3633	0.0001	0.9997
HKY	15	-1507.0682	3109.8637	0.0001	0.9998
GTRI	20	-1491.4092	3110.4548	0.0001	0.9999
K2P	12	-1516.9498	3110.4815	0.0001	1.0000
SYMG	17	-1502.8620	3114.2150	0.0000	1.0000
SYMI	17	-1503.3940	3115.2788	0.0000	1.0000
GTRIG	21	-1490.7281	3115.4743	0.0000	1.0000
SYMIG	18	-1502.7377	3120.3481	0.0000	1.0000
GTR	19	-1502.8672	3126.9889	0.0000	1.0000
SYM	16	-1514.2542	3130.6174	0.0000	1.0000
JC69G	12	-1529.6546	3135.8911	0.0000	1.0000
JC69I	12	-1529.9995	3136.5808	0.0000	1.0000
F81G	15	-1520.5860	3136.8992	0.0000	1.0000
F81I	15	-1520.9260	3137.5791	0.0000	1.0000
JC69IG	13	-1529.6066	3142.1768	0.0000	1.0000
F81IG	16	-1520.5395	3143.1880	0.0000	1.0000
JC69	11	-1539.5378	3149.2756	0.0000	1.0000
F81	14	-1530.4044	3150.1543	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/chr3_3116.phylip = 1 - 591;
 Partition Dummy = 1:591-tomodel-phylip/chr3_3116.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/chr3_3116.phylip = 1 - 591;
 Partition Dummy = 1:591-tomodel-phylip/chr3_3116.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/chr3_3116.phylip = 1 - 591;
 Partition Dummy = 1:591-tomodel-phylip/chr3_3116.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

