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

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

Minimum AIC  model: GTRG
Minimum AICc model: GTRG
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
GTRG	64	-2245.7816	4637.8489	0.5330	0.5330
SYMG	61	-2250.8821	4640.2794	0.1581	0.6911
GTRIG	65	-2245.6955	4640.2897	0.1573	0.8485
SYMIG	62	-2250.7604	4642.6149	0.0492	0.8976
K2PG	57	-2257.3023	4642.9163	0.0423	0.9399
HKYG	60	-2254.4555	4644.8586	0.0160	0.9560
GTRI	64	-2249.3141	4644.9140	0.0156	0.9715
K2PIG	58	-2257.0361	4644.9182	0.0155	0.9871
HKYIG	61	-2254.2250	4646.9653	0.0056	0.9927
SYMI	61	-2254.2687	4647.0526	0.0053	0.9980
K2PI	57	-2260.6599	4649.6316	0.0015	0.9995
HKYI	60	-2257.9173	4651.7824	0.0005	1.0000
GTR	63	-2272.1123	4687.9088	0.0000	1.0000
SYM	60	-2277.7739	4691.4955	0.0000	1.0000
HKY	59	-2283.7909	4700.9731	0.0000	1.0000
K2P	56	-2288.4057	4702.5998	0.0000	1.0000
JC69G	56	-2345.3316	4816.4516	0.0000	1.0000
F81G	59	-2341.6122	4816.6157	0.0000	1.0000
JC69IG	57	-2345.1846	4818.6809	0.0000	1.0000
F81IG	60	-2341.4595	4818.8667	0.0000	1.0000
JC69I	56	-2348.2415	4822.2713	0.0000	1.0000
F81I	59	-2344.5973	4822.5859	0.0000	1.0000
JC69	55	-2373.9892	4871.2542	0.0000	1.0000
F81	58	-2370.4913	4871.8285	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRG	64	-2245.7816	4619.5632	0.6105	0.6105
GTRIG	65	-2245.6955	4621.3910	0.2448	0.8553
SYMG	61	-2250.8821	4623.7641	0.0747	0.9300
SYMIG	62	-2250.7604	4625.5208	0.0310	0.9610
GTRI	64	-2249.3141	4626.6283	0.0178	0.9789
K2PG	57	-2257.3023	4628.6046	0.0066	0.9855
HKYG	60	-2254.4555	4628.9109	0.0057	0.9912
K2PIG	58	-2257.0361	4630.0722	0.0032	0.9944
HKYIG	61	-2254.2250	4630.4500	0.0026	0.9971
SYMI	61	-2254.2687	4630.5374	0.0025	0.9996
K2PI	57	-2260.6599	4635.3199	0.0002	0.9998
HKYI	60	-2257.9173	4635.8346	0.0002	1.0000
GTR	63	-2272.1123	4670.2245	0.0000	1.0000
SYM	60	-2277.7739	4675.5478	0.0000	1.0000
HKY	59	-2283.7909	4685.5818	0.0000	1.0000
K2P	56	-2288.4057	4688.8114	0.0000	1.0000
F81G	59	-2341.6122	4801.2244	0.0000	1.0000
JC69G	56	-2345.3316	4802.6633	0.0000	1.0000
F81IG	60	-2341.4595	4802.9190	0.0000	1.0000
JC69IG	57	-2345.1846	4804.3693	0.0000	1.0000
F81I	59	-2344.5973	4807.1946	0.0000	1.0000
JC69I	56	-2348.2415	4808.4829	0.0000	1.0000
F81	58	-2370.4913	4856.9825	0.0000	1.0000
JC69	55	-2373.9892	4857.9784	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
K2PG	57	-2257.3023	4871.0728	0.9123	0.9123
K2PIG	58	-2257.0361	4876.7943	0.0522	0.9645
K2PI	57	-2260.6599	4877.7881	0.0318	0.9963
SYMG	61	-2250.8821	4883.2477	0.0021	0.9984
HKYG	60	-2254.4555	4884.1406	0.0013	0.9997
SYMIG	62	-2250.7604	4889.2581	0.0001	0.9998
HKYIG	61	-2254.2250	4889.9336	0.0001	0.9999
SYMI	61	-2254.2687	4890.0209	0.0001	0.9999
HKYI	60	-2257.9173	4891.0644	0.0000	1.0000
GTRG	64	-2245.7816	4891.8082	0.0000	1.0000
GTRIG	65	-2245.6955	4897.8899	0.0000	1.0000
GTRI	64	-2249.3141	4898.8733	0.0000	1.0000
K2P	56	-2288.4057	4927.0259	0.0000	1.0000
SYM	60	-2277.7739	4930.7775	0.0000	1.0000
HKY	59	-2283.7909	4936.5577	0.0000	1.0000
GTR	63	-2272.1123	4938.2158	0.0000	1.0000
JC69G	56	-2345.3316	5040.8777	0.0000	1.0000
JC69I	56	-2348.2415	5046.6973	0.0000	1.0000
JC69IG	57	-2345.1846	5046.8375	0.0000	1.0000
F81G	59	-2341.6122	5052.2003	0.0000	1.0000
F81IG	60	-2341.4595	5058.1487	0.0000	1.0000
F81I	59	-2344.5973	5058.1705	0.0000	1.0000
JC69	55	-2373.9892	5091.9390	0.0000	1.0000
F81	58	-2370.4913	5103.7046	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/2250_hg18.phylip = 1 - 520;
 Partition Dummy = 1:261-aggressive-phylip/2250_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/2250_hg18.phylip = 1 - 520;
 Partition Dummy = 1:261-aggressive-phylip/2250_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 (K2PG)]

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

