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

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

Minimum AIC  model: GTRIG
Minimum AICc model: GTRIG
Minimum BIC  model: HKYG


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

Model	df	lnL		AIC		wAIC	cwAIC
GTRIG	21	-2741.3284	5524.6569	0.6636	0.6636
GTRG	20	-2743.5999	5527.1998	0.1861	0.8496
GTRI	20	-2744.1424	5528.2848	0.1082	0.9578
HKYIG	17	-2748.6329	5531.2658	0.0244	0.9822
HKYG	16	-2750.4704	5532.9408	0.0105	0.9927
HKYI	16	-2750.8415	5533.6830	0.0073	1.0000
K2PIG	14	-2789.2797	5606.5594	0.0000	1.0000
K2PG	13	-2790.8498	5607.6995	0.0000	1.0000
K2PI	13	-2791.1528	5608.3055	0.0000	1.0000
SYMIG	18	-2786.7345	5609.4689	0.0000	1.0000
SYMG	17	-2788.2205	5610.4410	0.0000	1.0000
SYMI	17	-2788.5142	5611.0283	0.0000	1.0000
GTR	19	-2792.6433	5623.2865	0.0000	1.0000
HKY	15	-2799.4262	5628.8524	0.0000	1.0000
F81IG	16	-2802.2886	5636.5772	0.0000	1.0000
F81G	15	-2803.3233	5636.6466	0.0000	1.0000
F81I	15	-2803.7036	5637.4071	0.0000	1.0000
JC69G	12	-2837.9806	5699.9611	0.0000	1.0000
JC69IG	13	-2837.0122	5700.0243	0.0000	1.0000
K2P	12	-2838.2209	5700.4418	0.0000	1.0000
JC69I	12	-2838.3497	5700.6993	0.0000	1.0000
SYM	16	-2835.1867	5702.3734	0.0000	1.0000
F81	14	-2846.4422	5720.8844	0.0000	1.0000
JC69	11	-2880.9231	5783.8462	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	16	-2750.4704	5613.5333	0.5326	0.5326
HKYI	16	-2750.8415	5614.2754	0.3675	0.9001
HKYIG	17	-2748.6329	5616.8952	0.0992	0.9993
GTRG	20	-2743.5999	5627.9404	0.0004	0.9997
GTRI	20	-2744.1424	5629.0254	0.0002	0.9999
GTRIG	21	-2741.3284	5630.4345	0.0001	1.0000
K2PG	13	-2790.8498	5673.1809	0.0000	1.0000
K2PI	13	-2791.1528	5673.7869	0.0000	1.0000
K2PIG	14	-2789.2797	5677.0778	0.0000	1.0000
SYMG	17	-2788.2205	5696.0704	0.0000	1.0000
SYMI	17	-2788.5142	5696.6578	0.0000	1.0000
SYMIG	18	-2786.7345	5700.1354	0.0000	1.0000
HKY	15	-2799.4262	5704.4079	0.0000	1.0000
F81G	15	-2803.3233	5712.2020	0.0000	1.0000
F81I	15	-2803.7036	5712.9625	0.0000	1.0000
F81IG	16	-2802.2886	5717.1696	0.0000	1.0000
GTR	19	-2792.6433	5718.9900	0.0000	1.0000
JC69G	12	-2837.9806	5760.4054	0.0000	1.0000
K2P	12	-2838.2209	5760.8862	0.0000	1.0000
JC69I	12	-2838.3497	5761.1436	0.0000	1.0000
JC69IG	13	-2837.0122	5765.5057	0.0000	1.0000
SYM	16	-2835.1867	5782.9659	0.0000	1.0000
F81	14	-2846.4422	5791.4028	0.0000	1.0000
JC69	11	-2880.9231	5839.2535	0.0000	1.0000



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

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

BEGIN MRBAYES;
 Charset 683-tomodel-phylip/chr3_576.phylip = 1 - 1138;
 Partition Dummy = 1:683-tomodel-phylip/chr3_576.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=6 rates=invgamma;
 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) shapepr=Uniform(0.1,50.0) pinvarpr=Uniform(0.0,1.0);
END;


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

BEGIN MRBAYES;
 Charset 683-tomodel-phylip/chr3_576.phylip = 1 - 1138;
 Partition Dummy = 1:683-tomodel-phylip/chr3_576.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=6 rates=invgamma;
 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) shapepr=Uniform(0.1,50.0) pinvarpr=Uniform(0.0,1.0);
END;


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

BEGIN MRBAYES;
 Charset 683-tomodel-phylip/chr3_576.phylip = 1 - 1138;
 Partition Dummy = 1:683-tomodel-phylip/chr3_576.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

