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

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

Minimum AIC  model: SYMI
Minimum AICc model: SYMI
Minimum BIC  model: K2PI


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
SYMI	17	-816.6701	1669.4433	0.2224	0.2224
K2PI	13	-821.1181	1669.4702	0.2194	0.4418
SYMG	17	-816.8709	1669.8450	0.1819	0.6237
K2PG	13	-821.3369	1669.9077	0.1763	0.8001
K2PIG	14	-821.3369	1672.1024	0.0588	0.8589
SYMIG	18	-816.9185	1672.1956	0.0562	0.9151
HKYI	16	-820.0160	1673.8951	0.0240	0.9391
HKYG	16	-820.2138	1674.2906	0.0197	0.9588
GTRI	20	-815.9328	1674.7822	0.0154	0.9742
GTRG	20	-816.1161	1675.1489	0.0128	0.9870
HKYIG	17	-820.2659	1676.6350	0.0061	0.9931
GTRIG	21	-816.1593	1677.5382	0.0039	0.9970
SYM	16	-822.7997	1679.4623	0.0015	0.9985
K2P	12	-827.3758	1679.8057	0.0013	0.9997
HKY	15	-826.2988	1684.2358	0.0001	0.9999
GTR	19	-822.0083	1684.6463	0.0001	1.0000
JC69I	12	-833.2800	1691.6140	0.0000	1.0000
JC69G	12	-833.5189	1692.0919	0.0000	1.0000
JC69IG	13	-833.5189	1694.2717	0.0000	1.0000
F81I	15	-832.2917	1696.2217	0.0000	1.0000
F81G	15	-832.5193	1696.6767	0.0000	1.0000
F81IG	16	-832.5192	1698.9015	0.0000	1.0000
JC69	11	-839.4722	1701.8333	0.0000	1.0000
F81	14	-838.4953	1706.4192	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
SYMI	17	-816.6701	1667.3402	0.2578	0.2578
SYMG	17	-816.8709	1667.7419	0.2109	0.4686
K2PI	13	-821.1181	1668.2363	0.1647	0.6333
K2PG	13	-821.3369	1668.6738	0.1323	0.7656
SYMIG	18	-816.9185	1669.8369	0.0740	0.8396
K2PIG	14	-821.3369	1670.6738	0.0487	0.8883
GTRI	20	-815.9328	1671.8656	0.0268	0.9151
HKYI	16	-820.0160	1672.0320	0.0247	0.9398
GTRG	20	-816.1161	1672.2322	0.0223	0.9621
HKYG	16	-820.2138	1672.4276	0.0203	0.9824
GTRIG	21	-816.1593	1674.3187	0.0079	0.9902
HKYIG	17	-820.2659	1674.5319	0.0071	0.9973
SYM	16	-822.7997	1677.5993	0.0015	0.9988
K2P	12	-827.3758	1678.7516	0.0009	0.9997
GTR	19	-822.0083	1682.0166	0.0002	0.9999
HKY	15	-826.2988	1682.5975	0.0001	1.0000
JC69I	12	-833.2800	1690.5599	0.0000	1.0000
JC69G	12	-833.5189	1691.0378	0.0000	1.0000
JC69IG	13	-833.5189	1693.0378	0.0000	1.0000
F81I	15	-832.2917	1694.5834	0.0000	1.0000
F81G	15	-832.5193	1695.0385	0.0000	1.0000
F81IG	16	-832.5192	1697.0384	0.0000	1.0000
JC69	11	-839.4722	1700.9444	0.0000	1.0000
F81	14	-838.4953	1704.9907	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
K2PI	13	-821.1181	1716.7697	0.5303	0.5303
K2PG	13	-821.3369	1717.2072	0.4261	0.9564
K2PIG	14	-821.3369	1722.9406	0.0242	0.9806
K2P	12	-827.3758	1723.5517	0.0179	0.9984
SYMI	17	-816.6701	1730.8070	0.0005	0.9989
SYMG	17	-816.8709	1731.2087	0.0004	0.9993
HKYI	16	-820.0160	1731.7655	0.0003	0.9996
HKYG	16	-820.2138	1732.1610	0.0002	0.9998
JC69I	12	-833.2800	1735.3600	0.0000	0.9999
JC69G	12	-833.5189	1735.8379	0.0000	0.9999
SYMIG	18	-816.9185	1737.0371	0.0000	1.0000
SYM	16	-822.7997	1737.3328	0.0000	1.0000
HKYIG	17	-820.2659	1737.9987	0.0000	1.0000
HKY	15	-826.2988	1738.5977	0.0000	1.0000
JC69IG	13	-833.5189	1741.5713	0.0000	1.0000
JC69	11	-839.4722	1742.0111	0.0000	1.0000
GTRI	20	-815.9328	1746.5324	0.0000	1.0000
GTRG	20	-816.1161	1746.8990	0.0000	1.0000
F81I	15	-832.2917	1750.5836	0.0000	1.0000
F81G	15	-832.5193	1751.0386	0.0000	1.0000
GTRIG	21	-816.1593	1752.7188	0.0000	1.0000
GTR	19	-822.0083	1752.9500	0.0000	1.0000
F81IG	16	-832.5192	1756.7719	0.0000	1.0000
F81	14	-838.4953	1757.2574	0.0000	1.0000



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

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

BEGIN MRBAYES;
 Charset 591-tomodel-phylip/chr5_3302.phylip = 1 - 309;
 Partition Dummy = 1:591-tomodel-phylip/chr5_3302.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=Fixed(Equal) pinvarpr=Uniform(0.0,1.0);
END;


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

BEGIN MRBAYES;
 Charset 591-tomodel-phylip/chr5_3302.phylip = 1 - 309;
 Partition Dummy = 1:591-tomodel-phylip/chr5_3302.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=Fixed(Equal) pinvarpr=Uniform(0.0,1.0);
END;


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

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


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

