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

Input data from file "261-aggressive-phylip/10976_hg18.phylip" (29 taxa, 584 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	-2900.6621	5945.3551	0.8458	0.8458
GTRIG	65	-2901.8724	5950.3085	0.0711	0.9168
HKYG	60	-2908.2946	5950.5854	0.0619	0.9787
HKYIG	61	-2908.1141	5952.7187	0.0213	1.0000
SYMG	61	-2921.6085	5979.7075	0.0000	1.0000
K2PIG	58	-2927.5461	5984.1285	0.0000	1.0000
SYMIG	62	-2923.4205	5985.8353	0.0000	1.0000
K2PG	57	-2930.3403	5987.2508	0.0000	1.0000
HKYI	60	-2930.8131	5995.6224	0.0000	1.0000
GTRI	64	-2926.8585	5997.7477	0.0000	1.0000
SYMI	61	-2945.2823	6027.0550	0.0000	1.0000
K2PI	57	-2950.2752	6027.1208	0.0000	1.0000
GTR	63	-2989.1224	6119.7525	0.0000	1.0000
HKY	59	-2995.3988	6122.3091	0.0000	1.0000
F81G	59	-2996.3857	6124.2829	0.0000	1.0000
F81IG	60	-2999.2808	6132.5578	0.0000	1.0000
SYM	60	-3003.4301	6140.8565	0.0000	1.0000
JC69G	56	-3012.9377	6149.9893	0.0000	1.0000
K2P	56	-3014.7264	6153.5667	0.0000	1.0000
JC69IG	57	-3013.6115	6153.7934	0.0000	1.0000
F81I	59	-3018.6169	6168.7452	0.0000	1.0000
JC69I	56	-3035.5327	6195.1792	0.0000	1.0000
F81	58	-3075.5308	6280.0979	0.0000	1.0000
JC69	55	-3089.8652	6301.3972	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRG	64	-2900.6621	5929.3242	0.8713	0.8713
GTRIG	65	-2901.8724	5933.7448	0.0956	0.9668
HKYG	60	-2908.2946	5936.5893	0.0230	0.9898
HKYIG	61	-2908.1141	5938.2283	0.0102	1.0000
SYMG	61	-2921.6085	5965.2171	0.0000	1.0000
SYMIG	62	-2923.4205	5970.8411	0.0000	1.0000
K2PIG	58	-2927.5461	5971.0923	0.0000	1.0000
K2PG	57	-2930.3403	5974.6805	0.0000	1.0000
HKYI	60	-2930.8131	5981.6263	0.0000	1.0000
GTRI	64	-2926.8585	5981.7169	0.0000	1.0000
SYMI	61	-2945.2823	6012.5646	0.0000	1.0000
K2PI	57	-2950.2752	6014.5504	0.0000	1.0000
GTR	63	-2989.1224	6104.2448	0.0000	1.0000
HKY	59	-2995.3988	6108.7976	0.0000	1.0000
F81G	59	-2996.3857	6110.7714	0.0000	1.0000
F81IG	60	-2999.2808	6118.5616	0.0000	1.0000
SYM	60	-3003.4301	6126.8603	0.0000	1.0000
JC69G	56	-3012.9377	6137.8754	0.0000	1.0000
JC69IG	57	-3013.6115	6141.2230	0.0000	1.0000
K2P	56	-3014.7264	6141.4529	0.0000	1.0000
F81I	59	-3018.6169	6155.2338	0.0000	1.0000
JC69I	56	-3035.5327	6183.0654	0.0000	1.0000
F81	58	-3075.5308	6267.0617	0.0000	1.0000
JC69	55	-3089.8652	6289.7305	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	60	-2908.2946	6198.7833	0.9472	0.9472
HKYIG	61	-2908.1141	6204.7922	0.0470	0.9942
GTRG	64	-2900.6621	6208.9979	0.0057	0.9999
GTRIG	65	-2901.8724	6217.7883	0.0001	1.0000
K2PG	57	-2930.3403	6223.7649	0.0000	1.0000
K2PIG	58	-2927.5461	6224.5465	0.0000	1.0000
SYMG	61	-2921.6085	6231.7810	0.0000	1.0000
SYMIG	62	-2923.4205	6241.7749	0.0000	1.0000
HKYI	60	-2930.8131	6243.8203	0.0000	1.0000
GTRI	64	-2926.8585	6261.3906	0.0000	1.0000
K2PI	57	-2950.2752	6263.6348	0.0000	1.0000
SYMI	61	-2945.2823	6279.1286	0.0000	1.0000
HKY	59	-2995.3988	6366.6218	0.0000	1.0000
F81G	59	-2996.3857	6368.5956	0.0000	1.0000
GTR	63	-2989.1224	6379.5485	0.0000	1.0000
F81IG	60	-2999.2808	6380.7557	0.0000	1.0000
JC69G	56	-3012.9377	6382.5899	0.0000	1.0000
K2P	56	-3014.7264	6386.1673	0.0000	1.0000
SYM	60	-3003.4301	6389.0543	0.0000	1.0000
JC69IG	57	-3013.6115	6390.3074	0.0000	1.0000
F81I	59	-3018.6169	6413.0579	0.0000	1.0000
JC69I	56	-3035.5327	6427.7798	0.0000	1.0000
F81	58	-3075.5308	6520.5159	0.0000	1.0000
JC69	55	-3089.8652	6530.0751	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/10976_hg18.phylip = 1 - 584;
 Partition Dummy = 1:261-aggressive-phylip/10976_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/10976_hg18.phylip = 1 - 584;
 Partition Dummy = 1:261-aggressive-phylip/10976_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/10976_hg18.phylip = 1 - 584;
 Partition Dummy = 1:261-aggressive-phylip/10976_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

