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

Input data from file "591-tomodel-phylip/chr7_10552.phylip" (7 taxa, 946 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	-3058.4277	6156.8554	0.5537	0.5537
GTRI	20	-3059.2936	6158.5872	0.2329	0.7866
GTRIG	21	-3058.4407	6158.8814	0.2011	0.9877
HKYG	16	-3066.8072	6165.6143	0.0069	0.9946
HKYI	16	-3067.6697	6167.3395	0.0029	0.9975
HKYIG	17	-3066.8416	6167.6831	0.0025	1.0000
K2PG	13	-3096.0141	6218.0281	0.0000	1.0000
GTR	19	-3090.4264	6218.8527	0.0000	1.0000
K2PIG	14	-3096.0347	6220.0695	0.0000	1.0000
K2PI	13	-3097.0391	6220.0783	0.0000	1.0000
SYMG	17	-3093.6394	6221.2789	0.0000	1.0000
SYMI	17	-3094.5035	6223.0070	0.0000	1.0000
SYMIG	18	-3093.6638	6223.3275	0.0000	1.0000
HKY	15	-3097.8232	6225.6463	0.0000	1.0000
F81G	15	-3118.2477	6266.4953	0.0000	1.0000
F81I	15	-3118.8583	6267.7165	0.0000	1.0000
F81IG	16	-3118.2899	6268.5799	0.0000	1.0000
K2P	12	-3127.2939	6278.5879	0.0000	1.0000
SYM	16	-3124.6116	6281.2231	0.0000	1.0000
JC69G	12	-3142.3292	6308.6584	0.0000	1.0000
JC69I	12	-3143.1026	6310.2051	0.0000	1.0000
JC69IG	13	-3142.3748	6310.7497	0.0000	1.0000
F81	14	-3147.2067	6322.4134	0.0000	1.0000
JC69	11	-3171.7243	6365.4486	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	16	-3066.8072	6243.2502	0.6847	0.6847
HKYI	16	-3067.6697	6244.9754	0.2890	0.9737
HKYIG	17	-3066.8416	6250.1713	0.0215	0.9952
GTRG	20	-3058.4277	6253.9003	0.0033	0.9985
GTRI	20	-3059.2936	6255.6321	0.0014	0.9999
GTRIG	21	-3058.4407	6260.7785	0.0001	1.0000
K2PG	13	-3096.0141	6281.1073	0.0000	1.0000
K2PI	13	-3097.0391	6283.1574	0.0000	1.0000
K2PIG	14	-3096.0347	6288.0009	0.0000	1.0000
HKY	15	-3097.8232	6298.4299	0.0000	1.0000
SYMG	17	-3093.6394	6303.7670	0.0000	1.0000
SYMI	17	-3094.5035	6305.4951	0.0000	1.0000
SYMIG	18	-3093.6638	6310.6679	0.0000	1.0000
GTR	19	-3090.4264	6311.0453	0.0000	1.0000
K2P	12	-3127.2939	6336.8148	0.0000	1.0000
F81G	15	-3118.2477	6339.2790	0.0000	1.0000
F81I	15	-3118.8583	6340.5002	0.0000	1.0000
F81IG	16	-3118.2899	6346.2158	0.0000	1.0000
SYM	16	-3124.6116	6358.8590	0.0000	1.0000
JC69G	12	-3142.3292	6366.8853	0.0000	1.0000
JC69I	12	-3143.1026	6368.4320	0.0000	1.0000
JC69IG	13	-3142.3748	6373.8288	0.0000	1.0000
F81	14	-3147.2067	6390.3448	0.0000	1.0000
JC69	11	-3171.7243	6418.8233	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 591-tomodel-phylip/chr7_10552.phylip = 1 - 946;
 Partition Dummy = 1:591-tomodel-phylip/chr7_10552.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 591-tomodel-phylip/chr7_10552.phylip = 1 - 946;
 Partition Dummy = 1:591-tomodel-phylip/chr7_10552.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 591-tomodel-phylip/chr7_10552.phylip = 1 - 946;
 Partition Dummy = 1:591-tomodel-phylip/chr7_10552.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

