
Output from mraic.pl version 1.4.4 by Johan Nylander
-------------------------------------------------------------
00:44 02/17/2011

Input data from file "183-phylip/chr3_18256.phylip" (29 taxa, 118 characters)

Minimum AIC  model: GTR
Minimum AICc model: K2P
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
K2P	56	-386.4313	989.5183	0.9006	0.9006
JC69	55	-392.9083	995.1714	0.0533	0.9539
K2PI	57	-386.3750	996.9500	0.0219	0.9758
K2PG	57	-386.4277	997.0555	0.0208	0.9966
JC69I	56	-392.8581	1002.3720	0.0015	0.9981
JC69G	56	-392.9069	1002.4695	0.0014	0.9995
K2PIG	58	-386.3800	1004.7601	0.0004	0.9999
HKY	59	-384.3752	1008.8194	0.0001	1.0000
JC69IG	57	-392.8639	1009.9278	0.0000	1.0000
SYM	60	-382.6295	1013.6801	0.0000	1.0000
F81	58	-391.1494	1014.2988	0.0000	1.0000
HKYI	60	-384.3027	1017.0264	0.0000	1.0000
HKYG	60	-384.3696	1017.1602	0.0000	1.0000
SYMI	61	-382.5612	1022.1938	0.0000	1.0000
F81I	59	-391.0890	1022.2469	0.0000	1.0000
SYMG	61	-382.6285	1022.3285	0.0000	1.0000
F81G	59	-391.1472	1022.3633	0.0000	1.0000
HKYIG	61	-384.3085	1025.6884	0.0000	1.0000
F81IG	60	-391.0949	1030.6109	0.0000	1.0000
SYMIG	62	-382.5682	1031.1728	0.0000	1.0000
GTR	63	-379.3236	1033.9806	0.0000	1.0000
GTRI	64	-379.0553	1043.0917	0.0000	1.0000
GTRG	64	-379.2691	1043.5193	0.0000	1.0000
GTRIG	65	-379.0660	1053.1321	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTR	63	-379.3236	884.6473	0.1715	0.1715
K2P	56	-386.4313	884.8625	0.1540	0.3254
SYM	60	-382.6295	885.2590	0.1263	0.4517
GTRI	64	-379.0553	886.1106	0.0825	0.5342
GTRG	64	-379.2691	886.5381	0.0666	0.6008
K2PI	57	-386.3750	886.7500	0.0599	0.6608
HKY	59	-384.3752	886.7504	0.0599	0.7207
K2PG	57	-386.4277	886.8555	0.0568	0.7775
SYMI	61	-382.5612	887.1224	0.0497	0.8273
SYMG	61	-382.6285	887.2571	0.0465	0.8738
GTRIG	65	-379.0660	888.1321	0.0300	0.9038
HKYI	60	-384.3027	888.6054	0.0237	0.9275
HKYG	60	-384.3696	888.7392	0.0222	0.9497
K2PIG	58	-386.3800	888.7601	0.0219	0.9716
SYMIG	62	-382.5682	889.1364	0.0182	0.9898
HKYIG	61	-384.3085	890.6170	0.0087	0.9984
JC69	55	-392.9083	895.8166	0.0006	0.9991
JC69I	56	-392.8581	897.7163	0.0002	0.9993
JC69G	56	-392.9069	897.8137	0.0002	0.9996
F81	58	-391.1494	898.2988	0.0002	0.9997
JC69IG	57	-392.8639	899.7278	0.0001	0.9998
F81I	59	-391.0890	900.1780	0.0001	0.9999
F81G	59	-391.1472	900.2943	0.0001	1.0000
F81IG	60	-391.0949	902.1898	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
K2P	56	-386.4313	1040.0209	0.8130	0.8130
K2PI	57	-386.3750	1044.6791	0.0792	0.8921
K2PG	57	-386.4277	1044.7845	0.0751	0.9672
JC69	55	-392.9083	1048.2042	0.0136	0.9808
K2PIG	58	-386.3800	1049.4598	0.0073	0.9881
HKY	59	-384.3752	1050.2208	0.0050	0.9930
SYM	60	-382.6295	1051.5001	0.0026	0.9957
JC69I	56	-392.8581	1052.8746	0.0013	0.9970
JC69G	56	-392.9069	1052.9721	0.0013	0.9982
HKYI	60	-384.3027	1054.8465	0.0005	0.9987
HKYG	60	-384.3696	1054.9803	0.0005	0.9992
SYMI	61	-382.5612	1056.1342	0.0003	0.9994
SYMG	61	-382.6285	1056.2688	0.0002	0.9997
JC69IG	57	-392.8639	1057.6568	0.0001	0.9998
F81	58	-391.1494	1058.9985	0.0001	0.9999
GTR	63	-379.3236	1059.2004	0.0001	0.9999
HKYIG	61	-384.3085	1059.6287	0.0000	1.0000
SYMIG	62	-382.5682	1060.9189	0.0000	1.0000
GTRI	64	-379.0553	1063.4344	0.0000	1.0000
F81I	59	-391.0890	1063.6484	0.0000	1.0000
F81G	59	-391.1472	1063.7647	0.0000	1.0000
GTRG	64	-379.2691	1063.8619	0.0000	1.0000
GTRIG	65	-379.0660	1068.2266	0.0000	1.0000
F81IG	60	-391.0949	1068.4309	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 183-phylip/chr3_18256.phylip = 1 - 118;
 Partition Dummy = 1:183-phylip/chr3_18256.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 (K2P)]

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


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

BEGIN MRBAYES;
 Charset 183-phylip/chr3_18256.phylip = 1 - 118;
 Partition Dummy = 1:183-phylip/chr3_18256.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

