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

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

Minimum AIC  model: HKYI
Minimum AICc model: HKYI
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
HKYI	16	-1310.0842	2652.9926	0.4056	0.4056
HKYG	16	-1310.1193	2653.0629	0.3916	0.7972
HKYIG	17	-1310.1378	2655.2042	0.1342	0.9314
HKY	15	-1313.3740	2657.4743	0.0431	0.9746
GTRI	20	-1309.5070	2660.2944	0.0105	0.9851
GTRG	20	-1309.5318	2660.3440	0.0103	0.9954
GTRIG	21	-1309.5454	2662.5014	0.0035	0.9989
GTR	19	-1312.8337	2664.8242	0.0011	1.0000
F81I	15	-1321.9565	2674.6391	0.0000	1.0000
F81G	15	-1321.9934	2674.7130	0.0000	1.0000
F81IG	16	-1321.9933	2676.8109	0.0000	1.0000
F81	14	-1325.1804	2678.9953	0.0000	1.0000
K2PI	13	-1358.6562	2743.8614	0.0000	1.0000
K2PG	13	-1358.7131	2743.9752	0.0000	1.0000
K2PIG	14	-1358.7401	2746.1147	0.0000	1.0000
K2P	12	-1361.8202	2748.1104	0.0000	1.0000
SYMI	17	-1357.9678	2750.8642	0.0000	1.0000
SYMG	17	-1358.0240	2750.9768	0.0000	1.0000
SYMIG	18	-1358.0256	2753.0907	0.0000	1.0000
SYM	16	-1361.0798	2754.9838	0.0000	1.0000
JC69I	12	-1367.8535	2760.1768	0.0000	1.0000
JC69G	12	-1367.9010	2760.2718	0.0000	1.0000
JC69IG	13	-1367.9245	2762.3980	0.0000	1.0000
JC69	11	-1371.0388	2764.4746	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
HKYI	16	-1310.0842	2652.1683	0.4009	0.4009
HKYG	16	-1310.1193	2652.2387	0.3870	0.7879
HKYIG	17	-1310.1378	2654.2756	0.1398	0.9276
HKY	15	-1313.3740	2656.7481	0.0406	0.9682
GTRI	20	-1309.5070	2659.0139	0.0131	0.9813
GTRG	20	-1309.5318	2659.0635	0.0128	0.9941
GTRIG	21	-1309.5454	2661.0907	0.0046	0.9987
GTR	19	-1312.8337	2663.6674	0.0013	1.0000
F81I	15	-1321.9565	2673.9130	0.0000	1.0000
F81G	15	-1321.9934	2673.9868	0.0000	1.0000
F81IG	16	-1321.9933	2675.9867	0.0000	1.0000
F81	14	-1325.1804	2678.3608	0.0000	1.0000
K2PI	13	-1358.6562	2743.3123	0.0000	1.0000
K2PG	13	-1358.7131	2743.4262	0.0000	1.0000
K2PIG	14	-1358.7401	2745.4803	0.0000	1.0000
K2P	12	-1361.8202	2747.6405	0.0000	1.0000
SYMI	17	-1357.9678	2749.9355	0.0000	1.0000
SYMG	17	-1358.0240	2750.0481	0.0000	1.0000
SYMIG	18	-1358.0256	2752.0512	0.0000	1.0000
SYM	16	-1361.0798	2754.1595	0.0000	1.0000
JC69I	12	-1367.8535	2759.7069	0.0000	1.0000
JC69G	12	-1367.9010	2759.8020	0.0000	1.0000
JC69IG	13	-1367.9245	2761.8490	0.0000	1.0000
JC69	11	-1371.0388	2764.0776	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYI	16	-1310.0842	2724.4511	0.3365	0.3365
HKY	15	-1313.3740	2724.5131	0.3262	0.6627
HKYG	16	-1310.1193	2724.5214	0.3249	0.9876
HKYIG	17	-1310.1378	2731.0760	0.0123	0.9998
F81	14	-1325.1804	2741.6082	0.0001	0.9999
F81I	15	-1321.9565	2741.6780	0.0001	0.9999
F81G	15	-1321.9934	2741.7518	0.0001	1.0000
F81IG	16	-1321.9933	2748.2694	0.0000	1.0000
GTRI	20	-1309.5070	2749.3673	0.0000	1.0000
GTRG	20	-1309.5318	2749.4169	0.0000	1.0000
GTR	19	-1312.8337	2749.5032	0.0000	1.0000
GTRIG	21	-1309.5454	2755.9618	0.0000	1.0000
K2P	12	-1361.8202	2801.8526	0.0000	1.0000
K2PI	13	-1358.6562	2802.0421	0.0000	1.0000
K2PG	13	-1358.7131	2802.1559	0.0000	1.0000
K2PIG	14	-1358.7401	2808.7277	0.0000	1.0000
JC69	11	-1371.0388	2813.7720	0.0000	1.0000
JC69I	12	-1367.8535	2813.9190	0.0000	1.0000
JC69G	12	-1367.9010	2814.0140	0.0000	1.0000
JC69IG	13	-1367.9245	2820.5787	0.0000	1.0000
SYM	16	-1361.0798	2826.4423	0.0000	1.0000
SYMI	17	-1357.9678	2826.7360	0.0000	1.0000
SYMG	17	-1358.0240	2826.8485	0.0000	1.0000
SYMIG	18	-1358.0256	2833.3693	0.0000	1.0000



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

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

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


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

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


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

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

