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

Input data from file "591-tomodel-phylip/chr1_32301.phylip" (7 taxa, 816 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	-3038.7337	6118.5239	0.7227	0.7227
GTRIG	21	-3038.7337	6120.6311	0.2520	0.9747
GTRI	20	-3042.1036	6125.2639	0.0249	0.9995
HKYG	16	-3050.5805	6133.8419	0.0003	0.9999
HKYIG	17	-3050.5830	6135.9330	0.0001	1.0000
HKYI	16	-3053.0431	6138.7670	0.0000	1.0000
GTR	19	-3067.7192	6174.3932	0.0000	1.0000
HKY	15	-3079.4632	6189.5265	0.0000	1.0000
F81G	15	-3119.1987	6268.9974	0.0000	1.0000
F81IG	16	-3119.2105	6271.1020	0.0000	1.0000
F81I	15	-3120.3999	6271.3998	0.0000	1.0000
K2PG	13	-3124.9367	6276.3273	0.0000	1.0000
K2PIG	14	-3124.9413	6278.4070	0.0000	1.0000
K2PI	13	-3126.6983	6279.8505	0.0000	1.0000
SYMG	17	-3124.4632	6283.6933	0.0000	1.0000
SYMIG	18	-3124.4944	6285.8471	0.0000	1.0000
SYMI	17	-3125.9267	6286.6203	0.0000	1.0000
F81	14	-3146.9231	6322.3706	0.0000	1.0000
K2P	12	-3158.8535	6342.0956	0.0000	1.0000
SYM	16	-3158.4734	6349.6277	0.0000	1.0000
JC69G	12	-3178.4014	6381.1914	0.0000	1.0000
JC69I	12	-3179.4150	6383.2185	0.0000	1.0000
JC69IG	13	-3178.4110	6383.2758	0.0000	1.0000
JC69	11	-3209.0203	6440.3689	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRG	20	-3038.7337	6117.4673	0.7128	0.7128
GTRIG	21	-3038.7337	6119.4673	0.2622	0.9751
GTRI	20	-3042.1036	6124.2073	0.0245	0.9996
HKYG	16	-3050.5805	6133.1610	0.0003	0.9999
HKYIG	17	-3050.5830	6135.1661	0.0001	1.0000
HKYI	16	-3053.0431	6138.0861	0.0000	1.0000
GTR	19	-3067.7192	6173.4384	0.0000	1.0000
HKY	15	-3079.4632	6188.9265	0.0000	1.0000
F81G	15	-3119.1987	6268.3974	0.0000	1.0000
F81IG	16	-3119.2105	6270.4211	0.0000	1.0000
F81I	15	-3120.3999	6270.7998	0.0000	1.0000
K2PG	13	-3124.9367	6275.8734	0.0000	1.0000
K2PIG	14	-3124.9413	6277.8826	0.0000	1.0000
K2PI	13	-3126.6983	6279.3966	0.0000	1.0000
SYMG	17	-3124.4632	6282.9264	0.0000	1.0000
SYMIG	18	-3124.4944	6284.9888	0.0000	1.0000
SYMI	17	-3125.9267	6285.8533	0.0000	1.0000
F81	14	-3146.9231	6321.8462	0.0000	1.0000
K2P	12	-3158.8535	6341.7071	0.0000	1.0000
SYM	16	-3158.4734	6348.9468	0.0000	1.0000
JC69G	12	-3178.4014	6380.8028	0.0000	1.0000
JC69IG	13	-3178.4110	6382.8219	0.0000	1.0000
JC69I	12	-3179.4150	6382.8299	0.0000	1.0000
JC69	11	-3209.0203	6440.0405	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	16	-3050.5805	6208.4317	0.7438	0.7438
GTRG	20	-3038.7337	6211.5556	0.1560	0.8998
HKYI	16	-3053.0431	6213.3567	0.0634	0.9632
HKYIG	17	-3050.5830	6215.1411	0.0260	0.9892
GTRIG	21	-3038.7337	6218.2600	0.0055	0.9946
GTRI	20	-3042.1036	6218.2955	0.0054	1.0000
HKY	15	-3079.4632	6259.4927	0.0000	1.0000
GTR	19	-3067.7192	6262.8223	0.0000	1.0000
K2PG	13	-3124.9367	6337.0308	0.0000	1.0000
F81G	15	-3119.1987	6338.9637	0.0000	1.0000
K2PI	13	-3126.6983	6340.5540	0.0000	1.0000
F81I	15	-3120.3999	6341.3660	0.0000	1.0000
K2PIG	14	-3124.9413	6343.7444	0.0000	1.0000
F81IG	16	-3119.2105	6345.6917	0.0000	1.0000
SYMG	17	-3124.4632	6362.9015	0.0000	1.0000
SYMI	17	-3125.9267	6365.8284	0.0000	1.0000
SYMIG	18	-3124.4944	6369.6683	0.0000	1.0000
F81	14	-3146.9231	6387.7080	0.0000	1.0000
K2P	12	-3158.8535	6398.1601	0.0000	1.0000
SYM	16	-3158.4734	6424.2174	0.0000	1.0000
JC69G	12	-3178.4014	6437.2558	0.0000	1.0000
JC69I	12	-3179.4150	6439.2829	0.0000	1.0000
JC69IG	13	-3178.4110	6443.9793	0.0000	1.0000
JC69	11	-3209.0203	6491.7891	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/chr1_32301.phylip = 1 - 816;
 Partition Dummy = 1:591-tomodel-phylip/chr1_32301.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/chr1_32301.phylip = 1 - 816;
 Partition Dummy = 1:591-tomodel-phylip/chr1_32301.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/chr1_32301.phylip = 1 - 816;
 Partition Dummy = 1:591-tomodel-phylip/chr1_32301.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

