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

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

Minimum AIC  model: HKYIG
Minimum AICc model: HKYIG
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
HKYIG	61	-1138.0737	2421.4930	0.3834	0.3834
HKYG	60	-1139.6235	2421.7700	0.3338	0.7171
HKYI	60	-1140.2062	2422.9355	0.1864	0.9035
GTRG	64	-1135.7041	2425.3272	0.0564	0.9599
GTRIG	65	-1134.9887	2426.7899	0.0271	0.9870
GTRI	64	-1137.4345	2428.7881	0.0100	0.9970
K2PIG	58	-1147.2668	2431.4632	0.0026	0.9996
K2PG	57	-1150.6638	2435.4862	0.0004	0.9999
K2PI	57	-1153.0315	2440.2216	0.0000	1.0000
SYMG	61	-1148.4292	2442.2042	0.0000	1.0000
SYMI	61	-1149.4732	2444.2921	0.0000	1.0000
SYMIG	62	-1148.2092	2444.6042	0.0000	1.0000
HKY	59	-1168.5008	2476.7193	0.0000	1.0000
GTR	63	-1165.5876	2482.2187	0.0000	1.0000
K2P	56	-1178.9714	2489.3471	0.0000	1.0000
F81G	59	-1177.6103	2494.9384	0.0000	1.0000
F81IG	60	-1176.2257	2494.9744	0.0000	1.0000
F81I	59	-1178.1290	2495.9757	0.0000	1.0000
SYM	60	-1177.9194	2498.3620	0.0000	1.0000
JC69IG	57	-1185.3620	2504.8825	0.0000	1.0000
JC69G	56	-1187.0387	2505.4817	0.0000	1.0000
JC69I	56	-1189.5989	2510.6020	0.0000	1.0000
F81	58	-1204.9708	2546.8712	0.0000	1.0000
JC69	55	-1213.7200	2556.1067	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
HKYIG	61	-1138.0737	2398.1474	0.3414	0.3414
HKYG	60	-1139.6235	2399.2469	0.1970	0.5385
GTRG	64	-1135.7041	2399.4082	0.1818	0.7202
GTRIG	65	-1134.9887	2399.9774	0.1368	0.8570
HKYI	60	-1140.2062	2400.4124	0.1100	0.9670
GTRI	64	-1137.4345	2402.8691	0.0322	0.9992
K2PIG	58	-1147.2668	2410.5336	0.0007	0.9999
K2PG	57	-1150.6638	2415.3277	0.0001	1.0000
SYMG	61	-1148.4292	2418.8585	0.0000	1.0000
K2PI	57	-1153.0315	2420.0630	0.0000	1.0000
SYMIG	62	-1148.2092	2420.4185	0.0000	1.0000
SYMI	61	-1149.4732	2420.9464	0.0000	1.0000
HKY	59	-1168.5008	2455.0015	0.0000	1.0000
GTR	63	-1165.5876	2457.1752	0.0000	1.0000
K2P	56	-1178.9714	2469.9429	0.0000	1.0000
F81IG	60	-1176.2257	2472.4513	0.0000	1.0000
F81G	59	-1177.6103	2473.2206	0.0000	1.0000
F81I	59	-1178.1290	2474.2579	0.0000	1.0000
SYM	60	-1177.9194	2475.8389	0.0000	1.0000
JC69IG	57	-1185.3620	2484.7240	0.0000	1.0000
JC69G	56	-1187.0387	2486.0774	0.0000	1.0000
JC69I	56	-1189.5989	2491.1977	0.0000	1.0000
F81	58	-1204.9708	2525.9415	0.0000	1.0000
JC69	55	-1213.7200	2537.4401	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	60	-1139.6235	2636.5972	0.4725	0.4725
HKYI	60	-1140.2062	2637.7626	0.2638	0.7363
HKYIG	61	-1138.0737	2639.4534	0.1133	0.8495
K2PIG	58	-1147.2668	2639.9721	0.0874	0.9369
K2PG	57	-1150.6638	2640.8104	0.0575	0.9944
K2PI	57	-1153.0315	2645.5458	0.0054	0.9998
GTRG	64	-1135.7041	2652.5818	0.0002	0.9999
GTRI	64	-1137.4345	2656.0427	0.0000	1.0000
GTRIG	65	-1134.9887	2657.1068	0.0000	1.0000
SYMG	61	-1148.4292	2660.1646	0.0000	1.0000
SYMI	61	-1149.4732	2662.2525	0.0000	1.0000
SYMIG	62	-1148.2092	2665.6804	0.0000	1.0000
HKY	59	-1168.5008	2688.3959	0.0000	1.0000
K2P	56	-1178.9714	2691.4698	0.0000	1.0000
GTR	63	-1165.5876	2706.3930	0.0000	1.0000
F81G	59	-1177.6103	2706.6150	0.0000	1.0000
JC69G	56	-1187.0387	2707.6043	0.0000	1.0000
F81I	59	-1178.1290	2707.6523	0.0000	1.0000
F81IG	60	-1176.2257	2709.8016	0.0000	1.0000
JC69IG	57	-1185.3620	2710.2067	0.0000	1.0000
JC69I	56	-1189.5989	2712.7246	0.0000	1.0000
SYM	60	-1177.9194	2713.1891	0.0000	1.0000
JC69	55	-1213.7200	2755.0111	0.0000	1.0000
F81	58	-1204.9708	2755.3801	0.0000	1.0000



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

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

BEGIN MRBAYES;
 Charset 183-phylip/chr12_5895.phylip = 1 - 386;
 Partition Dummy = 1:183-phylip/chr12_5895.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=2 rates=invgamma;
 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) pinvarpr=Uniform(0.0,1.0);
END;


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

BEGIN MRBAYES;
 Charset 183-phylip/chr12_5895.phylip = 1 - 386;
 Partition Dummy = 1:183-phylip/chr12_5895.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=2 rates=invgamma;
 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) pinvarpr=Uniform(0.0,1.0);
END;


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

BEGIN MRBAYES;
 Charset 183-phylip/chr12_5895.phylip = 1 - 386;
 Partition Dummy = 1:183-phylip/chr12_5895.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

