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

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

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


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
GTRI	20	-2599.3878	5240.1175	0.5181	0.5181
GTRG	20	-2600.0726	5241.4871	0.2612	0.7794
GTRIG	21	-2599.2226	5241.9235	0.2100	0.9894
HKYI	16	-2608.1310	5249.1255	0.0057	0.9951
HKYG	16	-2608.9337	5250.7309	0.0026	0.9977
HKYIG	17	-2607.9872	5250.9474	0.0023	1.0000
GTR	19	-2623.5204	5286.2530	0.0000	1.0000
HKY	15	-2632.9606	5296.6819	0.0000	1.0000
SYMG	17	-2666.3632	5367.6993	0.0000	1.0000
SYMI	17	-2666.5547	5368.0824	0.0000	1.0000
K2PI	13	-2671.2018	5368.9787	0.0000	1.0000
SYMIG	18	-2666.0493	5369.1879	0.0000	1.0000
K2PG	13	-2671.6773	5369.9296	0.0000	1.0000
K2PIG	14	-2671.0096	5370.6838	0.0000	1.0000
F81I	15	-2670.4838	5371.7282	0.0000	1.0000
F81G	15	-2670.9726	5372.7059	0.0000	1.0000
F81IG	16	-2670.3278	5373.5191	0.0000	1.0000
F81	14	-2693.1415	5414.9476	0.0000	1.0000
SYM	16	-2692.1927	5417.2489	0.0000	1.0000
K2P	12	-2698.9026	5422.2974	0.0000	1.0000
JC69G	12	-2721.2742	5467.0406	0.0000	1.0000
JC69I	12	-2721.3329	5467.1579	0.0000	1.0000
JC69IG	13	-2720.9277	5468.4305	0.0000	1.0000
JC69	11	-2746.4801	5515.3760	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRI	20	-2599.3878	5238.7756	0.5116	0.5116
GTRG	20	-2600.0726	5240.1452	0.2580	0.7696
GTRIG	21	-2599.2226	5240.4451	0.2220	0.9917
HKYI	16	-2608.1310	5248.2620	0.0045	0.9961
HKYG	16	-2608.9337	5249.8674	0.0020	0.9981
HKYIG	17	-2607.9872	5249.9744	0.0019	1.0000
GTR	19	-2623.5204	5285.0409	0.0000	1.0000
HKY	15	-2632.9606	5295.9212	0.0000	1.0000
SYMG	17	-2666.3632	5366.7263	0.0000	1.0000
SYMI	17	-2666.5547	5367.1094	0.0000	1.0000
SYMIG	18	-2666.0493	5368.0987	0.0000	1.0000
K2PI	13	-2671.2018	5368.4037	0.0000	1.0000
K2PG	13	-2671.6773	5369.3545	0.0000	1.0000
K2PIG	14	-2671.0096	5370.0193	0.0000	1.0000
F81I	15	-2670.4838	5370.9675	0.0000	1.0000
F81G	15	-2670.9726	5371.9452	0.0000	1.0000
F81IG	16	-2670.3278	5372.6556	0.0000	1.0000
F81	14	-2693.1415	5414.2831	0.0000	1.0000
SYM	16	-2692.1927	5416.3854	0.0000	1.0000
K2P	12	-2698.9026	5421.8053	0.0000	1.0000
JC69G	12	-2721.2742	5466.5485	0.0000	1.0000
JC69I	12	-2721.3329	5466.6658	0.0000	1.0000
JC69IG	13	-2720.9277	5467.8555	0.0000	1.0000
JC69	11	-2746.4801	5514.9603	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYI	16	-2608.1310	5319.8195	0.6593	0.6593
HKYG	16	-2608.9337	5321.4250	0.2954	0.9548
HKYIG	17	-2607.9872	5326.0043	0.0299	0.9847
GTRI	20	-2599.3878	5328.2225	0.0099	0.9946
GTRG	20	-2600.0726	5329.5922	0.0050	0.9995
GTRIG	21	-2599.2226	5334.3644	0.0005	1.0000
HKY	15	-2632.9606	5363.0064	0.0000	1.0000
GTR	19	-2623.5204	5370.0155	0.0000	1.0000
K2PI	13	-2671.2018	5426.5442	0.0000	1.0000
K2PG	13	-2671.6773	5427.4950	0.0000	1.0000
K2PIG	14	-2671.0096	5432.6321	0.0000	1.0000
F81I	15	-2670.4838	5438.0527	0.0000	1.0000
F81G	15	-2670.9726	5439.0304	0.0000	1.0000
SYMG	17	-2666.3632	5442.7562	0.0000	1.0000
SYMI	17	-2666.5547	5443.1393	0.0000	1.0000
F81IG	16	-2670.3278	5444.2131	0.0000	1.0000
SYMIG	18	-2666.0493	5448.6009	0.0000	1.0000
K2P	12	-2698.9026	5475.4734	0.0000	1.0000
F81	14	-2693.1415	5476.8959	0.0000	1.0000
SYM	16	-2692.1927	5487.9429	0.0000	1.0000
JC69G	12	-2721.2742	5520.2167	0.0000	1.0000
JC69I	12	-2721.3329	5520.3340	0.0000	1.0000
JC69IG	13	-2720.9277	5525.9960	0.0000	1.0000
JC69	11	-2746.4801	5564.1561	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/chr8_851.phylip = 1 - 647;
 Partition Dummy = 1:591-tomodel-phylip/chr8_851.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/chr8_851.phylip = 1 - 647;
 Partition Dummy = 1:591-tomodel-phylip/chr8_851.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/chr8_851.phylip = 1 - 647;
 Partition Dummy = 1:591-tomodel-phylip/chr8_851.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

