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

Input data from file "261-aggressive-phylip/11841_hg18.phylip" (29 taxa, 556 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	-1843.8257	3832.5964	0.5999	0.5999
GTRIG	65	-1843.7175	3834.9452	0.1854	0.7853
GTRI	64	-1845.2290	3835.4030	0.1474	0.9327
HKYG	60	-1851.5348	3837.8575	0.0432	0.9759
HKYIG	61	-1851.3944	3840.1006	0.0141	0.9900
HKYI	60	-1852.9978	3840.7835	0.0100	1.0000
GTR	63	-1861.0458	3864.4818	0.0000	1.0000
HKY	59	-1869.5781	3871.4305	0.0000	1.0000
K2PG	57	-1876.5354	3880.3479	0.0000	1.0000
K2PIG	58	-1876.3352	3882.4411	0.0000	1.0000
K2PI	57	-1878.1062	3883.4896	0.0000	1.0000
SYMG	61	-1874.1990	3885.7098	0.0000	1.0000
SYMIG	62	-1874.0320	3887.9098	0.0000	1.0000
SYMI	61	-1875.6982	3888.7081	0.0000	1.0000
K2P	56	-1895.7332	3916.2599	0.0000	1.0000
F81G	59	-1893.1663	3918.6068	0.0000	1.0000
SYM	60	-1892.2261	3919.2400	0.0000	1.0000
F81IG	60	-1893.0286	3920.8452	0.0000	1.0000
F81I	59	-1894.5229	3921.3201	0.0000	1.0000
F81	58	-1911.0081	3951.7868	0.0000	1.0000
JC69G	56	-1914.0313	3952.8561	0.0000	1.0000
JC69IG	57	-1913.8588	3954.9947	0.0000	1.0000
JC69I	56	-1915.4957	3955.7850	0.0000	1.0000
JC69	55	-1932.8441	3988.0082	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRG	64	-1843.8257	3815.6514	0.5895	0.5895
GTRIG	65	-1843.7175	3817.4350	0.2417	0.8312
GTRI	64	-1845.2290	3818.4580	0.1449	0.9761
HKYG	60	-1851.5348	3823.0696	0.0144	0.9905
HKYIG	61	-1851.3944	3824.7888	0.0061	0.9967
HKYI	60	-1852.9978	3825.9956	0.0033	1.0000
GTR	63	-1861.0458	3848.0916	0.0000	1.0000
HKY	59	-1869.5781	3857.1563	0.0000	1.0000
K2PG	57	-1876.5354	3867.0708	0.0000	1.0000
K2PIG	58	-1876.3352	3868.6705	0.0000	1.0000
K2PI	57	-1878.1062	3870.2125	0.0000	1.0000
SYMG	61	-1874.1990	3870.3981	0.0000	1.0000
SYMIG	62	-1874.0320	3872.0640	0.0000	1.0000
SYMI	61	-1875.6982	3873.3964	0.0000	1.0000
K2P	56	-1895.7332	3903.4663	0.0000	1.0000
F81G	59	-1893.1663	3904.3326	0.0000	1.0000
SYM	60	-1892.2261	3904.4521	0.0000	1.0000
F81IG	60	-1893.0286	3906.0573	0.0000	1.0000
F81I	59	-1894.5229	3907.0459	0.0000	1.0000
F81	58	-1911.0081	3938.0162	0.0000	1.0000
JC69G	56	-1914.0313	3940.0625	0.0000	1.0000
JC69IG	57	-1913.8588	3941.7176	0.0000	1.0000
JC69I	56	-1915.4957	3942.9914	0.0000	1.0000
JC69	55	-1932.8441	3975.6882	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	60	-1851.5348	4082.3157	0.7754	0.7754
HKYI	60	-1852.9978	4085.2417	0.1795	0.9549
HKYIG	61	-1851.3944	4088.3557	0.0378	0.9928
GTRG	64	-1843.8257	4092.1805	0.0056	0.9984
GTRI	64	-1845.2290	4094.9872	0.0014	0.9997
GTRIG	65	-1843.7175	4098.2850	0.0003	1.0000
HKY	59	-1869.5781	4112.0816	0.0000	1.0000
K2PG	57	-1876.5354	4113.3546	0.0000	1.0000
K2PI	57	-1878.1062	4116.4963	0.0000	1.0000
K2PIG	58	-1876.3352	4119.2750	0.0000	1.0000
GTR	63	-1861.0458	4120.3000	0.0000	1.0000
SYMG	61	-1874.1990	4133.9650	0.0000	1.0000
SYMI	61	-1875.6982	4136.9632	0.0000	1.0000
SYMIG	62	-1874.0320	4139.9516	0.0000	1.0000
K2P	56	-1895.7332	4145.4293	0.0000	1.0000
F81G	59	-1893.1663	4159.2579	0.0000	1.0000
F81I	59	-1894.5229	4161.9712	0.0000	1.0000
SYM	60	-1892.2261	4163.6982	0.0000	1.0000
F81IG	60	-1893.0286	4165.3034	0.0000	1.0000
JC69G	56	-1914.0313	4182.0255	0.0000	1.0000
JC69I	56	-1915.4957	4184.9544	0.0000	1.0000
JC69IG	57	-1913.8588	4188.0014	0.0000	1.0000
F81	58	-1911.0081	4188.6208	0.0000	1.0000
JC69	55	-1932.8441	4213.3305	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/11841_hg18.phylip = 1 - 556;
 Partition Dummy = 1:261-aggressive-phylip/11841_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/11841_hg18.phylip = 1 - 556;
 Partition Dummy = 1:261-aggressive-phylip/11841_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/11841_hg18.phylip = 1 - 556;
 Partition Dummy = 1:261-aggressive-phylip/11841_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

