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

Input data from file "591-tomodel-phylip/chr2_24827.phylip" (7 taxa, 513 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	-1587.4841	3216.6755	0.4384	0.4384
GTRIG	21	-1587.4734	3218.8287	0.1494	0.5878
GTRI	20	-1588.6869	3219.0811	0.1317	0.7194
SYMG	17	-1591.9730	3219.1824	0.1252	0.8446
SYMI	17	-1592.8693	3220.9749	0.0511	0.8957
SYMIG	18	-1591.9720	3221.3286	0.0428	0.9385
K2PG	13	-1598.0441	3222.8178	0.0203	0.9588
K2PI	13	-1598.3822	3223.4938	0.0145	0.9733
HKYG	16	-1595.5424	3224.1816	0.0103	0.9836
K2PIG	14	-1598.0457	3224.9348	0.0071	0.9906
HKYI	16	-1596.1124	3225.3216	0.0058	0.9964
HKYIG	17	-1595.5362	3226.3088	0.0035	1.0000
GTR	19	-1601.2434	3242.0284	0.0000	1.0000
SYM	16	-1605.2180	3243.5328	0.0000	1.0000
K2P	12	-1611.3252	3247.2745	0.0000	1.0000
HKY	15	-1608.6831	3248.3320	0.0000	1.0000
JC69G	12	-1611.9006	3248.4253	0.0000	1.0000
JC69I	12	-1612.5229	3249.6698	0.0000	1.0000
F81G	15	-1609.7328	3250.4313	0.0000	1.0000
JC69IG	13	-1611.9007	3250.5309	0.0000	1.0000
F81I	15	-1610.4132	3251.7922	0.0000	1.0000
F81IG	16	-1609.7298	3252.5564	0.0000	1.0000
JC69	11	-1624.3863	3271.2994	0.0000	1.0000
F81	14	-1622.2012	3273.2458	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRG	20	-1587.4841	3214.9682	0.4618	0.4618
GTRIG	21	-1587.4734	3216.9468	0.1717	0.6335
GTRI	20	-1588.6869	3217.3738	0.1387	0.7722
SYMG	17	-1591.9730	3217.9461	0.1042	0.8764
SYMI	17	-1592.8693	3219.7386	0.0425	0.9189
SYMIG	18	-1591.9720	3219.9440	0.0384	0.9572
K2PG	13	-1598.0441	3222.0883	0.0131	0.9704
K2PI	13	-1598.3822	3222.7644	0.0094	0.9797
HKYG	16	-1595.5424	3223.0849	0.0080	0.9877
K2PIG	14	-1598.0457	3224.0914	0.0048	0.9925
HKYI	16	-1596.1124	3224.2249	0.0045	0.9970
HKYIG	17	-1595.5362	3225.0724	0.0030	1.0000
GTR	19	-1601.2434	3240.4868	0.0000	1.0000
SYM	16	-1605.2180	3242.4360	0.0000	1.0000
K2P	12	-1611.3252	3246.6505	0.0000	1.0000
HKY	15	-1608.6831	3247.3662	0.0000	1.0000
JC69G	12	-1611.9006	3247.8013	0.0000	1.0000
JC69I	12	-1612.5229	3249.0458	0.0000	1.0000
F81G	15	-1609.7328	3249.4655	0.0000	1.0000
JC69IG	13	-1611.9007	3249.8014	0.0000	1.0000
F81I	15	-1610.4132	3250.8264	0.0000	1.0000
F81IG	16	-1609.7298	3251.4596	0.0000	1.0000
JC69	11	-1624.3863	3270.7725	0.0000	1.0000
F81	14	-1622.2012	3272.4024	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
K2PG	13	-1598.0441	3277.2119	0.5677	0.5677
K2PI	13	-1598.3822	3277.8880	0.4049	0.9726
K2PIG	14	-1598.0457	3283.4552	0.0250	0.9976
SYMG	17	-1591.9730	3290.0307	0.0009	0.9986
HKYG	16	-1595.5424	3290.9293	0.0006	0.9992
SYMI	17	-1592.8693	3291.8233	0.0004	0.9995
HKYI	16	-1596.1124	3292.0693	0.0003	0.9999
SYMIG	18	-1591.9720	3296.2690	0.0000	0.9999
HKYIG	17	-1595.5362	3297.1571	0.0000	0.9999
K2P	12	-1611.3252	3297.5338	0.0000	1.0000
JC69G	12	-1611.9006	3298.6846	0.0000	1.0000
GTRG	20	-1587.4841	3299.7737	0.0000	1.0000
JC69I	12	-1612.5229	3299.9291	0.0000	1.0000
GTRI	20	-1588.6869	3302.1793	0.0000	1.0000
JC69IG	13	-1611.9007	3304.9250	0.0000	1.0000
GTRIG	21	-1587.4734	3305.9926	0.0000	1.0000
SYM	16	-1605.2180	3310.2805	0.0000	1.0000
HKY	15	-1608.6831	3310.9703	0.0000	1.0000
F81G	15	-1609.7328	3313.0697	0.0000	1.0000
F81I	15	-1610.4132	3314.4306	0.0000	1.0000
JC69	11	-1624.3863	3317.4155	0.0000	1.0000
F81IG	16	-1609.7298	3319.3040	0.0000	1.0000
GTR	19	-1601.2434	3321.0521	0.0000	1.0000
F81	14	-1622.2012	3331.7662	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_24827.phylip = 1 - 513;
 Partition Dummy = 1:591-tomodel-phylip/chr2_24827.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_24827.phylip = 1 - 513;
 Partition Dummy = 1:591-tomodel-phylip/chr2_24827.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 591-tomodel-phylip/chr2_24827.phylip = 1 - 513;
 Partition Dummy = 1:591-tomodel-phylip/chr2_24827.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

