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

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

Minimum AIC  model: HKYG
Minimum AICc model: HKYG
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
HKYG	16	-701.0765	1436.2058	0.3195	0.3195
HKYI	16	-701.3874	1436.8276	0.2341	0.5536
K2PG	13	-705.3258	1438.0099	0.1296	0.6832
HKYIG	17	-701.0042	1438.3265	0.1106	0.7939
K2PI	13	-705.6276	1438.6134	0.0959	0.8897
K2PIG	14	-705.2476	1440.0681	0.0463	0.9361
GTRG	20	-699.5204	1442.2593	0.0155	0.9515
HKY	15	-705.5085	1442.8216	0.0117	0.9632
GTRI	20	-699.8081	1442.8345	0.0116	0.9749
SYMG	17	-703.7427	1443.8035	0.0072	0.9820
SYMI	17	-704.0444	1444.4069	0.0053	0.9873
GTRIG	21	-699.4595	1444.4729	0.0051	0.9924
K2P	12	-709.8351	1444.8300	0.0043	0.9967
SYMIG	18	-703.6661	1445.9330	0.0025	0.9992
GTR	19	-703.8575	1448.6157	0.0006	0.9998
SYM	16	-708.9140	1451.8809	0.0001	0.9999
F81G	15	-712.2063	1456.2172	0.0000	1.0000
JC69G	12	-715.7001	1456.5600	0.0000	1.0000
F81I	15	-712.4859	1456.7763	0.0000	1.0000
JC69I	12	-715.9826	1457.1251	0.0000	1.0000
F81IG	16	-712.1421	1458.3369	0.0000	1.0000
JC69IG	13	-715.6336	1458.6254	0.0000	1.0000
F81	14	-716.4183	1462.4096	0.0000	1.0000
JC69	11	-719.9479	1462.8735	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
HKYG	16	-701.0765	1434.1529	0.3305	0.3305
HKYI	16	-701.3874	1434.7748	0.2422	0.5726
HKYIG	17	-701.0042	1436.0083	0.1307	0.7033
K2PG	13	-705.3258	1436.6517	0.0947	0.7981
K2PI	13	-705.6276	1437.2552	0.0701	0.8681
K2PIG	14	-705.2476	1438.4951	0.0377	0.9058
GTRG	20	-699.5204	1439.0409	0.0287	0.9345
GTRI	20	-699.8081	1439.6161	0.0215	0.9560
GTRIG	21	-699.4595	1440.9190	0.0112	0.9672
HKY	15	-705.5085	1441.0170	0.0107	0.9779
SYMG	17	-703.7427	1441.4853	0.0085	0.9864
SYMI	17	-704.0444	1442.0887	0.0063	0.9926
SYMIG	18	-703.6661	1443.3323	0.0034	0.9960
K2P	12	-709.8351	1443.6701	0.0028	0.9988
GTR	19	-703.8575	1445.7149	0.0010	0.9998
SYM	16	-708.9140	1449.8280	0.0001	1.0000
F81G	15	-712.2063	1454.4126	0.0000	1.0000
F81I	15	-712.4859	1454.9718	0.0000	1.0000
JC69G	12	-715.7001	1455.4001	0.0000	1.0000
JC69I	12	-715.9826	1455.9653	0.0000	1.0000
F81IG	16	-712.1421	1456.2841	0.0000	1.0000
JC69IG	13	-715.6336	1457.2672	0.0000	1.0000
F81	14	-716.4183	1460.8365	0.0000	1.0000
JC69	11	-719.9479	1461.8958	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
K2PG	13	-705.3258	1483.9965	0.4952	0.4952
K2PI	13	-705.6276	1484.6000	0.3662	0.8614
K2P	12	-709.8351	1487.3730	0.0915	0.9529
K2PIG	14	-705.2476	1489.4818	0.0319	0.9848
HKYG	16	-701.0765	1492.4234	0.0073	0.9921
HKYI	16	-701.3874	1493.0453	0.0054	0.9975
HKY	15	-705.5085	1495.6456	0.0015	0.9989
HKYIG	17	-701.0042	1497.9208	0.0005	0.9994
JC69G	12	-715.7001	1499.1030	0.0003	0.9997
JC69I	12	-715.9826	1499.6681	0.0002	0.9999
JC69	11	-719.9479	1501.9567	0.0001	0.9999
SYMG	17	-703.7427	1503.3977	0.0000	1.0000
SYMI	17	-704.0444	1504.0012	0.0000	1.0000
JC69IG	13	-715.6336	1504.6120	0.0000	1.0000
SYM	16	-708.9140	1508.0985	0.0000	1.0000
SYMIG	18	-703.6661	1508.8866	0.0000	1.0000
F81G	15	-712.2063	1509.0412	0.0000	1.0000
F81I	15	-712.4859	1509.6004	0.0000	1.0000
F81	14	-716.4183	1511.8232	0.0000	1.0000
GTRG	20	-699.5204	1511.8790	0.0000	1.0000
GTRI	20	-699.8081	1512.4543	0.0000	1.0000
F81IG	16	-712.1421	1514.5546	0.0000	1.0000
GTR	19	-703.8575	1514.9112	0.0000	1.0000
GTRIG	21	-699.4595	1517.3991	0.0000	1.0000



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

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

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


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

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


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

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

