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

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

Minimum AIC  model: HKYG
Minimum AICc model: HKYG
Minimum BIC  model: HKYG


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	-1334.8339	2702.9724	0.5335	0.5335
HKYIG	17	-1334.7217	2704.9146	0.2020	0.7355
HKYI	16	-1336.4504	2706.2053	0.1060	0.8415
GTRG	20	-1332.1159	2706.2657	0.1028	0.9443
GTRIG	21	-1331.9244	2708.0916	0.0413	0.9855
GTRI	20	-1334.2135	2710.4610	0.0126	0.9982
K2PG	13	-1344.2349	2715.3364	0.0011	0.9993
K2PIG	14	-1344.1451	2717.2926	0.0004	0.9997
K2PI	13	-1345.6311	2718.1288	0.0003	0.9999
SYMG	17	-1343.4600	2722.3912	0.0000	1.0000
SYMIG	18	-1343.3503	2724.3487	0.0000	1.0000
SYMI	17	-1345.0122	2725.4956	0.0000	1.0000
HKY	15	-1351.8153	2734.7789	0.0000	1.0000
GTR	19	-1349.6485	2739.1327	0.0000	1.0000
K2P	12	-1361.2939	2747.3289	0.0000	1.0000
SYM	16	-1360.4487	2754.2019	0.0000	1.0000
F81G	15	-1362.1921	2755.5324	0.0000	1.0000
F81IG	16	-1362.1463	2757.5972	0.0000	1.0000
F81I	15	-1363.4973	2758.1429	0.0000	1.0000
JC69G	12	-1370.5378	2765.8168	0.0000	1.0000
JC69IG	13	-1370.5018	2767.8702	0.0000	1.0000
JC69I	12	-1371.7760	2768.2930	0.0000	1.0000
F81	14	-1377.2956	2783.5935	0.0000	1.0000
JC69	11	-1385.5628	2793.7512	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
HKYG	16	-1334.8339	2701.6678	0.4882	0.4882
HKYIG	17	-1334.7217	2703.4435	0.2009	0.6892
GTRG	20	-1332.1159	2704.2318	0.1355	0.8247
HKYI	16	-1336.4504	2704.9007	0.0970	0.9216
GTRIG	21	-1331.9244	2705.8489	0.0604	0.9820
GTRI	20	-1334.2135	2708.4271	0.0166	0.9986
K2PG	13	-1344.2349	2714.4698	0.0008	0.9994
K2PIG	14	-1344.1451	2716.2902	0.0003	0.9997
K2PI	13	-1345.6311	2717.2622	0.0002	0.9999
SYMG	17	-1343.4600	2720.9201	0.0000	1.0000
SYMIG	18	-1343.3503	2722.7005	0.0000	1.0000
SYMI	17	-1345.0122	2724.0245	0.0000	1.0000
HKY	15	-1351.8153	2733.6306	0.0000	1.0000
GTR	19	-1349.6485	2737.2970	0.0000	1.0000
K2P	12	-1361.2939	2746.5878	0.0000	1.0000
SYM	16	-1360.4487	2752.8973	0.0000	1.0000
F81G	15	-1362.1921	2754.3841	0.0000	1.0000
F81IG	16	-1362.1463	2756.2927	0.0000	1.0000
F81I	15	-1363.4973	2756.9946	0.0000	1.0000
JC69G	12	-1370.5378	2765.0757	0.0000	1.0000
JC69IG	13	-1370.5018	2767.0036	0.0000	1.0000
JC69I	12	-1371.7760	2767.5519	0.0000	1.0000
F81	14	-1377.2956	2782.5911	0.0000	1.0000
JC69	11	-1385.5628	2793.1256	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	16	-1334.8339	2766.8365	0.4497	0.4497
K2PG	13	-1344.2349	2767.4194	0.3360	0.7857
HKYI	16	-1336.4504	2770.0694	0.0893	0.8750
K2PI	13	-1345.6311	2770.2118	0.0832	0.9582
HKYIG	17	-1334.7217	2772.6852	0.0241	0.9823
K2PIG	14	-1344.1451	2773.3129	0.0176	1.0000
GTRG	20	-1332.1159	2785.6927	0.0000	1.0000
GTRI	20	-1334.2135	2789.8880	0.0000	1.0000
SYMG	17	-1343.4600	2790.1618	0.0000	1.0000
GTRIG	21	-1331.9244	2791.3828	0.0000	1.0000
SYMI	17	-1345.0122	2793.2662	0.0000	1.0000
HKY	15	-1351.8153	2794.7263	0.0000	1.0000
K2P	12	-1361.2939	2795.4644	0.0000	1.0000
SYMIG	18	-1343.3503	2796.0153	0.0000	1.0000
JC69G	12	-1370.5378	2813.9522	0.0000	1.0000
GTR	19	-1349.6485	2814.6848	0.0000	1.0000
F81G	15	-1362.1921	2815.4798	0.0000	1.0000
JC69I	12	-1371.7760	2816.4285	0.0000	1.0000
SYM	16	-1360.4487	2818.0660	0.0000	1.0000
F81I	15	-1363.4973	2818.0903	0.0000	1.0000
JC69IG	13	-1370.5018	2819.9531	0.0000	1.0000
F81IG	16	-1362.1463	2821.4614	0.0000	1.0000
JC69	11	-1385.5628	2837.9291	0.0000	1.0000
F81	14	-1377.2956	2839.6138	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 591-tomodel-phylip/chr2_20206.phylip = 1 - 434;
 Partition Dummy = 1:591-tomodel-phylip/chr2_20206.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 591-tomodel-phylip/chr2_20206.phylip = 1 - 434;
 Partition Dummy = 1:591-tomodel-phylip/chr2_20206.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 (HKYG)]

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


-------------------------------------------------------------
End of Output

