
Output from mraic.pl version 1.4.4 by Johan Nylander
-------------------------------------------------------------
16:19 04/04/2011

Input data from file "683-tomodel-phylip/chrZ_10758.phylip" (7 taxa, 579 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	-1866.7071	3766.3822	0.2716	0.2716
K2PG	13	-1869.9857	3766.6157	0.2417	0.5133
GTRG	20	-1863.3050	3768.1153	0.1142	0.6275
HKYIG	17	-1866.6139	3768.3187	0.1032	0.7307
K2PIG	14	-1869.9055	3768.5557	0.0916	0.8223
GTRIG	21	-1863.1905	3770.0399	0.0436	0.8659
HKYI	16	-1868.5582	3770.0844	0.0427	0.9086
K2PI	13	-1871.8150	3770.2742	0.0388	0.9474
GTRI	20	-1864.9507	3771.4069	0.0220	0.9694
SYMG	17	-1868.2744	3771.6397	0.0196	0.9890
SYMIG	18	-1868.1909	3773.6031	0.0073	0.9964
SYMI	17	-1869.9619	3775.0148	0.0036	1.0000
K2P	12	-1894.4712	3813.4937	0.0000	1.0000
HKY	15	-1891.6243	3814.1013	0.0000	1.0000
GTR	19	-1889.1527	3817.6650	0.0000	1.0000
SYM	16	-1893.2422	3819.4523	0.0000	1.0000
F81G	15	-1900.7616	3832.3758	0.0000	1.0000
JC69G	12	-1904.0322	3832.6156	0.0000	1.0000
F81IG	16	-1900.7126	3834.3931	0.0000	1.0000
JC69IG	13	-1903.9902	3834.6246	0.0000	1.0000
F81I	15	-1902.3164	3835.4854	0.0000	1.0000
JC69I	12	-1905.5902	3835.7316	0.0000	1.0000
F81	14	-1923.3728	3875.4903	0.0000	1.0000
JC69	11	-1926.7619	3875.9893	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
HKYG	16	-1866.7071	3765.4142	0.2670	0.2670
K2PG	13	-1869.9857	3765.9714	0.2021	0.4691
GTRG	20	-1863.3050	3766.6099	0.1468	0.6159
HKYIG	17	-1866.6139	3767.2278	0.1078	0.7237
K2PIG	14	-1869.9055	3767.8110	0.0805	0.8043
GTRIG	21	-1863.1905	3768.3811	0.0606	0.8648
HKYI	16	-1868.5582	3769.1164	0.0419	0.9068
K2PI	13	-1871.8150	3769.6300	0.0324	0.9392
GTRI	20	-1864.9507	3769.9015	0.0283	0.9675
SYMG	17	-1868.2744	3770.5488	0.0205	0.9880
SYMIG	18	-1868.1909	3772.3817	0.0082	0.9962
SYMI	17	-1869.9619	3773.9239	0.0038	1.0000
K2P	12	-1894.4712	3812.9425	0.0000	1.0000
HKY	15	-1891.6243	3813.2487	0.0000	1.0000
GTR	19	-1889.1527	3816.3054	0.0000	1.0000
SYM	16	-1893.2422	3818.4843	0.0000	1.0000
F81G	15	-1900.7616	3831.5232	0.0000	1.0000
JC69G	12	-1904.0322	3832.0644	0.0000	1.0000
F81IG	16	-1900.7126	3833.4251	0.0000	1.0000
JC69IG	13	-1903.9902	3833.9804	0.0000	1.0000
F81I	15	-1902.3164	3834.6329	0.0000	1.0000
JC69I	12	-1905.5902	3835.1804	0.0000	1.0000
F81	14	-1923.3728	3874.7456	0.0000	1.0000
JC69	11	-1926.7619	3875.5237	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
K2PG	13	-1869.9857	3822.6684	0.8279	0.8279
K2PI	13	-1871.8150	3826.3269	0.1329	0.9608
K2PIG	14	-1869.9055	3828.8693	0.0373	0.9981
HKYG	16	-1866.7071	3835.1951	0.0016	0.9997
HKYI	16	-1868.5582	3838.8972	0.0002	0.9999
HKYIG	17	-1866.6139	3841.3699	0.0001	1.0000
SYMG	17	-1868.2744	3844.6910	0.0000	1.0000
SYMI	17	-1869.9619	3848.0660	0.0000	1.0000
SYMIG	18	-1868.1909	3850.8852	0.0000	1.0000
GTRG	20	-1863.3050	3853.8360	0.0000	1.0000
GTRI	20	-1864.9507	3857.1275	0.0000	1.0000
GTRIG	21	-1863.1905	3859.9684	0.0000	1.0000
K2P	12	-1894.4712	3865.2781	0.0000	1.0000
HKY	15	-1891.6243	3878.6682	0.0000	1.0000
JC69G	12	-1904.0322	3884.4000	0.0000	1.0000
JC69I	12	-1905.5902	3887.5160	0.0000	1.0000
SYM	16	-1893.2422	3888.2651	0.0000	1.0000
JC69IG	13	-1903.9902	3890.6773	0.0000	1.0000
F81G	15	-1900.7616	3896.9427	0.0000	1.0000
GTR	19	-1889.1527	3899.1701	0.0000	1.0000
F81I	15	-1902.3164	3900.0524	0.0000	1.0000
F81IG	16	-1900.7126	3903.2060	0.0000	1.0000
JC69	11	-1926.7619	3923.4980	0.0000	1.0000
F81	14	-1923.3728	3935.8038	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 683-tomodel-phylip/chrZ_10758.phylip = 1 - 579;
 Partition Dummy = 1:683-tomodel-phylip/chrZ_10758.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 683-tomodel-phylip/chrZ_10758.phylip = 1 - 579;
 Partition Dummy = 1:683-tomodel-phylip/chrZ_10758.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 683-tomodel-phylip/chrZ_10758.phylip = 1 - 579;
 Partition Dummy = 1:683-tomodel-phylip/chrZ_10758.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

