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

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

Minimum AIC  model: GTRI
Minimum AICc model: GTRI
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
GTRI	20	-1169.5159	2380.4459	0.2447	0.2447
GTRG	20	-1169.5958	2380.6058	0.2259	0.4707
HKYI	16	-1174.1471	2381.2038	0.1675	0.6382
HKYG	16	-1174.2285	2381.3667	0.1544	0.7926
GTRIG	21	-1169.1000	2381.7582	0.1270	0.9196
HKYIG	17	-1173.8298	2382.6846	0.0799	0.9995
GTR	19	-1177.4406	2394.1585	0.0003	0.9998
HKY	15	-1181.7981	2394.3976	0.0002	1.0000
F81I	15	-1189.1870	2409.1753	0.0000	1.0000
F81G	15	-1189.2315	2409.2643	0.0000	1.0000
F81IG	16	-1188.8988	2410.7072	0.0000	1.0000
F81	14	-1196.4724	2421.6447	0.0000	1.0000
K2PI	13	-1203.1719	2432.9495	0.0000	1.0000
K2PG	13	-1203.2667	2433.1391	0.0000	1.0000
K2PIG	14	-1202.8756	2434.4512	0.0000	1.0000
SYMI	17	-1200.9603	2436.9457	0.0000	1.0000
SYMG	17	-1201.0511	2437.1272	0.0000	1.0000
SYMIG	18	-1200.6280	2438.4037	0.0000	1.0000
K2P	12	-1210.8124	2446.1430	0.0000	1.0000
SYM	16	-1208.6358	2450.1814	0.0000	1.0000
JC69I	12	-1216.2906	2457.0996	0.0000	1.0000
JC69G	12	-1216.3461	2457.2105	0.0000	1.0000
JC69IG	13	-1216.0081	2458.6218	0.0000	1.0000
JC69	11	-1223.6488	2469.7354	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRI	20	-1169.5159	2379.0317	0.2650	0.2650
GTRG	20	-1169.5958	2379.1917	0.2446	0.5096
GTRIG	21	-1169.1000	2380.2000	0.1478	0.6574
HKYI	16	-1174.1471	2380.2941	0.1410	0.7984
HKYG	16	-1174.2285	2380.4570	0.1299	0.9283
HKYIG	17	-1173.8298	2381.6595	0.0712	0.9996
GTR	19	-1177.4406	2392.8812	0.0003	0.9998
HKY	15	-1181.7981	2393.5963	0.0002	1.0000
F81I	15	-1189.1870	2408.3739	0.0000	1.0000
F81G	15	-1189.2315	2408.4630	0.0000	1.0000
F81IG	16	-1188.8988	2409.7975	0.0000	1.0000
F81	14	-1196.4724	2420.9447	0.0000	1.0000
K2PI	13	-1203.1719	2432.3438	0.0000	1.0000
K2PG	13	-1203.2667	2432.5334	0.0000	1.0000
K2PIG	14	-1202.8756	2433.7512	0.0000	1.0000
SYMI	17	-1200.9603	2435.9206	0.0000	1.0000
SYMG	17	-1201.0511	2436.1021	0.0000	1.0000
SYMIG	18	-1200.6280	2437.2560	0.0000	1.0000
K2P	12	-1210.8124	2445.6248	0.0000	1.0000
SYM	16	-1208.6358	2449.2717	0.0000	1.0000
JC69I	12	-1216.2906	2456.5813	0.0000	1.0000
JC69G	12	-1216.3461	2456.6922	0.0000	1.0000
JC69IG	13	-1216.0081	2458.0162	0.0000	1.0000
JC69	11	-1223.6488	2469.2976	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYI	16	-1174.1471	2451.0401	0.5026	0.5026
HKYG	16	-1174.2285	2451.2029	0.4633	0.9660
HKYIG	17	-1173.8298	2456.8271	0.0278	0.9938
HKY	15	-1181.7981	2459.9206	0.0059	0.9997
GTRI	20	-1169.5159	2467.4642	0.0001	0.9999
GTRG	20	-1169.5958	2467.6241	0.0001	1.0000
GTRIG	21	-1169.1000	2473.0541	0.0000	1.0000
F81I	15	-1189.1870	2474.6983	0.0000	1.0000
F81G	15	-1189.2315	2474.7873	0.0000	1.0000
GTR	19	-1177.4406	2476.8920	0.0000	1.0000
F81IG	16	-1188.8988	2480.5435	0.0000	1.0000
F81	14	-1196.4724	2482.8474	0.0000	1.0000
K2PI	13	-1203.1719	2489.8249	0.0000	1.0000
K2PG	13	-1203.2667	2490.0145	0.0000	1.0000
K2PIG	14	-1202.8756	2495.6539	0.0000	1.0000
K2P	12	-1210.8124	2498.6842	0.0000	1.0000
JC69I	12	-1216.2906	2509.6407	0.0000	1.0000
JC69G	12	-1216.3461	2509.7516	0.0000	1.0000
SYMI	17	-1200.9603	2511.0882	0.0000	1.0000
SYMG	17	-1201.0511	2511.2697	0.0000	1.0000
JC69IG	13	-1216.0081	2515.4972	0.0000	1.0000
SYMIG	18	-1200.6280	2516.8452	0.0000	1.0000
JC69	11	-1223.6488	2517.9354	0.0000	1.0000
SYM	16	-1208.6358	2520.0176	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_11797.phylip = 1 - 615;
 Partition Dummy = 1:591-tomodel-phylip/chr5_11797.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_11797.phylip = 1 - 615;
 Partition Dummy = 1:591-tomodel-phylip/chr5_11797.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 (HKYI)]

BEGIN MRBAYES;
 Charset 591-tomodel-phylip/chr5_11797.phylip = 1 - 615;
 Partition Dummy = 1:591-tomodel-phylip/chr5_11797.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

