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

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

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


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

Model	df	lnL		AIC		wAIC	cwAIC
HKYI	16	-2862.7241	5757.4482	0.4220	0.4220
HKYG	16	-2863.0202	5758.0404	0.3138	0.7358
HKYIG	17	-2863.0203	5760.0406	0.1154	0.8512
GTRI	20	-2860.4799	5760.9599	0.0729	0.9241
GTRG	20	-2860.7532	5761.5065	0.0555	0.9796
GTRIG	21	-2860.7532	5763.5064	0.0204	1.0000
K2PI	13	-2905.8328	5837.6655	0.0000	1.0000
K2PG	13	-2906.1157	5838.2314	0.0000	1.0000
K2PIG	14	-2906.1158	5840.2316	0.0000	1.0000
HKY	15	-2905.6420	5841.2841	0.0000	1.0000
SYMI	17	-2904.4555	5842.9110	0.0000	1.0000
SYMG	17	-2904.7373	5843.4746	0.0000	1.0000
GTR	19	-2903.1290	5844.2580	0.0000	1.0000
SYMIG	18	-2904.7373	5845.4746	0.0000	1.0000
F81I	15	-2908.2584	5846.5169	0.0000	1.0000
F81G	15	-2909.8386	5849.6771	0.0000	1.0000
F81IG	16	-2909.8386	5851.6773	0.0000	1.0000
JC69I	12	-2948.2896	5920.5792	0.0000	1.0000
K2P	12	-2948.7933	5921.5867	0.0000	1.0000
JC69G	12	-2950.0427	5924.0855	0.0000	1.0000
JC69IG	13	-2950.0428	5926.0855	0.0000	1.0000
F81	14	-2949.2678	5926.5356	0.0000	1.0000
SYM	16	-2947.4908	5926.9817	0.0000	1.0000
JC69	11	-2990.0205	6002.0411	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYI	16	-2862.7241	5837.3511	0.5662	0.5662
HKYG	16	-2863.0202	5837.9433	0.4211	0.9872
HKYIG	17	-2863.0203	5844.9374	0.0128	1.0000
GTRI	20	-2860.4799	5860.8386	0.0000	1.0000
GTRG	20	-2860.7532	5861.3851	0.0000	1.0000
GTRIG	21	-2860.7532	5868.3790	0.0000	1.0000
K2PI	13	-2905.8328	5902.5867	0.0000	1.0000
K2PG	13	-2906.1157	5903.1525	0.0000	1.0000
K2PIG	14	-2906.1158	5910.1467	0.0000	1.0000
HKY	15	-2905.6420	5916.1931	0.0000	1.0000
F81I	15	-2908.2584	5921.4259	0.0000	1.0000
F81G	15	-2909.8386	5924.5861	0.0000	1.0000
SYMI	17	-2904.4555	5927.8079	0.0000	1.0000
SYMG	17	-2904.7373	5928.3714	0.0000	1.0000
F81IG	16	-2909.8386	5931.5802	0.0000	1.0000
SYMIG	18	-2904.7373	5935.3654	0.0000	1.0000
GTR	19	-2903.1290	5939.1427	0.0000	1.0000
JC69I	12	-2948.2896	5980.5064	0.0000	1.0000
K2P	12	-2948.7933	5981.5139	0.0000	1.0000
JC69G	12	-2950.0427	5984.0127	0.0000	1.0000
JC69IG	13	-2950.0428	5991.0067	0.0000	1.0000
F81	14	-2949.2678	5996.4507	0.0000	1.0000
SYM	16	-2947.4908	6006.8846	0.0000	1.0000
JC69	11	-2990.0205	6056.9743	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 591-tomodel-phylip/chr5_10201.phylip = 1 - 1090;
 Partition Dummy = 1:591-tomodel-phylip/chr5_10201.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 (HKYI)]

BEGIN MRBAYES;
 Charset 591-tomodel-phylip/chr5_10201.phylip = 1 - 1090;
 Partition Dummy = 1:591-tomodel-phylip/chr5_10201.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 BIC model (HKYI)]

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

