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

Input data from file "591-tomodel-phylip/chr4_8796.phylip" (7 taxa, 578 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	-2217.5453	4476.5987	0.5722	0.5722
GTRIG	21	-2217.5490	4478.7599	0.1942	0.7664
GTRG	20	-2219.1738	4479.8556	0.1123	0.8787
HKYI	16	-2223.7257	4480.4212	0.0846	0.9633
HKYIG	17	-2223.8110	4482.7149	0.0269	0.9902
HKYG	16	-2225.8878	4484.7453	0.0097	1.0000
GTR	19	-2228.0803	4495.5226	0.0000	1.0000
HKY	15	-2237.3565	4505.5670	0.0000	1.0000
SYMI	17	-2255.8853	4546.8634	0.0000	1.0000
SYMIG	18	-2255.8279	4548.8794	0.0000	1.0000
SYMG	17	-2257.6979	4550.4886	0.0000	1.0000
K2PI	13	-2265.4923	4557.6301	0.0000	1.0000
K2PIG	14	-2265.5313	4559.8086	0.0000	1.0000
K2PG	13	-2267.6170	4561.8794	0.0000	1.0000
SYM	16	-2266.8059	4566.5814	0.0000	1.0000
K2P	12	-2277.3957	4579.3435	0.0000	1.0000
F81I	15	-2279.4214	4589.6970	0.0000	1.0000
F81IG	16	-2279.4554	4591.8805	0.0000	1.0000
F81G	15	-2281.1358	4593.1256	0.0000	1.0000
F81	14	-2290.1649	4609.0758	0.0000	1.0000
JC69I	12	-2312.2306	4649.0134	0.0000	1.0000
JC69IG	13	-2312.2708	4651.1870	0.0000	1.0000
JC69G	12	-2314.0631	4652.6783	0.0000	1.0000
JC69	11	-2322.8953	4668.2571	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRI	20	-2217.5453	4475.0906	0.5790	0.5790
GTRIG	21	-2217.5490	4477.0980	0.2122	0.7913
GTRG	20	-2219.1738	4478.3476	0.1136	0.9049
HKYI	16	-2223.7257	4479.4515	0.0654	0.9703
HKYIG	17	-2223.8110	4481.6221	0.0221	0.9924
HKYG	16	-2225.8878	4483.7756	0.0075	1.0000
GTR	19	-2228.0803	4494.1606	0.0000	1.0000
HKY	15	-2237.3565	4504.7129	0.0000	1.0000
SYMI	17	-2255.8853	4545.7705	0.0000	1.0000
SYMIG	18	-2255.8279	4547.6557	0.0000	1.0000
SYMG	17	-2257.6979	4549.3958	0.0000	1.0000
K2PI	13	-2265.4923	4556.9847	0.0000	1.0000
K2PIG	14	-2265.5313	4559.0626	0.0000	1.0000
K2PG	13	-2267.6170	4561.2340	0.0000	1.0000
SYM	16	-2266.8059	4565.6117	0.0000	1.0000
K2P	12	-2277.3957	4578.7913	0.0000	1.0000
F81I	15	-2279.4214	4588.8429	0.0000	1.0000
F81IG	16	-2279.4554	4590.9108	0.0000	1.0000
F81G	15	-2281.1358	4592.2715	0.0000	1.0000
F81	14	-2290.1649	4608.3298	0.0000	1.0000
JC69I	12	-2312.2306	4648.4611	0.0000	1.0000
JC69IG	13	-2312.2708	4650.5416	0.0000	1.0000
JC69G	12	-2314.0631	4652.1261	0.0000	1.0000
JC69	11	-2322.8953	4667.7906	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYI	16	-2223.7257	4549.2046	0.8657	0.8657
HKYG	16	-2225.8878	4553.5287	0.0996	0.9654
HKYIG	17	-2223.8110	4555.7348	0.0331	0.9984
GTRI	20	-2217.5453	4562.2821	0.0013	0.9997
GTRG	20	-2219.1738	4565.5390	0.0002	0.9999
GTRIG	21	-2217.5490	4568.6490	0.0001	1.0000
HKY	15	-2237.3565	4570.1065	0.0000	1.0000
GTR	19	-2228.0803	4576.9925	0.0000	1.0000
K2PI	13	-2265.4923	4613.6592	0.0000	1.0000
K2PG	13	-2267.6170	4617.9085	0.0000	1.0000
SYMI	17	-2255.8853	4619.8833	0.0000	1.0000
K2PIG	14	-2265.5313	4620.0966	0.0000	1.0000
SYMG	17	-2257.6979	4623.5085	0.0000	1.0000
SYMIG	18	-2255.8279	4626.1281	0.0000	1.0000
K2P	12	-2277.3957	4631.1062	0.0000	1.0000
SYM	16	-2266.8059	4635.3649	0.0000	1.0000
F81I	15	-2279.4214	4654.2365	0.0000	1.0000
F81G	15	-2281.1358	4657.6651	0.0000	1.0000
F81IG	16	-2279.4554	4660.6640	0.0000	1.0000
F81	14	-2290.1649	4669.3639	0.0000	1.0000
JC69I	12	-2312.2306	4700.7760	0.0000	1.0000
JC69G	12	-2314.0631	4704.4410	0.0000	1.0000
JC69IG	13	-2312.2708	4707.2160	0.0000	1.0000
JC69	11	-2322.8953	4715.7459	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_8796.phylip = 1 - 578;
 Partition Dummy = 1:591-tomodel-phylip/chr4_8796.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_8796.phylip = 1 - 578;
 Partition Dummy = 1:591-tomodel-phylip/chr4_8796.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_8796.phylip = 1 - 578;
 Partition Dummy = 1:591-tomodel-phylip/chr4_8796.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

