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

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

Minimum AIC  model: GTRG
Minimum AICc model: GTRG
Minimum BIC  model: HKYI


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	-1200.8816	2443.5736	0.5787	0.5787
GTRIG	21	-1200.8815	2445.7587	0.1941	0.7727
GTRI	20	-1202.0361	2445.8826	0.1824	0.9551
HKYI	16	-1208.6935	2450.5493	0.0177	0.9728
HKYG	16	-1208.7032	2450.5688	0.0175	0.9903
HKYIG	17	-1208.2249	2451.7603	0.0097	1.0000
GTR	19	-1214.7402	2469.1148	0.0000	1.0000
K2PI	13	-1222.1945	2471.1618	0.0000	1.0000
K2PG	13	-1222.2630	2471.2987	0.0000	1.0000
SYMI	17	-1218.1735	2471.6575	0.0000	1.0000
SYMG	17	-1218.2626	2471.8358	0.0000	1.0000
HKY	15	-1220.7913	2472.6061	0.0000	1.0000
SYMIG	18	-1217.7457	2472.9592	0.0000	1.0000
K2PIG	14	-1222.0518	2472.9972	0.0000	1.0000
K2P	12	-1234.1996	2493.0603	0.0000	1.0000
SYM	16	-1230.1075	2493.3774	0.0000	1.0000
F81G	15	-1231.7774	2494.5782	0.0000	1.0000
F81I	15	-1231.8115	2494.6464	0.0000	1.0000
F81IG	16	-1231.7794	2496.7212	0.0000	1.0000
JC69G	12	-1243.5939	2511.8488	0.0000	1.0000
JC69I	12	-1243.6001	2511.8613	0.0000	1.0000
JC69IG	13	-1243.5985	2513.9697	0.0000	1.0000
F81	14	-1243.0834	2515.0605	0.0000	1.0000
JC69	11	-1254.8777	2532.3136	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRG	20	-1200.8816	2441.7633	0.5747	0.5747
GTRIG	21	-1200.8815	2443.7630	0.2114	0.7861
GTRI	20	-1202.0361	2444.0723	0.1811	0.9672
HKYI	16	-1208.6935	2449.3869	0.0127	0.9799
HKYG	16	-1208.7032	2449.4064	0.0126	0.9925
HKYIG	17	-1208.2249	2450.4498	0.0075	1.0000
GTR	19	-1214.7402	2467.4804	0.0000	1.0000
SYMI	17	-1218.1735	2470.3470	0.0000	1.0000
K2PI	13	-1222.1945	2470.3890	0.0000	1.0000
SYMG	17	-1218.2626	2470.5253	0.0000	1.0000
K2PG	13	-1222.2630	2470.5259	0.0000	1.0000
SYMIG	18	-1217.7457	2471.4913	0.0000	1.0000
HKY	15	-1220.7913	2471.5826	0.0000	1.0000
K2PIG	14	-1222.0518	2472.1036	0.0000	1.0000
SYM	16	-1230.1075	2492.2150	0.0000	1.0000
K2P	12	-1234.1996	2492.3993	0.0000	1.0000
F81G	15	-1231.7774	2493.5547	0.0000	1.0000
F81I	15	-1231.8115	2493.6229	0.0000	1.0000
F81IG	16	-1231.7794	2495.5588	0.0000	1.0000
JC69G	12	-1243.5939	2511.1878	0.0000	1.0000
JC69I	12	-1243.6001	2511.2002	0.0000	1.0000
JC69IG	13	-1243.5985	2513.1969	0.0000	1.0000
F81	14	-1243.0834	2514.1669	0.0000	1.0000
JC69	11	-1254.8777	2531.7555	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYI	16	-1208.6935	2516.3333	0.4748	0.4748
HKYG	16	-1208.7032	2516.3528	0.4702	0.9449
HKYIG	17	-1208.2249	2521.5803	0.0344	0.9794
K2PI	13	-1222.1945	2524.7829	0.0069	0.9863
K2PG	13	-1222.2630	2524.9198	0.0065	0.9928
GTRG	20	-1200.8816	2525.4462	0.0050	0.9978
GTRI	20	-1202.0361	2527.7552	0.0016	0.9993
K2PIG	14	-1222.0518	2530.6817	0.0004	0.9997
GTRIG	21	-1200.8815	2531.6301	0.0002	0.9999
HKY	15	-1220.7913	2534.3449	0.0001	1.0000
SYMI	17	-1218.1735	2541.4775	0.0000	1.0000
SYMG	17	-1218.2626	2541.6558	0.0000	1.0000
K2P	12	-1234.1996	2542.6090	0.0000	1.0000
SYMIG	18	-1217.7457	2546.8060	0.0000	1.0000
GTR	19	-1214.7402	2546.9792	0.0000	1.0000
F81G	15	-1231.7774	2556.3170	0.0000	1.0000
F81I	15	-1231.8115	2556.3852	0.0000	1.0000
SYM	16	-1230.1075	2559.1614	0.0000	1.0000
JC69G	12	-1243.5939	2561.3976	0.0000	1.0000
JC69I	12	-1243.6001	2561.4100	0.0000	1.0000
F81IG	16	-1231.7794	2562.5052	0.0000	1.0000
JC69IG	13	-1243.5985	2567.5909	0.0000	1.0000
F81	14	-1243.0834	2572.7449	0.0000	1.0000
JC69	11	-1254.8777	2577.7811	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 591-tomodel-phylip/chr3_9749.phylip = 1 - 485;
 Partition Dummy = 1:591-tomodel-phylip/chr3_9749.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 591-tomodel-phylip/chr3_9749.phylip = 1 - 485;
 Partition Dummy = 1:591-tomodel-phylip/chr3_9749.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 (HKYI)]

BEGIN MRBAYES;
 Charset 591-tomodel-phylip/chr3_9749.phylip = 1 - 485;
 Partition Dummy = 1:591-tomodel-phylip/chr3_9749.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=2 rates=propinv;
 Prset applyto=(1) statefreqpr=Dirichlet(1.0,1.0,1.0,1.0) tratiopr=Beta(1.0,1.0) pinvarpr=Uniform(0.0,1.0);
END;


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

