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

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

Minimum AIC  model: GTRG
Minimum AICc model: GTRG
Minimum BIC  model: K2PI


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

Model	df	lnL		AIC		wAIC	cwAIC
GTRG	20	-3289.9888	6619.9775	0.3273	0.3273
GTRI	20	-3290.0883	6620.1766	0.2963	0.6236
HKYI	16	-3294.7054	6621.4108	0.1598	0.7834
HKYG	16	-3294.7541	6621.5081	0.1523	0.9357
HKYIG	17	-3294.7642	6623.5284	0.0554	0.9912
GTRIG	21	-3292.7844	6627.5688	0.0074	0.9985
K2PI	13	-3303.3935	6632.7870	0.0005	0.9990
K2PG	13	-3303.4894	6632.9789	0.0005	0.9995
K2PIG	14	-3303.4986	6634.9972	0.0002	0.9997
SYMG	17	-3300.7235	6635.4471	0.0001	0.9999
SYMI	17	-3300.7764	6635.5528	0.0001	1.0000
SYMIG	18	-3303.6680	6643.3360	0.0000	1.0000
HKY	15	-3310.7175	6651.4350	0.0000	1.0000
GTR	19	-3307.1478	6652.2957	0.0000	1.0000
K2P	12	-3318.7314	6661.4629	0.0000	1.0000
SYM	16	-3316.3198	6664.6395	0.0000	1.0000
F81G	15	-3342.3216	6714.6433	0.0000	1.0000
F81I	15	-3344.8024	6719.6047	0.0000	1.0000
F81IG	16	-3344.8060	6721.6119	0.0000	1.0000
JC69G	12	-3349.1796	6722.3593	0.0000	1.0000
JC69I	12	-3351.6114	6727.2227	0.0000	1.0000
JC69IG	13	-3351.6145	6729.2290	0.0000	1.0000
F81	14	-3357.0493	6742.0987	0.0000	1.0000
JC69	11	-3363.7110	6749.4220	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
K2PI	13	-3303.3935	6698.1190	0.4451	0.4451
K2PG	13	-3303.4894	6698.3109	0.4044	0.8494
HKYI	16	-3294.7054	6701.8195	0.0700	0.9194
HKYG	16	-3294.7541	6701.9167	0.0666	0.9861
K2PIG	14	-3303.4986	6705.3547	0.0119	0.9980
HKYIG	17	-3294.7642	6708.9625	0.0020	1.0000
GTRG	20	-3289.9888	6720.4883	0.0000	1.0000
GTRI	20	-3290.0883	6720.6874	0.0000	1.0000
SYMG	17	-3300.7235	6720.8812	0.0000	1.0000
SYMI	17	-3300.7764	6720.9869	0.0000	1.0000
K2P	12	-3318.7314	6721.7693	0.0000	1.0000
HKY	15	-3310.7175	6726.8181	0.0000	1.0000
GTRIG	21	-3292.7844	6733.1051	0.0000	1.0000
SYMIG	18	-3303.6680	6733.7957	0.0000	1.0000
SYM	16	-3316.3198	6745.0481	0.0000	1.0000
GTR	19	-3307.1478	6747.7809	0.0000	1.0000
JC69G	12	-3349.1796	6782.6657	0.0000	1.0000
JC69I	12	-3351.6114	6787.5292	0.0000	1.0000
F81G	15	-3342.3216	6790.0264	0.0000	1.0000
JC69IG	13	-3351.6145	6794.5610	0.0000	1.0000
F81I	15	-3344.8024	6794.9878	0.0000	1.0000
F81IG	16	-3344.8060	6802.0206	0.0000	1.0000
JC69	11	-3363.7110	6804.7029	0.0000	1.0000
F81	14	-3357.0493	6812.4562	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/chr7_10505.phylip = 1 - 1125;
 Partition Dummy = 1:683-tomodel-phylip/chr7_10505.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/chr7_10505.phylip = 1 - 1125;
 Partition Dummy = 1:683-tomodel-phylip/chr7_10505.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 (K2PI)]

BEGIN MRBAYES;
 Charset 683-tomodel-phylip/chr7_10505.phylip = 1 - 1125;
 Partition Dummy = 1:683-tomodel-phylip/chr7_10505.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=2 rates=propinv;
 Prset applyto=(1) statefreqpr=Fixed(Equal) tratiopr=Beta(1.0,1.0) pinvarpr=Uniform(0.0,1.0);
END;


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

