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

Input data from file "683-tomodel-phylip/chr1_32429.phylip" (7 taxa, 815 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	-1648.2546	3329.1909	0.3210	0.3210
HKYI	16	-1648.2784	3329.2385	0.3135	0.6345
HKYIG	17	-1648.2529	3331.2736	0.1133	0.7478
GTRG	20	-1645.2060	3331.4700	0.1027	0.8505
GTRI	20	-1645.2405	3331.5390	0.0992	0.9498
GTRIG	21	-1645.2002	3333.5656	0.0360	0.9858
HKY	15	-1652.7258	3336.0523	0.0104	0.9962
GTR	19	-1649.5445	3338.0450	0.0038	1.0000
F81G	15	-1663.9759	3358.5526	0.0000	1.0000
F81I	15	-1664.0257	3358.6522	0.0000	1.0000
F81IG	16	-1663.9689	3360.6195	0.0000	1.0000
F81	14	-1667.9195	3364.3639	0.0000	1.0000
K2PG	13	-1677.8428	3382.1400	0.0000	1.0000
K2PI	13	-1677.9019	3382.2583	0.0000	1.0000
K2PIG	14	-1677.8298	3384.1846	0.0000	1.0000
SYMG	17	-1676.5485	3387.8648	0.0000	1.0000
SYMI	17	-1676.6433	3388.0544	0.0000	1.0000
K2P	12	-1682.5654	3389.5197	0.0000	1.0000
SYMIG	18	-1676.5250	3389.9094	0.0000	1.0000
SYM	16	-1680.5604	3393.8026	0.0000	1.0000
JC69G	12	-1692.0346	3408.4581	0.0000	1.0000
JC69I	12	-1692.1047	3408.5985	0.0000	1.0000
JC69IG	13	-1692.0215	3410.4974	0.0000	1.0000
JC69	11	-1695.9666	3414.2619	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
HKYG	16	-1648.2546	3328.5092	0.3038	0.3038
HKYI	16	-1648.2784	3328.5568	0.2966	0.6004
GTRG	20	-1645.2060	3330.4120	0.1173	0.7177
GTRI	20	-1645.2405	3330.4810	0.1133	0.8310
HKYIG	17	-1648.2529	3330.5058	0.1119	0.9430
GTRIG	21	-1645.2002	3332.4004	0.0434	0.9864
HKY	15	-1652.7258	3335.4516	0.0094	0.9958
GTR	19	-1649.5445	3337.0890	0.0042	1.0000
F81G	15	-1663.9759	3357.9518	0.0000	1.0000
F81I	15	-1664.0257	3358.0514	0.0000	1.0000
F81IG	16	-1663.9689	3359.9378	0.0000	1.0000
F81	14	-1667.9195	3363.8389	0.0000	1.0000
K2PG	13	-1677.8428	3381.6855	0.0000	1.0000
K2PI	13	-1677.9019	3381.8039	0.0000	1.0000
K2PIG	14	-1677.8298	3383.6596	0.0000	1.0000
SYMG	17	-1676.5485	3387.0969	0.0000	1.0000
SYMI	17	-1676.6433	3387.2865	0.0000	1.0000
SYMIG	18	-1676.5250	3389.0501	0.0000	1.0000
K2P	12	-1682.5654	3389.1307	0.0000	1.0000
SYM	16	-1680.5604	3393.1209	0.0000	1.0000
JC69G	12	-1692.0346	3408.0691	0.0000	1.0000
JC69I	12	-1692.1047	3408.2095	0.0000	1.0000
JC69IG	13	-1692.0215	3410.0430	0.0000	1.0000
JC69	11	-1695.9666	3413.9331	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	16	-1648.2546	3403.7602	0.4277	0.4277
HKYI	16	-1648.2784	3403.8078	0.4176	0.8453
HKY	15	-1652.7258	3405.9994	0.1396	0.9850
HKYIG	17	-1648.2529	3410.4600	0.0150	1.0000
GTRG	20	-1645.2060	3424.4758	0.0000	1.0000
GTRI	20	-1645.2405	3424.5448	0.0000	1.0000
GTR	19	-1649.5445	3426.4496	0.0000	1.0000
F81G	15	-1663.9759	3428.4996	0.0000	1.0000
F81I	15	-1664.0257	3428.5992	0.0000	1.0000
F81	14	-1667.9195	3429.6836	0.0000	1.0000
GTRIG	21	-1645.2002	3431.1674	0.0000	1.0000
F81IG	16	-1663.9689	3435.1888	0.0000	1.0000
K2PG	13	-1677.8428	3442.8270	0.0000	1.0000
K2PI	13	-1677.9019	3442.9453	0.0000	1.0000
K2P	12	-1682.5654	3445.5690	0.0000	1.0000
K2PIG	14	-1677.8298	3449.5043	0.0000	1.0000
JC69G	12	-1692.0346	3464.5074	0.0000	1.0000
JC69I	12	-1692.1047	3464.6477	0.0000	1.0000
JC69	11	-1695.9666	3465.6682	0.0000	1.0000
SYMG	17	-1676.5485	3467.0511	0.0000	1.0000
SYMI	17	-1676.6433	3467.2407	0.0000	1.0000
SYM	16	-1680.5604	3468.3719	0.0000	1.0000
JC69IG	13	-1692.0215	3471.1844	0.0000	1.0000
SYMIG	18	-1676.5250	3473.7074	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/chr1_32429.phylip = 1 - 815;
 Partition Dummy = 1:683-tomodel-phylip/chr1_32429.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/chr1_32429.phylip = 1 - 815;
 Partition Dummy = 1:683-tomodel-phylip/chr1_32429.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 683-tomodel-phylip/chr1_32429.phylip = 1 - 815;
 Partition Dummy = 1:683-tomodel-phylip/chr1_32429.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

