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

Input data from file "261-aggressive-phylip/7355_hg18.phylip" (29 taxa, 582 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	-2045.7024	4235.4977	0.5500	0.5500
GTRIG	65	-2044.7558	4236.1396	0.3990	0.9490
HKYG	60	-2053.8809	4241.8116	0.0234	0.9725
HKYIG	61	-2052.9645	4242.4752	0.0168	0.9892
GTRI	64	-2049.6678	4243.4285	0.0104	0.9997
HKYI	60	-2058.1593	4250.3684	0.0003	1.0000
SYMG	61	-2078.3917	4293.3296	0.0000	1.0000
K2PG	57	-2083.8598	4294.3380	0.0000	1.0000
SYMIG	62	-2079.2352	4297.5225	0.0000	1.0000
K2PIG	58	-2084.6680	4298.4220	0.0000	1.0000
SYMI	61	-2081.6111	4299.7683	0.0000	1.0000
K2PI	57	-2086.9999	4300.6181	0.0000	1.0000
GTR	63	-2085.6118	4312.7911	0.0000	1.0000
HKY	59	-2094.6722	4320.9077	0.0000	1.0000
F81G	59	-2106.4131	4344.3894	0.0000	1.0000
F81IG	60	-2105.5068	4345.0635	0.0000	1.0000
F81I	59	-2109.5575	4350.6782	0.0000	1.0000
SYM	60	-2120.4875	4375.0250	0.0000	1.0000
K2P	56	-2127.4966	4379.1533	0.0000	1.0000
JC69G	56	-2132.9662	4390.0925	0.0000	1.0000
JC69IG	57	-2134.0294	4394.6771	0.0000	1.0000
JC69I	56	-2136.1513	4396.4626	0.0000	1.0000
F81	58	-2146.9321	4422.9502	0.0000	1.0000
JC69	55	-2174.3627	4470.4365	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRG	64	-2045.7024	4219.4048	0.5009	0.5009
GTRIG	65	-2044.7558	4219.5117	0.4748	0.9757
GTRI	64	-2049.6678	4227.3357	0.0095	0.9852
HKYG	60	-2053.8809	4227.7617	0.0077	0.9928
HKYIG	61	-2052.9645	4227.9291	0.0071	0.9999
HKYI	60	-2058.1593	4236.3185	0.0001	1.0000
SYMG	61	-2078.3917	4278.7834	0.0000	1.0000
K2PG	57	-2083.8598	4281.7196	0.0000	1.0000
SYMIG	62	-2079.2352	4282.4705	0.0000	1.0000
SYMI	61	-2081.6111	4285.2221	0.0000	1.0000
K2PIG	58	-2084.6680	4285.3359	0.0000	1.0000
K2PI	57	-2086.9999	4287.9998	0.0000	1.0000
GTR	63	-2085.6118	4297.2235	0.0000	1.0000
HKY	59	-2094.6722	4307.3445	0.0000	1.0000
F81G	59	-2106.4131	4330.8262	0.0000	1.0000
F81IG	60	-2105.5068	4331.0136	0.0000	1.0000
F81I	59	-2109.5575	4337.1150	0.0000	1.0000
SYM	60	-2120.4875	4360.9751	0.0000	1.0000
K2P	56	-2127.4966	4366.9933	0.0000	1.0000
JC69G	56	-2132.9662	4377.9325	0.0000	1.0000
JC69IG	57	-2134.0294	4382.0588	0.0000	1.0000
JC69I	56	-2136.1513	4384.3026	0.0000	1.0000
F81	58	-2146.9321	4409.8641	0.0000	1.0000
JC69	55	-2174.3627	4458.7255	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	60	-2053.8809	4489.7499	0.8855	0.8855
HKYIG	61	-2052.9645	4494.2838	0.0918	0.9772
HKYI	60	-2058.1593	4498.3067	0.0123	0.9895
GTRG	64	-2045.7024	4498.8589	0.0093	0.9988
GTRIG	65	-2044.7558	4503.3322	0.0010	0.9998
GTRI	64	-2049.6678	4506.7898	0.0002	1.0000
K2PG	57	-2083.8598	4530.6085	0.0000	1.0000
K2PI	57	-2086.9999	4536.8886	0.0000	1.0000
K2PIG	58	-2084.6680	4538.5912	0.0000	1.0000
SYMG	61	-2078.3917	4545.1381	0.0000	1.0000
SYMI	61	-2081.6111	4551.5768	0.0000	1.0000
SYMIG	62	-2079.2352	4553.1916	0.0000	1.0000
HKY	59	-2094.6722	4564.9662	0.0000	1.0000
GTR	63	-2085.6118	4572.3112	0.0000	1.0000
F81G	59	-2106.4131	4588.4480	0.0000	1.0000
F81IG	60	-2105.5068	4593.0018	0.0000	1.0000
F81I	59	-2109.5575	4594.7367	0.0000	1.0000
K2P	56	-2127.4966	4611.5156	0.0000	1.0000
JC69G	56	-2132.9662	4622.4548	0.0000	1.0000
SYM	60	-2120.4875	4622.9633	0.0000	1.0000
JC69I	56	-2136.1513	4628.8249	0.0000	1.0000
JC69IG	57	-2134.0294	4630.9476	0.0000	1.0000
F81	58	-2146.9321	4663.1194	0.0000	1.0000
JC69	55	-2174.3627	4698.8814	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/7355_hg18.phylip = 1 - 582;
 Partition Dummy = 1:261-aggressive-phylip/7355_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/7355_hg18.phylip = 1 - 582;
 Partition Dummy = 1:261-aggressive-phylip/7355_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/7355_hg18.phylip = 1 - 582;
 Partition Dummy = 1:261-aggressive-phylip/7355_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

