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

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

Minimum AIC  model: GTRG
Minimum AICc model: GTRG
Minimum BIC  model: K2PG


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
GTRG	64	-1159.7312	2470.5094	0.2654	0.2654
HKYG	60	-1165.2601	2470.5750	0.2568	0.5223
HKYI	60	-1165.8130	2471.6808	0.1478	0.6700
GTRI	64	-1160.4962	2472.0394	0.1235	0.7935
GTRIG	65	-1159.6851	2473.2036	0.0690	0.8625
HKYIG	61	-1165.2279	2473.2361	0.0679	0.9304
K2PG	57	-1171.5371	2475.0416	0.0275	0.9580
K2PI	57	-1172.0345	2476.0363	0.0167	0.9747
K2PIG	58	-1171.5072	2477.6630	0.0074	0.9821
HKY	59	-1170.5524	2478.4491	0.0050	0.9871
SYMG	61	-1167.8645	2478.5092	0.0049	0.9920
GTR	63	-1165.4454	2479.1671	0.0035	0.9955
SYMI	61	-1168.4945	2479.7692	0.0026	0.9981
SYMIG	62	-1167.8300	2481.1806	0.0013	0.9994
K2P	56	-1176.7723	2482.8454	0.0006	0.9999
SYM	60	-1173.2243	2486.5034	0.0001	1.0000
F81G	59	-1184.9524	2507.2490	0.0000	1.0000
F81I	59	-1185.5078	2508.3599	0.0000	1.0000
F81IG	60	-1184.9260	2509.9068	0.0000	1.0000
JC69G	56	-1190.3246	2509.9501	0.0000	1.0000
JC69I	56	-1190.8538	2511.0083	0.0000	1.0000
JC69IG	57	-1190.2975	2512.5624	0.0000	1.0000
F81	58	-1190.1082	2514.8649	0.0000	1.0000
JC69	55	-1195.4993	2517.6473	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRG	64	-1159.7312	2447.4623	0.4322	0.4322
GTRI	64	-1160.4962	2448.9923	0.2011	0.6332
GTRIG	65	-1159.6851	2449.3703	0.1665	0.7997
HKYG	60	-1165.2601	2450.5202	0.0937	0.8934
HKYI	60	-1165.8130	2451.6260	0.0539	0.9473
HKYIG	61	-1165.2279	2452.4559	0.0356	0.9829
GTR	63	-1165.4454	2456.8908	0.0039	0.9867
K2PG	57	-1171.5371	2457.0742	0.0035	0.9903
SYMG	61	-1167.8645	2457.7290	0.0025	0.9928
K2PI	57	-1172.0345	2458.0689	0.0022	0.9950
SYMI	61	-1168.4945	2458.9889	0.0014	0.9963
K2PIG	58	-1171.5072	2459.0145	0.0013	0.9977
HKY	59	-1170.5524	2459.1049	0.0013	0.9989
SYMIG	62	-1167.8300	2459.6599	0.0010	0.9999
K2P	56	-1176.7723	2465.5446	0.0001	1.0000
SYM	60	-1173.2243	2466.4486	0.0000	1.0000
F81G	59	-1184.9524	2487.9047	0.0000	1.0000
F81I	59	-1185.5078	2489.0157	0.0000	1.0000
F81IG	60	-1184.9260	2489.8520	0.0000	1.0000
JC69G	56	-1190.3246	2492.6493	0.0000	1.0000
JC69I	56	-1190.8538	2493.7075	0.0000	1.0000
JC69IG	57	-1190.2975	2494.5950	0.0000	1.0000
F81	58	-1190.1082	2496.2164	0.0000	1.0000
JC69	55	-1195.4993	2500.9987	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
K2PG	57	-1171.5371	2688.1772	0.5339	0.5339
K2PI	57	-1172.0345	2689.1719	0.3247	0.8585
K2P	56	-1176.7723	2692.5932	0.0587	0.9172
HKYG	60	-1165.2601	2693.7866	0.0323	0.9495
K2PIG	58	-1171.5072	2694.1720	0.0266	0.9762
HKYI	60	-1165.8130	2694.8923	0.0186	0.9948
HKY	59	-1170.5524	2698.3168	0.0034	0.9981
HKYIG	61	-1165.2279	2699.7767	0.0016	0.9997
SYMG	61	-1167.8645	2705.0498	0.0001	0.9999
SYMI	61	-1168.4945	2706.3097	0.0001	0.9999
GTRG	64	-1159.7312	2706.9464	0.0000	1.0000
GTRI	64	-1160.4962	2708.4765	0.0000	1.0000
SYM	60	-1173.2243	2709.7150	0.0000	1.0000
SYMIG	62	-1167.8300	2711.0352	0.0000	1.0000
GTR	63	-1165.4454	2712.3205	0.0000	1.0000
GTRIG	65	-1159.6851	2712.9088	0.0000	1.0000
JC69G	56	-1190.3246	2719.6979	0.0000	1.0000
JC69I	56	-1190.8538	2720.7561	0.0000	1.0000
JC69	55	-1195.4993	2723.9928	0.0000	1.0000
JC69IG	57	-1190.2975	2725.6981	0.0000	1.0000
F81G	59	-1184.9524	2727.1167	0.0000	1.0000
F81I	59	-1185.5078	2728.2276	0.0000	1.0000
F81	58	-1190.1082	2731.3738	0.0000	1.0000
F81IG	60	-1184.9260	2733.1183	0.0000	1.0000



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

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

BEGIN MRBAYES;
 Charset 183-phylip/chr7_11565.phylip = 1 - 426;
 Partition Dummy = 1:183-phylip/chr7_11565.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=6 rates=gamma;
 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) shapepr=Uniform(0.1,50.0);
END;


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

BEGIN MRBAYES;
 Charset 183-phylip/chr7_11565.phylip = 1 - 426;
 Partition Dummy = 1:183-phylip/chr7_11565.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=6 rates=gamma;
 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) shapepr=Uniform(0.1,50.0);
END;


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

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


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

