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

Input data from file "591-tomodel-phylip/chr5_3310.phylip" (7 taxa, 921 characters)

Minimum AIC  model: GTRI
Minimum AICc model: GTRI
Minimum BIC  model: HKYG


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

Model	df	lnL		AIC		wAIC	cwAIC
GTRI	20	-2306.9900	4653.9800	0.4459	0.4459
GTRG	20	-2307.0877	4654.1753	0.4044	0.8504
GTRIG	21	-2307.0876	4656.1753	0.1488	0.9992
HKYG	16	-2318.0482	4668.0963	0.0004	0.9995
HKYI	16	-2318.2624	4668.5248	0.0003	0.9998
HKYIG	17	-2317.9542	4669.9084	0.0002	1.0000
GTR	19	-2331.7162	4701.4324	0.0000	1.0000
HKY	15	-2342.9594	4715.9187	0.0000	1.0000
SYMG	17	-2345.1859	4724.3718	0.0000	1.0000
SYMI	17	-2345.2234	4724.4469	0.0000	1.0000
SYMIG	18	-2345.2412	4726.4825	0.0000	1.0000
K2PG	13	-2351.4819	4728.9638	0.0000	1.0000
K2PI	13	-2351.7493	4729.4986	0.0000	1.0000
K2PIG	14	-2351.3253	4730.6506	0.0000	1.0000
F81G	15	-2365.5925	4761.1850	0.0000	1.0000
F81I	15	-2365.6760	4761.3520	0.0000	1.0000
F81IG	16	-2365.5933	4763.1865	0.0000	1.0000
SYM	16	-2369.6333	4771.2667	0.0000	1.0000
K2P	12	-2376.8973	4777.7946	0.0000	1.0000
F81	14	-2388.8423	4805.6845	0.0000	1.0000
JC69G	12	-2394.1560	4812.3119	0.0000	1.0000
JC69I	12	-2394.2845	4812.5689	0.0000	1.0000
JC69IG	13	-2394.1304	4814.2608	0.0000	1.0000
JC69	11	-2417.8005	4857.6009	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	16	-2318.0482	4745.3037	0.5030	0.5030
HKYI	16	-2318.2624	4745.7322	0.4060	0.9089
GTRI	20	-2306.9900	4750.4892	0.0376	0.9465
GTRG	20	-2307.0877	4750.6845	0.0341	0.9807
HKYIG	17	-2317.9542	4751.9412	0.0182	0.9989
GTRIG	21	-2307.0876	4757.5099	0.0011	1.0000
HKY	15	-2342.9594	4788.3006	0.0000	1.0000
K2PG	13	-2351.4819	4791.6947	0.0000	1.0000
K2PI	13	-2351.7493	4792.2296	0.0000	1.0000
GTR	19	-2331.7162	4793.1162	0.0000	1.0000
K2PIG	14	-2351.3253	4798.2071	0.0000	1.0000
SYMG	17	-2345.1859	4806.4046	0.0000	1.0000
SYMI	17	-2345.2234	4806.4797	0.0000	1.0000
SYMIG	18	-2345.2412	4813.3408	0.0000	1.0000
F81G	15	-2365.5925	4833.5669	0.0000	1.0000
F81I	15	-2365.6760	4833.7339	0.0000	1.0000
K2P	12	-2376.8973	4835.7001	0.0000	1.0000
F81IG	16	-2365.5933	4840.3939	0.0000	1.0000
SYM	16	-2369.6333	4848.4740	0.0000	1.0000
JC69G	12	-2394.1560	4870.2175	0.0000	1.0000
JC69I	12	-2394.2845	4870.4745	0.0000	1.0000
F81	14	-2388.8423	4873.2410	0.0000	1.0000
JC69IG	13	-2394.1304	4876.9918	0.0000	1.0000
JC69	11	-2417.8005	4910.6810	0.0000	1.0000



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

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

BEGIN MRBAYES;
 Charset 591-tomodel-phylip/chr5_3310.phylip = 1 - 921;
 Partition Dummy = 1:591-tomodel-phylip/chr5_3310.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=6 rates=propinv;
 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) pinvarpr=Uniform(0.0,1.0);
END;


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

BEGIN MRBAYES;
 Charset 591-tomodel-phylip/chr5_3310.phylip = 1 - 921;
 Partition Dummy = 1:591-tomodel-phylip/chr5_3310.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=6 rates=propinv;
 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) pinvarpr=Uniform(0.0,1.0);
END;


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

BEGIN MRBAYES;
 Charset 591-tomodel-phylip/chr5_3310.phylip = 1 - 921;
 Partition Dummy = 1:591-tomodel-phylip/chr5_3310.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

