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

Input data from file "683-tomodel-phylip/chr2_20477.phylip" (7 taxa, 701 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	-1605.0251	3251.2855	0.4477	0.4477
GTRG	20	-1605.2806	3251.7965	0.3467	0.7944
GTRIG	21	-1605.3122	3253.9851	0.1161	0.9105
HKYI	16	-1611.6691	3256.1336	0.0396	0.9502
HKYG	16	-1611.7766	3256.3486	0.0356	0.9858
HKYIG	17	-1611.8088	3258.5136	0.0121	0.9978
GTR	19	-1611.4680	3262.0519	0.0021	0.9999
HKY	15	-1618.5785	3267.8577	0.0001	1.0000
K2PG	13	-1649.7077	3325.9452	0.0000	1.0000
SYMG	17	-1645.5407	3325.9774	0.0000	1.0000
K2PI	13	-1649.7351	3326.0001	0.0000	1.0000
SYMI	17	-1645.5551	3326.0063	0.0000	1.0000
K2PIG	14	-1649.7181	3328.0484	0.0000	1.0000
SYMIG	18	-1645.5491	3328.1011	0.0000	1.0000
SYM	16	-1652.2661	3337.3276	0.0000	1.0000
K2P	12	-1657.1081	3338.6697	0.0000	1.0000
F81I	15	-1659.0377	3348.7761	0.0000	1.0000
F81G	15	-1659.2146	3349.1299	0.0000	1.0000
F81IG	16	-1659.2354	3351.2660	0.0000	1.0000
F81	14	-1664.5624	3357.7370	0.0000	1.0000
JC69I	12	-1692.3605	3409.1746	0.0000	1.0000
JC69G	12	-1692.4513	3409.3561	0.0000	1.0000
JC69IG	13	-1692.4682	3411.4663	0.0000	1.0000
JC69	11	-1698.1991	3418.7815	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRI	20	-1605.0251	3250.0502	0.4519	0.4519
GTRG	20	-1605.2806	3250.5612	0.3500	0.8019
GTRIG	21	-1605.3122	3252.6243	0.1248	0.9267
HKYI	16	-1611.6691	3255.3382	0.0321	0.9588
HKYG	16	-1611.7766	3255.5533	0.0288	0.9877
HKYIG	17	-1611.8088	3257.6175	0.0103	0.9980
GTR	19	-1611.4680	3260.9359	0.0020	0.9999
HKY	15	-1618.5785	3267.1569	0.0001	1.0000
SYMG	17	-1645.5407	3325.0814	0.0000	1.0000
SYMI	17	-1645.5551	3325.1103	0.0000	1.0000
K2PG	13	-1649.7077	3325.4153	0.0000	1.0000
K2PI	13	-1649.7351	3325.4702	0.0000	1.0000
SYMIG	18	-1645.5491	3327.0981	0.0000	1.0000
K2PIG	14	-1649.7181	3327.4362	0.0000	1.0000
SYM	16	-1652.2661	3336.5322	0.0000	1.0000
K2P	12	-1657.1081	3338.2162	0.0000	1.0000
F81I	15	-1659.0377	3348.0754	0.0000	1.0000
F81G	15	-1659.2146	3348.4291	0.0000	1.0000
F81IG	16	-1659.2354	3350.4707	0.0000	1.0000
F81	14	-1664.5624	3357.1248	0.0000	1.0000
JC69I	12	-1692.3605	3408.7211	0.0000	1.0000
JC69G	12	-1692.4513	3408.9026	0.0000	1.0000
JC69IG	13	-1692.4682	3410.9365	0.0000	1.0000
JC69	11	-1698.1991	3418.3983	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYI	16	-1611.6691	3328.1784	0.5101	0.5101
HKYG	16	-1611.7766	3328.3934	0.4581	0.9683
HKYIG	17	-1611.8088	3335.0102	0.0168	0.9850
HKY	15	-1618.5785	3335.4446	0.0135	0.9985
GTRI	20	-1605.0251	3341.1004	0.0008	0.9993
GTRG	20	-1605.2806	3341.6114	0.0006	0.9999
GTR	19	-1611.4680	3347.4336	0.0000	1.0000
GTRIG	21	-1605.3122	3348.2270	0.0000	1.0000
K2PG	13	-1649.7077	3384.5979	0.0000	1.0000
K2PI	13	-1649.7351	3384.6528	0.0000	1.0000
K2PIG	14	-1649.7181	3391.1713	0.0000	1.0000
K2P	12	-1657.1081	3392.8463	0.0000	1.0000
SYMG	17	-1645.5407	3402.4740	0.0000	1.0000
SYMI	17	-1645.5551	3402.5029	0.0000	1.0000
SYMIG	18	-1645.5491	3409.0433	0.0000	1.0000
SYM	16	-1652.2661	3409.3724	0.0000	1.0000
F81I	15	-1659.0377	3416.3630	0.0000	1.0000
F81G	15	-1659.2146	3416.7168	0.0000	1.0000
F81	14	-1664.5624	3420.8599	0.0000	1.0000
F81IG	16	-1659.2354	3423.3108	0.0000	1.0000
JC69I	12	-1692.3605	3463.3512	0.0000	1.0000
JC69G	12	-1692.4513	3463.5327	0.0000	1.0000
JC69	11	-1698.1991	3468.4759	0.0000	1.0000
JC69IG	13	-1692.4682	3470.1191	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 683-tomodel-phylip/chr2_20477.phylip = 1 - 701;
 Partition Dummy = 1:683-tomodel-phylip/chr2_20477.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 683-tomodel-phylip/chr2_20477.phylip = 1 - 701;
 Partition Dummy = 1:683-tomodel-phylip/chr2_20477.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 683-tomodel-phylip/chr2_20477.phylip = 1 - 701;
 Partition Dummy = 1:683-tomodel-phylip/chr2_20477.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

