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

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

Minimum AIC  model: HKYI
Minimum AICc model: HKYI
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
HKYI	60	-1391.8081	2919.0918	0.3893	0.3893
HKYG	60	-1391.8348	2919.1453	0.3791	0.7684
HKYIG	61	-1391.8083	2921.6421	0.1088	0.8771
GTRI	64	-1388.6610	2923.0618	0.0535	0.9306
GTRG	64	-1388.6739	2923.0877	0.0528	0.9834
GTRIG	65	-1388.6454	2925.6240	0.0149	0.9983
K2PI	57	-1401.9232	2931.7372	0.0007	0.9990
K2PG	57	-1401.9620	2931.8148	0.0007	0.9997
K2PIG	58	-1401.9236	2934.2555	0.0002	0.9999
F81I	59	-1402.3075	2937.5516	0.0000	0.9999
F81G	59	-1402.3407	2937.6181	0.0000	0.9999
SYMI	61	-1400.0884	2938.2021	0.0000	1.0000
SYMG	61	-1400.1130	2938.2514	0.0000	1.0000
F81IG	60	-1402.3078	2940.0913	0.0000	1.0000
SYMIG	62	-1400.0887	2940.7634	0.0000	1.0000
JC69I	56	-1411.6778	2948.7392	0.0000	1.0000
JC69G	56	-1411.7160	2948.8156	0.0000	1.0000
JC69IG	57	-1411.6782	2951.2471	0.0000	1.0000
HKY	59	-1418.6930	2970.3227	0.0000	1.0000
GTR	63	-1415.9011	2974.9596	0.0000	1.0000
K2P	56	-1428.9522	2983.2881	0.0000	1.0000
F81	58	-1429.2536	2988.9157	0.0000	1.0000
SYM	60	-1427.0206	2989.5169	0.0000	1.0000
JC69	55	-1438.6883	3000.2636	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
HKYI	60	-1391.8081	2903.6161	0.2985	0.2985
HKYG	60	-1391.8348	2903.6696	0.2907	0.5892
GTRI	64	-1388.6610	2905.3220	0.1272	0.7164
GTRG	64	-1388.6739	2905.3478	0.1256	0.8420
HKYIG	61	-1391.8083	2905.6167	0.1098	0.9518
GTRIG	65	-1388.6454	2907.2907	0.0475	0.9993
K2PI	57	-1401.9232	2917.8465	0.0002	0.9996
K2PG	57	-1401.9620	2917.9240	0.0002	0.9998
K2PIG	58	-1401.9236	2919.8471	0.0001	0.9999
SYMI	61	-1400.0884	2922.1767	0.0000	0.9999
SYMG	61	-1400.1130	2922.2260	0.0000	0.9999
F81I	59	-1402.3075	2922.6149	0.0000	1.0000
F81G	59	-1402.3407	2922.6814	0.0000	1.0000
SYMIG	62	-1400.0887	2924.1774	0.0000	1.0000
F81IG	60	-1402.3078	2924.6156	0.0000	1.0000
JC69I	56	-1411.6778	2935.3556	0.0000	1.0000
JC69G	56	-1411.7160	2935.4320	0.0000	1.0000
JC69IG	57	-1411.6782	2937.3563	0.0000	1.0000
HKY	59	-1418.6930	2955.3860	0.0000	1.0000
GTR	63	-1415.9011	2957.8022	0.0000	1.0000
K2P	56	-1428.9522	2969.9044	0.0000	1.0000
SYM	60	-1427.0206	2974.0412	0.0000	1.0000
F81	58	-1429.2536	2974.5072	0.0000	1.0000
JC69	55	-1438.6883	2987.3766	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYI	60	-1391.8081	3160.4399	0.3310	0.3310
HKYG	60	-1391.8348	3160.4933	0.3223	0.6533
K2PI	57	-1401.9232	3161.8290	0.1653	0.8186
K2PG	57	-1401.9620	3161.9066	0.1590	0.9776
HKYIG	61	-1391.8083	3166.7208	0.0143	0.9919
K2PIG	58	-1401.9236	3168.1101	0.0071	0.9991
JC69I	56	-1411.6778	3175.0578	0.0002	0.9993
JC69G	56	-1411.7160	3175.1341	0.0002	0.9995
F81I	59	-1402.3075	3175.1583	0.0002	0.9997
F81G	59	-1402.3407	3175.2248	0.0002	0.9999
GTRI	64	-1388.6610	3179.2673	0.0000	0.9999
GTRG	64	-1388.6739	3179.2931	0.0000	1.0000
JC69IG	57	-1411.6782	3181.3389	0.0000	1.0000
F81IG	60	-1402.3078	3181.4394	0.0000	1.0000
SYMI	61	-1400.0884	3183.2809	0.0000	1.0000
SYMG	61	-1400.1130	3183.3302	0.0000	1.0000
GTRIG	65	-1388.6454	3185.5164	0.0000	1.0000
SYMIG	62	-1400.0887	3189.5619	0.0000	1.0000
HKY	59	-1418.6930	3207.9293	0.0000	1.0000
K2P	56	-1428.9522	3209.6066	0.0000	1.0000
F81	58	-1429.2536	3222.7702	0.0000	1.0000
JC69	55	-1438.6883	3222.7984	0.0000	1.0000
GTR	63	-1415.9011	3227.4671	0.0000	1.0000
SYM	60	-1427.0206	3230.8650	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/1203_hg18.phylip = 1 - 534;
 Partition Dummy = 1:261-aggressive-phylip/1203_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/1203_hg18.phylip = 1 - 534;
 Partition Dummy = 1:261-aggressive-phylip/1203_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 (HKYI)]

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


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

