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

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

Minimum AIC  model: GTRI
Minimum AICc model: GTRI
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
GTRI	20	-693.2311	1429.0390	0.2360	0.2360
HKYI	16	-697.6956	1429.0398	0.2359	0.4720
GTRG	20	-693.4260	1429.4287	0.1942	0.6662
HKYG	16	-698.0325	1429.7134	0.1685	0.8346
GTRIG	21	-693.0939	1431.0308	0.0872	0.9218
HKYIG	17	-697.6957	1431.2516	0.0781	0.9999
GTR	19	-702.7279	1445.7800	0.0001	1.0000
HKY	15	-707.4842	1446.4186	0.0000	1.0000
F81I	15	-710.2936	1452.0375	0.0000	1.0000
F81G	15	-710.7666	1452.9833	0.0000	1.0000
F81IG	16	-709.8730	1453.3944	0.0000	1.0000
SYMI	17	-713.3238	1462.5077	0.0000	1.0000
K2PI	13	-717.7278	1462.5487	0.0000	1.0000
K2PG	13	-718.4924	1464.0780	0.0000	1.0000
SYMG	17	-714.2270	1464.3142	0.0000	1.0000
K2PIG	14	-717.7278	1464.7206	0.0000	1.0000
SYMIG	18	-713.3237	1464.7328	0.0000	1.0000
F81	14	-719.8813	1469.0276	0.0000	1.0000
K2P	12	-728.7455	1482.4251	0.0000	1.0000
SYM	16	-724.7663	1483.1810	0.0000	1.0000
JC69I	12	-729.5953	1484.1248	0.0000	1.0000
JC69IG	13	-729.2655	1485.6240	0.0000	1.0000
JC69G	12	-730.3958	1485.7257	0.0000	1.0000
JC69	11	-739.9815	1502.7510	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRI	20	-693.2311	1426.4623	0.2814	0.2814
GTRG	20	-693.4260	1426.8520	0.2316	0.5130
HKYI	16	-697.6956	1427.3913	0.1768	0.6898
HKYG	16	-698.0325	1428.0650	0.1263	0.8161
GTRIG	21	-693.0939	1428.1877	0.1188	0.9349
HKYIG	17	-697.6957	1429.3915	0.0651	0.9999
GTR	19	-702.7279	1443.4558	0.0001	1.0000
HKY	15	-707.4842	1444.9685	0.0000	1.0000
F81I	15	-710.2936	1450.5873	0.0000	1.0000
F81G	15	-710.7666	1451.5332	0.0000	1.0000
F81IG	16	-709.8730	1451.7459	0.0000	1.0000
SYMI	17	-713.3238	1460.6476	0.0000	1.0000
K2PI	13	-717.7278	1461.4556	0.0000	1.0000
SYMG	17	-714.2270	1462.4540	0.0000	1.0000
SYMIG	18	-713.3237	1462.6474	0.0000	1.0000
K2PG	13	-718.4924	1462.9849	0.0000	1.0000
K2PIG	14	-717.7278	1463.4556	0.0000	1.0000
F81	14	-719.8813	1467.7626	0.0000	1.0000
K2P	12	-728.7455	1481.4910	0.0000	1.0000
SYM	16	-724.7663	1481.5325	0.0000	1.0000
JC69I	12	-729.5953	1483.1907	0.0000	1.0000
JC69IG	13	-729.2655	1484.5309	0.0000	1.0000
JC69G	12	-730.3958	1484.7916	0.0000	1.0000
JC69	11	-739.9815	1501.9629	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYI	16	-697.6956	1488.9805	0.5649	0.5649
HKYG	16	-698.0325	1489.6542	0.4033	0.9682
HKYIG	17	-697.6957	1494.8300	0.0303	0.9986
HKY	15	-707.4842	1502.7084	0.0006	0.9992
GTRI	20	-693.2311	1503.4488	0.0004	0.9996
GTRG	20	-693.4260	1503.8385	0.0003	0.9999
F81I	15	-710.2936	1508.3272	0.0000	0.9999
GTRIG	21	-693.0939	1509.0236	0.0000	1.0000
F81G	15	-710.7666	1509.2730	0.0000	1.0000
K2PI	13	-717.7278	1511.4968	0.0000	1.0000
K2PG	13	-718.4924	1513.0261	0.0000	1.0000
F81IG	16	-709.8730	1513.3351	0.0000	1.0000
GTR	19	-702.7279	1516.5930	0.0000	1.0000
K2PIG	14	-717.7278	1517.3461	0.0000	1.0000
F81	14	-719.8813	1521.6531	0.0000	1.0000
SYMI	17	-713.3238	1526.0861	0.0000	1.0000
K2P	12	-728.7455	1527.6829	0.0000	1.0000
SYMG	17	-714.2270	1527.8926	0.0000	1.0000
JC69I	12	-729.5953	1529.3826	0.0000	1.0000
JC69G	12	-730.3958	1530.9835	0.0000	1.0000
SYMIG	18	-713.3237	1531.9352	0.0000	1.0000
JC69IG	13	-729.2655	1534.5722	0.0000	1.0000
SYM	16	-724.7663	1543.1217	0.0000	1.0000
JC69	11	-739.9815	1544.3055	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/chr4_10598.phylip = 1 - 347;
 Partition Dummy = 1:591-tomodel-phylip/chr4_10598.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/chr4_10598.phylip = 1 - 347;
 Partition Dummy = 1:591-tomodel-phylip/chr4_10598.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 (HKYI)]

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

