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

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

Minimum AIC  model: HKYIG
Minimum AICc model: HKYG
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
HKYG	60	-1134.8643	2411.2581	0.4204	0.4204
HKYIG	61	-1133.7861	2411.8848	0.3073	0.7278
HKYI	60	-1135.5471	2412.6237	0.2124	0.9402
GTRI	64	-1132.1446	2417.0511	0.0232	0.9634
GTRG	64	-1132.1885	2417.1390	0.0222	0.9856
GTRIG	65	-1131.3776	2418.3671	0.0120	0.9976
K2PG	57	-1144.9940	2423.2650	0.0010	0.9987
K2PIG	58	-1143.9975	2424.0066	0.0007	0.9994
K2PI	57	-1145.6554	2424.5878	0.0005	0.9999
SYMG	61	-1143.0562	2430.4251	0.0000	1.0000
SYMI	61	-1143.4775	2431.2676	0.0000	1.0000
SYMIG	62	-1142.4181	2431.9487	0.0000	1.0000
F81G	59	-1154.2689	2447.3003	0.0000	1.0000
HKY	59	-1154.4720	2447.7064	0.0000	1.0000
F81IG	60	-1153.3282	2448.1857	0.0000	1.0000
F81I	59	-1155.0012	2448.7648	0.0000	1.0000
JC69G	56	-1162.2467	2455.0516	0.0000	1.0000
GTR	63	-1152.8588	2455.6464	0.0000	1.0000
JC69IG	57	-1161.3222	2455.9213	0.0000	1.0000
JC69I	56	-1163.0633	2456.6847	0.0000	1.0000
K2P	56	-1163.1234	2456.8049	0.0000	1.0000
SYM	60	-1161.5625	2464.6544	0.0000	1.0000
F81	58	-1172.7822	2481.5762	0.0000	1.0000
JC69	55	-1180.1522	2488.1595	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
HKYIG	61	-1133.7861	2389.5721	0.3225	0.3225
HKYG	60	-1134.8643	2389.7287	0.2982	0.6208
HKYI	60	-1135.5471	2391.0943	0.1507	0.7714
GTRI	64	-1132.1446	2392.2892	0.0829	0.8543
GTRG	64	-1132.1885	2392.3771	0.0793	0.9337
GTRIG	65	-1131.3776	2392.7552	0.0657	0.9993
K2PG	57	-1144.9940	2403.9881	0.0002	0.9996
K2PIG	58	-1143.9975	2403.9949	0.0002	0.9998
K2PI	57	-1145.6554	2405.3109	0.0001	0.9999
SYMG	61	-1143.0562	2408.1124	0.0000	1.0000
SYMIG	62	-1142.4181	2408.8363	0.0000	1.0000
SYMI	61	-1143.4775	2408.9550	0.0000	1.0000
F81G	59	-1154.2689	2426.5379	0.0000	1.0000
F81IG	60	-1153.3282	2426.6563	0.0000	1.0000
HKY	59	-1154.4720	2426.9439	0.0000	1.0000
F81I	59	-1155.0012	2428.0024	0.0000	1.0000
GTR	63	-1152.8588	2431.7176	0.0000	1.0000
JC69G	56	-1162.2467	2436.4935	0.0000	1.0000
JC69IG	57	-1161.3222	2436.6444	0.0000	1.0000
JC69I	56	-1163.0633	2438.1266	0.0000	1.0000
K2P	56	-1163.1234	2438.2468	0.0000	1.0000
SYM	60	-1161.5625	2443.1250	0.0000	1.0000
F81	58	-1172.7822	2461.5645	0.0000	1.0000
JC69	55	-1180.1522	2470.3044	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	60	-1134.8643	2629.3664	0.4585	0.4585
HKYI	60	-1135.5471	2630.7320	0.2316	0.6901
K2PG	57	-1144.9940	2631.6439	0.1468	0.8370
K2PI	57	-1145.6554	2632.9667	0.0758	0.9127
HKYIG	61	-1133.7861	2633.2038	0.0673	0.9800
K2PIG	58	-1143.9975	2635.6447	0.0199	0.9999
GTRI	64	-1132.1446	2647.9027	0.0000	0.9999
GTRG	64	-1132.1885	2647.9906	0.0000	1.0000
SYMG	61	-1143.0562	2651.7441	0.0000	1.0000
GTRIG	65	-1131.3776	2652.3627	0.0000	1.0000
SYMI	61	-1143.4775	2652.5866	0.0000	1.0000
SYMIG	62	-1142.4181	2656.4619	0.0000	1.0000
JC69G	56	-1162.2467	2660.1553	0.0000	1.0000
JC69I	56	-1163.0633	2661.7884	0.0000	1.0000
K2P	56	-1163.1234	2661.9086	0.0000	1.0000
F81G	59	-1154.2689	2662.1816	0.0000	1.0000
HKY	59	-1154.4720	2662.5876	0.0000	1.0000
F81I	59	-1155.0012	2663.6461	0.0000	1.0000
JC69IG	57	-1161.3222	2664.3002	0.0000	1.0000
F81IG	60	-1153.3282	2666.2940	0.0000	1.0000
SYM	60	-1161.5625	2682.7626	0.0000	1.0000
GTR	63	-1152.8588	2683.3372	0.0000	1.0000
JC69	55	-1180.1522	2689.9723	0.0000	1.0000
F81	58	-1172.7822	2693.2142	0.0000	1.0000



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

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

BEGIN MRBAYES;
 Charset 183-phylip/chr6_8057.phylip = 1 - 401;
 Partition Dummy = 1:183-phylip/chr6_8057.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=2 rates=invgamma;
 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) pinvarpr=Uniform(0.0,1.0);
END;


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

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


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

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

