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

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

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


Note: number of parameters compared to sample size
is low (Nchar/Nparams < 40). Consider the use of AICc.


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

Model	df	lnL		AICc		wAICc	cwAICc
GTRG	20	-3519.5639	7080.1713	0.6421	0.6421
GTRIG	21	-3519.5642	7082.2777	0.2240	0.8661
SYMG	17	-3524.6068	7083.9709	0.0961	0.9621
SYMIG	18	-3524.6069	7086.0613	0.0338	0.9959
K2PG	13	-3532.5619	7091.5720	0.0021	0.9980
K2PIG	14	-3532.5621	7093.6422	0.0008	0.9988
HKYG	16	-3530.9494	7094.5712	0.0005	0.9993
GTRI	20	-3526.8646	7094.7728	0.0004	0.9997
HKYIG	17	-3530.9497	7096.6567	0.0002	0.9999
SYMI	17	-3531.3090	7097.3755	0.0001	1.0000
K2PI	13	-3538.8252	7104.0988	0.0000	1.0000
HKYI	16	-3537.2378	7107.1480	0.0000	1.0000
GTR	19	-3552.4838	7143.9104	0.0000	1.0000
SYM	16	-3555.9086	7144.4896	0.0000	1.0000
K2P	12	-3563.5125	7151.4088	0.0000	1.0000
HKY	15	-3562.0987	7154.7900	0.0000	1.0000
JC69G	12	-3606.0957	7236.5752	0.0000	1.0000
JC69IG	13	-3606.0960	7238.6402	0.0000	1.0000
F81G	15	-3605.5438	7241.6802	0.0000	1.0000
F81IG	16	-3605.5441	7243.7606	0.0000	1.0000
JC69I	12	-3611.7488	7247.8814	0.0000	1.0000
F81I	15	-3611.2200	7253.0326	0.0000	1.0000
JC69	11	-3634.6393	7291.6029	0.0000	1.0000
F81	14	-3634.0439	7296.6057	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRG	20	-3519.5639	7079.1278	0.6451	0.6451
GTRIG	21	-3519.5642	7081.1284	0.2372	0.8823
SYMG	17	-3524.6068	7083.2135	0.0836	0.9660
SYMIG	18	-3524.6069	7085.2138	0.0308	0.9967
K2PG	13	-3532.5619	7091.1237	0.0016	0.9983
K2PIG	14	-3532.5621	7093.1243	0.0006	0.9989
GTRI	20	-3526.8646	7093.7293	0.0004	0.9993
HKYG	16	-3530.9494	7093.8988	0.0004	0.9997
HKYIG	17	-3530.9497	7095.8993	0.0001	0.9999
SYMI	17	-3531.3090	7096.6181	0.0001	1.0000
K2PI	13	-3538.8252	7103.6505	0.0000	1.0000
HKYI	16	-3537.2378	7106.4756	0.0000	1.0000
GTR	19	-3552.4838	7142.9675	0.0000	1.0000
SYM	16	-3555.9086	7143.8171	0.0000	1.0000
K2P	12	-3563.5125	7151.0250	0.0000	1.0000
HKY	15	-3562.0987	7154.1974	0.0000	1.0000
JC69G	12	-3606.0957	7236.1914	0.0000	1.0000
JC69IG	13	-3606.0960	7238.1919	0.0000	1.0000
F81G	15	-3605.5438	7241.0876	0.0000	1.0000
F81IG	16	-3605.5441	7243.0881	0.0000	1.0000
JC69I	12	-3611.7488	7247.4976	0.0000	1.0000
F81I	15	-3611.2200	7252.4400	0.0000	1.0000
JC69	11	-3634.6393	7291.2786	0.0000	1.0000
F81	14	-3634.0439	7296.0878	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
K2PG	13	-3532.5619	7152.4395	0.9604	0.9604
K2PIG	14	-3532.5621	7159.1566	0.0334	0.9938
SYMG	17	-3524.6068	7163.3956	0.0040	0.9978
K2PI	13	-3538.8252	7164.9662	0.0018	0.9996
HKYG	16	-3530.9494	7169.3643	0.0002	0.9998
SYMIG	18	-3524.6069	7170.1125	0.0001	1.0000
GTRG	20	-3519.5639	7173.4597	0.0000	1.0000
HKYIG	17	-3530.9497	7176.0814	0.0000	1.0000
SYMI	17	-3531.3090	7176.8002	0.0000	1.0000
GTRIG	21	-3519.5642	7180.1769	0.0000	1.0000
HKYI	16	-3537.2378	7181.9411	0.0000	1.0000
GTRI	20	-3526.8646	7188.0612	0.0000	1.0000
K2P	12	-3563.5125	7207.6242	0.0000	1.0000
SYM	16	-3555.9086	7219.2827	0.0000	1.0000
HKY	15	-3562.0987	7224.9463	0.0000	1.0000
GTR	19	-3552.4838	7232.5828	0.0000	1.0000
JC69G	12	-3606.0957	7292.7905	0.0000	1.0000
JC69IG	13	-3606.0960	7299.5077	0.0000	1.0000
JC69I	12	-3611.7488	7304.0968	0.0000	1.0000
F81G	15	-3605.5438	7311.8365	0.0000	1.0000
F81IG	16	-3605.5441	7318.5537	0.0000	1.0000
F81I	15	-3611.2200	7323.1889	0.0000	1.0000
JC69	11	-3634.6393	7343.1611	0.0000	1.0000
F81	14	-3634.0439	7362.1201	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/chr2_4153.phylip = 1 - 826;
 Partition Dummy = 1:683-tomodel-phylip/chr2_4153.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/chr2_4153.phylip = 1 - 826;
 Partition Dummy = 1:683-tomodel-phylip/chr2_4153.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/chr2_4153.phylip = 1 - 826;
 Partition Dummy = 1:683-tomodel-phylip/chr2_4153.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

