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

Input data from file "261-aggressive-phylip/8242_hg18.phylip" (29 taxa, 400 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	-1334.1711	2809.9351	0.6105	0.6105
HKYIG	61	-1333.5364	2811.4516	0.2860	0.8966
HKYI	60	-1336.1360	2813.8649	0.0856	0.9821
GTRG	64	-1332.7180	2818.2718	0.0094	0.9916
GTRIG	65	-1331.5550	2818.7987	0.0073	0.9989
GTRI	64	-1334.8295	2822.4947	0.0011	1.0000
K2PG	57	-1355.7239	2844.7812	0.0000	1.0000
K2PIG	58	-1355.2310	2846.5324	0.0000	1.0000
K2PI	57	-1356.7895	2846.9123	0.0000	1.0000
SYMG	61	-1354.3865	2853.1516	0.0000	1.0000
SYMIG	62	-1353.9167	2855.0144	0.0000	1.0000
F81G	59	-1358.2121	2855.2477	0.0000	1.0000
SYMI	61	-1355.5988	2855.5764	0.0000	1.0000
F81IG	60	-1357.6600	2856.9129	0.0000	1.0000
F81I	59	-1359.4484	2857.7203	0.0000	1.0000
HKY	59	-1374.7228	2888.2692	0.0000	1.0000
JC69G	56	-1380.4717	2891.5556	0.0000	1.0000
GTR	63	-1371.6318	2893.2635	0.0000	1.0000
JC69IG	57	-1380.0559	2893.4451	0.0000	1.0000
JC69I	56	-1381.4454	2893.5031	0.0000	1.0000
K2P	56	-1394.5849	2919.7820	0.0000	1.0000
SYM	60	-1390.0338	2921.6606	0.0000	1.0000
F81	58	-1399.6750	2935.4204	0.0000	1.0000
JC69	55	-1418.8890	2965.6849	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
HKYG	60	-1334.1711	2788.3421	0.4965	0.4965
HKYIG	61	-1333.5364	2789.0729	0.3445	0.8410
HKYI	60	-1336.1360	2792.2720	0.0696	0.9106
GTRIG	65	-1331.5550	2793.1101	0.0458	0.9564
GTRG	64	-1332.7180	2793.4360	0.0389	0.9953
GTRI	64	-1334.8295	2797.6589	0.0047	1.0000
K2PG	57	-1355.7239	2825.4478	0.0000	1.0000
K2PIG	58	-1355.2310	2826.4620	0.0000	1.0000
K2PI	57	-1356.7895	2827.5790	0.0000	1.0000
SYMG	61	-1354.3865	2830.7729	0.0000	1.0000
SYMIG	62	-1353.9167	2831.8334	0.0000	1.0000
SYMI	61	-1355.5988	2833.1977	0.0000	1.0000
F81G	59	-1358.2121	2834.4242	0.0000	1.0000
F81IG	60	-1357.6600	2835.3200	0.0000	1.0000
F81I	59	-1359.4484	2836.8968	0.0000	1.0000
HKY	59	-1374.7228	2867.4457	0.0000	1.0000
GTR	63	-1371.6318	2869.2635	0.0000	1.0000
JC69G	56	-1380.4717	2872.9433	0.0000	1.0000
JC69IG	57	-1380.0559	2874.1118	0.0000	1.0000
JC69I	56	-1381.4454	2874.8909	0.0000	1.0000
SYM	60	-1390.0338	2900.0677	0.0000	1.0000
K2P	56	-1394.5849	2901.1697	0.0000	1.0000
F81	58	-1399.6750	2915.3500	0.0000	1.0000
JC69	55	-1418.8890	2947.7779	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	60	-1334.1711	3027.8300	0.8100	0.8100
HKYI	60	-1336.1360	3031.7598	0.1135	0.9236
HKYIG	61	-1333.5364	3032.5522	0.0764	1.0000
GTRG	64	-1332.7180	3048.8897	0.0000	1.0000
GTRIG	65	-1331.5550	3052.5553	0.0000	1.0000
K2PG	57	-1355.7239	3052.9613	0.0000	1.0000
GTRI	64	-1334.8295	3053.1127	0.0000	1.0000
K2PI	57	-1356.7895	3055.0925	0.0000	1.0000
K2PIG	58	-1355.2310	3057.9670	0.0000	1.0000
F81G	59	-1358.2121	3069.9206	0.0000	1.0000
F81I	59	-1359.4484	3072.3932	0.0000	1.0000
SYMG	61	-1354.3865	3074.2523	0.0000	1.0000
F81IG	60	-1357.6600	3074.8079	0.0000	1.0000
SYMI	61	-1355.5988	3076.6770	0.0000	1.0000
SYMIG	62	-1353.9167	3079.3042	0.0000	1.0000
JC69G	56	-1380.4717	3096.4654	0.0000	1.0000
JC69I	56	-1381.4454	3098.4129	0.0000	1.0000
JC69IG	57	-1380.0559	3101.6252	0.0000	1.0000
HKY	59	-1374.7228	3102.9421	0.0000	1.0000
GTR	63	-1371.6318	3120.7258	0.0000	1.0000
K2P	56	-1394.5849	3124.6917	0.0000	1.0000
SYM	60	-1390.0338	3139.5556	0.0000	1.0000
F81	58	-1399.6750	3146.8549	0.0000	1.0000
JC69	55	-1418.8890	3167.3085	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 261-aggressive-phylip/8242_hg18.phylip = 1 - 400;
 Partition Dummy = 1:261-aggressive-phylip/8242_hg18.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 261-aggressive-phylip/8242_hg18.phylip = 1 - 400;
 Partition Dummy = 1:261-aggressive-phylip/8242_hg18.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 261-aggressive-phylip/8242_hg18.phylip = 1 - 400;
 Partition Dummy = 1:261-aggressive-phylip/8242_hg18.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

