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

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

Minimum AIC  model: GTRG
Minimum AICc model: GTRG
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
GTRG	64	-1393.9968	2933.5834	0.4265	0.4265
HKYG	60	-1399.2415	2933.8289	0.3772	0.8037
HKYIG	61	-1399.2506	2936.3919	0.1047	0.9085
GTRIG	65	-1394.2416	2936.6611	0.0915	1.0000
K2PG	57	-1413.7700	2955.3149	0.0000	1.0000
K2PIG	58	-1414.5991	2959.4862	0.0000	1.0000
SYMG	61	-1411.5956	2961.0819	0.0000	1.0000
SYMIG	62	-1412.7074	2965.8611	0.0000	1.0000
GTR	63	-1424.5699	2992.1524	0.0000	1.0000
HKY	59	-1432.5772	2997.9661	0.0000	1.0000
GTRI	64	-1427.0318	2999.6535	0.0000	1.0000
HKYI	60	-1433.2868	3001.9195	0.0000	1.0000
F81G	59	-1438.5088	3009.8294	0.0000	1.0000
F81IG	60	-1440.7032	3016.7523	0.0000	1.0000
K2P	56	-1451.0082	3027.2887	0.0000	1.0000
SYM	60	-1446.7896	3028.9251	0.0000	1.0000
JC69G	56	-1452.2349	3029.7421	0.0000	1.0000
K2PI	57	-1451.8581	3031.4913	0.0000	1.0000
SYMI	61	-1447.8902	3033.6711	0.0000	1.0000
JC69IG	57	-1454.2871	3036.3493	0.0000	1.0000
F81	58	-1472.8861	3076.0602	0.0000	1.0000
F81I	59	-1475.1519	3083.1155	0.0000	1.0000
JC69	55	-1488.4245	3099.6292	0.0000	1.0000
JC69I	56	-1490.5186	3106.3096	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRG	64	-1393.9968	2915.9936	0.5949	0.5949
HKYG	60	-1399.2415	2918.4830	0.1713	0.7662
GTRIG	65	-1394.2416	2918.4832	0.1713	0.9375
HKYIG	61	-1399.2506	2920.5012	0.0625	1.0000
K2PG	57	-1413.7700	2941.5399	0.0000	1.0000
SYMG	61	-1411.5956	2945.1911	0.0000	1.0000
K2PIG	58	-1414.5991	2945.1981	0.0000	1.0000
SYMIG	62	-1412.7074	2949.4148	0.0000	1.0000
GTR	63	-1424.5699	2975.1398	0.0000	1.0000
GTRI	64	-1427.0318	2982.0636	0.0000	1.0000
HKY	59	-1432.5772	2983.1544	0.0000	1.0000
HKYI	60	-1433.2868	2986.5736	0.0000	1.0000
F81G	59	-1438.5088	2995.0177	0.0000	1.0000
F81IG	60	-1440.7032	3001.4064	0.0000	1.0000
SYM	60	-1446.7896	3013.5792	0.0000	1.0000
K2P	56	-1451.0082	3014.0164	0.0000	1.0000
JC69G	56	-1452.2349	3016.4697	0.0000	1.0000
K2PI	57	-1451.8581	3017.7163	0.0000	1.0000
SYMI	61	-1447.8902	3017.7803	0.0000	1.0000
JC69IG	57	-1454.2871	3022.5743	0.0000	1.0000
F81	58	-1472.8861	3061.7721	0.0000	1.0000
F81I	59	-1475.1519	3068.3038	0.0000	1.0000
JC69	55	-1488.4245	3086.8491	0.0000	1.0000
JC69I	56	-1490.5186	3093.0372	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	60	-1399.2415	3175.7545	0.9527	0.9527
HKYIG	61	-1399.2506	3182.0605	0.0407	0.9934
K2PG	57	-1413.7700	3185.9479	0.0058	0.9992
GTRG	64	-1393.9968	3190.4165	0.0006	0.9999
K2PIG	58	-1414.5991	3193.8939	0.0001	1.0000
GTRIG	65	-1394.2416	3197.1940	0.0000	1.0000
SYMG	61	-1411.5956	3206.7505	0.0000	1.0000
SYMIG	62	-1412.7074	3215.2620	0.0000	1.0000
HKY	59	-1432.5772	3236.1380	0.0000	1.0000
HKYI	60	-1433.2868	3243.8451	0.0000	1.0000
GTR	63	-1424.5699	3245.2748	0.0000	1.0000
F81G	59	-1438.5088	3248.0013	0.0000	1.0000
K2P	56	-1451.0082	3254.1365	0.0000	1.0000
GTRI	64	-1427.0318	3256.4866	0.0000	1.0000
JC69G	56	-1452.2349	3256.5898	0.0000	1.0000
F81IG	60	-1440.7032	3258.6779	0.0000	1.0000
K2PI	57	-1451.8581	3262.1242	0.0000	1.0000
JC69IG	57	-1454.2871	3266.9822	0.0000	1.0000
SYM	60	-1446.7896	3270.8507	0.0000	1.0000
SYMI	61	-1447.8902	3279.3397	0.0000	1.0000
F81	58	-1472.8861	3310.4679	0.0000	1.0000
F81I	59	-1475.1519	3321.2875	0.0000	1.0000
JC69	55	-1488.4245	3322.6813	0.0000	1.0000
JC69I	56	-1490.5186	3333.1573	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 261-aggressive-phylip/12228_hg18.phylip = 1 - 538;
 Partition Dummy = 1:261-aggressive-phylip/12228_hg18.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 261-aggressive-phylip/12228_hg18.phylip = 1 - 538;
 Partition Dummy = 1:261-aggressive-phylip/12228_hg18.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 (HKYG)]

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

