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

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

Minimum AIC  model: HKYI
Minimum AICc model: HKYI
Minimum BIC  model: HKYI


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
HKYI	16	-719.4633	1472.6109	0.3638	0.3638
HKYG	16	-719.4646	1472.6135	0.3634	0.7272
HKYIG	17	-719.4645	1474.8296	0.1200	0.8472
GTRG	20	-716.7369	1476.1070	0.0633	0.9105
GTRI	20	-716.7496	1476.1323	0.0626	0.9731
GTRIG	21	-716.6967	1478.2991	0.0212	0.9943
F81I	15	-726.2346	1483.9507	0.0013	0.9955
F81G	15	-726.2614	1484.0043	0.0012	0.9967
HKY	15	-726.3718	1484.2251	0.0011	0.9978
GTR	19	-722.2221	1484.8193	0.0008	0.9987
F81IG	16	-726.2615	1486.2071	0.0004	0.9991
K2PI	13	-729.6453	1486.4071	0.0004	0.9994
K2PG	13	-729.6517	1486.4200	0.0004	0.9998
K2PIG	14	-729.6516	1488.5956	0.0001	0.9999
SYMI	17	-727.8376	1491.5758	0.0000	0.9999
SYMG	17	-727.8649	1491.6304	0.0000	1.0000
F81	14	-731.8283	1492.9489	0.0000	1.0000
SYMIG	18	-727.8648	1493.8605	0.0000	1.0000
JC69I	12	-735.5969	1496.1480	0.0000	1.0000
JC69G	12	-735.6284	1496.2109	0.0000	1.0000
K2P	12	-736.3291	1497.6124	0.0000	1.0000
JC69IG	13	-735.6284	1498.3733	0.0000	1.0000
SYM	16	-733.0563	1499.7968	0.0000	1.0000
JC69	11	-741.1469	1505.0987	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
HKYI	16	-719.4633	1470.9267	0.3284	0.3284
HKYG	16	-719.4646	1470.9292	0.3280	0.6565
HKYIG	17	-719.4645	1472.9290	0.1207	0.7772
GTRG	20	-716.7369	1473.4738	0.0919	0.8691
GTRI	20	-716.7496	1473.4991	0.0908	0.9598
GTRIG	21	-716.6967	1475.3935	0.0352	0.9950
GTR	19	-722.2221	1482.4443	0.0010	0.9961
F81I	15	-726.2346	1482.4693	0.0010	0.9971
F81G	15	-726.2614	1482.5228	0.0010	0.9981
HKY	15	-726.3718	1482.7436	0.0009	0.9990
F81IG	16	-726.2615	1484.5229	0.0004	0.9993
K2PI	13	-729.6453	1485.2905	0.0002	0.9996
K2PG	13	-729.6517	1485.3034	0.0002	0.9998
K2PIG	14	-729.6516	1487.3033	0.0001	0.9999
SYMI	17	-727.8376	1489.6752	0.0000	0.9999
SYMG	17	-727.8649	1489.7298	0.0000	1.0000
F81	14	-731.8283	1491.6566	0.0000	1.0000
SYMIG	18	-727.8648	1491.7297	0.0000	1.0000
JC69I	12	-735.5969	1495.1939	0.0000	1.0000
JC69G	12	-735.6284	1495.2568	0.0000	1.0000
K2P	12	-736.3291	1496.6582	0.0000	1.0000
JC69IG	13	-735.6284	1497.2567	0.0000	1.0000
SYM	16	-733.0563	1498.1126	0.0000	1.0000
JC69	11	-741.1469	1504.2938	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYI	16	-719.4633	1532.1898	0.3801	0.3801
HKYG	16	-719.4646	1532.1924	0.3796	0.7597
K2PI	13	-729.6453	1535.0668	0.0902	0.8499
K2PG	13	-729.6517	1535.0797	0.0896	0.9395
HKYIG	17	-719.4645	1538.0211	0.0206	0.9601
F81I	15	-726.2346	1539.9034	0.0080	0.9681
F81G	15	-726.2614	1539.9570	0.0078	0.9759
HKY	15	-726.3718	1540.1778	0.0070	0.9829
K2PIG	14	-729.6516	1540.9085	0.0049	0.9878
JC69I	12	-735.5969	1541.1412	0.0043	0.9921
JC69G	12	-735.6284	1541.2041	0.0042	0.9963
K2P	12	-736.3291	1542.6056	0.0021	0.9984
F81	14	-731.8283	1545.2618	0.0006	0.9989
F81IG	16	-726.2615	1545.7860	0.0004	0.9993
JC69	11	-741.1469	1546.4122	0.0003	0.9997
JC69IG	13	-735.6284	1547.0330	0.0002	0.9999
GTRG	20	-716.7369	1550.0527	0.0001	0.9999
GTRI	20	-716.7496	1550.0780	0.0000	1.0000
SYMI	17	-727.8376	1554.7672	0.0000	1.0000
SYMG	17	-727.8649	1554.8219	0.0000	1.0000
GTR	19	-722.2221	1555.1942	0.0000	1.0000
GTRIG	21	-716.6967	1555.8013	0.0000	1.0000
SYM	16	-733.0563	1559.3757	0.0000	1.0000
SYMIG	18	-727.8648	1560.6507	0.0000	1.0000



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

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

BEGIN MRBAYES;
 Charset 591-tomodel-phylip/chr7_9143.phylip = 1 - 340;
 Partition Dummy = 1:591-tomodel-phylip/chr7_9143.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=2 rates=propinv;
 Prset applyto=(1) statefreqpr=Dirichlet(1.0,1.0,1.0,1.0) tratiopr=Beta(1.0,1.0) pinvarpr=Uniform(0.0,1.0);
END;


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

BEGIN MRBAYES;
 Charset 591-tomodel-phylip/chr7_9143.phylip = 1 - 340;
 Partition Dummy = 1:591-tomodel-phylip/chr7_9143.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=2 rates=propinv;
 Prset applyto=(1) statefreqpr=Dirichlet(1.0,1.0,1.0,1.0) tratiopr=Beta(1.0,1.0) pinvarpr=Uniform(0.0,1.0);
END;


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

BEGIN MRBAYES;
 Charset 591-tomodel-phylip/chr7_9143.phylip = 1 - 340;
 Partition Dummy = 1:591-tomodel-phylip/chr7_9143.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=2 rates=propinv;
 Prset applyto=(1) statefreqpr=Dirichlet(1.0,1.0,1.0,1.0) tratiopr=Beta(1.0,1.0) pinvarpr=Uniform(0.0,1.0);
END;


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

