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

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

Minimum AIC  model: K2PIG
Minimum AICc model: K2PG
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
K2PG	57	-1268.5163	2669.8167	0.3711	0.3711
K2PIG	58	-1267.4449	2670.3883	0.2788	0.6499
HKYG	60	-1265.8423	2672.6589	0.0896	0.7395
SYMG	61	-1264.6487	2673.0330	0.0743	0.8138
HKYIG	61	-1264.8141	2673.3638	0.0630	0.8768
SYMIG	62	-1263.5540	2673.6209	0.0554	0.9322
K2PI	57	-1270.6734	2674.1309	0.0429	0.9751
HKYI	60	-1267.9226	2676.8194	0.0112	0.9863
GTRG	64	-1263.0085	2678.1329	0.0058	0.9921
GTRIG	65	-1261.9787	2678.8992	0.0040	0.9960
SYMI	61	-1267.7051	2679.1459	0.0035	0.9995
GTRI	64	-1265.5290	2683.1740	0.0005	1.0000
K2P	56	-1293.8652	2717.8153	0.0000	1.0000
SYM	60	-1288.9166	2718.8074	0.0000	1.0000
JC69G	56	-1294.8664	2719.8179	0.0000	1.0000
HKY	59	-1290.8374	2719.9035	0.0000	1.0000
F81G	59	-1291.3401	2720.9088	0.0000	1.0000
JC69IG	57	-1294.1670	2721.1181	0.0000	1.0000
F81IG	60	-1290.5921	2722.1585	0.0000	1.0000
GTR	63	-1287.2066	2723.7196	0.0000	1.0000
JC69I	56	-1296.9289	2723.9428	0.0000	1.0000
F81I	59	-1293.5022	2725.2330	0.0000	1.0000
JC69	55	-1316.6977	2760.7966	0.0000	1.0000
F81	58	-1313.6902	2762.8790	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
K2PIG	58	-1267.4449	2650.8897	0.1761	0.1761
K2PG	57	-1268.5163	2651.0326	0.1640	0.3401
SYMIG	62	-1263.5540	2651.1080	0.1579	0.4980
SYMG	61	-1264.6487	2651.2973	0.1436	0.6416
HKYIG	61	-1264.8141	2651.6282	0.1217	0.7633
HKYG	60	-1265.8423	2651.6847	0.1183	0.8817
GTRIG	65	-1261.9787	2653.9573	0.0380	0.9197
GTRG	64	-1263.0085	2654.0169	0.0369	0.9565
K2PI	57	-1270.6734	2655.3468	0.0190	0.9755
HKYI	60	-1267.9226	2655.8451	0.0148	0.9903
SYMI	61	-1267.7051	2657.4103	0.0068	0.9970
GTRI	64	-1265.5290	2659.0581	0.0030	1.0000
SYM	60	-1288.9166	2697.8332	0.0000	1.0000
HKY	59	-1290.8374	2699.6749	0.0000	1.0000
K2P	56	-1293.8652	2699.7303	0.0000	1.0000
GTR	63	-1287.2066	2700.4132	0.0000	1.0000
F81G	59	-1291.3401	2700.6803	0.0000	1.0000
F81IG	60	-1290.5921	2701.1842	0.0000	1.0000
JC69G	56	-1294.8664	2701.7329	0.0000	1.0000
JC69IG	57	-1294.1670	2702.3340	0.0000	1.0000
F81I	59	-1293.5022	2705.0044	0.0000	1.0000
JC69I	56	-1296.9289	2705.8578	0.0000	1.0000
F81	58	-1313.6902	2743.3804	0.0000	1.0000
JC69	55	-1316.6977	2743.3954	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
K2PG	57	-1268.5163	2879.9535	0.7922	0.7922
K2PIG	58	-1267.4449	2883.8269	0.1142	0.9064
K2PI	57	-1270.6734	2884.2678	0.0916	0.9980
HKYG	60	-1265.8423	2892.6541	0.0014	0.9994
SYMG	61	-1264.6487	2896.2829	0.0002	0.9996
HKYIG	61	-1264.8141	2896.6138	0.0002	0.9998
HKYI	60	-1267.9226	2896.8146	0.0002	1.0000
SYMIG	62	-1263.5540	2900.1097	0.0000	1.0000
SYMI	61	-1267.7051	2902.3959	0.0000	1.0000
GTRG	64	-1263.0085	2911.0510	0.0000	1.0000
GTRIG	65	-1261.9787	2915.0075	0.0000	1.0000
GTRI	64	-1265.5290	2916.0921	0.0000	1.0000
K2P	56	-1293.8652	2924.6351	0.0000	1.0000
JC69G	56	-1294.8664	2926.6377	0.0000	1.0000
JC69I	56	-1296.9289	2930.7626	0.0000	1.0000
JC69IG	57	-1294.1670	2931.2550	0.0000	1.0000
HKY	59	-1290.8374	2936.6282	0.0000	1.0000
F81G	59	-1291.3401	2937.6335	0.0000	1.0000
SYM	60	-1288.9166	2938.8026	0.0000	1.0000
F81I	59	-1293.5022	2941.9577	0.0000	1.0000
F81IG	60	-1290.5921	2942.1537	0.0000	1.0000
GTR	63	-1287.2066	2953.4311	0.0000	1.0000
JC69	55	-1316.6977	2964.2841	0.0000	1.0000
F81	58	-1313.6902	2976.3175	0.0000	1.0000



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

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

BEGIN MRBAYES;
 Charset 261-aggressive-phylip/10965_hg18.phylip = 1 - 410;
 Partition Dummy = 1:261-aggressive-phylip/10965_hg18.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=2 rates=invgamma;
 Prset applyto=(1) statefreqpr=Fixed(Equal) tratiopr=Beta(1.0,1.0) shapepr=Uniform(0.1,50.0) pinvarpr=Uniform(0.0,1.0);
END;


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

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


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

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

