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

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

Minimum AIC  model: HKY
Minimum AICc model: HKY
Minimum BIC  model: HKY


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
HKY	15	-717.3634	1466.0527	0.3036	0.3036
HKYI	16	-716.4069	1466.3208	0.2655	0.5691
HKYG	16	-716.4328	1466.3725	0.2587	0.8278
HKYIG	17	-716.4376	1468.5753	0.0860	0.9138
GTR	19	-715.2722	1470.6673	0.0302	0.9440
GTRI	20	-714.4116	1471.1762	0.0234	0.9674
GTRG	20	-714.4306	1471.2141	0.0230	0.9904
GTRIG	21	-714.4331	1473.4617	0.0075	0.9979
F81	14	-724.5580	1478.2731	0.0007	0.9986
F81I	15	-723.5444	1478.4148	0.0006	0.9992
F81G	15	-723.5759	1478.4777	0.0006	0.9998
F81IG	16	-723.5829	1480.6727	0.0002	1.0000
K2PI	13	-733.9670	1494.9341	0.0000	1.0000
K2PG	13	-733.9903	1494.9807	0.0000	1.0000
K2P	12	-735.2184	1495.2916	0.0000	1.0000
K2PIG	14	-733.9996	1497.1563	0.0000	1.0000
SYMI	17	-733.1377	1501.9753	0.0000	1.0000
SYM	16	-734.2500	1502.0069	0.0000	1.0000
SYMG	17	-733.1625	1502.0251	0.0000	1.0000
SYMIG	18	-733.1625	1504.2303	0.0000	1.0000
JC69I	12	-740.3978	1505.6504	0.0000	1.0000
JC69G	12	-740.4278	1505.7104	0.0000	1.0000
JC69	11	-741.5693	1505.8598	0.0000	1.0000
JC69IG	13	-740.4375	1507.8751	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
HKY	15	-717.3634	1464.7267	0.2711	0.2711
HKYI	16	-716.4069	1464.8139	0.2595	0.5306
HKYG	16	-716.4328	1464.8656	0.2529	0.7835
HKYIG	17	-716.4376	1466.8753	0.0926	0.8761
GTR	19	-715.2722	1468.5444	0.0402	0.9163
GTRI	20	-714.4116	1468.8233	0.0350	0.9513
GTRG	20	-714.4306	1468.8611	0.0343	0.9856
GTRIG	21	-714.4331	1470.8662	0.0126	0.9981
F81I	15	-723.5444	1477.0888	0.0006	0.9987
F81	14	-724.5580	1477.1161	0.0006	0.9993
F81G	15	-723.5759	1477.1518	0.0005	0.9998
F81IG	16	-723.5829	1479.1658	0.0002	1.0000
K2PI	13	-733.9670	1493.9341	0.0000	1.0000
K2PG	13	-733.9903	1493.9807	0.0000	1.0000
K2P	12	-735.2184	1494.4368	0.0000	1.0000
K2PIG	14	-733.9996	1495.9993	0.0000	1.0000
SYMI	17	-733.1377	1500.2753	0.0000	1.0000
SYMG	17	-733.1625	1500.3251	0.0000	1.0000
SYM	16	-734.2500	1500.5000	0.0000	1.0000
SYMIG	18	-733.1625	1502.3250	0.0000	1.0000
JC69I	12	-740.3978	1504.7956	0.0000	1.0000
JC69G	12	-740.4278	1504.8556	0.0000	1.0000
JC69	11	-741.5693	1505.1385	0.0000	1.0000
JC69IG	13	-740.4375	1506.8751	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKY	15	-717.3634	1523.7502	0.7752	0.7752
HKYI	16	-716.4069	1527.7722	0.1038	0.8790
HKYG	16	-716.4328	1527.8239	0.1011	0.9801
F81	14	-724.5580	1532.2046	0.0113	0.9914
HKYIG	17	-716.4376	1533.7685	0.0052	0.9966
F81I	15	-723.5444	1536.1122	0.0016	0.9982
F81G	15	-723.5759	1536.1752	0.0016	0.9997
K2P	12	-735.2184	1541.6556	0.0001	0.9998
F81IG	16	-723.5829	1542.1241	0.0001	0.9999
GTR	19	-715.2722	1543.3073	0.0000	0.9999
K2PI	13	-733.9670	1545.0877	0.0000	1.0000
K2PG	13	-733.9903	1545.1343	0.0000	1.0000
GTRI	20	-714.4116	1547.5212	0.0000	1.0000
GTRG	20	-714.4306	1547.5590	0.0000	1.0000
JC69	11	-741.5693	1548.4223	0.0000	1.0000
K2PIG	14	-733.9996	1551.0878	0.0000	1.0000
JC69I	12	-740.3978	1552.0143	0.0000	1.0000
JC69G	12	-740.4278	1552.0743	0.0000	1.0000
GTRIG	21	-714.4331	1553.4990	0.0000	1.0000
JC69IG	13	-740.4375	1558.0287	0.0000	1.0000
SYM	16	-734.2500	1563.4583	0.0000	1.0000
SYMI	17	-733.1377	1567.1685	0.0000	1.0000
SYMG	17	-733.1625	1567.2183	0.0000	1.0000
SYMIG	18	-733.1625	1573.1531	0.0000	1.0000



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

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

BEGIN MRBAYES;
 Charset 591-tomodel-phylip/chr3_9719.phylip = 1 - 378;
 Partition Dummy = 1:591-tomodel-phylip/chr3_9719.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=2 rates=equal;
 Prset applyto=(1) statefreqpr=Dirichlet(1.0,1.0,1.0,1.0) tratiopr=Beta(1.0,1.0);
END;


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

BEGIN MRBAYES;
 Charset 591-tomodel-phylip/chr3_9719.phylip = 1 - 378;
 Partition Dummy = 1:591-tomodel-phylip/chr3_9719.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=2 rates=equal;
 Prset applyto=(1) statefreqpr=Dirichlet(1.0,1.0,1.0,1.0) tratiopr=Beta(1.0,1.0);
END;


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

BEGIN MRBAYES;
 Charset 591-tomodel-phylip/chr3_9719.phylip = 1 - 378;
 Partition Dummy = 1:591-tomodel-phylip/chr3_9719.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=2 rates=equal;
 Prset applyto=(1) statefreqpr=Dirichlet(1.0,1.0,1.0,1.0) tratiopr=Beta(1.0,1.0);
END;


-------------------------------------------------------------
End of Output

