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

Input data from file "591-tomodel-phylip/chr3_24775.phylip" (7 taxa, 509 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	-1583.9603	3209.6418	0.4866	0.4866
GTRI	20	-1584.7185	3211.1583	0.2280	0.7146
GTRIG	21	-1583.9437	3211.7847	0.1667	0.8813
HKYG	16	-1590.4262	3213.9581	0.0562	0.9375
HKYI	16	-1590.6893	3214.4843	0.0432	0.9807
HKYIG	17	-1590.4258	3216.0981	0.0193	1.0000
GTR	19	-1597.5869	3234.7280	0.0000	1.0000
HKY	15	-1604.7375	3240.4487	0.0000	1.0000
SYMG	17	-1605.2298	3245.7061	0.0000	1.0000
SYMIG	18	-1605.2222	3247.8403	0.0000	1.0000
SYMI	17	-1606.3094	3247.8652	0.0000	1.0000
K2PG	13	-1610.8323	3248.3999	0.0000	1.0000
K2PI	13	-1611.3586	3249.4527	0.0000	1.0000
K2PIG	14	-1610.8351	3250.5204	0.0000	1.0000
SYM	16	-1619.2640	3271.6337	0.0000	1.0000
F81G	15	-1620.9900	3272.9537	0.0000	1.0000
F81I	15	-1621.2225	3273.4187	0.0000	1.0000
F81IG	16	-1620.9757	3275.0572	0.0000	1.0000
K2P	12	-1625.9194	3276.4679	0.0000	1.0000
F81	14	-1634.4144	3297.6790	0.0000	1.0000
JC69G	12	-1638.9129	3302.4549	0.0000	1.0000
JC69I	12	-1639.2990	3303.2269	0.0000	1.0000
JC69IG	13	-1638.9086	3304.5526	0.0000	1.0000
JC69	11	-1652.7314	3327.9939	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRG	20	-1583.9603	3207.9205	0.4941	0.4941
GTRI	20	-1584.7185	3209.4370	0.2315	0.7255
GTRIG	21	-1583.9437	3209.8874	0.1848	0.9103
HKYG	16	-1590.4262	3212.8524	0.0420	0.9523
HKYI	16	-1590.6893	3213.3786	0.0323	0.9846
HKYIG	17	-1590.4258	3214.8517	0.0154	1.0000
GTR	19	-1597.5869	3233.1738	0.0000	1.0000
HKY	15	-1604.7375	3239.4751	0.0000	1.0000
SYMG	17	-1605.2298	3244.4596	0.0000	1.0000
SYMIG	18	-1605.2222	3246.4443	0.0000	1.0000
SYMI	17	-1606.3094	3246.6187	0.0000	1.0000
K2PG	13	-1610.8323	3247.6646	0.0000	1.0000
K2PI	13	-1611.3586	3248.7173	0.0000	1.0000
K2PIG	14	-1610.8351	3249.6702	0.0000	1.0000
SYM	16	-1619.2640	3270.5280	0.0000	1.0000
F81G	15	-1620.9900	3271.9800	0.0000	1.0000
F81I	15	-1621.2225	3272.4451	0.0000	1.0000
F81IG	16	-1620.9757	3273.9515	0.0000	1.0000
K2P	12	-1625.9194	3275.8389	0.0000	1.0000
F81	14	-1634.4144	3296.8288	0.0000	1.0000
JC69G	12	-1638.9129	3301.8259	0.0000	1.0000
JC69I	12	-1639.2990	3302.5979	0.0000	1.0000
JC69IG	13	-1638.9086	3303.8172	0.0000	1.0000
JC69	11	-1652.7314	3327.4627	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	16	-1590.4262	3280.5715	0.5504	0.5504
HKYI	16	-1590.6893	3281.0978	0.4231	0.9735
HKYIG	17	-1590.4258	3286.8033	0.0244	0.9979
GTRG	20	-1583.9603	3292.5695	0.0014	0.9993
GTRI	20	-1584.7185	3294.0859	0.0006	0.9999
GTRIG	21	-1583.9437	3298.7688	0.0001	1.0000
K2PG	13	-1610.8323	3302.6864	0.0000	1.0000
HKY	15	-1604.7375	3302.9618	0.0000	1.0000
K2PI	13	-1611.3586	3303.7391	0.0000	1.0000
K2PIG	14	-1610.8351	3308.9245	0.0000	1.0000
GTR	19	-1597.5869	3313.5903	0.0000	1.0000
SYMG	17	-1605.2298	3316.4112	0.0000	1.0000
SYMI	17	-1606.3094	3318.5703	0.0000	1.0000
SYMIG	18	-1605.2222	3322.6284	0.0000	1.0000
K2P	12	-1625.9194	3326.6282	0.0000	1.0000
F81G	15	-1620.9900	3335.4668	0.0000	1.0000
F81I	15	-1621.2225	3335.9318	0.0000	1.0000
SYM	16	-1619.2640	3338.2472	0.0000	1.0000
F81IG	16	-1620.9757	3341.6706	0.0000	1.0000
JC69G	12	-1638.9129	3352.6153	0.0000	1.0000
JC69I	12	-1639.2990	3353.3873	0.0000	1.0000
F81	14	-1634.4144	3356.0831	0.0000	1.0000
JC69IG	13	-1638.9086	3358.8390	0.0000	1.0000
JC69	11	-1652.7314	3374.0196	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/chr3_24775.phylip = 1 - 509;
 Partition Dummy = 1:591-tomodel-phylip/chr3_24775.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/chr3_24775.phylip = 1 - 509;
 Partition Dummy = 1:591-tomodel-phylip/chr3_24775.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/chr3_24775.phylip = 1 - 509;
 Partition Dummy = 1:591-tomodel-phylip/chr3_24775.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

