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

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

Minimum AIC  model: GTRI
Minimum AICc model: GTRI
Minimum BIC  model: K2P


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	-1551.8965	3145.6311	0.3600	0.3600
GTRG	20	-1552.4615	3146.7611	0.2046	0.5647
GTRIG	21	-1551.9148	3147.8560	0.1184	0.6830
GTR	19	-1554.1137	3147.8867	0.1166	0.7996
HKYI	16	-1557.6703	3148.5207	0.0849	0.8845
HKYG	16	-1558.1144	3149.4089	0.0545	0.9389
HKYIG	17	-1557.6843	3150.6991	0.0286	0.9675
HKY	15	-1559.8736	3150.7861	0.0273	0.9949
SYMI	17	-1560.5998	3156.5300	0.0015	0.9964
SYMG	17	-1561.1856	3157.7016	0.0009	0.9973
K2PI	13	-1565.5682	3157.9209	0.0008	0.9980
SYMIG	18	-1560.6189	3158.7279	0.0005	0.9986
K2PG	13	-1565.9917	3158.7680	0.0005	0.9991
SYM	16	-1562.9342	3159.0484	0.0004	0.9995
K2PIG	14	-1565.5819	3160.0709	0.0003	0.9998
K2P	12	-1567.8215	3160.3140	0.0002	1.0000
F81I	15	-1581.7488	3194.5365	0.0000	1.0000
F81G	15	-1582.1126	3195.2642	0.0000	1.0000
F81	14	-1583.2987	3195.5046	0.0000	1.0000
F81IG	16	-1581.7621	3196.7043	0.0000	1.0000
JC69I	12	-1588.5713	3201.8136	0.0000	1.0000
JC69G	12	-1588.9623	3202.5956	0.0000	1.0000
JC69	11	-1590.2064	3202.9793	0.0000	1.0000
JC69IG	13	-1588.5854	3203.9553	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRI	20	-1551.8965	3143.7930	0.3807	0.3807
GTRG	20	-1552.4615	3144.9230	0.2164	0.5971
GTRIG	21	-1551.9148	3145.8296	0.1375	0.7346
GTR	19	-1554.1137	3146.2273	0.1127	0.8473
HKYI	16	-1557.6703	3147.3406	0.0646	0.9119
HKYG	16	-1558.1144	3148.2289	0.0414	0.9533
HKYIG	17	-1557.6843	3149.3686	0.0234	0.9767
HKY	15	-1559.8736	3149.7471	0.0194	0.9961
SYMI	17	-1560.5998	3155.1995	0.0013	0.9974
SYMG	17	-1561.1856	3156.3712	0.0007	0.9981
K2PI	13	-1565.5682	3157.1364	0.0005	0.9986
SYMIG	18	-1560.6189	3157.2377	0.0005	0.9990
SYM	16	-1562.9342	3157.8684	0.0003	0.9994
K2PG	13	-1565.9917	3157.9835	0.0003	0.9997
K2PIG	14	-1565.5819	3159.1638	0.0002	0.9999
K2P	12	-1567.8215	3159.6430	0.0001	1.0000
F81I	15	-1581.7488	3193.4976	0.0000	1.0000
F81G	15	-1582.1126	3194.2252	0.0000	1.0000
F81	14	-1583.2987	3194.5975	0.0000	1.0000
F81IG	16	-1581.7621	3195.5242	0.0000	1.0000
JC69I	12	-1588.5713	3201.1426	0.0000	1.0000
JC69G	12	-1588.9623	3201.9246	0.0000	1.0000
JC69	11	-1590.2064	3202.4127	0.0000	1.0000
JC69IG	13	-1588.5854	3203.1708	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
K2P	12	-1567.8215	3209.6784	0.4542	0.4542
K2PI	13	-1565.5682	3211.3414	0.1978	0.6520
K2PG	13	-1565.9917	3212.1884	0.1295	0.7814
HKY	15	-1559.8736	3212.2913	0.1230	0.9044
HKYI	16	-1557.6703	3214.0544	0.0509	0.9554
HKYG	16	-1558.1144	3214.9427	0.0327	0.9880
K2PIG	14	-1565.5819	3217.5383	0.0089	0.9970
HKYIG	17	-1557.6843	3220.2520	0.0023	0.9993
SYM	16	-1562.9342	3224.5822	0.0003	0.9995
GTR	19	-1554.1137	3225.4499	0.0002	0.9997
SYMI	17	-1560.5998	3226.0829	0.0001	0.9998
GTRI	20	-1551.8965	3227.1852	0.0001	0.9999
SYMG	17	-1561.1856	3227.2545	0.0001	1.0000
GTRG	20	-1552.4615	3228.3152	0.0000	1.0000
SYMIG	18	-1560.6189	3232.2907	0.0000	1.0000
GTRIG	21	-1551.9148	3233.3915	0.0000	1.0000
JC69	11	-1590.2064	3248.2785	0.0000	1.0000
JC69I	12	-1588.5713	3251.1779	0.0000	1.0000
JC69G	12	-1588.9623	3251.9599	0.0000	1.0000
F81	14	-1583.2987	3252.9720	0.0000	1.0000
F81I	15	-1581.7488	3256.0417	0.0000	1.0000
F81G	15	-1582.1126	3256.7694	0.0000	1.0000
JC69IG	13	-1588.5854	3257.3757	0.0000	1.0000
F81IG	16	-1581.7621	3262.2380	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 591-tomodel-phylip/chr5_10051.phylip = 1 - 478;
 Partition Dummy = 1:591-tomodel-phylip/chr5_10051.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 591-tomodel-phylip/chr5_10051.phylip = 1 - 478;
 Partition Dummy = 1:591-tomodel-phylip/chr5_10051.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 (K2P)]

BEGIN MRBAYES;
 Charset 591-tomodel-phylip/chr5_10051.phylip = 1 - 478;
 Partition Dummy = 1:591-tomodel-phylip/chr5_10051.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=2 rates=equal;
 Prset applyto=(1) statefreqpr=Fixed(Equal) tratiopr=Beta(1.0,1.0);
END;


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

