
Output from mraic.pl version 1.4.4 by Johan Nylander
-------------------------------------------------------------
16:17 04/04/2011

Input data from file "683-tomodel-phylip/chr8_4342.phylip" (7 taxa, 611 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	-1163.2633	2365.8125	0.3501	0.3501
GTRI	20	-1162.4688	2366.3612	0.2661	0.6162
GTRG	20	-1162.5131	2366.4499	0.2546	0.8708
GTRIG	21	-1162.5164	2368.6015	0.0868	0.9576
HKY	15	-1170.4725	2371.7517	0.0180	0.9756
HKYI	16	-1169.9567	2372.8293	0.0105	0.9861
HKYG	16	-1169.9812	2372.8783	0.0102	0.9963
HKYIG	17	-1169.9798	2374.9916	0.0036	0.9999
SYM	16	-1175.2949	2383.5057	0.0001	0.9999
SYMI	17	-1174.6564	2384.3449	0.0000	0.9999
SYMG	17	-1174.6907	2384.4133	0.0000	1.0000
SYMIG	18	-1174.6907	2386.5368	0.0000	1.0000
K2P	12	-1182.0482	2388.6181	0.0000	1.0000
K2PI	13	-1181.5800	2389.7696	0.0000	1.0000
K2PG	13	-1181.5991	2389.8079	0.0000	1.0000
F81	14	-1181.4453	2391.5953	0.0000	1.0000
K2PIG	14	-1181.5974	2391.8994	0.0000	1.0000
F81I	15	-1181.0475	2392.9016	0.0000	1.0000
F81G	15	-1181.0673	2392.9412	0.0000	1.0000
F81IG	16	-1181.0644	2395.0445	0.0000	1.0000
JC69	11	-1192.6269	2407.6946	0.0000	1.0000
JC69I	12	-1192.2431	2409.0079	0.0000	1.0000
JC69G	12	-1192.2623	2409.0464	0.0000	1.0000
JC69IG	13	-1192.2594	2411.1284	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTR	19	-1163.2633	2364.5265	0.3358	0.3358
GTRI	20	-1162.4688	2364.9375	0.2735	0.6093
GTRG	20	-1162.5131	2365.0261	0.2616	0.8709
GTRIG	21	-1162.5164	2367.0327	0.0959	0.9668
HKY	15	-1170.4725	2370.9450	0.0136	0.9804
HKYI	16	-1169.9567	2371.9135	0.0084	0.9887
HKYG	16	-1169.9812	2371.9625	0.0082	0.9969
HKYIG	17	-1169.9798	2373.9596	0.0030	0.9999
SYM	16	-1175.2949	2382.5899	0.0000	0.9999
SYMI	17	-1174.6564	2383.3128	0.0000	1.0000
SYMG	17	-1174.6907	2383.3813	0.0000	1.0000
SYMIG	18	-1174.6907	2385.3814	0.0000	1.0000
K2P	12	-1182.0482	2388.0964	0.0000	1.0000
K2PI	13	-1181.5800	2389.1599	0.0000	1.0000
K2PG	13	-1181.5991	2389.1982	0.0000	1.0000
F81	14	-1181.4453	2390.8906	0.0000	1.0000
K2PIG	14	-1181.5974	2391.1947	0.0000	1.0000
F81I	15	-1181.0475	2392.0949	0.0000	1.0000
F81G	15	-1181.0673	2392.1345	0.0000	1.0000
F81IG	16	-1181.0644	2394.1287	0.0000	1.0000
JC69	11	-1192.6269	2407.2538	0.0000	1.0000
JC69I	12	-1192.2431	2408.4861	0.0000	1.0000
JC69G	12	-1192.2623	2408.5247	0.0000	1.0000
JC69IG	13	-1192.2594	2410.5187	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKY	15	-1170.4725	2437.1715	0.7680	0.7680
K2P	12	-1182.0482	2441.0775	0.1089	0.8769
HKYI	16	-1169.9567	2442.5550	0.0520	0.9289
HKYG	16	-1169.9812	2442.6040	0.0508	0.9797
K2PI	13	-1181.5800	2446.5562	0.0070	0.9867
K2PG	13	-1181.5991	2446.5944	0.0069	0.9937
GTR	19	-1163.2633	2448.4134	0.0028	0.9964
HKYIG	17	-1169.9798	2449.0162	0.0021	0.9985
F81	14	-1181.4453	2452.7019	0.0003	0.9988
K2PIG	14	-1181.5974	2453.0061	0.0003	0.9991
SYM	16	-1175.2949	2453.2315	0.0003	0.9993
GTRI	20	-1162.4688	2453.2395	0.0002	0.9996
GTRG	20	-1162.5131	2453.3281	0.0002	0.9998
JC69	11	-1192.6269	2455.8199	0.0001	0.9999
F81I	15	-1181.0475	2458.3214	0.0000	0.9999
F81G	15	-1181.0673	2458.3610	0.0000	0.9999
SYMI	17	-1174.6564	2458.3695	0.0000	1.0000
SYMG	17	-1174.6907	2458.4379	0.0000	1.0000
GTRIG	21	-1162.5164	2459.7497	0.0000	1.0000
JC69I	12	-1192.2431	2461.4673	0.0000	1.0000
JC69G	12	-1192.2623	2461.5059	0.0000	1.0000
F81IG	16	-1181.0644	2464.7703	0.0000	1.0000
SYMIG	18	-1174.6907	2464.8531	0.0000	1.0000
JC69IG	13	-1192.2594	2467.9150	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 683-tomodel-phylip/chr8_4342.phylip = 1 - 611;
 Partition Dummy = 1:683-tomodel-phylip/chr8_4342.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 683-tomodel-phylip/chr8_4342.phylip = 1 - 611;
 Partition Dummy = 1:683-tomodel-phylip/chr8_4342.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 683-tomodel-phylip/chr8_4342.phylip = 1 - 611;
 Partition Dummy = 1:683-tomodel-phylip/chr8_4342.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

