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

Input data from file "591-tomodel-phylip/chr5_1737.phylip" (7 taxa, 530 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.6511	2654.3626	0.4590	0.4590
HKYG	16	-1311.1492	2655.3589	0.2789	0.7379
HKYIG	17	-1310.6668	2656.5288	0.1554	0.8933
GTRI	20	-1308.6225	2658.8952	0.0476	0.9409
GTRG	20	-1309.1899	2660.0300	0.0270	0.9679
GTRIG	21	-1308.6406	2661.1002	0.0158	0.9837
HKY	15	-1315.1565	2661.2468	0.0147	0.9984
GTR	19	-1313.0787	2665.6477	0.0016	1.0000
SYMI	17	-1327.0088	2689.2129	0.0000	1.0000
SYMG	17	-1327.5372	2690.2696	0.0000	1.0000
K2PI	13	-1331.8700	2690.4454	0.0000	1.0000
F81I	15	-1330.1157	2691.1654	0.0000	1.0000
K2PG	13	-1332.3283	2691.3621	0.0000	1.0000
SYMIG	18	-1327.0257	2691.3900	0.0000	1.0000
F81G	15	-1330.5718	2692.0774	0.0000	1.0000
K2PIG	14	-1331.8844	2692.5844	0.0000	1.0000
F81IG	16	-1330.1301	2693.3206	0.0000	1.0000
SYM	16	-1331.5873	2696.2351	0.0000	1.0000
F81	14	-1333.8799	2696.5754	0.0000	1.0000
K2P	12	-1336.1805	2696.9645	0.0000	1.0000
JC69I	12	-1349.3116	2723.2266	0.0000	1.0000
JC69G	12	-1349.7402	2724.0839	0.0000	1.0000
JC69IG	13	-1349.3249	2725.3552	0.0000	1.0000
JC69	11	-1353.0666	2728.6428	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
HKYI	16	-1310.6511	2653.3022	0.4400	0.4400
HKYG	16	-1311.1492	2654.2985	0.2674	0.7073
HKYIG	17	-1310.6668	2655.3335	0.1593	0.8667
GTRI	20	-1308.6225	2657.2449	0.0613	0.9280
GTRG	20	-1309.1899	2658.3797	0.0347	0.9627
GTRIG	21	-1308.6406	2659.2813	0.0221	0.9849
HKY	15	-1315.1565	2660.3130	0.0132	0.9981
GTR	19	-1313.0787	2664.1575	0.0019	1.0000
SYMI	17	-1327.0088	2688.0176	0.0000	1.0000
SYMG	17	-1327.5372	2689.0743	0.0000	1.0000
K2PI	13	-1331.8700	2689.7400	0.0000	1.0000
SYMIG	18	-1327.0257	2690.0515	0.0000	1.0000
F81I	15	-1330.1157	2690.2315	0.0000	1.0000
K2PG	13	-1332.3283	2690.6567	0.0000	1.0000
F81G	15	-1330.5718	2691.1435	0.0000	1.0000
K2PIG	14	-1331.8844	2691.7689	0.0000	1.0000
F81IG	16	-1330.1301	2692.2602	0.0000	1.0000
SYM	16	-1331.5873	2695.1746	0.0000	1.0000
F81	14	-1333.8799	2695.7598	0.0000	1.0000
K2P	12	-1336.1805	2696.3610	0.0000	1.0000
JC69I	12	-1349.3116	2722.6232	0.0000	1.0000
JC69G	12	-1349.7402	2723.4804	0.0000	1.0000
JC69IG	13	-1349.3249	2724.6497	0.0000	1.0000
JC69	11	-1353.0666	2728.1332	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYI	16	-1310.6511	2721.6682	0.5250	0.5250
HKYG	16	-1311.1492	2722.6645	0.3190	0.8440
HKY	15	-1315.1565	2724.4061	0.1335	0.9775
HKYIG	17	-1310.6668	2727.9724	0.0224	1.0000
GTRI	20	-1308.6225	2742.7025	0.0000	1.0000
GTRG	20	-1309.1899	2743.8373	0.0000	1.0000
K2PI	13	-1331.8700	2745.2874	0.0000	1.0000
GTR	19	-1313.0787	2745.3422	0.0000	1.0000
K2PG	13	-1332.3283	2746.2041	0.0000	1.0000
K2P	12	-1336.1805	2747.6355	0.0000	1.0000
GTRIG	21	-1308.6406	2749.0117	0.0000	1.0000
K2PIG	14	-1331.8844	2751.5891	0.0000	1.0000
F81I	15	-1330.1157	2754.3247	0.0000	1.0000
F81G	15	-1330.5718	2755.2367	0.0000	1.0000
F81	14	-1333.8799	2755.5801	0.0000	1.0000
F81IG	16	-1330.1301	2760.6263	0.0000	1.0000
SYMI	17	-1327.0088	2760.6565	0.0000	1.0000
SYMG	17	-1327.5372	2761.7132	0.0000	1.0000
SYM	16	-1331.5873	2763.5407	0.0000	1.0000
SYMIG	18	-1327.0257	2766.9632	0.0000	1.0000
JC69I	12	-1349.3116	2773.8977	0.0000	1.0000
JC69G	12	-1349.7402	2774.7550	0.0000	1.0000
JC69	11	-1353.0666	2775.1348	0.0000	1.0000
JC69IG	13	-1349.3249	2780.1971	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/chr5_1737.phylip = 1 - 530;
 Partition Dummy = 1:591-tomodel-phylip/chr5_1737.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/chr5_1737.phylip = 1 - 530;
 Partition Dummy = 1:591-tomodel-phylip/chr5_1737.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/chr5_1737.phylip = 1 - 530;
 Partition Dummy = 1:591-tomodel-phylip/chr5_1737.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

