
Output from mraic.pl version 1.4.4 by Johan Nylander
-------------------------------------------------------------
12:08 08/01/2011

Input data from file "261-aggressive-phylip/8724_hg18.phylip" (29 taxa, 214 characters)

Minimum AIC  model: HKYI
Minimum AICc model: HKYI
Minimum BIC  model: K2PI


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
HKYI	60	-392.5419	952.9270	0.4072	0.4072
HKYG	60	-392.8959	953.6350	0.2858	0.6930
HKY	59	-395.9929	955.9599	0.0894	0.7823
K2PI	57	-400.1597	956.7041	0.0616	0.8439
HKYIG	61	-392.5420	956.8472	0.0573	0.9013
K2PG	57	-400.4837	957.3520	0.0446	0.9459
K2P	56	-403.5199	959.7022	0.0138	0.9596
K2PIG	58	-400.1602	960.4753	0.0093	0.9690
GTRI	64	-388.5647	960.9683	0.0073	0.9763
GTRG	64	-388.8542	961.5474	0.0055	0.9817
F81I	59	-398.8037	961.5813	0.0054	0.9871
F81G	59	-399.1446	962.2633	0.0038	0.9909
GTR	63	-391.9335	963.6270	0.0019	0.9929
JC69I	56	-405.8562	964.3748	0.0013	0.9942
F81	58	-402.2065	964.5679	0.0012	0.9954
GTRIG	65	-388.5181	965.0092	0.0010	0.9964
JC69G	56	-406.1899	965.0422	0.0010	0.9973
F81IG	60	-398.8041	965.4514	0.0008	0.9981
SYMI	61	-397.0608	965.8848	0.0006	0.9987
SYMG	61	-397.2804	966.3241	0.0005	0.9992
JC69	55	-409.2384	967.4641	0.0003	0.9995
JC69IG	57	-405.8565	968.0976	0.0002	0.9997
SYM	60	-400.1797	968.2025	0.0002	0.9999
SYMIG	62	-397.0607	969.8565	0.0001	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
HKYI	60	-392.5419	905.0839	0.2285	0.2285
GTRI	64	-388.5647	905.1294	0.2234	0.4519
GTRG	64	-388.8542	905.7084	0.1672	0.6192
HKYG	60	-392.8959	905.7919	0.1604	0.7796
GTRIG	65	-388.5181	907.0362	0.0861	0.8657
HKYIG	61	-392.5420	907.0840	0.0841	0.9497
GTR	63	-391.9335	909.8670	0.0209	0.9707
HKY	59	-395.9929	909.9859	0.0197	0.9904
K2PI	57	-400.1597	914.3195	0.0023	0.9926
K2PG	57	-400.4837	914.9674	0.0016	0.9942
F81I	59	-398.8037	915.6073	0.0012	0.9954
SYMI	61	-397.0608	916.1216	0.0009	0.9963
F81G	59	-399.1446	916.2892	0.0008	0.9972
K2PIG	58	-400.1602	916.3204	0.0008	0.9980
SYMG	61	-397.2804	916.5609	0.0007	0.9988
F81IG	60	-398.8041	917.6083	0.0004	0.9992
SYMIG	62	-397.0607	918.1214	0.0003	0.9995
K2P	56	-403.5199	919.0398	0.0002	0.9997
SYM	60	-400.1797	920.3593	0.0001	0.9998
F81	58	-402.2065	920.4130	0.0001	1.0000
JC69I	56	-405.8562	923.7124	0.0000	1.0000
JC69G	56	-406.1899	924.3798	0.0000	1.0000
JC69IG	57	-405.8565	925.7130	0.0000	1.0000
JC69	55	-409.2384	928.4767	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
K2PI	57	-400.1597	1106.1801	0.2559	0.2559
K2PG	57	-400.4837	1106.8280	0.1851	0.4410
HKYI	60	-392.5419	1107.0424	0.1663	0.6073
K2P	56	-403.5199	1107.5345	0.1300	0.7373
HKYG	60	-392.8959	1107.7504	0.1167	0.8540
HKY	59	-395.9929	1108.5785	0.0771	0.9312
K2PIG	58	-400.1602	1111.5470	0.0175	0.9487
JC69I	56	-405.8562	1112.2070	0.0126	0.9612
HKYIG	61	-392.5420	1112.4086	0.0114	0.9726
JC69G	56	-406.1899	1112.8745	0.0090	0.9816
JC69	55	-409.2384	1113.6054	0.0062	0.9879
F81I	59	-398.8037	1114.1999	0.0046	0.9925
F81G	59	-399.1446	1114.8818	0.0033	0.9958
F81	58	-402.2065	1115.6396	0.0023	0.9981
JC69IG	57	-405.8565	1117.5736	0.0009	0.9989
F81IG	60	-398.8041	1119.5668	0.0003	0.9992
GTRI	64	-388.5647	1120.5518	0.0002	0.9994
GTRG	64	-388.8542	1121.1309	0.0001	0.9996
SYMI	61	-397.0608	1121.4462	0.0001	0.9997
SYMG	61	-397.2804	1121.8854	0.0001	0.9998
GTR	63	-391.9335	1121.9235	0.0001	0.9999
SYM	60	-400.1797	1122.3179	0.0001	1.0000
GTRIG	65	-388.5181	1125.8246	0.0000	1.0000
SYMIG	62	-397.0607	1126.8119	0.0000	1.0000



To specify default values for best AIC, AICc, or BIC models in MrBayes:

[Mrbayes block for the best AIC model (HKYI)]

BEGIN MRBAYES;
 Charset 261-aggressive-phylip/8724_hg18.phylip = 1 - 214;
 Partition Dummy = 1:261-aggressive-phylip/8724_hg18.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;


[Mrbayes block for the best AICc model (HKYI)]

BEGIN MRBAYES;
 Charset 261-aggressive-phylip/8724_hg18.phylip = 1 - 214;
 Partition Dummy = 1:261-aggressive-phylip/8724_hg18.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;


[Mrbayes block for the best BIC model (K2PI)]

BEGIN MRBAYES;
 Charset 261-aggressive-phylip/8724_hg18.phylip = 1 - 214;
 Partition Dummy = 1:261-aggressive-phylip/8724_hg18.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=2 rates=propinv;
 Prset applyto=(1) statefreqpr=Fixed(Equal) tratiopr=Beta(1.0,1.0) pinvarpr=Uniform(0.0,1.0);
END;


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

