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

Input data from file "591-tomodel-phylip/chr6_3291.phylip" (7 taxa, 540 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	-1115.6598	2272.9381	0.4126	0.4126
GTRI	20	-1115.6870	2272.9926	0.4015	0.8141
GTRIG	21	-1115.5696	2274.9229	0.1529	0.9671
HKYG	16	-1123.3747	2279.7895	0.0134	0.9805
HKYI	16	-1123.3873	2279.8148	0.0133	0.9938
HKYIG	17	-1123.3375	2281.8474	0.0048	0.9986
GTR	19	-1122.9377	2285.3369	0.0008	0.9994
F81I	15	-1128.5147	2287.9455	0.0002	0.9996
F81G	15	-1128.5155	2287.9471	0.0002	0.9999
F81IG	16	-1128.5171	2290.0743	0.0001	0.9999
HKY	15	-1130.5553	2292.0267	0.0000	1.0000
SYMG	17	-1129.6148	2294.4021	0.0000	1.0000
SYMI	17	-1129.6331	2294.4386	0.0000	1.0000
K2PG	13	-1134.5846	2295.8611	0.0000	1.0000
K2PI	13	-1134.5932	2295.8784	0.0000	1.0000
SYMIG	18	-1129.5886	2296.4901	0.0000	1.0000
K2PIG	14	-1134.5659	2297.9319	0.0000	1.0000
F81	14	-1135.5971	2299.9942	0.0000	1.0000
JC69I	12	-1139.1310	2302.8541	0.0000	1.0000
JC69G	12	-1139.1334	2302.8589	0.0000	1.0000
JC69IG	13	-1139.1556	2305.0032	0.0000	1.0000
SYM	16	-1136.3555	2305.7511	0.0000	1.0000
K2P	12	-1141.6710	2307.9340	0.0000	1.0000
JC69	11	-1146.1323	2314.7646	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRG	20	-1115.6598	2271.3196	0.4104	0.4104
GTRI	20	-1115.6870	2271.3741	0.3994	0.8099
GTRIG	21	-1115.5696	2273.1391	0.1653	0.9751
HKYG	16	-1123.3747	2278.7494	0.0100	0.9851
HKYI	16	-1123.3873	2278.7747	0.0099	0.9950
HKYIG	17	-1123.3375	2280.6749	0.0038	0.9988
GTR	19	-1122.9377	2283.8754	0.0008	0.9996
F81I	15	-1128.5147	2287.0294	0.0002	0.9997
F81G	15	-1128.5155	2287.0310	0.0002	0.9999
F81IG	16	-1128.5171	2289.0341	0.0001	1.0000
HKY	15	-1130.5553	2291.1107	0.0000	1.0000
SYMG	17	-1129.6148	2293.2297	0.0000	1.0000
SYMI	17	-1129.6331	2293.2662	0.0000	1.0000
K2PG	13	-1134.5846	2295.1691	0.0000	1.0000
SYMIG	18	-1129.5886	2295.1773	0.0000	1.0000
K2PI	13	-1134.5932	2295.1864	0.0000	1.0000
K2PIG	14	-1134.5659	2297.1319	0.0000	1.0000
F81	14	-1135.5971	2299.1942	0.0000	1.0000
JC69I	12	-1139.1310	2302.2621	0.0000	1.0000
JC69G	12	-1139.1334	2302.2669	0.0000	1.0000
JC69IG	13	-1139.1556	2304.3112	0.0000	1.0000
SYM	16	-1136.3555	2304.7109	0.0000	1.0000
K2P	12	-1141.6710	2307.3420	0.0000	1.0000
JC69	11	-1146.1323	2314.2646	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	16	-1123.3747	2347.4145	0.3595	0.3595
HKYI	16	-1123.3873	2347.4398	0.3550	0.7145
K2PG	13	-1134.5846	2350.9595	0.0611	0.7756
K2PI	13	-1134.5932	2350.9768	0.0606	0.8361
F81I	15	-1128.5147	2351.4030	0.0489	0.8851
F81G	15	-1128.5155	2351.4046	0.0489	0.9340
HKYIG	17	-1123.3375	2353.6316	0.0161	0.9500
JC69I	12	-1139.1310	2353.7609	0.0151	0.9651
JC69G	12	-1139.1334	2353.7657	0.0150	0.9801
HKY	15	-1130.5553	2355.4842	0.0064	0.9864
GTRG	20	-1115.6598	2357.1510	0.0028	0.9892
GTRI	20	-1115.6870	2357.2054	0.0027	0.9919
K2PIG	14	-1134.5659	2357.2139	0.0027	0.9946
F81IG	16	-1128.5171	2357.6992	0.0021	0.9967
K2P	12	-1141.6710	2358.8408	0.0012	0.9979
F81	14	-1135.5971	2359.2761	0.0010	0.9988
JC69IG	13	-1139.1556	2360.1016	0.0006	0.9994
JC69	11	-1146.1323	2361.4719	0.0003	0.9998
GTRIG	21	-1115.5696	2363.2621	0.0001	0.9999
GTR	19	-1122.9377	2365.4152	0.0000	0.9999
SYMG	17	-1129.6148	2366.1863	0.0000	1.0000
SYMI	17	-1129.6331	2366.2229	0.0000	1.0000
SYMIG	18	-1129.5886	2372.4255	0.0000	1.0000
SYM	16	-1136.3555	2373.3760	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/chr6_3291.phylip = 1 - 540;
 Partition Dummy = 1:591-tomodel-phylip/chr6_3291.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/chr6_3291.phylip = 1 - 540;
 Partition Dummy = 1:591-tomodel-phylip/chr6_3291.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/chr6_3291.phylip = 1 - 540;
 Partition Dummy = 1:591-tomodel-phylip/chr6_3291.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

