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

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

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


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	-1403.2454	2951.6455	0.7062	0.7062
GTRIG	65	-1402.9211	2953.5695	0.2699	0.9760
GTRI	64	-1406.6286	2958.4118	0.0240	1.0000
HKYG	60	-1421.8620	2978.6932	0.0000	1.0000
HKYIG	61	-1420.7555	2979.0109	0.0000	1.0000
HKYI	60	-1422.9235	2980.8163	0.0000	1.0000
GTR	63	-1433.3064	3009.2055	0.0000	1.0000
F81G	59	-1446.4074	3025.2637	0.0000	1.0000
SYMIG	62	-1442.6131	3025.2672	0.0000	1.0000
F81IG	60	-1445.4660	3025.9013	0.0000	1.0000
F81I	59	-1446.9164	3026.2817	0.0000	1.0000
SYMG	61	-1444.6569	3026.8138	0.0000	1.0000
SYMI	61	-1444.7824	3027.0649	0.0000	1.0000
HKY	59	-1452.5449	3037.5387	0.0000	1.0000
K2PIG	58	-1455.0993	3040.1375	0.0000	1.0000
K2PI	57	-1456.9085	3041.2560	0.0000	1.0000
K2PG	57	-1456.9900	3041.4189	0.0000	1.0000
JC69IG	57	-1473.1076	3073.6543	0.0000	1.0000
JC69G	56	-1476.7841	3078.5176	0.0000	1.0000
JC69I	56	-1478.0714	3081.0921	0.0000	1.0000
SYM	60	-1474.3760	3083.7212	0.0000	1.0000
F81	58	-1476.9523	3083.8435	0.0000	1.0000
K2P	56	-1488.9658	3102.8808	0.0000	1.0000
JC69	55	-1508.0610	3138.5916	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRG	64	-1403.2454	2934.4909	0.6482	0.6482
GTRIG	65	-1402.9211	2935.8422	0.3298	0.9780
GTRI	64	-1406.6286	2941.2572	0.0220	1.0000
HKYIG	61	-1420.7555	2963.5109	0.0000	1.0000
HKYG	60	-1421.8620	2963.7239	0.0000	1.0000
HKYI	60	-1422.9235	2965.8470	0.0000	1.0000
GTR	63	-1433.3064	2992.6129	0.0000	1.0000
SYMIG	62	-1442.6131	3009.2262	0.0000	1.0000
F81G	59	-1446.4074	3010.8147	0.0000	1.0000
F81IG	60	-1445.4660	3010.9319	0.0000	1.0000
SYMG	61	-1444.6569	3011.3138	0.0000	1.0000
SYMI	61	-1444.7824	3011.5649	0.0000	1.0000
F81I	59	-1446.9164	3011.8328	0.0000	1.0000
HKY	59	-1452.5449	3023.0898	0.0000	1.0000
K2PIG	58	-1455.0993	3026.1986	0.0000	1.0000
K2PI	57	-1456.9085	3027.8170	0.0000	1.0000
K2PG	57	-1456.9900	3027.9799	0.0000	1.0000
JC69IG	57	-1473.1076	3060.2153	0.0000	1.0000
JC69G	56	-1476.7841	3065.5683	0.0000	1.0000
JC69I	56	-1478.0714	3068.1428	0.0000	1.0000
SYM	60	-1474.3760	3068.7519	0.0000	1.0000
F81	58	-1476.9523	3069.9046	0.0000	1.0000
K2P	56	-1488.9658	3089.9315	0.0000	1.0000
JC69	55	-1508.0610	3126.1219	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
GTRG	64	-1403.2454	3210.3257	0.9119	0.9119
GTRIG	65	-1402.9211	3215.9869	0.0538	0.9657
GTRI	64	-1406.6286	3217.0919	0.0310	0.9967
HKYG	60	-1421.8620	3222.3190	0.0023	0.9989
HKYI	60	-1422.9235	3224.4421	0.0008	0.9997
HKYIG	61	-1420.7555	3226.4160	0.0003	1.0000
GTR	63	-1433.3064	3264.1377	0.0000	1.0000
F81G	59	-1446.4074	3265.0999	0.0000	1.0000
F81I	59	-1446.9164	3266.1179	0.0000	1.0000
F81IG	60	-1445.4660	3269.5270	0.0000	1.0000
K2PI	57	-1456.9085	3273.4824	0.0000	1.0000
K2PG	57	-1456.9900	3273.6453	0.0000	1.0000
SYMG	61	-1444.6569	3274.2188	0.0000	1.0000
SYMI	61	-1444.7824	3274.4699	0.0000	1.0000
K2PIG	58	-1455.0993	3276.1739	0.0000	1.0000
SYMIG	62	-1442.6131	3276.4411	0.0000	1.0000
HKY	59	-1452.5449	3277.3749	0.0000	1.0000
JC69IG	57	-1473.1076	3305.8806	0.0000	1.0000
JC69G	56	-1476.7841	3306.9237	0.0000	1.0000
JC69I	56	-1478.0714	3309.4982	0.0000	1.0000
F81	58	-1476.9523	3319.8798	0.0000	1.0000
SYM	60	-1474.3760	3327.3470	0.0000	1.0000
K2P	56	-1488.9658	3331.2870	0.0000	1.0000
JC69	55	-1508.0610	3363.1674	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/11793_hg18.phylip = 1 - 550;
 Partition Dummy = 1:261-aggressive-phylip/11793_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/11793_hg18.phylip = 1 - 550;
 Partition Dummy = 1:261-aggressive-phylip/11793_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 (GTRG)]

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


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

