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

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

Minimum AIC  model: GTRI
Minimum AICc model: GTRI
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
GTRI	20	-646.5457	1335.9007	0.4284	0.4284
GTRG	20	-646.5930	1335.9953	0.4086	0.8370
GTRIG	21	-646.5808	1338.2623	0.1315	0.9686
GTR	19	-650.4315	1341.3963	0.0274	0.9960
HKYG	16	-656.6032	1347.0017	0.0017	0.9977
HKYI	16	-656.6443	1347.0840	0.0016	0.9993
HKYIG	17	-656.5715	1349.1694	0.0006	0.9999
HKY	15	-660.1764	1351.9317	0.0001	1.0000
F81G	15	-664.6680	1360.9149	0.0000	1.0000
F81I	15	-664.7157	1361.0103	0.0000	1.0000
F81IG	16	-664.6358	1363.0670	0.0000	1.0000
F81	14	-668.1237	1365.6245	0.0000	1.0000
SYMG	17	-671.4582	1378.9429	0.0000	1.0000
SYMI	17	-671.4875	1379.0014	0.0000	1.0000
SYMIG	18	-671.4410	1381.1544	0.0000	1.0000
K2PG	13	-677.5447	1382.2789	0.0000	1.0000
K2PI	13	-677.5948	1382.3792	0.0000	1.0000
SYM	16	-674.6952	1383.1858	0.0000	1.0000
K2PIG	14	-677.4981	1384.3733	0.0000	1.0000
K2P	12	-681.1190	1387.2543	0.0000	1.0000
JC69G	12	-684.4385	1393.8933	0.0000	1.0000
JC69I	12	-684.4995	1394.0153	0.0000	1.0000
JC69IG	13	-684.3835	1395.9566	0.0000	1.0000
JC69	11	-687.9413	1398.7397	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRI	20	-646.5457	1333.0914	0.4219	0.4219
GTRG	20	-646.5930	1333.1859	0.4024	0.8242
GTRIG	21	-646.5808	1335.1616	0.1498	0.9741
GTR	19	-650.4315	1338.8629	0.0235	0.9976
HKYG	16	-656.6032	1345.2063	0.0010	0.9986
HKYI	16	-656.6443	1345.2886	0.0009	0.9995
HKYIG	17	-656.5715	1347.1429	0.0004	0.9999
HKY	15	-660.1764	1350.3527	0.0001	1.0000
F81G	15	-664.6680	1359.3359	0.0000	1.0000
F81I	15	-664.7157	1359.4314	0.0000	1.0000
F81IG	16	-664.6358	1361.2716	0.0000	1.0000
F81	14	-668.1237	1364.2475	0.0000	1.0000
SYMG	17	-671.4582	1376.9164	0.0000	1.0000
SYMI	17	-671.4875	1376.9749	0.0000	1.0000
SYMIG	18	-671.4410	1378.8820	0.0000	1.0000
K2PG	13	-677.5447	1381.0894	0.0000	1.0000
K2PI	13	-677.5948	1381.1896	0.0000	1.0000
SYM	16	-674.6952	1381.3904	0.0000	1.0000
K2PIG	14	-677.4981	1382.9962	0.0000	1.0000
K2P	12	-681.1190	1386.2380	0.0000	1.0000
JC69G	12	-684.4385	1392.8770	0.0000	1.0000
JC69I	12	-684.4995	1392.9990	0.0000	1.0000
JC69IG	13	-684.3835	1394.7671	0.0000	1.0000
JC69	11	-687.9413	1397.8826	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	16	-656.6032	1405.4995	0.3252	0.3252
HKYI	16	-656.6443	1405.5818	0.3121	0.6373
HKY	15	-660.1764	1406.8775	0.1633	0.8006
GTRI	20	-646.5457	1408.4578	0.0741	0.8746
GTRG	20	-646.5930	1408.5523	0.0707	0.9453
GTR	19	-650.4315	1410.4610	0.0272	0.9725
HKYIG	17	-656.5715	1411.2044	0.0188	0.9913
GTRIG	21	-646.5808	1414.2963	0.0040	0.9953
F81G	15	-664.6680	1415.8608	0.0018	0.9971
F81I	15	-664.7157	1415.9562	0.0017	0.9989
F81	14	-668.1237	1417.0040	0.0010	0.9999
F81IG	16	-664.6358	1421.5648	0.0001	1.0000
K2PG	13	-677.5447	1430.0775	0.0000	1.0000
K2PI	13	-677.5948	1430.1778	0.0000	1.0000
K2P	12	-681.1190	1431.4579	0.0000	1.0000
K2PIG	14	-677.4981	1435.7527	0.0000	1.0000
JC69G	12	-684.4385	1438.0969	0.0000	1.0000
JC69I	12	-684.4995	1438.2189	0.0000	1.0000
JC69	11	-687.9413	1439.3341	0.0000	1.0000
SYMG	17	-671.4582	1440.9779	0.0000	1.0000
SYMI	17	-671.4875	1441.0364	0.0000	1.0000
SYM	16	-674.6952	1441.6836	0.0000	1.0000
JC69IG	13	-684.3835	1443.7552	0.0000	1.0000
SYMIG	18	-671.4410	1446.7118	0.0000	1.0000



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

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

BEGIN MRBAYES;
 Charset 591-tomodel-phylip/chr9_5190.phylip = 1 - 320;
 Partition Dummy = 1:591-tomodel-phylip/chr9_5190.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=6 rates=propinv;
 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) pinvarpr=Uniform(0.0,1.0);
END;


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

BEGIN MRBAYES;
 Charset 591-tomodel-phylip/chr9_5190.phylip = 1 - 320;
 Partition Dummy = 1:591-tomodel-phylip/chr9_5190.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=6 rates=propinv;
 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) pinvarpr=Uniform(0.0,1.0);
END;


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

BEGIN MRBAYES;
 Charset 591-tomodel-phylip/chr9_5190.phylip = 1 - 320;
 Partition Dummy = 1:591-tomodel-phylip/chr9_5190.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

