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

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

Minimum AIC  model: HKYG
Minimum AICc model: HKYG
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
HKYG	16	-2125.9203	4284.6465	0.5668	0.5668
HKYI	16	-2126.9254	4286.6567	0.2075	0.7743
HKYIG	17	-2125.9217	4286.7514	0.1979	0.9721
GTRG	20	-2125.2685	4291.7888	0.0159	0.9881
GTRI	20	-2126.1714	4293.5946	0.0065	0.9945
GTRIG	21	-2125.2754	4293.9299	0.0055	1.0000
HKY	15	-2142.9747	4316.6594	0.0000	1.0000
GTR	19	-2142.0017	4323.1345	0.0000	1.0000
F81G	15	-2159.2381	4349.1862	0.0000	1.0000
F81I	15	-2160.0105	4350.7311	0.0000	1.0000
F81IG	16	-2159.2442	4351.2943	0.0000	1.0000
K2PG	13	-2165.2195	4356.9759	0.0000	1.0000
SYMG	17	-2161.3157	4357.5394	0.0000	1.0000
K2PI	13	-2166.1438	4358.8245	0.0000	1.0000
K2PIG	14	-2165.2246	4359.0696	0.0000	1.0000
SYMI	17	-2162.2707	4359.4495	0.0000	1.0000
SYMIG	18	-2161.3133	4359.6429	0.0000	1.0000
F81	14	-2173.6602	4375.9408	0.0000	1.0000
K2P	12	-2180.6640	4385.7874	0.0000	1.0000
SYM	16	-2177.0111	4386.8281	0.0000	1.0000
JC69G	12	-2192.3887	4409.2369	0.0000	1.0000
JC69I	12	-2193.1190	4410.6975	0.0000	1.0000
JC69IG	13	-2192.3969	4411.3307	0.0000	1.0000
JC69	11	-2206.2424	4434.8730	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
HKYG	16	-2125.9203	4283.8406	0.5569	0.5569
HKYIG	17	-2125.9217	4285.8434	0.2046	0.7615
HKYI	16	-2126.9254	4285.8508	0.2038	0.9653
GTRG	20	-2125.2685	4290.5369	0.0196	0.9849
GTRI	20	-2126.1714	4292.3428	0.0079	0.9928
GTRIG	21	-2125.2754	4292.5508	0.0072	1.0000
HKY	15	-2142.9747	4315.9494	0.0000	1.0000
GTR	19	-2142.0017	4322.0035	0.0000	1.0000
F81G	15	-2159.2381	4348.4761	0.0000	1.0000
F81I	15	-2160.0105	4350.0210	0.0000	1.0000
F81IG	16	-2159.2442	4350.4884	0.0000	1.0000
K2PG	13	-2165.2195	4356.4390	0.0000	1.0000
SYMG	17	-2161.3157	4356.6314	0.0000	1.0000
K2PI	13	-2166.1438	4358.2876	0.0000	1.0000
K2PIG	14	-2165.2246	4358.4492	0.0000	1.0000
SYMI	17	-2162.2707	4358.5415	0.0000	1.0000
SYMIG	18	-2161.3133	4358.6266	0.0000	1.0000
F81	14	-2173.6602	4375.3205	0.0000	1.0000
K2P	12	-2180.6640	4385.3279	0.0000	1.0000
SYM	16	-2177.0111	4386.0222	0.0000	1.0000
JC69G	12	-2192.3887	4408.7774	0.0000	1.0000
JC69I	12	-2193.1190	4410.2380	0.0000	1.0000
JC69IG	13	-2192.3969	4410.7938	0.0000	1.0000
JC69	11	-2206.2424	4434.4847	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	16	-2125.9203	4356.4740	0.7123	0.7123
HKYI	16	-2126.9254	4358.4842	0.2607	0.9730
HKYIG	17	-2125.9217	4363.0163	0.0270	1.0000
GTRG	20	-2125.2685	4381.3286	0.0000	1.0000
GTRI	20	-2126.1714	4383.1345	0.0000	1.0000
HKY	15	-2142.9747	4384.0431	0.0000	1.0000
GTRIG	21	-2125.2754	4387.8821	0.0000	1.0000
GTR	19	-2142.0017	4408.2556	0.0000	1.0000
K2PG	13	-2165.2195	4415.4536	0.0000	1.0000
F81G	15	-2159.2381	4416.5699	0.0000	1.0000
K2PI	13	-2166.1438	4417.3023	0.0000	1.0000
F81I	15	-2160.0105	4418.1148	0.0000	1.0000
K2PIG	14	-2165.2246	4422.0034	0.0000	1.0000
F81IG	16	-2159.2442	4423.1218	0.0000	1.0000
SYMG	17	-2161.3157	4433.8044	0.0000	1.0000
SYMI	17	-2162.2707	4435.7144	0.0000	1.0000
F81	14	-2173.6602	4438.8747	0.0000	1.0000
K2P	12	-2180.6640	4439.8029	0.0000	1.0000
SYMIG	18	-2161.3133	4440.3391	0.0000	1.0000
SYM	16	-2177.0111	4458.6556	0.0000	1.0000
JC69G	12	-2192.3887	4463.2524	0.0000	1.0000
JC69I	12	-2193.1190	4464.7131	0.0000	1.0000
JC69IG	13	-2192.3969	4469.8084	0.0000	1.0000
JC69	11	-2206.2424	4484.4202	0.0000	1.0000



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

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

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


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

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


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

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

