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

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

Minimum AIC  model: GTRG
Minimum AICc model: GTRG
Minimum BIC  model: HKYG


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

Model	df	lnL		AIC		wAIC	cwAIC
GTRG	20	-3276.3944	6592.7888	0.4800	0.4800
GTRI	20	-3276.8341	6593.6682	0.3092	0.7892
GTRIG	21	-3276.4077	6594.8154	0.1743	0.9635
HKYG	16	-3283.6343	6599.2686	0.0188	0.9823
HKYI	16	-3284.1867	6600.3734	0.0108	0.9931
HKYIG	17	-3283.6361	6601.2721	0.0069	1.0000
GTR	19	-3309.2073	6656.4145	0.0000	1.0000
HKY	15	-3315.6868	6661.3736	0.0000	1.0000
K2PG	13	-3340.1210	6706.2419	0.0000	1.0000
SYMG	17	-3336.9969	6707.9938	0.0000	1.0000
K2PIG	14	-3340.0664	6708.1328	0.0000	1.0000
K2PI	13	-3341.1944	6708.3889	0.0000	1.0000
SYMIG	18	-3336.9406	6709.8812	0.0000	1.0000
SYMI	17	-3338.0307	6710.0615	0.0000	1.0000
K2P	12	-3375.2682	6774.5364	0.0000	1.0000
SYM	16	-3373.6086	6779.2172	0.0000	1.0000
F81G	15	-3380.4196	6790.8392	0.0000	1.0000
F81I	15	-3380.8740	6791.7481	0.0000	1.0000
F81IG	16	-3380.4436	6792.8871	0.0000	1.0000
F81	14	-3411.7964	6851.5929	0.0000	1.0000
JC69G	12	-3427.5809	6879.1618	0.0000	1.0000
JC69I	12	-3428.4591	6880.9183	0.0000	1.0000
JC69IG	13	-3427.5593	6881.1186	0.0000	1.0000
JC69	11	-3460.0543	6942.1087	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	16	-3283.6343	6679.5200	0.6221	0.6221
HKYI	16	-3284.1867	6680.6248	0.3581	0.9802
HKYIG	17	-3283.6361	6686.5393	0.0186	0.9988
GTRG	20	-3276.3944	6693.1031	0.0007	0.9995
GTRI	20	-3276.8341	6693.9824	0.0005	1.0000
GTRIG	21	-3276.4077	6700.1453	0.0000	1.0000
HKY	15	-3315.6868	6736.6093	0.0000	1.0000
GTR	19	-3309.2073	6751.7131	0.0000	1.0000
K2PG	13	-3340.1210	6771.4462	0.0000	1.0000
K2PI	13	-3341.1944	6773.5931	0.0000	1.0000
K2PIG	14	-3340.0664	6778.3528	0.0000	1.0000
SYMG	17	-3336.9969	6793.2610	0.0000	1.0000
SYMI	17	-3338.0307	6795.3286	0.0000	1.0000
SYMIG	18	-3336.9406	6800.1640	0.0000	1.0000
K2P	12	-3375.2682	6834.7249	0.0000	1.0000
SYM	16	-3373.6086	6859.4686	0.0000	1.0000
F81G	15	-3380.4196	6866.0749	0.0000	1.0000
F81I	15	-3380.8740	6866.9838	0.0000	1.0000
F81IG	16	-3380.4436	6873.1385	0.0000	1.0000
F81	14	-3411.7964	6921.8129	0.0000	1.0000
JC69G	12	-3427.5809	6939.3503	0.0000	1.0000
JC69I	12	-3428.4591	6941.1068	0.0000	1.0000
JC69IG	13	-3427.5593	6946.3229	0.0000	1.0000
JC69	11	-3460.0543	6997.2815	0.0000	1.0000



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

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

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


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

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


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

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

