
Output from mraic.pl version 1.4.4 by Johan Nylander
-------------------------------------------------------------
00:55 02/17/2011

Input data from file "183-phylip/chr5_11146.phylip" (29 taxa, 124 characters)

Minimum AIC  model: HKYI
Minimum AICc model: K2P
Minimum BIC  model: K2P


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
K2P	56	-371.3241	949.9318	0.7770	0.7770
K2PI	57	-369.6852	953.5521	0.1271	0.9042
K2PG	57	-370.0754	954.3326	0.0861	0.9902
HKY	59	-365.6557	959.9364	0.0052	0.9955
K2PIG	58	-369.6965	960.6852	0.0036	0.9991
HKYI	60	-364.2259	964.6422	0.0005	0.9996
HKYG	60	-364.5405	965.2715	0.0004	0.9999
JC69	55	-384.3922	969.3727	0.0000	1.0000
HKYIG	61	-364.2348	972.4697	0.0000	1.0000
JC69I	56	-382.7529	972.7895	0.0000	1.0000
JC69G	56	-383.1480	973.5796	0.0000	1.0000
SYM	60	-368.9001	973.9907	0.0000	1.0000
SYMI	61	-367.3120	978.6240	0.0000	1.0000
SYMG	61	-367.7254	979.4507	0.0000	1.0000
JC69IG	57	-382.7645	979.7108	0.0000	1.0000
F81	58	-379.4774	980.2471	0.0000	1.0000
F81I	59	-377.9411	984.5072	0.0000	1.0000
F81G	59	-378.3385	985.3019	0.0000	1.0000
GTR	63	-363.1000	986.6000	0.0000	1.0000
SYMIG	62	-367.3238	986.7132	0.0000	1.0000
F81IG	60	-377.9527	992.0959	0.0000	1.0000
GTRI	64	-361.7548	992.5266	0.0000	1.0000
GTRG	64	-362.1190	993.2550	0.0000	1.0000
GTRIG	65	-361.7650	1001.4611	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
HKYI	60	-364.2259	848.4517	0.2772	0.2772
HKYG	60	-364.5405	849.0810	0.2023	0.4795
HKY	59	-365.6557	849.3114	0.1803	0.6598
HKYIG	61	-364.2348	850.4697	0.1011	0.7609
GTRI	64	-361.7548	851.5097	0.0601	0.8210
GTR	63	-363.1000	852.2000	0.0425	0.8635
GTRG	64	-362.1190	852.2381	0.0417	0.9052
K2PI	57	-369.6852	853.3703	0.0237	0.9289
GTRIG	65	-361.7650	853.5300	0.0219	0.9508
K2PG	57	-370.0754	854.1508	0.0160	0.9669
K2P	56	-371.3241	854.6482	0.0125	0.9794
K2PIG	58	-369.6965	855.3929	0.0086	0.9880
SYMI	61	-367.3120	856.6240	0.0047	0.9926
SYMG	61	-367.7254	857.4507	0.0031	0.9957
SYM	60	-368.9001	857.8002	0.0026	0.9983
SYMIG	62	-367.3238	858.6477	0.0017	1.0000
F81I	59	-377.9411	873.8822	0.0000	1.0000
F81G	59	-378.3385	874.6769	0.0000	1.0000
F81	58	-379.4774	874.9548	0.0000	1.0000
F81IG	60	-377.9527	875.9055	0.0000	1.0000
JC69I	56	-382.7529	877.5059	0.0000	1.0000
JC69G	56	-383.1480	878.2960	0.0000	1.0000
JC69	55	-384.3922	878.7845	0.0000	1.0000
JC69IG	57	-382.7645	879.5290	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
K2P	56	-371.3241	1012.5839	0.4606	0.4606
K2PI	57	-369.6852	1014.1264	0.2130	0.6736
K2PG	57	-370.0754	1014.9068	0.1442	0.8178
HKY	59	-365.6557	1015.7080	0.0966	0.9144
HKYI	60	-364.2259	1017.6686	0.0362	0.9507
HKYG	60	-364.5405	1018.2979	0.0265	0.9771
K2PIG	58	-369.6965	1018.9693	0.0189	0.9960
HKYIG	61	-364.2348	1022.5069	0.0032	0.9992
SYM	60	-368.9001	1027.0171	0.0003	0.9996
SYMI	61	-367.3120	1028.6612	0.0001	0.9997
SYMG	61	-367.7254	1029.4879	0.0001	0.9998
GTR	63	-363.1000	1029.8778	0.0001	0.9999
GTRI	64	-361.7548	1032.0077	0.0000	0.9999
GTRG	64	-362.1190	1032.7361	0.0000	1.0000
SYMIG	62	-367.3238	1033.5051	0.0000	1.0000
JC69	55	-384.3922	1033.9000	0.0000	1.0000
JC69I	56	-382.7529	1035.4417	0.0000	1.0000
JC69G	56	-383.1480	1036.2318	0.0000	1.0000
GTRIG	65	-361.7650	1036.8483	0.0000	1.0000
F81	58	-379.4774	1038.5311	0.0000	1.0000
F81I	59	-377.9411	1040.2788	0.0000	1.0000
JC69IG	57	-382.7645	1040.2850	0.0000	1.0000
F81G	59	-378.3385	1041.0736	0.0000	1.0000
F81IG	60	-377.9527	1045.1224	0.0000	1.0000



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

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

BEGIN MRBAYES;
 Charset 183-phylip/chr5_11146.phylip = 1 - 124;
 Partition Dummy = 1:183-phylip/chr5_11146.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=2 rates=propinv;
 Prset applyto=(1) statefreqpr=Dirichlet(1.0,1.0,1.0,1.0) tratiopr=Beta(1.0,1.0) pinvarpr=Uniform(0.0,1.0);
END;


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

BEGIN MRBAYES;
 Charset 183-phylip/chr5_11146.phylip = 1 - 124;
 Partition Dummy = 1:183-phylip/chr5_11146.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=2 rates=equal;
 Prset applyto=(1) statefreqpr=Fixed(Equal) tratiopr=Beta(1.0,1.0);
END;


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

BEGIN MRBAYES;
 Charset 183-phylip/chr5_11146.phylip = 1 - 124;
 Partition Dummy = 1:183-phylip/chr5_11146.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=2 rates=equal;
 Prset applyto=(1) statefreqpr=Fixed(Equal) tratiopr=Beta(1.0,1.0);
END;


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

