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

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

Minimum AIC  model: SYMG
Minimum AICc model: SYMG
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
SYMG	61	-1335.4664	2817.2543	0.4035	0.4035
SYMIG	62	-1334.5167	2818.2333	0.2473	0.6507
K2PIG	58	-1341.0830	2819.9622	0.1042	0.7549
GTRG	64	-1332.8860	2820.7850	0.0690	0.8240
K2PG	57	-1342.9483	2820.8871	0.0656	0.8896
GTRIG	65	-1331.9546	2821.8570	0.0404	0.9300
SYMI	61	-1338.5450	2823.4115	0.0186	0.9485
K2PI	57	-1344.2723	2823.5351	0.0175	0.9660
HKYG	60	-1340.2504	2823.9624	0.0141	0.9801
HKYIG	61	-1338.8309	2823.9833	0.0140	0.9940
GTRI	64	-1335.9219	2826.8569	0.0033	0.9974
HKYI	60	-1341.9240	2827.3095	0.0026	1.0000
SYM	60	-1367.9497	2879.3609	0.0000	1.0000
K2P	56	-1374.4617	2881.1259	0.0000	1.0000
GTR	63	-1365.5856	2883.2683	0.0000	1.0000
HKY	59	-1372.3296	2885.2790	0.0000	1.0000
JC69G	56	-1384.7832	2901.7690	0.0000	1.0000
JC69IG	57	-1383.4272	2901.8449	0.0000	1.0000
JC69I	56	-1386.6236	2905.4497	0.0000	1.0000
F81IG	60	-1381.7170	2906.8956	0.0000	1.0000
F81G	59	-1383.1407	2906.9012	0.0000	1.0000
F81I	59	-1384.9561	2910.5321	0.0000	1.0000
JC69	55	-1414.3292	2958.0907	0.0000	1.0000
F81	58	-1412.6503	2963.0968	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
SYMG	61	-1335.4664	2792.9328	0.2879	0.2879
SYMIG	62	-1334.5167	2793.0333	0.2738	0.5617
GTRG	64	-1332.8860	2793.7720	0.1892	0.7510
GTRIG	65	-1331.9546	2793.9092	0.1767	0.9277
K2PIG	58	-1341.0830	2798.1660	0.0210	0.9487
SYMI	61	-1338.5450	2799.0900	0.0133	0.9620
HKYIG	61	-1338.8309	2799.6617	0.0100	0.9719
GTRI	64	-1335.9219	2799.8439	0.0091	0.9810
K2PG	57	-1342.9483	2799.8966	0.0089	0.9899
HKYG	60	-1340.2504	2800.5008	0.0065	0.9964
K2PI	57	-1344.2723	2802.5446	0.0024	0.9988
HKYI	60	-1341.9240	2803.8480	0.0012	1.0000
SYM	60	-1367.9497	2855.8994	0.0000	1.0000
GTR	63	-1365.5856	2857.1712	0.0000	1.0000
K2P	56	-1374.4617	2860.9233	0.0000	1.0000
HKY	59	-1372.3296	2862.6592	0.0000	1.0000
JC69IG	57	-1383.4272	2880.8544	0.0000	1.0000
JC69G	56	-1384.7832	2881.5665	0.0000	1.0000
F81IG	60	-1381.7170	2883.4341	0.0000	1.0000
F81G	59	-1383.1407	2884.2814	0.0000	1.0000
JC69I	56	-1386.6236	2885.2472	0.0000	1.0000
F81I	59	-1384.9561	2887.9123	0.0000	1.0000
JC69	55	-1414.3292	2938.6585	0.0000	1.0000
F81	58	-1412.6503	2941.3007	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
K2PG	57	-1342.9483	3023.4266	0.6179	0.6179
K2PIG	58	-1341.0830	3025.6175	0.2066	0.8245
K2PI	57	-1344.2723	3026.0746	0.1644	0.9889
SYMG	61	-1335.4664	3032.1491	0.0079	0.9968
HKYG	60	-1340.2504	3035.7955	0.0013	0.9981
SYMIG	62	-1334.5167	3036.1712	0.0011	0.9991
SYMI	61	-1338.5450	3038.3063	0.0004	0.9995
HKYIG	61	-1338.8309	3038.8780	0.0003	0.9997
HKYI	60	-1341.9240	3039.1427	0.0002	1.0000
GTRG	64	-1332.8860	3044.7531	0.0000	1.0000
GTRIG	65	-1331.9546	3048.8118	0.0000	1.0000
GTRI	64	-1335.9219	3050.8249	0.0000	1.0000
K2P	56	-1374.4617	3080.5317	0.0000	1.0000
SYM	60	-1367.9497	3091.1941	0.0000	1.0000
HKY	59	-1372.3296	3094.0323	0.0000	1.0000
JC69G	56	-1384.7832	3101.1749	0.0000	1.0000
GTR	63	-1365.5856	3104.2307	0.0000	1.0000
JC69IG	57	-1383.4272	3104.3844	0.0000	1.0000
JC69I	56	-1386.6236	3104.8556	0.0000	1.0000
F81G	59	-1383.1407	3115.6545	0.0000	1.0000
F81IG	60	-1381.7170	3118.7288	0.0000	1.0000
F81I	59	-1384.9561	3119.2854	0.0000	1.0000
JC69	55	-1414.3292	3154.3453	0.0000	1.0000
F81	58	-1412.6503	3168.7522	0.0000	1.0000



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

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

BEGIN MRBAYES;
 Charset 261-aggressive-phylip/13148_hg18.phylip = 1 - 373;
 Partition Dummy = 1:261-aggressive-phylip/13148_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=Fixed(Equal);
END;


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

BEGIN MRBAYES;
 Charset 261-aggressive-phylip/13148_hg18.phylip = 1 - 373;
 Partition Dummy = 1:261-aggressive-phylip/13148_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=Fixed(Equal);
END;


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

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

