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

Input data from file "591-tomodel-phylip/chr7_10536.phylip" (7 taxa, 683 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	-1359.7880	2760.8449	0.2490	0.2490
GTRG	20	-1359.8131	2760.8950	0.2428	0.4918
HKYI	16	-1364.3354	2761.4877	0.1805	0.6723
HKYG	16	-1364.3360	2761.4888	0.1804	0.8528
GTRIG	21	-1359.8131	2763.0240	0.0837	0.9365
HKYIG	17	-1364.3359	2763.5921	0.0630	0.9996
GTR	19	-1367.6009	2774.3482	0.0003	0.9999
HKY	15	-1372.5908	2775.9011	0.0001	1.0000
F81I	15	-1377.5679	2785.8554	0.0000	1.0000
F81G	15	-1377.5949	2785.9094	0.0000	1.0000
F81IG	16	-1377.5949	2788.0066	0.0000	1.0000
F81	14	-1385.2313	2799.0914	0.0000	1.0000
SYMI	17	-1385.5317	2805.9837	0.0000	1.0000
SYMG	17	-1385.5525	2806.0253	0.0000	1.0000
K2PI	13	-1390.2363	2807.0167	0.0000	1.0000
K2PG	13	-1390.2408	2807.0256	0.0000	1.0000
SYMIG	18	-1385.5525	2808.1351	0.0000	1.0000
K2PIG	14	-1390.2407	2809.1102	0.0000	1.0000
SYM	16	-1393.3118	2819.4405	0.0000	1.0000
K2P	12	-1398.2470	2820.9596	0.0000	1.0000
JC69I	12	-1402.4097	2829.2851	0.0000	1.0000
JC69G	12	-1402.4134	2829.2926	0.0000	1.0000
JC69IG	13	-1402.4167	2831.3774	0.0000	1.0000
JC69	11	-1410.0658	2842.5250	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRI	20	-1359.7880	2759.5760	0.2699	0.2699
GTRG	20	-1359.8131	2759.6261	0.2632	0.5332
HKYI	16	-1364.3354	2760.6708	0.1561	0.6893
HKYG	16	-1364.3360	2760.6720	0.1560	0.8453
GTRIG	21	-1359.8131	2761.6261	0.0968	0.9422
HKYIG	17	-1364.3359	2762.6718	0.0574	0.9996
GTR	19	-1367.6009	2773.2019	0.0003	0.9999
HKY	15	-1372.5908	2775.1815	0.0001	1.0000
F81I	15	-1377.5679	2785.1357	0.0000	1.0000
F81G	15	-1377.5949	2785.1897	0.0000	1.0000
F81IG	16	-1377.5949	2787.1898	0.0000	1.0000
F81	14	-1385.2313	2798.4626	0.0000	1.0000
SYMI	17	-1385.5317	2805.0634	0.0000	1.0000
SYMG	17	-1385.5525	2805.1050	0.0000	1.0000
K2PI	13	-1390.2363	2806.4726	0.0000	1.0000
K2PG	13	-1390.2408	2806.4815	0.0000	1.0000
SYMIG	18	-1385.5525	2807.1050	0.0000	1.0000
K2PIG	14	-1390.2407	2808.4815	0.0000	1.0000
SYM	16	-1393.3118	2818.6237	0.0000	1.0000
K2P	12	-1398.2470	2820.4940	0.0000	1.0000
JC69I	12	-1402.4097	2828.8194	0.0000	1.0000
JC69G	12	-1402.4134	2828.8269	0.0000	1.0000
JC69IG	13	-1402.4167	2830.8333	0.0000	1.0000
JC69	11	-1410.0658	2842.1316	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYI	16	-1364.3354	2833.0948	0.4890	0.4890
HKYG	16	-1364.3360	2833.0959	0.4887	0.9777
HKYIG	17	-1364.3359	2839.6223	0.0187	0.9964
HKY	15	-1372.5908	2843.0789	0.0033	0.9998
GTRI	20	-1359.7880	2850.1059	0.0001	0.9999
GTRG	20	-1359.8131	2850.1560	0.0001	0.9999
F81I	15	-1377.5679	2853.0331	0.0000	1.0000
F81G	15	-1377.5949	2853.0871	0.0000	1.0000
GTRIG	21	-1359.8131	2856.6825	0.0000	1.0000
GTR	19	-1367.6009	2859.2053	0.0000	1.0000
F81IG	16	-1377.5949	2859.6137	0.0000	1.0000
F81	14	-1385.2313	2861.8335	0.0000	1.0000
K2PI	13	-1390.2363	2865.3170	0.0000	1.0000
K2PG	13	-1390.2408	2865.3260	0.0000	1.0000
K2PIG	14	-1390.2407	2871.8524	0.0000	1.0000
K2P	12	-1398.2470	2874.8119	0.0000	1.0000
SYMI	17	-1385.5317	2882.0138	0.0000	1.0000
SYMG	17	-1385.5525	2882.0554	0.0000	1.0000
JC69I	12	-1402.4097	2883.1373	0.0000	1.0000
JC69G	12	-1402.4134	2883.1448	0.0000	1.0000
SYMIG	18	-1385.5525	2888.5819	0.0000	1.0000
JC69IG	13	-1402.4167	2889.6778	0.0000	1.0000
SYM	16	-1393.3118	2891.0476	0.0000	1.0000
JC69	11	-1410.0658	2891.9230	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/chr7_10536.phylip = 1 - 683;
 Partition Dummy = 1:591-tomodel-phylip/chr7_10536.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/chr7_10536.phylip = 1 - 683;
 Partition Dummy = 1:591-tomodel-phylip/chr7_10536.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/chr7_10536.phylip = 1 - 683;
 Partition Dummy = 1:591-tomodel-phylip/chr7_10536.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

