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

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

Minimum AIC  model: GTRIG
Minimum AICc model: GTRIG
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
GTRIG	65	-4928.9639	9998.9420	0.8039	0.8039
GTRG	64	-4931.5487	10001.7640	0.1961	1.0000
HKYIG	61	-4947.2600	10026.1803	0.0000	1.0000
SYMIG	62	-4946.3754	10026.7405	0.0000	1.0000
HKYG	60	-4949.8014	10028.9395	0.0000	1.0000
SYMG	61	-4948.9763	10029.6129	0.0000	1.0000
K2PIG	58	-4954.4988	10033.7050	0.0000	1.0000
K2PG	57	-4957.2361	10036.8737	0.0000	1.0000
GTRI	64	-4953.5160	10045.6986	0.0000	1.0000
SYMI	61	-4971.8275	10075.3153	0.0000	1.0000
HKYI	60	-4973.2537	10075.8441	0.0000	1.0000
K2PI	57	-4980.0001	10082.4016	0.0000	1.0000
F81IG	60	-5080.8803	10291.0972	0.0000	1.0000
F81G	59	-5083.7814	10294.5818	0.0000	1.0000
JC69IG	57	-5087.6565	10297.7145	0.0000	1.0000
JC69G	56	-5090.5363	10301.1742	0.0000	1.0000
F81I	59	-5103.5675	10334.1541	0.0000	1.0000
JC69I	56	-5110.2816	10340.6648	0.0000	1.0000
GTR	63	-5137.9315	10412.1882	0.0000	1.0000
HKY	59	-5153.2434	10433.5059	0.0000	1.0000
SYM	60	-5154.7992	10438.9351	0.0000	1.0000
K2P	56	-5167.8401	10455.7818	0.0000	1.0000
F81	58	-5282.5521	10689.8116	0.0000	1.0000
JC69	55	-5291.5818	10700.9709	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRIG	65	-4928.9639	9987.9279	0.8299	0.8299
GTRG	64	-4931.5487	9991.0974	0.1701	1.0000
HKYIG	61	-4947.2600	10016.5200	0.0000	1.0000
SYMIG	62	-4946.3754	10016.7507	0.0000	1.0000
HKYG	60	-4949.8014	10019.6028	0.0000	1.0000
SYMG	61	-4948.9763	10019.9527	0.0000	1.0000
K2PIG	58	-4954.4988	10024.9976	0.0000	1.0000
K2PG	57	-4957.2361	10028.4722	0.0000	1.0000
GTRI	64	-4953.5160	10035.0320	0.0000	1.0000
SYMI	61	-4971.8275	10065.6550	0.0000	1.0000
HKYI	60	-4973.2537	10066.5074	0.0000	1.0000
K2PI	57	-4980.0001	10074.0001	0.0000	1.0000
F81IG	60	-5080.8803	10281.7605	0.0000	1.0000
F81G	59	-5083.7814	10285.5627	0.0000	1.0000
JC69IG	57	-5087.6565	10289.3130	0.0000	1.0000
JC69G	56	-5090.5363	10293.0727	0.0000	1.0000
F81I	59	-5103.5675	10325.1350	0.0000	1.0000
JC69I	56	-5110.2816	10332.5633	0.0000	1.0000
GTR	63	-5137.9315	10401.8630	0.0000	1.0000
HKY	59	-5153.2434	10424.4868	0.0000	1.0000
SYM	60	-5154.7992	10429.5984	0.0000	1.0000
K2P	56	-5167.8401	10447.6802	0.0000	1.0000
F81	58	-5282.5521	10681.1042	0.0000	1.0000
JC69	55	-5291.5818	10693.1636	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
GTRG	64	-4931.5487	10294.4149	0.6036	0.6036
GTRIG	65	-4928.9639	10295.9847	0.2753	0.8789
K2PG	57	-4957.2361	10298.6143	0.0739	0.9528
K2PIG	58	-4954.4988	10299.8791	0.0393	0.9921
HKYG	60	-4949.8014	10303.9630	0.0051	0.9972
HKYIG	61	-4947.2600	10305.6196	0.0022	0.9994
SYMG	61	-4948.9763	10309.0522	0.0004	0.9998
SYMIG	62	-4946.3754	10310.5896	0.0002	1.0000
GTRI	64	-4953.5160	10338.3495	0.0000	1.0000
K2PI	57	-4980.0001	10344.1423	0.0000	1.0000
HKYI	60	-4973.2537	10350.8676	0.0000	1.0000
SYMI	61	-4971.8275	10354.7545	0.0000	1.0000
JC69G	56	-5090.5363	10558.4755	0.0000	1.0000
JC69IG	57	-5087.6565	10559.4552	0.0000	1.0000
F81G	59	-5083.7814	10565.1836	0.0000	1.0000
F81IG	60	-5080.8803	10566.1207	0.0000	1.0000
JC69I	56	-5110.2816	10597.9662	0.0000	1.0000
F81I	59	-5103.5675	10604.7558	0.0000	1.0000
GTR	63	-5137.9315	10700.4412	0.0000	1.0000
HKY	59	-5153.2434	10704.1076	0.0000	1.0000
K2P	56	-5167.8401	10713.0831	0.0000	1.0000
SYM	60	-5154.7992	10713.9586	0.0000	1.0000
JC69	55	-5291.5818	10953.8271	0.0000	1.0000
F81	58	-5282.5521	10955.9857	0.0000	1.0000



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

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

BEGIN MRBAYES;
 Charset 261-aggressive-phylip/11541_hg18.phylip = 1 - 845;
 Partition Dummy = 1:261-aggressive-phylip/11541_hg18.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=6 rates=invgamma;
 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) pinvarpr=Uniform(0.0,1.0);
END;


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

BEGIN MRBAYES;
 Charset 261-aggressive-phylip/11541_hg18.phylip = 1 - 845;
 Partition Dummy = 1:261-aggressive-phylip/11541_hg18.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=6 rates=invgamma;
 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) pinvarpr=Uniform(0.0,1.0);
END;


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

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

