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

Input data from file "683-tomodel-phylip/chr3_378.phylip" (7 taxa, 671 characters)

Minimum AIC  model: GTRG
Minimum AICc model: GTRG
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
GTRG	20	-1534.7559	3110.8041	0.3916	0.3916
GTRI	20	-1535.1487	3111.5896	0.2644	0.6560
GTRIG	21	-1534.4208	3112.2653	0.1886	0.8446
HKYG	16	-1540.7039	3114.2395	0.0703	0.9149
HKYI	16	-1541.0079	3114.8477	0.0519	0.9667
HKYIG	17	-1540.3990	3115.7352	0.0333	1.0000
GTR	19	-1547.4644	3134.0962	0.0000	1.0000
SYMG	17	-1550.7534	3136.4439	0.0000	1.0000
SYMI	17	-1551.1131	3137.1635	0.0000	1.0000
SYMIG	18	-1550.4316	3137.9122	0.0000	1.0000
HKY	15	-1553.6776	3138.0880	0.0000	1.0000
K2PG	13	-1556.4438	3139.4417	0.0000	1.0000
K2PI	13	-1556.7122	3139.9784	0.0000	1.0000
K2PIG	14	-1556.1188	3140.8779	0.0000	1.0000
F81G	15	-1564.4078	3159.5484	0.0000	1.0000
SYM	16	-1563.9116	3160.6549	0.0000	1.0000
F81IG	16	-1564.0237	3160.8791	0.0000	1.0000
F81I	15	-1565.9145	3162.5619	0.0000	1.0000
K2P	12	-1569.9328	3164.3398	0.0000	1.0000
F81	14	-1579.1252	3186.8906	0.0000	1.0000
JC69G	12	-1581.5657	3187.6055	0.0000	1.0000
JC69IG	13	-1581.1658	3188.8856	0.0000	1.0000
JC69I	12	-1582.9967	3190.4676	0.0000	1.0000
JC69	11	-1595.8455	3214.0917	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRG	20	-1534.7559	3109.5118	0.3987	0.3987
GTRI	20	-1535.1487	3110.2973	0.2692	0.6678
GTRIG	21	-1534.4208	3110.8416	0.2050	0.8729
HKYG	16	-1540.7039	3113.4077	0.0568	0.9297
HKYI	16	-1541.0079	3114.0159	0.0419	0.9716
HKYIG	17	-1540.3990	3114.7980	0.0284	1.0000
GTR	19	-1547.4644	3132.9288	0.0000	1.0000
SYMG	17	-1550.7534	3135.5067	0.0000	1.0000
SYMI	17	-1551.1131	3136.2262	0.0000	1.0000
SYMIG	18	-1550.4316	3136.8631	0.0000	1.0000
HKY	15	-1553.6776	3137.3552	0.0000	1.0000
K2PG	13	-1556.4438	3138.8876	0.0000	1.0000
K2PI	13	-1556.7122	3139.4244	0.0000	1.0000
K2PIG	14	-1556.1188	3140.2377	0.0000	1.0000
F81G	15	-1564.4078	3158.8156	0.0000	1.0000
SYM	16	-1563.9116	3159.8231	0.0000	1.0000
F81IG	16	-1564.0237	3160.0473	0.0000	1.0000
F81I	15	-1565.9145	3161.8291	0.0000	1.0000
K2P	12	-1569.9328	3163.8656	0.0000	1.0000
F81	14	-1579.1252	3186.2503	0.0000	1.0000
JC69G	12	-1581.5657	3187.1314	0.0000	1.0000
JC69IG	13	-1581.1658	3188.3316	0.0000	1.0000
JC69I	12	-1582.9967	3189.9935	0.0000	1.0000
JC69	11	-1595.8455	3213.6911	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	16	-1540.7039	3185.5480	0.5567	0.5567
HKYI	16	-1541.0079	3186.1562	0.4107	0.9674
HKYIG	17	-1540.3990	3191.4471	0.0292	0.9966
K2PG	13	-1556.4438	3197.5016	0.0014	0.9980
K2PI	13	-1556.7122	3198.0384	0.0011	0.9991
GTRG	20	-1534.7559	3199.6872	0.0005	0.9995
GTRI	20	-1535.1487	3200.4727	0.0003	0.9999
K2PIG	14	-1556.1188	3203.3604	0.0001	0.9999
HKY	15	-1553.6776	3204.9868	0.0000	1.0000
GTRIG	21	-1534.4208	3205.5257	0.0000	1.0000
SYMG	17	-1550.7534	3212.1558	0.0000	1.0000
SYMI	17	-1551.1131	3212.8753	0.0000	1.0000
K2P	12	-1569.9328	3217.9709	0.0000	1.0000
SYMIG	18	-1550.4316	3218.0210	0.0000	1.0000
GTR	19	-1547.4644	3218.5954	0.0000	1.0000
F81G	15	-1564.4078	3226.4471	0.0000	1.0000
F81I	15	-1565.9145	3229.4606	0.0000	1.0000
SYM	16	-1563.9116	3231.9634	0.0000	1.0000
F81IG	16	-1564.0237	3232.1876	0.0000	1.0000
JC69G	12	-1581.5657	3241.2366	0.0000	1.0000
JC69I	12	-1582.9967	3244.0987	0.0000	1.0000
JC69IG	13	-1581.1658	3246.9456	0.0000	1.0000
F81	14	-1579.1252	3249.3731	0.0000	1.0000
JC69	11	-1595.8455	3263.2875	0.0000	1.0000



To specify default values for best AIC, AICc, or BIC models in MrBayes:

[Mrbayes block for the best AIC model (GTRG)]

BEGIN MRBAYES;
 Charset 683-tomodel-phylip/chr3_378.phylip = 1 - 671;
 Partition Dummy = 1:683-tomodel-phylip/chr3_378.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=6 rates=gamma;
 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) shapepr=Uniform(0.1,50.0);
END;


[Mrbayes block for the best AICc model (GTRG)]

BEGIN MRBAYES;
 Charset 683-tomodel-phylip/chr3_378.phylip = 1 - 671;
 Partition Dummy = 1:683-tomodel-phylip/chr3_378.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=6 rates=gamma;
 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) shapepr=Uniform(0.1,50.0);
END;


[Mrbayes block for the best BIC model (HKYG)]

BEGIN MRBAYES;
 Charset 683-tomodel-phylip/chr3_378.phylip = 1 - 671;
 Partition Dummy = 1:683-tomodel-phylip/chr3_378.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

