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

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

Minimum AIC  model: GTRG
Minimum AICc model: GTRG
Minimum BIC  model: HKYG


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	-2447.6464	4936.4340	0.5086	0.5086
GTRI	20	-2448.5637	4938.2688	0.2032	0.7118
GTRIG	21	-2447.6581	4938.5733	0.1745	0.8863
HKYG	16	-2453.8592	4940.4536	0.0682	0.9545
HKYIG	17	-2453.8598	4942.5478	0.0239	0.9784
HKYI	16	-2455.0228	4942.7808	0.0213	0.9997
SYMI	17	-2458.9467	4952.7215	0.0001	0.9998
SYMG	17	-2459.1071	4953.0424	0.0001	1.0000
SYMIG	18	-2459.1227	4955.1722	0.0000	1.0000
K2PG	13	-2467.6317	4961.7534	0.0000	1.0000
K2PI	13	-2468.5150	4963.5198	0.0000	1.0000
K2PIG	14	-2467.6392	4963.8445	0.0000	1.0000
GTR	19	-2466.1009	4971.2330	0.0000	1.0000
HKY	15	-2473.1194	4976.8866	0.0000	1.0000
SYM	16	-2476.0133	4984.7617	0.0000	1.0000
K2P	12	-2485.2413	4994.9020	0.0000	1.0000
F81G	15	-2491.5510	5013.7498	0.0000	1.0000
F81I	15	-2492.3548	5015.3574	0.0000	1.0000
F81IG	16	-2491.5584	5015.8520	0.0000	1.0000
JC69G	12	-2502.9446	5030.3086	0.0000	1.0000
JC69I	12	-2503.6077	5031.6347	0.0000	1.0000
JC69IG	13	-2502.9592	5032.4083	0.0000	1.0000
F81	14	-2508.9519	5046.4698	0.0000	1.0000
JC69	11	-2519.4834	5061.3212	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRG	20	-2447.6464	4935.2927	0.5135	0.5135
GTRI	20	-2448.5637	4937.1275	0.2052	0.7187
GTRIG	21	-2447.6581	4937.3162	0.1867	0.9054
HKYG	16	-2453.8592	4939.7184	0.0562	0.9615
HKYIG	17	-2453.8598	4941.7197	0.0207	0.9822
HKYI	16	-2455.0228	4942.0457	0.0175	0.9997
SYMI	17	-2458.9467	4951.8933	0.0001	0.9998
SYMG	17	-2459.1071	4952.2142	0.0001	1.0000
SYMIG	18	-2459.1227	4954.2453	0.0000	1.0000
K2PG	13	-2467.6317	4961.2635	0.0000	1.0000
K2PI	13	-2468.5150	4963.0299	0.0000	1.0000
K2PIG	14	-2467.6392	4963.2785	0.0000	1.0000
GTR	19	-2466.1009	4970.2018	0.0000	1.0000
HKY	15	-2473.1194	4976.2388	0.0000	1.0000
SYM	16	-2476.0133	4984.0266	0.0000	1.0000
K2P	12	-2485.2413	4994.4826	0.0000	1.0000
F81G	15	-2491.5510	5013.1021	0.0000	1.0000
F81I	15	-2492.3548	5014.7096	0.0000	1.0000
F81IG	16	-2491.5584	5015.1169	0.0000	1.0000
JC69G	12	-2502.9446	5029.8893	0.0000	1.0000
JC69I	12	-2503.6077	5031.2154	0.0000	1.0000
JC69IG	13	-2502.9592	5031.9184	0.0000	1.0000
F81	14	-2508.9519	5045.9038	0.0000	1.0000
JC69	11	-2519.4834	5060.9668	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	16	-2453.8592	5013.7882	0.7237	0.7237
HKYI	16	-2455.0228	5016.1155	0.2260	0.9497
HKYIG	17	-2453.8598	5020.4189	0.0263	0.9760
K2PG	13	-2467.6317	5021.4452	0.0157	0.9917
K2PI	13	-2468.5150	5023.2116	0.0065	0.9982
GTRG	20	-2447.6464	5027.8800	0.0006	0.9989
K2PIG	14	-2467.6392	5028.0895	0.0006	0.9994
GTRI	20	-2448.5637	5029.7147	0.0003	0.9997
SYMI	17	-2458.9467	5030.5925	0.0002	0.9998
SYMG	17	-2459.1071	5030.9134	0.0001	1.0000
GTRIG	21	-2447.6581	5034.5328	0.0000	1.0000
SYMIG	18	-2459.1227	5037.5739	0.0000	1.0000
HKY	15	-2473.1194	5045.6792	0.0000	1.0000
K2P	12	-2485.2413	5050.0350	0.0000	1.0000
SYM	16	-2476.0133	5058.0964	0.0000	1.0000
GTR	19	-2466.1009	5058.1597	0.0000	1.0000
F81G	15	-2491.5510	5082.5425	0.0000	1.0000
F81I	15	-2492.3548	5084.1500	0.0000	1.0000
JC69G	12	-2502.9446	5085.4416	0.0000	1.0000
JC69I	12	-2503.6077	5086.7677	0.0000	1.0000
F81IG	16	-2491.5584	5089.1867	0.0000	1.0000
JC69IG	13	-2502.9592	5092.1001	0.0000	1.0000
F81	14	-2508.9519	5110.7148	0.0000	1.0000
JC69	11	-2519.4834	5111.8898	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/chr2_21271.phylip = 1 - 757;
 Partition Dummy = 1:591-tomodel-phylip/chr2_21271.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/chr2_21271.phylip = 1 - 757;
 Partition Dummy = 1:591-tomodel-phylip/chr2_21271.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/chr2_21271.phylip = 1 - 757;
 Partition Dummy = 1:591-tomodel-phylip/chr2_21271.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

