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

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

Minimum AIC  model: GTR
Minimum AICc model: GTR
Minimum BIC  model: HKY


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
GTR	19	-1165.0755	2369.7149	0.2876	0.2876
HKY	15	-1169.9573	2370.8943	0.1594	0.4470
GTRG	20	-1164.7025	2371.1370	0.1412	0.5882
HKYG	16	-1169.1634	2371.4392	0.1214	0.7097
HKYI	16	-1169.3318	2371.7762	0.1026	0.8123
GTRI	20	-1165.0611	2371.8541	0.0987	0.9109
GTRIG	21	-1164.7025	2373.3141	0.0476	0.9585
HKYIG	17	-1169.1656	2373.5852	0.0415	1.0000
K2P	12	-1189.0658	2402.7645	0.0000	1.0000
K2PG	13	-1188.3750	2403.4899	0.0000	1.0000
SYM	16	-1185.1946	2403.5016	0.0000	1.0000
K2PI	13	-1188.5376	2403.8150	0.0000	1.0000
SYMI	17	-1184.7965	2404.8471	0.0000	1.0000
K2PIG	14	-1188.3750	2405.6054	0.0000	1.0000
SYMG	17	-1185.1817	2405.6175	0.0000	1.0000
SYMIG	18	-1185.1817	2407.7680	0.0000	1.0000
F81	14	-1192.8083	2414.4720	0.0000	1.0000
F81G	15	-1192.3117	2415.6029	0.0000	1.0000
F81I	15	-1192.5384	2416.0565	0.0000	1.0000
F81IG	16	-1192.3117	2417.7358	0.0000	1.0000
JC69	11	-1210.0046	2442.5437	0.0000	1.0000
JC69G	12	-1209.5425	2443.7178	0.0000	1.0000
JC69I	12	-1209.7449	2444.1227	0.0000	1.0000
JC69IG	13	-1209.5425	2445.8248	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTR	19	-1165.0755	2368.1511	0.3064	0.3064
GTRG	20	-1164.7025	2369.4050	0.1637	0.4701
HKY	15	-1169.9573	2369.9147	0.1269	0.5970
GTRI	20	-1165.0611	2370.1221	0.1144	0.7114
HKYG	16	-1169.1634	2370.3267	0.1033	0.8146
HKYI	16	-1169.3318	2370.6637	0.0872	0.9019
GTRIG	21	-1164.7025	2371.4050	0.0602	0.9621
HKYIG	17	-1169.1656	2372.3311	0.0379	1.0000
K2P	12	-1189.0658	2402.1317	0.0000	1.0000
SYM	16	-1185.1946	2402.3892	0.0000	1.0000
K2PG	13	-1188.3750	2402.7501	0.0000	1.0000
K2PI	13	-1188.5376	2403.0751	0.0000	1.0000
SYMI	17	-1184.7965	2403.5930	0.0000	1.0000
SYMG	17	-1185.1817	2404.3634	0.0000	1.0000
K2PIG	14	-1188.3750	2404.7500	0.0000	1.0000
SYMIG	18	-1185.1817	2406.3635	0.0000	1.0000
F81	14	-1192.8083	2413.6166	0.0000	1.0000
F81G	15	-1192.3117	2414.6233	0.0000	1.0000
F81I	15	-1192.5384	2415.0769	0.0000	1.0000
F81IG	16	-1192.3117	2416.6234	0.0000	1.0000
JC69	11	-1210.0046	2442.0093	0.0000	1.0000
JC69G	12	-1209.5425	2443.0850	0.0000	1.0000
JC69I	12	-1209.7449	2443.4899	0.0000	1.0000
JC69IG	13	-1209.5425	2445.0850	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKY	15	-1169.9573	2433.3127	0.8429	0.8429
HKYG	16	-1169.1634	2437.9513	0.0829	0.9257
HKYI	16	-1169.3318	2438.2883	0.0700	0.9958
HKYIG	17	-1169.1656	2444.1823	0.0037	0.9995
GTR	19	-1165.0755	2448.4553	0.0004	0.9999
K2P	12	-1189.0658	2452.8501	0.0000	0.9999
GTRG	20	-1164.7025	2453.9358	0.0000	1.0000
GTRI	20	-1165.0611	2454.6529	0.0000	1.0000
K2PG	13	-1188.3750	2457.6951	0.0000	1.0000
K2PI	13	-1188.5376	2458.0201	0.0000	1.0000
GTRIG	21	-1164.7025	2460.1623	0.0000	1.0000
K2PIG	14	-1188.3750	2463.9215	0.0000	1.0000
SYM	16	-1185.1946	2470.0137	0.0000	1.0000
F81	14	-1192.8083	2472.7881	0.0000	1.0000
SYMI	17	-1184.7965	2475.4441	0.0000	1.0000
SYMG	17	-1185.1817	2476.2146	0.0000	1.0000
F81G	15	-1192.3117	2478.0214	0.0000	1.0000
F81I	15	-1192.5384	2478.4749	0.0000	1.0000
SYMIG	18	-1185.1817	2482.4411	0.0000	1.0000
F81IG	16	-1192.3117	2484.2479	0.0000	1.0000
JC69	11	-1210.0046	2488.5012	0.0000	1.0000
JC69G	12	-1209.5425	2493.8034	0.0000	1.0000
JC69I	12	-1209.7449	2494.2083	0.0000	1.0000
JC69IG	13	-1209.5425	2500.0300	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 591-tomodel-phylip/chr3_5766.phylip = 1 - 506;
 Partition Dummy = 1:591-tomodel-phylip/chr3_5766.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 (GTR)]

BEGIN MRBAYES;
 Charset 591-tomodel-phylip/chr3_5766.phylip = 1 - 506;
 Partition Dummy = 1:591-tomodel-phylip/chr3_5766.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 BIC model (HKY)]

BEGIN MRBAYES;
 Charset 591-tomodel-phylip/chr3_5766.phylip = 1 - 506;
 Partition Dummy = 1:591-tomodel-phylip/chr3_5766.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=2 rates=equal;
 Prset applyto=(1) statefreqpr=Dirichlet(1.0,1.0,1.0,1.0) tratiopr=Beta(1.0,1.0);
END;


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

