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

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

Minimum AIC  model: GTRG
Minimum AICc model: GTRG
Minimum BIC  model: K2PG


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

Model	df	lnL		AIC		wAIC	cwAIC
GTRG	20	-3805.1288	7650.2576	0.6791	0.6791
GTRIG	21	-3805.0809	7652.1618	0.2621	0.9411
GTRI	20	-3807.7982	7655.5965	0.0471	0.9882
HKYG	16	-3813.6900	7659.3801	0.0071	0.9953
HKYIG	17	-3813.6621	7661.3242	0.0027	0.9980
SYMG	17	-3814.6355	7663.2710	0.0010	0.9990
HKYI	16	-3816.3454	7664.6907	0.0005	0.9995
SYMIG	18	-3814.6066	7665.2133	0.0004	0.9999
SYMI	17	-3817.3616	7668.7232	0.0001	0.9999
K2PG	13	-3821.6741	7669.3482	0.0000	1.0000
K2PIG	14	-3821.6500	7671.3000	0.0000	1.0000
K2PI	13	-3824.2767	7674.5534	0.0000	1.0000
GTR	19	-3846.4512	7730.9025	0.0000	1.0000
HKY	15	-3853.5142	7737.0283	0.0000	1.0000
SYM	16	-3854.4236	7740.8472	0.0000	1.0000
K2P	12	-3861.3431	7746.6862	0.0000	1.0000
F81G	15	-3900.2827	7830.5654	0.0000	1.0000
F81IG	16	-3900.2743	7832.5486	0.0000	1.0000
F81I	15	-3901.7817	7833.5635	0.0000	1.0000
JC69G	12	-3910.8074	7845.6147	0.0000	1.0000
JC69IG	13	-3910.8062	7847.6124	0.0000	1.0000
JC69I	12	-3912.3832	7848.7664	0.0000	1.0000
F81	14	-3936.4139	7900.8277	0.0000	1.0000
JC69	11	-3946.1516	7914.3032	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
K2PG	13	-3821.6741	7734.0164	0.8337	0.8337
HKYG	16	-3813.6900	7738.9717	0.0700	0.9037
K2PI	13	-3824.2767	7739.2216	0.0618	0.9655
K2PIG	14	-3821.6500	7740.9427	0.0261	0.9916
HKYI	16	-3816.3454	7744.2824	0.0049	0.9965
HKYIG	17	-3813.6621	7745.8904	0.0022	0.9987
SYMG	17	-3814.6355	7747.8372	0.0008	0.9996
GTRG	20	-3805.1288	7749.7472	0.0003	0.9999
SYMI	17	-3817.3616	7753.2893	0.0001	0.9999
SYMIG	18	-3814.6066	7754.7539	0.0000	1.0000
GTRI	20	-3807.7982	7755.0861	0.0000	1.0000
GTRIG	21	-3805.0809	7756.6258	0.0000	1.0000
K2P	12	-3861.3431	7806.3799	0.0000	1.0000
HKY	15	-3853.5142	7811.6455	0.0000	1.0000
SYM	16	-3854.4236	7820.4389	0.0000	1.0000
GTR	19	-3846.4512	7825.4176	0.0000	1.0000
F81G	15	-3900.2827	7905.1826	0.0000	1.0000
JC69G	12	-3910.8074	7905.3084	0.0000	1.0000
F81I	15	-3901.7817	7908.1807	0.0000	1.0000
JC69I	12	-3912.3832	7908.4602	0.0000	1.0000
F81IG	16	-3900.2743	7912.1403	0.0000	1.0000
JC69IG	13	-3910.8062	7912.2807	0.0000	1.0000
JC69	11	-3946.1516	7969.0225	0.0000	1.0000
F81	14	-3936.4139	7970.4704	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/chr1_28671.phylip = 1 - 1069;
 Partition Dummy = 1:683-tomodel-phylip/chr1_28671.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/chr1_28671.phylip = 1 - 1069;
 Partition Dummy = 1:683-tomodel-phylip/chr1_28671.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 (K2PG)]

BEGIN MRBAYES;
 Charset 683-tomodel-phylip/chr1_28671.phylip = 1 - 1069;
 Partition Dummy = 1:683-tomodel-phylip/chr1_28671.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=2 rates=gamma;
 Prset applyto=(1) statefreqpr=Fixed(Equal) tratiopr=Beta(1.0,1.0) shapepr=Uniform(0.1,50.0);
END;


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

