
Output from mraic.pl version 1.4.4 by Johan Nylander
-------------------------------------------------------------
16:10 04/04/2011

Input data from file "683-tomodel-phylip/chr1_5316.phylip" (7 taxa, 635 characters)

Minimum AIC  model: GTR
Minimum AICc model: HKY
Minimum BIC  model: HKY


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
HKY	15	-1006.8374	2044.4503	0.2882	0.2882
GTR	19	-1002.7132	2044.6622	0.2592	0.5473
HKYI	16	-1006.8373	2046.5549	0.1006	0.6479
HKYG	16	-1006.8384	2046.5570	0.1005	0.7484
GTRG	20	-1002.7125	2046.7931	0.0893	0.8378
GTRI	20	-1002.7131	2046.7944	0.0893	0.9270
HKYIG	17	-1006.8397	2048.6714	0.0349	0.9619
GTRIG	21	-1002.7156	2048.9386	0.0305	0.9925
F81	14	-1012.1312	2052.9397	0.0041	0.9966
F81I	15	-1012.1311	2055.0376	0.0014	0.9980
F81G	15	-1012.1317	2055.0388	0.0014	0.9995
F81IG	16	-1012.1337	2057.1476	0.0005	1.0000
SYM	16	-1029.9640	2092.8082	0.0000	1.0000
K2P	12	-1034.6793	2093.8601	0.0000	1.0000
SYMI	17	-1029.9644	2094.9206	0.0000	1.0000
SYMG	17	-1029.9650	2094.9219	0.0000	1.0000
K2PI	13	-1034.6792	2095.9445	0.0000	1.0000
K2PG	13	-1034.6798	2095.9458	0.0000	1.0000
SYMIG	18	-1029.9674	2097.0451	0.0000	1.0000
K2PIG	14	-1034.6819	2098.0412	0.0000	1.0000
JC69	11	-1039.5375	2101.4988	0.0000	1.0000
JC69I	12	-1039.5373	2103.5763	0.0000	1.0000
JC69G	12	-1039.5379	2103.5774	0.0000	1.0000
JC69IG	13	-1039.5400	2105.6662	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTR	19	-1002.7132	2043.4264	0.2820	0.2820
HKY	15	-1006.8374	2043.6748	0.2491	0.5311
GTRG	20	-1002.7125	2045.4250	0.1038	0.6350
GTRI	20	-1002.7131	2045.4263	0.1038	0.7387
HKYI	16	-1006.8373	2045.6747	0.0916	0.8304
HKYG	16	-1006.8384	2045.6767	0.0916	0.9219
GTRIG	21	-1002.7156	2047.4313	0.0381	0.9600
HKYIG	17	-1006.8397	2047.6795	0.0336	0.9936
F81	14	-1012.1312	2052.2623	0.0034	0.9970
F81I	15	-1012.1311	2054.2621	0.0013	0.9983
F81G	15	-1012.1317	2054.2633	0.0013	0.9995
F81IG	16	-1012.1337	2056.2674	0.0005	1.0000
SYM	16	-1029.9640	2091.9279	0.0000	1.0000
K2P	12	-1034.6793	2093.3585	0.0000	1.0000
SYMI	17	-1029.9644	2093.9287	0.0000	1.0000
SYMG	17	-1029.9650	2093.9300	0.0000	1.0000
K2PI	13	-1034.6792	2095.3584	0.0000	1.0000
K2PG	13	-1034.6798	2095.3596	0.0000	1.0000
SYMIG	18	-1029.9674	2095.9347	0.0000	1.0000
K2PIG	14	-1034.6819	2097.3638	0.0000	1.0000
JC69	11	-1039.5375	2101.0750	0.0000	1.0000
JC69I	12	-1039.5373	2103.0747	0.0000	1.0000
JC69G	12	-1039.5379	2103.0758	0.0000	1.0000
JC69IG	13	-1039.5400	2105.0800	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKY	15	-1006.8374	2110.4792	0.8211	0.8211
F81	14	-1012.1312	2114.6130	0.1039	0.9250
HKYI	16	-1006.8373	2116.9327	0.0326	0.9576
HKYG	16	-1006.8384	2116.9347	0.0326	0.9902
F81I	15	-1012.1311	2121.0665	0.0041	0.9943
F81G	15	-1012.1317	2121.0677	0.0041	0.9984
HKYIG	17	-1006.8397	2123.3911	0.0013	0.9997
F81IG	16	-1012.1337	2127.5254	0.0002	0.9999
GTR	19	-1002.7132	2128.0453	0.0001	1.0000
GTRG	20	-1002.7125	2134.4975	0.0000	1.0000
GTRI	20	-1002.7131	2134.4988	0.0000	1.0000
GTRIG	21	-1002.7156	2140.9574	0.0000	1.0000
K2P	12	-1034.6793	2146.8020	0.0000	1.0000
JC69	11	-1039.5375	2150.0649	0.0000	1.0000
K2PI	13	-1034.6792	2153.2555	0.0000	1.0000
K2PG	13	-1034.6798	2153.2567	0.0000	1.0000
JC69I	12	-1039.5373	2156.5182	0.0000	1.0000
JC69G	12	-1039.5379	2156.5193	0.0000	1.0000
K2PIG	14	-1034.6819	2159.7145	0.0000	1.0000
JC69IG	13	-1039.5400	2162.9771	0.0000	1.0000
SYM	16	-1029.9640	2163.1859	0.0000	1.0000
SYMI	17	-1029.9644	2169.6403	0.0000	1.0000
SYMG	17	-1029.9650	2169.6417	0.0000	1.0000
SYMIG	18	-1029.9674	2176.1000	0.0000	1.0000



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

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

BEGIN MRBAYES;
 Charset 683-tomodel-phylip/chr1_5316.phylip = 1 - 635;
 Partition Dummy = 1:683-tomodel-phylip/chr1_5316.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=6 rates=equal;
 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);
END;


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

BEGIN MRBAYES;
 Charset 683-tomodel-phylip/chr1_5316.phylip = 1 - 635;
 Partition Dummy = 1:683-tomodel-phylip/chr1_5316.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=2 rates=equal;
 Prset applyto=(1) statefreqpr=Dirichlet(1.0,1.0,1.0,1.0) tratiopr=Beta(1.0,1.0);
END;


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

BEGIN MRBAYES;
 Charset 683-tomodel-phylip/chr1_5316.phylip = 1 - 635;
 Partition Dummy = 1:683-tomodel-phylip/chr1_5316.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=2 rates=equal;
 Prset applyto=(1) statefreqpr=Dirichlet(1.0,1.0,1.0,1.0) tratiopr=Beta(1.0,1.0);
END;


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

