
Output from mraic.pl version 1.4.4 by Johan Nylander
-------------------------------------------------------------
00:34 02/17/2011

Input data from file "183-phylip/chr11_1600.phylip" (29 taxa, 251 characters)

Minimum AIC  model: GTRI
Minimum AICc model: GTRI
Minimum BIC  model: K2PG


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	64	-678.0898	1528.9109	0.2710	0.2710
GTRG	64	-678.1793	1529.0897	0.2479	0.5189
HKYG	60	-685.7903	1530.1068	0.1491	0.6679
HKYI	60	-685.8007	1530.1278	0.1475	0.8154
K2PG	57	-692.1314	1532.5220	0.0446	0.8600
K2PI	57	-692.1593	1532.5777	0.0433	0.9033
GTRIG	65	-678.1787	1532.7358	0.0400	0.9434
HKYIG	61	-685.7942	1533.6096	0.0259	0.9692
SYMI	61	-686.6837	1535.3886	0.0106	0.9799
SYMG	61	-686.7820	1535.5851	0.0096	0.9895
K2PIG	58	-692.1353	1535.9163	0.0082	0.9976
SYMIG	62	-686.7814	1539.1160	0.0016	0.9993
GTR	63	-686.4213	1541.9655	0.0004	0.9997
HKY	59	-693.8627	1542.7934	0.0003	1.0000
K2P	56	-700.9385	1546.7842	0.0000	1.0000
SYM	60	-695.1857	1548.8976	0.0000	1.0000
F81I	59	-711.3319	1577.7319	0.0000	1.0000
F81G	59	-711.3456	1577.7592	0.0000	1.0000
JC69I	56	-716.8145	1578.5363	0.0000	1.0000
JC69G	56	-716.8226	1578.5525	0.0000	1.0000
F81IG	60	-711.3487	1581.2238	0.0000	1.0000
JC69IG	57	-716.8272	1581.9134	0.0000	1.0000
F81	58	-719.7121	1591.0701	0.0000	1.0000
JC69	55	-725.4583	1592.5064	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRI	64	-678.0898	1484.1797	0.4310	0.4310
GTRG	64	-678.1793	1484.3585	0.3941	0.8251
GTRIG	65	-678.1787	1486.3575	0.1451	0.9701
HKYG	60	-685.7903	1491.5805	0.0107	0.9808
HKYI	60	-685.8007	1491.6015	0.0105	0.9913
HKYIG	61	-685.7942	1493.5884	0.0039	0.9952
SYMI	61	-686.6837	1495.3675	0.0016	0.9968
SYMG	61	-686.7820	1495.5640	0.0015	0.9983
SYMIG	62	-686.7814	1497.5628	0.0005	0.9988
K2PG	57	-692.1314	1498.2629	0.0004	0.9992
K2PI	57	-692.1593	1498.3187	0.0004	0.9996
GTR	63	-686.4213	1498.8425	0.0003	0.9999
K2PIG	58	-692.1353	1500.2705	0.0001	1.0000
HKY	59	-693.8627	1505.7253	0.0000	1.0000
SYM	60	-695.1857	1510.3713	0.0000	1.0000
K2P	56	-700.9385	1513.8770	0.0000	1.0000
F81I	59	-711.3319	1540.6639	0.0000	1.0000
F81G	59	-711.3456	1540.6911	0.0000	1.0000
F81IG	60	-711.3487	1542.6975	0.0000	1.0000
JC69I	56	-716.8145	1545.6291	0.0000	1.0000
JC69G	56	-716.8226	1545.6453	0.0000	1.0000
JC69IG	57	-716.8272	1547.6544	0.0000	1.0000
F81	58	-719.7121	1555.4243	0.0000	1.0000
JC69	55	-725.4583	1560.9167	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
K2PG	57	-692.1314	1699.2137	0.4258	0.4258
K2PI	57	-692.1593	1699.2695	0.4141	0.8399
HKYG	60	-685.7903	1703.1077	0.0608	0.9007
HKYI	60	-685.8007	1703.1287	0.0601	0.9608
K2PIG	58	-692.1353	1704.7468	0.0268	0.9876
HKYIG	61	-685.7942	1708.6411	0.0038	0.9914
GTRI	64	-678.0898	1709.8087	0.0021	0.9935
GTRG	64	-678.1793	1709.9875	0.0019	0.9955
SYMI	61	-686.6837	1710.4201	0.0016	0.9970
SYMG	61	-686.7820	1710.6166	0.0014	0.9985
K2P	56	-700.9385	1711.3024	0.0010	0.9995
HKY	59	-693.8627	1713.7270	0.0003	0.9998
GTRIG	65	-678.1787	1715.5119	0.0001	0.9999
SYMIG	62	-686.7814	1716.1409	0.0001	1.0000
GTR	63	-686.4213	1720.9461	0.0000	1.0000
SYM	60	-695.1857	1721.8985	0.0000	1.0000
JC69I	56	-716.8145	1743.0544	0.0000	1.0000
JC69G	56	-716.8226	1743.0706	0.0000	1.0000
JC69IG	57	-716.8272	1748.6052	0.0000	1.0000
F81I	59	-711.3319	1748.6656	0.0000	1.0000
F81G	59	-711.3456	1748.6928	0.0000	1.0000
F81IG	60	-711.3487	1754.2247	0.0000	1.0000
JC69	55	-725.4583	1754.8166	0.0000	1.0000
F81	58	-719.7121	1759.9005	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 183-phylip/chr11_1600.phylip = 1 - 251;
 Partition Dummy = 1:183-phylip/chr11_1600.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 183-phylip/chr11_1600.phylip = 1 - 251;
 Partition Dummy = 1:183-phylip/chr11_1600.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 (K2PG)]

BEGIN MRBAYES;
 Charset 183-phylip/chr11_1600.phylip = 1 - 251;
 Partition Dummy = 1:183-phylip/chr11_1600.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=2 rates=gamma;
 Prset applyto=(1) statefreqpr=Fixed(Equal) tratiopr=Beta(1.0,1.0) shapepr=Uniform(0.1,50.0);
END;


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

