
Output from mraic.pl version 1.4.4 by Johan Nylander
-------------------------------------------------------------
16:12 04/04/2011

Input data from file "683-tomodel-phylip/chr2_8145.phylip" (7 taxa, 678 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	-1485.1739	3011.6264	0.3988	0.3988
GTRI	20	-1485.5541	3012.3868	0.2727	0.6715
GTRIG	21	-1485.0406	3013.4897	0.1571	0.8286
HKYG	16	-1491.0347	3014.8924	0.0779	0.9065
HKYI	16	-1491.3929	3015.6089	0.0545	0.9610
HKYIG	17	-1490.8218	3016.5708	0.0337	0.9947
GTR	19	-1491.0043	3021.1636	0.0034	0.9980
SYMG	17	-1494.7240	3024.3753	0.0007	0.9987
SYMI	17	-1495.0742	3025.0757	0.0005	0.9992
HKY	15	-1497.5969	3025.9190	0.0003	0.9995
SYMIG	18	-1494.5909	3026.2196	0.0003	0.9998
K2PG	13	-1500.8938	3028.3358	0.0001	0.9999
K2PI	13	-1501.2341	3029.0165	0.0001	1.0000
K2PIG	14	-1500.6841	3030.0016	0.0000	1.0000
SYM	16	-1500.5623	3033.9475	0.0000	1.0000
K2P	12	-1507.4572	3039.3836	0.0000	1.0000
F81G	15	-1520.8090	3072.3431	0.0000	1.0000
F81I	15	-1521.1196	3072.9642	0.0000	1.0000
F81IG	16	-1520.6828	3074.1887	0.0000	1.0000
F81	14	-1526.4839	3081.6012	0.0000	1.0000
JC69G	12	-1529.7254	3083.9199	0.0000	1.0000
JC69I	12	-1530.0271	3084.5234	0.0000	1.0000
JC69IG	13	-1529.6032	3085.7546	0.0000	1.0000
JC69	11	-1535.3690	3093.1344	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRG	20	-1485.1739	3010.3479	0.4080	0.4080
GTRI	20	-1485.5541	3011.1083	0.2789	0.6869
GTRIG	21	-1485.0406	3012.0812	0.1715	0.8584
HKYG	16	-1491.0347	3014.0694	0.0635	0.9219
HKYI	16	-1491.3929	3014.7859	0.0444	0.9662
HKYIG	17	-1490.8218	3015.6436	0.0289	0.9951
GTR	19	-1491.0043	3020.0086	0.0033	0.9984
SYMG	17	-1494.7240	3023.4480	0.0006	0.9990
SYMI	17	-1495.0742	3024.1484	0.0004	0.9994
SYMIG	18	-1494.5909	3025.1817	0.0002	0.9996
HKY	15	-1497.5969	3025.1939	0.0002	0.9999
K2PG	13	-1500.8938	3027.7876	0.0001	0.9999
K2PI	13	-1501.2341	3028.4683	0.0000	1.0000
K2PIG	14	-1500.6841	3029.3681	0.0000	1.0000
SYM	16	-1500.5623	3033.1245	0.0000	1.0000
K2P	12	-1507.4572	3038.9144	0.0000	1.0000
F81G	15	-1520.8090	3071.6180	0.0000	1.0000
F81I	15	-1521.1196	3072.2392	0.0000	1.0000
F81IG	16	-1520.6828	3073.3657	0.0000	1.0000
F81	14	-1526.4839	3080.9677	0.0000	1.0000
JC69G	12	-1529.7254	3083.4508	0.0000	1.0000
JC69I	12	-1530.0271	3084.0542	0.0000	1.0000
JC69IG	13	-1529.6032	3085.2064	0.0000	1.0000
JC69	11	-1535.3690	3092.7380	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	16	-1491.0347	3086.3757	0.2904	0.2904
K2PG	13	-1500.8938	3086.5365	0.2680	0.5584
HKYI	16	-1491.3929	3087.0922	0.2030	0.7614
K2PI	13	-1501.2341	3087.2172	0.1907	0.9521
HKYIG	17	-1490.8218	3092.4691	0.0138	0.9659
K2PIG	14	-1500.6841	3092.6362	0.0127	0.9786
HKY	15	-1497.5969	3092.9811	0.0107	0.9892
K2P	12	-1507.4572	3093.1442	0.0098	0.9991
SYMG	17	-1494.7240	3100.2735	0.0003	0.9994
GTRG	20	-1485.1739	3100.7308	0.0002	0.9996
SYMI	17	-1495.0742	3100.9739	0.0002	0.9998
GTRI	20	-1485.5541	3101.4912	0.0002	0.9999
SYM	16	-1500.5623	3105.4309	0.0000	1.0000
GTR	19	-1491.0043	3105.8724	0.0000	1.0000
SYMIG	18	-1494.5909	3106.5264	0.0000	1.0000
GTRIG	21	-1485.0406	3106.9833	0.0000	1.0000
JC69G	12	-1529.7254	3137.6805	0.0000	1.0000
JC69I	12	-1530.0271	3138.2840	0.0000	1.0000
F81G	15	-1520.8090	3139.4052	0.0000	1.0000
F81I	15	-1521.1196	3140.0264	0.0000	1.0000
JC69	11	-1535.3690	3142.4487	0.0000	1.0000
JC69IG	13	-1529.6032	3143.9553	0.0000	1.0000
F81	14	-1526.4839	3144.2358	0.0000	1.0000
F81IG	16	-1520.6828	3145.6721	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 683-tomodel-phylip/chr2_8145.phylip = 1 - 678;
 Partition Dummy = 1:683-tomodel-phylip/chr2_8145.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 683-tomodel-phylip/chr2_8145.phylip = 1 - 678;
 Partition Dummy = 1:683-tomodel-phylip/chr2_8145.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 683-tomodel-phylip/chr2_8145.phylip = 1 - 678;
 Partition Dummy = 1:683-tomodel-phylip/chr2_8145.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

