
Output from mraic.pl version 1.4.4 by Johan Nylander
-------------------------------------------------------------
00:52 02/17/2011

Input data from file "183-phylip/chr3_6248.phylip" (29 taxa, 356 characters)

Minimum AIC  model: HKYG
Minimum AICc model: HKYG
Minimum BIC  model: HKYG


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
HKYG	60	-1051.9550	2248.7236	0.7628	0.7628
HKYI	60	-1053.6740	2252.1616	0.1367	0.8995
HKYIG	61	-1053.1362	2254.0003	0.0545	0.9540
GTRG	64	-1049.5113	2255.6137	0.0243	0.9783
GTRI	64	-1049.9428	2256.4766	0.0158	0.9941
GTRIG	65	-1049.4402	2258.4667	0.0058	1.0000
HKY	59	-1064.9361	2271.7911	0.0000	1.0000
K2PG	57	-1070.2383	2276.6646	0.0000	1.0000
K2PI	57	-1070.4756	2277.1391	0.0000	1.0000
GTR	63	-1062.5188	2278.6541	0.0000	1.0000
K2PIG	58	-1070.1856	2279.4149	0.0000	1.0000
SYMG	61	-1068.0486	2283.8251	0.0000	1.0000
SYMI	61	-1068.3475	2284.4228	0.0000	1.0000
SYMIG	62	-1067.9769	2286.6160	0.0000	1.0000
K2P	56	-1082.3957	2298.1426	0.0000	1.0000
SYM	60	-1080.4824	2305.7783	0.0000	1.0000
F81G	59	-1089.9704	2321.8596	0.0000	1.0000
F81I	59	-1092.1559	2326.2307	0.0000	1.0000
F81IG	60	-1091.7738	2328.3612	0.0000	1.0000
JC69G	56	-1103.2541	2339.8594	0.0000	1.0000
JC69I	56	-1103.4318	2340.2148	0.0000	1.0000
JC69IG	57	-1103.3666	2342.9211	0.0000	1.0000
F81	58	-1102.1631	2343.3699	0.0000	1.0000
JC69	55	-1115.3725	2361.2783	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
HKYG	60	-1051.9550	2223.9101	0.5802	0.5802
GTRG	64	-1049.5113	2227.0226	0.1224	0.7026
HKYI	60	-1053.6740	2227.3480	0.1040	0.8066
GTRI	64	-1049.9428	2227.8855	0.0795	0.8861
HKYIG	61	-1053.1362	2228.2724	0.0655	0.9517
GTRIG	65	-1049.4402	2228.8805	0.0483	1.0000
HKY	59	-1064.9361	2247.8722	0.0000	1.0000
GTR	63	-1062.5188	2251.0377	0.0000	1.0000
K2PG	57	-1070.2383	2254.4766	0.0000	1.0000
K2PI	57	-1070.4756	2254.9512	0.0000	1.0000
K2PIG	58	-1070.1856	2256.3711	0.0000	1.0000
SYMG	61	-1068.0486	2258.0972	0.0000	1.0000
SYMI	61	-1068.3475	2258.6949	0.0000	1.0000
SYMIG	62	-1067.9769	2259.9539	0.0000	1.0000
K2P	56	-1082.3957	2276.7914	0.0000	1.0000
SYM	60	-1080.4824	2280.9647	0.0000	1.0000
F81G	59	-1089.9704	2297.9407	0.0000	1.0000
F81I	59	-1092.1559	2302.3118	0.0000	1.0000
F81IG	60	-1091.7738	2303.5476	0.0000	1.0000
JC69G	56	-1103.2541	2318.5082	0.0000	1.0000
JC69I	56	-1103.4318	2318.8636	0.0000	1.0000
F81	58	-1102.1631	2320.3261	0.0000	1.0000
JC69IG	57	-1103.3666	2320.7332	0.0000	1.0000
JC69	55	-1115.3725	2340.7450	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	60	-1051.9550	2456.4059	0.8362	0.8362
HKYI	60	-1053.6740	2459.8439	0.1499	0.9861
HKYIG	61	-1053.1362	2464.6432	0.0136	0.9997
GTRG	64	-1049.5113	2475.0182	0.0001	0.9998
K2PG	57	-1070.2383	2475.3477	0.0001	0.9999
K2PI	57	-1070.4756	2475.8222	0.0001	0.9999
GTRI	64	-1049.9428	2475.8811	0.0000	1.0000
HKY	59	-1064.9361	2476.4931	0.0000	1.0000
GTRIG	65	-1049.4402	2480.7510	0.0000	1.0000
K2PIG	58	-1070.1856	2481.1171	0.0000	1.0000
K2P	56	-1082.3957	2493.7875	0.0000	1.0000
SYMG	61	-1068.0486	2494.4680	0.0000	1.0000
SYMI	61	-1068.3475	2495.0657	0.0000	1.0000
GTR	63	-1062.5188	2495.1583	0.0000	1.0000
SYMIG	62	-1067.9769	2500.1996	0.0000	1.0000
SYM	60	-1080.4824	2513.4606	0.0000	1.0000
F81G	59	-1089.9704	2526.5616	0.0000	1.0000
F81I	59	-1092.1559	2530.9327	0.0000	1.0000
JC69G	56	-1103.2541	2535.5043	0.0000	1.0000
JC69I	56	-1103.4318	2535.8598	0.0000	1.0000
F81IG	60	-1091.7738	2536.0435	0.0000	1.0000
JC69IG	57	-1103.3666	2541.6043	0.0000	1.0000
F81	58	-1102.1631	2545.0721	0.0000	1.0000
JC69	55	-1115.3725	2553.8662	0.0000	1.0000



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

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

BEGIN MRBAYES;
 Charset 183-phylip/chr3_6248.phylip = 1 - 356;
 Partition Dummy = 1:183-phylip/chr3_6248.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=2 rates=gamma;
 Prset applyto=(1) statefreqpr=Dirichlet(1.0,1.0,1.0,1.0) tratiopr=Beta(1.0,1.0) shapepr=Uniform(0.1,50.0);
END;


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

BEGIN MRBAYES;
 Charset 183-phylip/chr3_6248.phylip = 1 - 356;
 Partition Dummy = 1:183-phylip/chr3_6248.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=2 rates=gamma;
 Prset applyto=(1) statefreqpr=Dirichlet(1.0,1.0,1.0,1.0) tratiopr=Beta(1.0,1.0) shapepr=Uniform(0.1,50.0);
END;


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

BEGIN MRBAYES;
 Charset 183-phylip/chr3_6248.phylip = 1 - 356;
 Partition Dummy = 1:183-phylip/chr3_6248.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=2 rates=gamma;
 Prset applyto=(1) statefreqpr=Dirichlet(1.0,1.0,1.0,1.0) tratiopr=Beta(1.0,1.0) shapepr=Uniform(0.1,50.0);
END;


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

