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

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

Minimum AIC  model: HKYG
Minimum AICc model: HKYG
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
HKYG	16	-1437.3956	2908.0593	0.3181	0.3181
HKYI	16	-1437.5825	2908.4331	0.2639	0.5821
GTRG	20	-1433.8794	2909.7353	0.1376	0.7197
GTRI	20	-1434.0264	2910.0293	0.1188	0.8385
HKYIG	17	-1437.4079	2910.2456	0.1066	0.9451
GTRIG	21	-1433.8792	2911.9376	0.0458	0.9909
HKY	15	-1442.3472	2915.8107	0.0066	0.9975
GTR	19	-1438.9680	2917.7201	0.0025	1.0000
K2PG	13	-1455.1696	2937.1818	0.0000	1.0000
K2PI	13	-1455.4471	2937.7367	0.0000	1.0000
K2PIG	14	-1455.1768	2939.3280	0.0000	1.0000
SYMG	17	-1452.9921	2941.4141	0.0000	1.0000
SYMI	17	-1453.2056	2941.8412	0.0000	1.0000
SYMIG	18	-1452.9990	2943.5998	0.0000	1.0000
K2P	12	-1460.5419	2945.8043	0.0000	1.0000
SYM	16	-1458.4688	2950.2057	0.0000	1.0000
F81G	15	-1472.6311	2976.3785	0.0000	1.0000
F81I	15	-1472.8531	2976.8225	0.0000	1.0000
F81IG	16	-1472.6375	2978.5431	0.0000	1.0000
F81	14	-1477.0966	2983.1678	0.0000	1.0000
JC69G	12	-1487.4908	2999.7021	0.0000	1.0000
JC69I	12	-1487.7647	3000.2500	0.0000	1.0000
JC69IG	13	-1487.4990	3001.8406	0.0000	1.0000
JC69	11	-1492.2705	3007.1494	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
HKYG	16	-1437.3956	2906.7912	0.2779	0.2779
HKYI	16	-1437.5825	2907.1650	0.2306	0.5085
GTRG	20	-1433.8794	2907.7588	0.1713	0.6798
GTRI	20	-1434.0264	2908.0529	0.1479	0.8277
HKYIG	17	-1437.4079	2908.8157	0.1010	0.9287
GTRIG	21	-1433.8792	2909.7584	0.0630	0.9918
HKY	15	-1442.3472	2914.6944	0.0053	0.9971
GTR	19	-1438.9680	2915.9361	0.0029	1.0000
K2PG	13	-1455.1696	2936.3392	0.0000	1.0000
K2PI	13	-1455.4471	2936.8941	0.0000	1.0000
K2PIG	14	-1455.1768	2938.3535	0.0000	1.0000
SYMG	17	-1452.9921	2939.9842	0.0000	1.0000
SYMI	17	-1453.2056	2940.4113	0.0000	1.0000
SYMIG	18	-1452.9990	2941.9979	0.0000	1.0000
K2P	12	-1460.5419	2945.0837	0.0000	1.0000
SYM	16	-1458.4688	2948.9377	0.0000	1.0000
F81G	15	-1472.6311	2975.2622	0.0000	1.0000
F81I	15	-1472.8531	2975.7062	0.0000	1.0000
F81IG	16	-1472.6375	2977.2751	0.0000	1.0000
F81	14	-1477.0966	2982.1933	0.0000	1.0000
JC69G	12	-1487.4908	2998.9816	0.0000	1.0000
JC69I	12	-1487.7647	2999.5294	0.0000	1.0000
JC69IG	13	-1487.4990	3000.9980	0.0000	1.0000
JC69	11	-1492.2705	3006.5411	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	16	-1437.3956	2972.3963	0.4935	0.4935
HKYI	16	-1437.5825	2972.7701	0.4094	0.9029
HKY	15	-1442.3472	2976.1992	0.0737	0.9766
HKYIG	17	-1437.4079	2978.5212	0.0231	0.9997
K2PG	13	-1455.1696	2989.6434	0.0001	0.9997
GTRG	20	-1433.8794	2989.7652	0.0001	0.9998
GTRI	20	-1434.0264	2990.0592	0.0001	0.9999
K2PI	13	-1455.4471	2990.1983	0.0001	1.0000
GTR	19	-1438.9680	2993.8421	0.0000	1.0000
K2P	12	-1460.5419	2994.2875	0.0000	1.0000
K2PIG	14	-1455.1768	2995.7580	0.0000	1.0000
GTRIG	21	-1433.8792	2995.8651	0.0000	1.0000
SYMG	17	-1452.9921	3009.6896	0.0000	1.0000
SYMI	17	-1453.2056	3010.1167	0.0000	1.0000
SYM	16	-1458.4688	3014.5428	0.0000	1.0000
SYMIG	18	-1452.9990	3015.8037	0.0000	1.0000
F81G	15	-1472.6311	3036.7670	0.0000	1.0000
F81I	15	-1472.8531	3037.2110	0.0000	1.0000
F81	14	-1477.0966	3039.5977	0.0000	1.0000
F81IG	16	-1472.6375	3042.8802	0.0000	1.0000
JC69G	12	-1487.4908	3048.1854	0.0000	1.0000
JC69I	12	-1487.7647	3048.7332	0.0000	1.0000
JC69	11	-1492.2705	3051.6446	0.0000	1.0000
JC69IG	13	-1487.4990	3054.3022	0.0000	1.0000



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

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

BEGIN MRBAYES;
 Charset 591-tomodel-phylip/chr2_24734.phylip = 1 - 446;
 Partition Dummy = 1:591-tomodel-phylip/chr2_24734.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;


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

BEGIN MRBAYES;
 Charset 591-tomodel-phylip/chr2_24734.phylip = 1 - 446;
 Partition Dummy = 1:591-tomodel-phylip/chr2_24734.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;


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

BEGIN MRBAYES;
 Charset 591-tomodel-phylip/chr2_24734.phylip = 1 - 446;
 Partition Dummy = 1:591-tomodel-phylip/chr2_24734.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

