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

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

Minimum AIC  model: GTRI
Minimum AICc model: GTRI
Minimum BIC  model: HKYI


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

Model	df	lnL		AIC		wAIC	cwAIC
GTRI	20	-2979.7663	5999.5327	0.7047	0.7047
GTRG	20	-2981.2635	6002.5270	0.1577	0.8623
GTRIG	21	-2980.5011	6003.0022	0.1243	0.9867
HKYI	16	-2988.2412	6008.4824	0.0080	0.9947
HKYIG	17	-2988.2559	6010.5119	0.0029	0.9976
HKYG	16	-2989.4474	6010.8948	0.0024	1.0000
K2PI	13	-3016.8505	6059.7010	0.0000	1.0000
SYMI	17	-3013.5926	6061.1853	0.0000	1.0000
K2PG	13	-3018.2616	6062.5232	0.0000	1.0000
SYMIG	18	-3014.4996	6064.9993	0.0000	1.0000
K2PIG	14	-3018.5276	6065.0553	0.0000	1.0000
SYMG	17	-3015.5630	6065.1259	0.0000	1.0000
F81I	15	-3028.8468	6087.6937	0.0000	1.0000
F81IG	16	-3028.1036	6088.2072	0.0000	1.0000
F81G	15	-3030.0136	6090.0271	0.0000	1.0000
GTR	19	-3032.9719	6103.9438	0.0000	1.0000
HKY	15	-3039.0041	6108.0082	0.0000	1.0000
JC69I	12	-3052.4738	6128.9476	0.0000	1.0000
JC69IG	13	-3051.8197	6129.6393	0.0000	1.0000
JC69G	12	-3053.7570	6131.5139	0.0000	1.0000
K2P	12	-3065.6554	6155.3108	0.0000	1.0000
SYM	16	-3062.7366	6157.4733	0.0000	1.0000
F81	14	-3076.1911	6180.3822	0.0000	1.0000
JC69	11	-3098.7649	6219.5298	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYI	16	-2988.2412	6086.0335	0.7473	0.7473
HKYG	16	-2989.4474	6088.4459	0.2237	0.9710
HKYIG	17	-2988.2559	6092.9099	0.0240	0.9950
GTRI	20	-2979.7663	6096.4715	0.0040	0.9990
GTRG	20	-2981.2635	6099.4659	0.0009	0.9999
GTRIG	21	-2980.5011	6104.7880	0.0001	1.0000
K2PI	13	-3016.8505	6122.7113	0.0000	1.0000
K2PG	13	-3018.2616	6125.5335	0.0000	1.0000
K2PIG	14	-3018.5276	6132.9125	0.0000	1.0000
SYMI	17	-3013.5926	6143.5833	0.0000	1.0000
SYMG	17	-3015.5630	6147.5239	0.0000	1.0000
SYMIG	18	-3014.4996	6152.2442	0.0000	1.0000
F81I	15	-3028.8468	6160.3978	0.0000	1.0000
F81G	15	-3030.0136	6162.7312	0.0000	1.0000
F81IG	16	-3028.1036	6165.7583	0.0000	1.0000
HKY	15	-3039.0041	6180.7123	0.0000	1.0000
JC69I	12	-3052.4738	6187.1109	0.0000	1.0000
JC69G	12	-3053.7570	6189.6773	0.0000	1.0000
JC69IG	13	-3051.8197	6192.6496	0.0000	1.0000
GTR	19	-3032.9719	6196.0357	0.0000	1.0000
K2P	12	-3065.6554	6213.4741	0.0000	1.0000
SYM	16	-3062.7366	6235.0244	0.0000	1.0000
F81	14	-3076.1911	6248.2394	0.0000	1.0000
JC69	11	-3098.7649	6272.8462	0.0000	1.0000



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

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

BEGIN MRBAYES;
 Charset 591-tomodel-phylip/chr6_7113.phylip = 1 - 941;
 Partition Dummy = 1:591-tomodel-phylip/chr6_7113.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=6 rates=propinv;
 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) pinvarpr=Uniform(0.0,1.0);
END;


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

BEGIN MRBAYES;
 Charset 591-tomodel-phylip/chr6_7113.phylip = 1 - 941;
 Partition Dummy = 1:591-tomodel-phylip/chr6_7113.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=6 rates=propinv;
 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) pinvarpr=Uniform(0.0,1.0);
END;


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

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


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

