
Output from mraic.pl version 1.4.4 by Johan Nylander
-------------------------------------------------------------
21:34 04/04/2011

Input data from file "591-tomodel-phylip/chr3_24921.phylip" (7 taxa, 753 characters)

Minimum AIC  model: GTRI
Minimum AICc model: GTRI
Minimum BIC  model: K2PI


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
GTRI	20	-2405.7681	4852.6838	0.4135	0.4135
HKYI	16	-2410.7625	4854.2641	0.1876	0.6011
GTRIG	21	-2405.7933	4854.8505	0.1399	0.7411
GTRG	20	-2407.1244	4855.3963	0.1065	0.8476
HKYIG	17	-2410.7881	4856.4089	0.0642	0.9118
HKYG	16	-2412.0937	4856.9265	0.0496	0.9614
K2PI	13	-2416.2268	4858.9462	0.0181	0.9794
K2PIG	14	-2416.2467	4861.0626	0.0063	0.9857
K2PG	13	-2417.3497	4861.1920	0.0059	0.9916
SYMI	17	-2413.3558	4861.5444	0.0049	0.9965
SYMG	17	-2414.3509	4863.5344	0.0018	0.9983
SYMIG	18	-2413.3717	4863.6752	0.0017	1.0000
GTR	19	-2429.9688	4898.9744	0.0000	1.0000
HKY	15	-2435.1023	4900.8558	0.0000	1.0000
K2P	12	-2440.3005	4905.0227	0.0000	1.0000
SYM	16	-2436.9423	4906.6238	0.0000	1.0000
F81I	15	-2463.0275	4956.7062	0.0000	1.0000
F81IG	16	-2463.0446	4958.8283	0.0000	1.0000
F81G	15	-2464.0917	4958.8346	0.0000	1.0000
JC69G	12	-2468.2498	4960.9211	0.0000	1.0000
JC69IG	13	-2467.2931	4961.0788	0.0000	1.0000
JC69I	12	-2472.7922	4970.0059	0.0000	1.0000
F81	14	-2484.3782	4997.3256	0.0000	1.0000
JC69	11	-2488.5980	4999.5523	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRI	20	-2405.7681	4851.5363	0.4372	0.4372
HKYI	16	-2410.7625	4853.5249	0.1617	0.5989
GTRIG	21	-2405.7933	4853.5865	0.1568	0.7558
GTRG	20	-2407.1244	4854.2488	0.1126	0.8684
HKYIG	17	-2410.7881	4855.5762	0.0580	0.9264
HKYG	16	-2412.0937	4856.1873	0.0427	0.9691
K2PI	13	-2416.2268	4858.4536	0.0138	0.9829
K2PIG	14	-2416.2467	4860.4935	0.0050	0.9878
K2PG	13	-2417.3497	4860.6995	0.0045	0.9923
SYMI	17	-2413.3558	4860.7117	0.0044	0.9967
SYMG	17	-2414.3509	4862.7017	0.0016	0.9984
SYMIG	18	-2413.3717	4862.7434	0.0016	1.0000
GTR	19	-2429.9688	4897.9376	0.0000	1.0000
HKY	15	-2435.1023	4900.2046	0.0000	1.0000
K2P	12	-2440.3005	4904.6011	0.0000	1.0000
SYM	16	-2436.9423	4905.8846	0.0000	1.0000
F81I	15	-2463.0275	4956.0549	0.0000	1.0000
F81IG	16	-2463.0446	4958.0892	0.0000	1.0000
F81G	15	-2464.0917	4958.1833	0.0000	1.0000
JC69G	12	-2468.2498	4960.4995	0.0000	1.0000
JC69IG	13	-2467.2931	4960.5863	0.0000	1.0000
JC69I	12	-2472.7922	4969.5843	0.0000	1.0000
F81	14	-2484.3782	4996.7565	0.0000	1.0000
JC69	11	-2488.5980	4999.1960	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
K2PI	13	-2416.2268	4918.5665	0.7268	0.7268
K2PG	13	-2417.3497	4920.8123	0.2364	0.9632
K2PIG	14	-2416.2467	4925.2304	0.0260	0.9892
HKYI	16	-2410.7625	4927.5100	0.0083	0.9975
HKYG	16	-2412.0937	4930.1724	0.0022	0.9997
HKYIG	17	-2410.7881	4934.1853	0.0003	1.0000
SYMI	17	-2413.3558	4939.3208	0.0000	1.0000
SYMG	17	-2414.3509	4941.3108	0.0000	1.0000
GTRI	20	-2405.7681	4944.0176	0.0000	1.0000
SYMIG	18	-2413.3717	4945.9765	0.0000	1.0000
GTRG	20	-2407.1244	4946.7301	0.0000	1.0000
GTRIG	21	-2405.7933	4950.6919	0.0000	1.0000
K2P	12	-2440.3005	4960.0898	0.0000	1.0000
HKY	15	-2435.1023	4969.5655	0.0000	1.0000
SYM	16	-2436.9423	4979.8697	0.0000	1.0000
GTR	19	-2429.9688	4985.7948	0.0000	1.0000
JC69G	12	-2468.2498	5015.9883	0.0000	1.0000
JC69IG	13	-2467.2931	5020.6991	0.0000	1.0000
JC69I	12	-2472.7922	5025.0731	0.0000	1.0000
F81I	15	-2463.0275	5025.4159	0.0000	1.0000
F81G	15	-2464.0917	5027.5443	0.0000	1.0000
F81IG	16	-2463.0446	5032.0742	0.0000	1.0000
JC69	11	-2488.5980	5050.0607	0.0000	1.0000
F81	14	-2484.3782	5061.4934	0.0000	1.0000



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

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

BEGIN MRBAYES;
 Charset 591-tomodel-phylip/chr3_24921.phylip = 1 - 753;
 Partition Dummy = 1:591-tomodel-phylip/chr3_24921.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=6 rates=propinv;
 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) pinvarpr=Uniform(0.0,1.0);
END;


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

BEGIN MRBAYES;
 Charset 591-tomodel-phylip/chr3_24921.phylip = 1 - 753;
 Partition Dummy = 1:591-tomodel-phylip/chr3_24921.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=6 rates=propinv;
 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) pinvarpr=Uniform(0.0,1.0);
END;


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

BEGIN MRBAYES;
 Charset 591-tomodel-phylip/chr3_24921.phylip = 1 - 753;
 Partition Dummy = 1:591-tomodel-phylip/chr3_24921.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=2 rates=propinv;
 Prset applyto=(1) statefreqpr=Fixed(Equal) tratiopr=Beta(1.0,1.0) pinvarpr=Uniform(0.0,1.0);
END;


-------------------------------------------------------------
End of Output

