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

Input data from file "591-tomodel-phylip/chr2_6767.phylip" (7 taxa, 530 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	-1330.6913	2703.0329	0.4696	0.4696
GTRG	20	-1331.1864	2704.0231	0.2862	0.7558
GTRIG	21	-1330.7044	2705.2278	0.1567	0.9125
GTR	19	-1333.4852	2706.4607	0.0846	0.9971
HKYI	16	-1340.7994	2714.6592	0.0014	0.9985
HKYG	16	-1341.2573	2715.5750	0.0009	0.9994
HKYIG	17	-1340.8112	2716.8178	0.0005	0.9999
HKY	15	-1344.5774	2720.0888	0.0001	1.0000
SYMI	17	-1357.4782	2750.1516	0.0000	1.0000
SYMG	17	-1357.7485	2750.6924	0.0000	1.0000
SYMIG	18	-1357.4850	2752.3085	0.0000	1.0000
K2PI	13	-1363.0752	2752.8558	0.0000	1.0000
K2PG	13	-1363.3687	2753.4428	0.0000	1.0000
SYM	16	-1360.3475	2753.7555	0.0000	1.0000
K2PIG	14	-1363.0818	2754.9791	0.0000	1.0000
K2P	12	-1367.1741	2758.9516	0.0000	1.0000
F81I	15	-1364.4534	2759.8406	0.0000	1.0000
F81G	15	-1364.7385	2760.4109	0.0000	1.0000
F81IG	16	-1364.4608	2761.9821	0.0000	1.0000
F81	14	-1367.4031	2763.6217	0.0000	1.0000
JC69I	12	-1384.0172	2792.6379	0.0000	1.0000
JC69G	12	-1384.1891	2792.9816	0.0000	1.0000
JC69IG	13	-1384.0210	2794.7473	0.0000	1.0000
JC69	11	-1387.1665	2796.8426	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRI	20	-1330.6913	2701.3826	0.4665	0.4665
GTRG	20	-1331.1864	2702.3728	0.2844	0.7509
GTRIG	21	-1330.7044	2703.4089	0.1694	0.9203
GTR	19	-1333.4852	2704.9705	0.0776	0.9979
HKYI	16	-1340.7994	2713.5987	0.0010	0.9989
HKYG	16	-1341.2573	2714.5146	0.0007	0.9996
HKYIG	17	-1340.8112	2715.6225	0.0004	0.9999
HKY	15	-1344.5774	2719.1549	0.0001	1.0000
SYMI	17	-1357.4782	2748.9563	0.0000	1.0000
SYMG	17	-1357.7485	2749.4971	0.0000	1.0000
SYMIG	18	-1357.4850	2750.9700	0.0000	1.0000
K2PI	13	-1363.0752	2752.1504	0.0000	1.0000
SYM	16	-1360.3475	2752.6951	0.0000	1.0000
K2PG	13	-1363.3687	2752.7373	0.0000	1.0000
K2PIG	14	-1363.0818	2754.1635	0.0000	1.0000
K2P	12	-1367.1741	2758.3482	0.0000	1.0000
F81I	15	-1364.4534	2758.9068	0.0000	1.0000
F81G	15	-1364.7385	2759.4771	0.0000	1.0000
F81IG	16	-1364.4608	2760.9217	0.0000	1.0000
F81	14	-1367.4031	2762.8062	0.0000	1.0000
JC69I	12	-1384.0172	2792.0344	0.0000	1.0000
JC69G	12	-1384.1891	2792.3781	0.0000	1.0000
JC69IG	13	-1384.0210	2794.0419	0.0000	1.0000
JC69	11	-1387.1665	2796.3330	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYI	16	-1340.7994	2781.9648	0.4050	0.4050
HKYG	16	-1341.2573	2782.8806	0.2562	0.6611
HKY	15	-1344.5774	2783.2481	0.2132	0.8743
GTR	19	-1333.4852	2786.1551	0.0498	0.9242
GTRI	20	-1330.6913	2786.8402	0.0354	0.9595
GTRG	20	-1331.1864	2787.8304	0.0216	0.9811
HKYIG	17	-1340.8112	2788.2614	0.0174	0.9985
GTRIG	21	-1330.7044	2793.1393	0.0015	1.0000
K2PI	13	-1363.0752	2807.6978	0.0000	1.0000
K2PG	13	-1363.3687	2808.2847	0.0000	1.0000
K2P	12	-1367.1741	2809.6227	0.0000	1.0000
K2PIG	14	-1363.0818	2813.9838	0.0000	1.0000
SYM	16	-1360.3475	2821.0611	0.0000	1.0000
SYMI	17	-1357.4782	2821.5952	0.0000	1.0000
SYMG	17	-1357.7485	2822.1360	0.0000	1.0000
F81	14	-1367.4031	2822.6265	0.0000	1.0000
F81I	15	-1364.4534	2822.9999	0.0000	1.0000
F81G	15	-1364.7385	2823.5702	0.0000	1.0000
SYMIG	18	-1357.4850	2827.8818	0.0000	1.0000
F81IG	16	-1364.4608	2829.2877	0.0000	1.0000
JC69I	12	-1384.0172	2843.3089	0.0000	1.0000
JC69	11	-1387.1665	2843.3346	0.0000	1.0000
JC69G	12	-1384.1891	2843.6526	0.0000	1.0000
JC69IG	13	-1384.0210	2849.5893	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/chr2_6767.phylip = 1 - 530;
 Partition Dummy = 1:591-tomodel-phylip/chr2_6767.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/chr2_6767.phylip = 1 - 530;
 Partition Dummy = 1:591-tomodel-phylip/chr2_6767.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/chr2_6767.phylip = 1 - 530;
 Partition Dummy = 1:591-tomodel-phylip/chr2_6767.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

