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

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

Minimum AIC  model: GTR
Minimum AICc model: GTR
Minimum BIC  model: F81


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
GTR	19	-408.4699	859.1620	0.2414	0.2414
HKY	15	-413.2928	859.1944	0.2376	0.4790
GTRG	20	-408.2060	861.1047	0.0914	0.5704
HKYG	16	-413.0671	861.1070	0.0913	0.6617
HKYI	16	-413.0729	861.1184	0.0908	0.7525
GTRI	20	-408.2435	861.1798	0.0880	0.8405
F81	14	-415.9090	862.0884	0.0559	0.8964
HKYIG	17	-413.0678	863.4983	0.0276	0.9240
GTRIG	21	-408.2059	863.6028	0.0262	0.9502
F81I	15	-415.6854	863.9795	0.0217	0.9719
F81G	15	-415.7006	864.0098	0.0214	0.9933
F81IG	16	-415.6973	866.3673	0.0066	0.9999
SYM	16	-421.6626	878.2979	0.0000	0.9999
SYMI	17	-420.6345	878.6317	0.0000	0.9999
K2P	12	-426.4928	878.6540	0.0000	1.0000
SYMG	17	-420.6726	878.7079	0.0000	1.0000
K2PI	13	-426.0919	880.1407	0.0000	1.0000
K2PG	13	-426.0949	880.1468	0.0000	1.0000
JC69	11	-428.8580	881.1203	0.0000	1.0000
SYMIG	18	-420.6878	881.1546	0.0000	1.0000
K2PIG	14	-426.0949	882.4601	0.0000	1.0000
JC69I	12	-428.5180	882.7045	0.0000	1.0000
JC69G	12	-428.5344	882.7372	0.0000	1.0000
JC69IG	13	-428.5331	885.0232	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTR	19	-408.4699	854.9398	0.3110	0.3110
GTRG	20	-408.2060	856.4120	0.1490	0.4599
GTRI	20	-408.2435	856.4870	0.1435	0.6034
HKY	15	-413.2928	856.5857	0.1366	0.7400
HKYG	16	-413.0671	858.1343	0.0630	0.8030
HKYI	16	-413.0729	858.1457	0.0626	0.8656
GTRIG	21	-408.2059	858.4118	0.0548	0.9204
F81	14	-415.9090	859.8181	0.0271	0.9475
HKYIG	17	-413.0678	860.1357	0.0231	0.9706
F81I	15	-415.6854	861.3708	0.0125	0.9831
F81G	15	-415.7006	861.4011	0.0123	0.9954
F81IG	16	-415.6973	863.3946	0.0045	0.9999
SYMI	17	-420.6345	875.2691	0.0000	1.0000
SYM	16	-421.6626	875.3252	0.0000	1.0000
SYMG	17	-420.6726	875.3453	0.0000	1.0000
K2P	12	-426.4928	876.9855	0.0000	1.0000
SYMIG	18	-420.6878	877.3756	0.0000	1.0000
K2PI	13	-426.0919	878.1837	0.0000	1.0000
K2PG	13	-426.0949	878.1898	0.0000	1.0000
JC69	11	-428.8580	879.7160	0.0000	1.0000
K2PIG	14	-426.0949	880.1898	0.0000	1.0000
JC69I	12	-428.5180	881.0361	0.0000	1.0000
JC69G	12	-428.5344	881.0687	0.0000	1.0000
JC69IG	13	-428.5331	883.0662	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
F81	14	-415.9090	905.9945	0.4266	0.4266
HKY	15	-413.2928	906.0604	0.4128	0.8393
F81I	15	-415.6854	910.8456	0.0377	0.8771
F81G	15	-415.7006	910.8759	0.0372	0.9142
HKYG	16	-413.0671	910.9074	0.0366	0.9508
HKYI	16	-413.0729	910.9188	0.0364	0.9872
JC69	11	-428.8580	915.9975	0.0029	0.9900
F81IG	16	-415.6973	916.1677	0.0026	0.9927
HKYIG	17	-413.0678	916.2071	0.0026	0.9952
K2P	12	-426.4928	916.5653	0.0022	0.9974
GTR	19	-408.4699	917.6078	0.0013	0.9987
JC69I	12	-428.5180	920.6159	0.0003	0.9990
JC69G	12	-428.5344	920.6485	0.0003	0.9993
K2PI	13	-426.0919	921.0618	0.0002	0.9995
K2PG	13	-426.0949	921.0679	0.0002	0.9997
GTRG	20	-408.2060	922.3783	0.0001	0.9998
GTRI	20	-408.2435	922.4534	0.0001	0.9999
JC69IG	13	-428.5331	925.9444	0.0000	1.0000
K2PIG	14	-426.0949	926.3663	0.0000	1.0000
GTRIG	21	-408.2059	927.6764	0.0000	1.0000
SYM	16	-421.6626	928.0983	0.0000	1.0000
SYMI	17	-420.6345	931.3405	0.0000	1.0000
SYMG	17	-420.6726	931.4167	0.0000	1.0000
SYMIG	18	-420.6878	936.7453	0.0000	1.0000



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

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

BEGIN MRBAYES;
 Charset 591-tomodel-phylip/chr3_3158.phylip = 1 - 200;
 Partition Dummy = 1:591-tomodel-phylip/chr3_3158.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=6 rates=equal;
 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);
END;


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

BEGIN MRBAYES;
 Charset 591-tomodel-phylip/chr3_3158.phylip = 1 - 200;
 Partition Dummy = 1:591-tomodel-phylip/chr3_3158.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=6 rates=equal;
 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);
END;


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

BEGIN MRBAYES;
 Charset 591-tomodel-phylip/chr3_3158.phylip = 1 - 200;
 Partition Dummy = 1:591-tomodel-phylip/chr3_3158.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=1 rates=equal;
 Prset applyto=(1) statefreqpr=Dirichlet(1.0,1.0,1.0,1.0);
END;


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

