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

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

Minimum AIC  model: GTRG
Minimum AICc model: GTRG
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
GTRG	20	-2586.6643	5214.4239	0.6231	0.6231
GTRIG	21	-2586.6884	5216.5830	0.2117	0.8348
HKYG	16	-2592.7614	5218.2283	0.0930	0.9278
GTRI	20	-2589.5143	5220.1238	0.0360	0.9638
HKYIG	17	-2592.8574	5220.5095	0.0297	0.9935
HKYI	16	-2595.4271	5223.5597	0.0065	1.0000
SYMG	17	-2618.1227	5271.0402	0.0000	1.0000
K2PG	13	-2623.0580	5272.5862	0.0000	1.0000
SYMIG	18	-2618.1405	5273.1704	0.0000	1.0000
K2PIG	14	-2623.0584	5274.6601	0.0000	1.0000
SYMI	17	-2621.0982	5276.9911	0.0000	1.0000
K2PI	13	-2625.8295	5278.1293	0.0000	1.0000
GTR	19	-2644.2865	5327.5626	0.0000	1.0000
HKY	15	-2650.0158	5330.6534	0.0000	1.0000
F81G	15	-2673.4196	5377.4609	0.0000	1.0000
F81IG	16	-2673.4252	5379.5560	0.0000	1.0000
F81I	15	-2675.2296	5381.0809	0.0000	1.0000
K2P	12	-2681.8858	5388.1741	0.0000	1.0000
SYM	16	-2678.2717	5389.2489	0.0000	1.0000
JC69G	12	-2697.4520	5419.3066	0.0000	1.0000
JC69IG	13	-2696.4807	5419.4318	0.0000	1.0000
JC69I	12	-2698.3761	5421.1547	0.0000	1.0000
F81	14	-2723.3037	5475.1507	0.0000	1.0000
JC69	11	-2747.6555	5517.6512	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRG	20	-2586.6643	5213.3287	0.6292	0.6292
GTRIG	21	-2586.6884	5215.3768	0.2260	0.8551
HKYG	16	-2592.7614	5217.5227	0.0773	0.9324
GTRI	20	-2589.5143	5219.0286	0.0364	0.9688
HKYIG	17	-2592.8574	5219.7147	0.0258	0.9946
HKYI	16	-2595.4271	5222.8542	0.0054	1.0000
SYMG	17	-2618.1227	5270.2454	0.0000	1.0000
K2PG	13	-2623.0580	5272.1160	0.0000	1.0000
SYMIG	18	-2618.1405	5272.2810	0.0000	1.0000
K2PIG	14	-2623.0584	5274.1167	0.0000	1.0000
SYMI	17	-2621.0982	5276.1963	0.0000	1.0000
K2PI	13	-2625.8295	5277.6590	0.0000	1.0000
GTR	19	-2644.2865	5326.5730	0.0000	1.0000
HKY	15	-2650.0158	5330.0316	0.0000	1.0000
F81G	15	-2673.4196	5376.8392	0.0000	1.0000
F81IG	16	-2673.4252	5378.8504	0.0000	1.0000
F81I	15	-2675.2296	5380.4592	0.0000	1.0000
K2P	12	-2681.8858	5387.7715	0.0000	1.0000
SYM	16	-2678.2717	5388.5433	0.0000	1.0000
JC69G	12	-2697.4520	5418.9040	0.0000	1.0000
JC69IG	13	-2696.4807	5418.9615	0.0000	1.0000
JC69I	12	-2698.3761	5420.7521	0.0000	1.0000
F81	14	-2723.3037	5474.6074	0.0000	1.0000
JC69	11	-2747.6555	5517.3110	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	16	-2592.7614	5292.2347	0.9069	0.9069
HKYI	16	-2595.4271	5297.5661	0.0631	0.9699
HKYIG	17	-2592.8574	5299.0962	0.0293	0.9993
GTRG	20	-2586.6643	5306.7186	0.0006	0.9999
GTRI	20	-2589.5143	5312.4186	0.0000	1.0000
GTRIG	21	-2586.6884	5313.4362	0.0000	1.0000
K2PG	13	-2623.0580	5332.8194	0.0000	1.0000
K2PI	13	-2625.8295	5338.3625	0.0000	1.0000
K2PIG	14	-2623.0584	5339.4897	0.0000	1.0000
SYMG	17	-2618.1227	5349.6269	0.0000	1.0000
SYMI	17	-2621.0982	5355.5778	0.0000	1.0000
SYMIG	18	-2618.1405	5356.3319	0.0000	1.0000
HKY	15	-2650.0158	5400.0741	0.0000	1.0000
GTR	19	-2644.2865	5415.2935	0.0000	1.0000
K2P	12	-2681.8858	5443.8055	0.0000	1.0000
F81G	15	-2673.4196	5446.8817	0.0000	1.0000
F81I	15	-2675.2296	5450.5016	0.0000	1.0000
F81IG	16	-2673.4252	5453.5624	0.0000	1.0000
SYM	16	-2678.2717	5463.2553	0.0000	1.0000
JC69G	12	-2697.4520	5474.9380	0.0000	1.0000
JC69I	12	-2698.3761	5476.7861	0.0000	1.0000
JC69IG	13	-2696.4807	5479.6650	0.0000	1.0000
F81	14	-2723.3037	5539.9804	0.0000	1.0000
JC69	11	-2747.6555	5568.6754	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/chr7_10377.phylip = 1 - 788;
 Partition Dummy = 1:591-tomodel-phylip/chr7_10377.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/chr7_10377.phylip = 1 - 788;
 Partition Dummy = 1:591-tomodel-phylip/chr7_10377.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 (HKYG)]

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

