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

Input data from file "591-tomodel-phylip/chr2_23521.phylip" (7 taxa, 776 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	-3004.5663	6050.2452	0.6958	0.6958
GTRIG	21	-3004.5663	6052.3581	0.2419	0.9377
GTRI	20	-3007.3756	6055.8638	0.0419	0.9797
HKYG	16	-3012.7643	6058.2454	0.0127	0.9924
HKYIG	17	-3012.7726	6060.3527	0.0044	0.9968
HKYI	16	-3014.8382	6062.3931	0.0016	0.9984
SYMG	17	-3014.5308	6063.8690	0.0008	0.9992
K2PG	13	-3019.7042	6065.8861	0.0003	0.9995
SYMIG	18	-3014.5332	6065.9700	0.0003	0.9998
K2PIG	14	-3019.7164	6067.9847	0.0001	0.9999
SYMI	17	-3016.5985	6068.0044	0.0001	1.0000
K2PI	13	-3021.5212	6069.5202	0.0000	1.0000
GTR	19	-3035.1031	6109.2115	0.0000	1.0000
HKY	15	-3041.9319	6114.4954	0.0000	1.0000
SYM	16	-3043.5745	6119.8657	0.0000	1.0000
K2P	12	-3048.7279	6121.8646	0.0000	1.0000
F81G	15	-3071.6514	6173.9344	0.0000	1.0000
F81IG	16	-3071.6453	6176.0074	0.0000	1.0000
F81I	15	-3073.0432	6176.7180	0.0000	1.0000
JC69G	12	-3079.0053	6182.4196	0.0000	1.0000
JC69IG	13	-3078.9785	6184.4346	0.0000	1.0000
JC69I	12	-3080.1774	6184.7637	0.0000	1.0000
F81	14	-3098.2934	6225.1387	0.0000	1.0000
JC69	11	-3105.1020	6232.5495	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRG	20	-3004.5663	6049.1327	0.6885	0.6885
GTRIG	21	-3004.5663	6051.1327	0.2533	0.9418
GTRI	20	-3007.3756	6054.7513	0.0415	0.9833
HKYG	16	-3012.7643	6057.5287	0.0103	0.9936
HKYIG	17	-3012.7726	6059.5453	0.0038	0.9974
HKYI	16	-3014.8382	6061.6763	0.0013	0.9987
SYMG	17	-3014.5308	6063.0617	0.0007	0.9994
SYMIG	18	-3014.5332	6065.0664	0.0002	0.9996
K2PG	13	-3019.7042	6065.4084	0.0002	0.9998
SYMI	17	-3016.5985	6067.1970	0.0001	0.9999
K2PIG	14	-3019.7164	6067.4328	0.0001	1.0000
K2PI	13	-3021.5212	6069.0425	0.0000	1.0000
GTR	19	-3035.1031	6108.2062	0.0000	1.0000
HKY	15	-3041.9319	6113.8638	0.0000	1.0000
SYM	16	-3043.5745	6119.1490	0.0000	1.0000
K2P	12	-3048.7279	6121.4557	0.0000	1.0000
F81G	15	-3071.6514	6173.3029	0.0000	1.0000
F81IG	16	-3071.6453	6175.2907	0.0000	1.0000
F81I	15	-3073.0432	6176.0864	0.0000	1.0000
JC69G	12	-3079.0053	6182.0107	0.0000	1.0000
JC69IG	13	-3078.9785	6183.9569	0.0000	1.0000
JC69I	12	-3080.1774	6184.3548	0.0000	1.0000
F81	14	-3098.2934	6224.5868	0.0000	1.0000
JC69	11	-3105.1020	6232.2039	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
K2PG	13	-3019.7042	6125.9123	0.7974	0.7974
K2PI	13	-3021.5212	6129.5465	0.1296	0.9270
HKYG	16	-3012.7643	6131.9951	0.0381	0.9651
K2PIG	14	-3019.7164	6132.5910	0.0283	0.9933
HKYI	16	-3014.8382	6136.1428	0.0048	0.9981
HKYIG	17	-3012.7726	6138.6659	0.0014	0.9995
SYMG	17	-3014.5308	6142.1823	0.0002	0.9997
GTRG	20	-3004.5663	6142.2157	0.0002	0.9999
SYMI	17	-3016.5985	6146.3176	0.0000	1.0000
GTRI	20	-3007.3756	6147.8343	0.0000	1.0000
SYMIG	18	-3014.5332	6148.8412	0.0000	1.0000
GTRIG	21	-3004.5663	6148.8699	0.0000	1.0000
K2P	12	-3048.7279	6177.3056	0.0000	1.0000
HKY	15	-3041.9319	6183.6761	0.0000	1.0000
SYM	16	-3043.5745	6193.6154	0.0000	1.0000
GTR	19	-3035.1031	6196.6351	0.0000	1.0000
JC69G	12	-3079.0053	6237.8605	0.0000	1.0000
JC69I	12	-3080.1774	6240.2046	0.0000	1.0000
F81G	15	-3071.6514	6243.1151	0.0000	1.0000
JC69IG	13	-3078.9785	6244.4609	0.0000	1.0000
F81I	15	-3073.0432	6245.8987	0.0000	1.0000
F81IG	16	-3071.6453	6249.7571	0.0000	1.0000
JC69	11	-3105.1020	6283.3996	0.0000	1.0000
F81	14	-3098.2934	6289.7450	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_23521.phylip = 1 - 776;
 Partition Dummy = 1:591-tomodel-phylip/chr2_23521.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_23521.phylip = 1 - 776;
 Partition Dummy = 1:591-tomodel-phylip/chr2_23521.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_23521.phylip = 1 - 776;
 Partition Dummy = 1:591-tomodel-phylip/chr2_23521.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

