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

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

Minimum AIC  model: GTRI
Minimum AICc model: GTRI
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
GTRI	64	-1153.0645	2459.3412	0.5005	0.5005
GTRG	64	-1153.2993	2459.8107	0.3958	0.8963
GTRIG	65	-1153.3152	2462.7094	0.0929	0.9892
HKYI	60	-1163.8125	2469.5412	0.0031	0.9923
HKYG	60	-1163.9210	2469.7582	0.0027	0.9950
SYMI	61	-1162.9960	2470.7067	0.0017	0.9967
SYMG	61	-1163.1488	2471.0123	0.0015	0.9982
HKYIG	61	-1163.9315	2472.5777	0.0007	0.9988
SYMIG	62	-1163.1581	2473.8463	0.0004	0.9992
K2PI	57	-1170.1491	2473.9183	0.0003	0.9995
K2PG	57	-1170.2052	2474.0307	0.0003	0.9999
K2PIG	58	-1170.2235	2476.8161	0.0001	0.9999
GTR	63	-1163.6637	2477.6900	0.0001	1.0000
HKY	59	-1173.3297	2485.7937	0.0000	1.0000
SYM	60	-1172.3629	2486.6421	0.0000	1.0000
K2P	56	-1179.4589	2489.8053	0.0000	1.0000
F81I	59	-1183.5843	2506.3029	0.0000	1.0000
F81G	59	-1183.7141	2506.5624	0.0000	1.0000
JC69I	56	-1188.5403	2507.9682	0.0000	1.0000
JC69G	56	-1188.6387	2508.1650	0.0000	1.0000
F81IG	60	-1183.7158	2509.3478	0.0000	1.0000
JC69IG	57	-1188.6397	2510.8995	0.0000	1.0000
F81	58	-1192.8197	2522.0084	0.0000	1.0000
JC69	55	-1197.7131	2523.5972	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRI	64	-1153.0645	2434.1290	0.4803	0.4803
GTRG	64	-1153.2993	2434.5986	0.3798	0.8601
GTRIG	65	-1153.3152	2436.6304	0.1375	0.9976
HKYI	60	-1163.8125	2447.6250	0.0006	0.9982
HKYG	60	-1163.9210	2447.8421	0.0005	0.9987
SYMI	61	-1162.9960	2447.9920	0.0005	0.9992
SYMG	61	-1163.1488	2448.2976	0.0004	0.9996
HKYIG	61	-1163.9315	2449.8630	0.0002	0.9998
SYMIG	62	-1163.1581	2450.3162	0.0001	0.9999
GTR	63	-1163.6637	2453.3274	0.0000	1.0000
K2PI	57	-1170.1491	2454.2981	0.0000	1.0000
K2PG	57	-1170.2052	2454.4105	0.0000	1.0000
K2PIG	58	-1170.2235	2456.4471	0.0000	1.0000
HKY	59	-1173.3297	2464.6594	0.0000	1.0000
SYM	60	-1172.3629	2464.7259	0.0000	1.0000
K2P	56	-1179.4589	2470.9177	0.0000	1.0000
F81I	59	-1183.5843	2485.1685	0.0000	1.0000
F81G	59	-1183.7141	2485.4281	0.0000	1.0000
F81IG	60	-1183.7158	2487.4316	0.0000	1.0000
JC69I	56	-1188.5403	2489.0806	0.0000	1.0000
JC69G	56	-1188.6387	2489.2775	0.0000	1.0000
JC69IG	57	-1188.6397	2491.2793	0.0000	1.0000
F81	58	-1192.8197	2501.6393	0.0000	1.0000
JC69	55	-1197.7131	2505.4261	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
K2PI	57	-1170.1491	2681.0946	0.4569	0.4569
K2PG	57	-1170.2052	2681.2070	0.4320	0.8889
HKYI	60	-1163.8125	2686.3582	0.0329	0.9218
HKYG	60	-1163.9210	2686.5752	0.0295	0.9513
K2PIG	58	-1170.2235	2687.2224	0.0213	0.9726
GTRI	64	-1153.0645	2688.7777	0.0098	0.9824
GTRG	64	-1153.2993	2689.2472	0.0078	0.9902
SYMI	61	-1162.9960	2690.7040	0.0037	0.9939
SYMG	61	-1163.1488	2691.0097	0.0032	0.9971
HKYIG	61	-1163.9315	2692.5750	0.0015	0.9986
K2P	56	-1179.4589	2693.7353	0.0008	0.9994
GTRIG	65	-1153.3152	2695.2580	0.0004	0.9998
SYMIG	62	-1163.1581	2697.0071	0.0002	0.9999
HKY	59	-1173.3297	2699.4137	0.0000	1.0000
SYM	60	-1172.3629	2703.4590	0.0000	1.0000
GTR	63	-1163.6637	2703.9972	0.0000	1.0000
JC69I	56	-1188.5403	2711.8982	0.0000	1.0000
JC69G	56	-1188.6387	2712.0951	0.0000	1.0000
JC69IG	57	-1188.6397	2718.0758	0.0000	1.0000
F81I	59	-1183.5843	2719.9228	0.0000	1.0000
F81G	59	-1183.7141	2720.1824	0.0000	1.0000
JC69	55	-1197.7131	2724.2648	0.0000	1.0000
F81IG	60	-1183.7158	2726.1648	0.0000	1.0000
F81	58	-1192.8197	2732.4147	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 261-aggressive-phylip/3830_hg18.phylip = 1 - 395;
 Partition Dummy = 1:261-aggressive-phylip/3830_hg18.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 261-aggressive-phylip/3830_hg18.phylip = 1 - 395;
 Partition Dummy = 1:261-aggressive-phylip/3830_hg18.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 (K2PI)]

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

