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

Input data from file "591-tomodel-phylip/chr4_7557.phylip" (7 taxa, 668 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	-2328.5824	4698.4631	0.3514	0.3514
GTRG	20	-2329.0570	4699.4122	0.2186	0.5700
HKYI	16	-2333.4780	4699.7916	0.1808	0.7509
GTRIG	21	-2328.5460	4700.5224	0.1255	0.8763
HKYIG	17	-2333.4867	4701.9149	0.0626	0.9389
HKYG	16	-2334.5631	4701.9619	0.0611	1.0000
HKY	15	-2351.7014	4734.1390	0.0000	1.0000
GTR	19	-2347.4910	4734.1548	0.0000	1.0000
K2PI	13	-2372.6397	4771.8360	0.0000	1.0000
K2PIG	14	-2372.6547	4773.9525	0.0000	1.0000
K2PG	13	-2373.8734	4774.3034	0.0000	1.0000
SYMI	17	-2372.3735	4779.6886	0.0000	1.0000
SYMIG	18	-2372.3880	4781.8299	0.0000	1.0000
SYMG	17	-2373.5787	4782.0990	0.0000	1.0000
F81I	15	-2384.7035	4800.1432	0.0000	1.0000
F81IG	16	-2384.7247	4802.2850	0.0000	1.0000
F81G	15	-2386.1546	4803.0455	0.0000	1.0000
K2P	12	-2390.7088	4805.8940	0.0000	1.0000
SYM	16	-2390.4895	4813.8147	0.0000	1.0000
F81	14	-2401.5940	4831.8312	0.0000	1.0000
JC69I	12	-2415.6239	4855.7241	0.0000	1.0000
JC69IG	13	-2415.6445	4857.8456	0.0000	1.0000
JC69G	12	-2417.0786	4858.6335	0.0000	1.0000
JC69	11	-2432.7257	4887.8538	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRI	20	-2328.5824	4697.1648	0.3705	0.3705
GTRG	20	-2329.0570	4698.1139	0.2305	0.6010
HKYI	16	-2333.4780	4698.9560	0.1513	0.7523
GTRIG	21	-2328.5460	4699.0920	0.1414	0.8937
HKYIG	17	-2333.4867	4700.9734	0.0552	0.9489
HKYG	16	-2334.5631	4701.1262	0.0511	1.0000
GTR	19	-2347.4910	4732.9819	0.0000	1.0000
HKY	15	-2351.7014	4733.4028	0.0000	1.0000
K2PI	13	-2372.6397	4771.2794	0.0000	1.0000
K2PIG	14	-2372.6547	4773.3094	0.0000	1.0000
K2PG	13	-2373.8734	4773.7468	0.0000	1.0000
SYMI	17	-2372.3735	4778.7470	0.0000	1.0000
SYMIG	18	-2372.3880	4780.7759	0.0000	1.0000
SYMG	17	-2373.5787	4781.1574	0.0000	1.0000
F81I	15	-2384.7035	4799.4070	0.0000	1.0000
F81IG	16	-2384.7247	4801.4494	0.0000	1.0000
F81G	15	-2386.1546	4802.3093	0.0000	1.0000
K2P	12	-2390.7088	4805.4176	0.0000	1.0000
SYM	16	-2390.4895	4812.9791	0.0000	1.0000
F81	14	-2401.5940	4831.1880	0.0000	1.0000
JC69I	12	-2415.6239	4855.2478	0.0000	1.0000
JC69IG	13	-2415.6445	4857.2891	0.0000	1.0000
JC69G	12	-2417.0786	4858.1571	0.0000	1.0000
JC69	11	-2432.7257	4887.4514	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYI	16	-2333.4780	4771.0246	0.7264	0.7264
HKYG	16	-2334.5631	4773.1949	0.2454	0.9718
HKYIG	17	-2333.4867	4777.5463	0.0279	0.9996
GTRI	20	-2328.5824	4787.2505	0.0002	0.9999
GTRG	20	-2329.0570	4788.1997	0.0001	1.0000
GTRIG	21	-2328.5460	4793.6821	0.0000	1.0000
HKY	15	-2351.7014	4800.9671	0.0000	1.0000
GTR	19	-2347.4910	4818.5634	0.0000	1.0000
K2PI	13	-2372.6397	4829.8352	0.0000	1.0000
K2PG	13	-2373.8734	4832.3025	0.0000	1.0000
K2PIG	14	-2372.6547	4836.3694	0.0000	1.0000
SYMI	17	-2372.3735	4855.3199	0.0000	1.0000
SYMG	17	-2373.5787	4857.7303	0.0000	1.0000
K2P	12	-2390.7088	4859.4691	0.0000	1.0000
SYMIG	18	-2372.3880	4861.8531	0.0000	1.0000
F81I	15	-2384.7035	4866.9714	0.0000	1.0000
F81G	15	-2386.1546	4869.8736	0.0000	1.0000
F81IG	16	-2384.7247	4873.5180	0.0000	1.0000
SYM	16	-2390.4895	4885.0477	0.0000	1.0000
F81	14	-2401.5940	4894.2481	0.0000	1.0000
JC69I	12	-2415.6239	4909.2993	0.0000	1.0000
JC69G	12	-2417.0786	4912.2086	0.0000	1.0000
JC69IG	13	-2415.6445	4915.8448	0.0000	1.0000
JC69	11	-2432.7257	4936.9985	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/chr4_7557.phylip = 1 - 668;
 Partition Dummy = 1:591-tomodel-phylip/chr4_7557.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/chr4_7557.phylip = 1 - 668;
 Partition Dummy = 1:591-tomodel-phylip/chr4_7557.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/chr4_7557.phylip = 1 - 668;
 Partition Dummy = 1:591-tomodel-phylip/chr4_7557.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

