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

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

Minimum AIC  model: K2PI
Minimum AICc model: K2PI
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
K2PI	13	-1690.9484	3408.3751	0.2156	0.2156
HKYI	16	-1688.0775	3408.8727	0.1681	0.3838
K2PIG	14	-1690.2245	3409.0017	0.1576	0.5414
K2PG	13	-1691.3461	3409.1706	0.1449	0.6863
HKYIG	17	-1687.4043	3409.6171	0.1159	0.8022
HKYG	16	-1688.5164	3409.7505	0.1084	0.9106
SYMI	17	-1688.8658	3412.5401	0.0269	0.9374
SYMG	17	-1689.3146	3413.4376	0.0172	0.9546
SYMIG	18	-1688.2950	3413.4948	0.0167	0.9713
GTRI	20	-1686.4284	3413.9709	0.0131	0.9844
GTRG	20	-1686.9433	3415.0007	0.0079	0.9923
GTRIG	21	-1685.9012	3415.0295	0.0077	1.0000
JC69I	12	-1710.4236	3445.2567	0.0000	1.0000
JC69IG	13	-1709.4871	3445.4525	0.0000	1.0000
JC69G	12	-1710.7931	3445.9957	0.0000	1.0000
F81I	15	-1707.7094	3446.0511	0.0000	1.0000
F81IG	16	-1706.7627	3446.2431	0.0000	1.0000
F81G	15	-1708.0963	3446.8251	0.0000	1.0000
K2P	12	-1715.4172	3455.2439	0.0000	1.0000
HKY	15	-1712.7987	3456.2299	0.0000	1.0000
SYM	16	-1713.7763	3460.2703	0.0000	1.0000
GTR	19	-1711.6544	3462.3154	0.0000	1.0000
JC69	11	-1734.6740	3491.6939	0.0000	1.0000
F81	14	-1732.0762	3492.7050	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
K2PI	13	-1690.9484	3407.8968	0.1987	0.1987
HKYI	16	-1688.0775	3408.1550	0.1746	0.3734
K2PIG	14	-1690.2245	3408.4491	0.1508	0.5241
K2PG	13	-1691.3461	3408.6922	0.1335	0.6576
HKYIG	17	-1687.4043	3408.8086	0.1260	0.7836
HKYG	16	-1688.5164	3409.0328	0.1126	0.8962
SYMI	17	-1688.8658	3411.7317	0.0292	0.9254
SYMIG	18	-1688.2950	3412.5900	0.0190	0.9444
SYMG	17	-1689.3146	3412.6291	0.0186	0.9630
GTRI	20	-1686.4284	3412.8569	0.0166	0.9797
GTRIG	21	-1685.9012	3413.8024	0.0104	0.9901
GTRG	20	-1686.9433	3413.8866	0.0099	1.0000
JC69I	12	-1710.4236	3444.8472	0.0000	1.0000
JC69IG	13	-1709.4871	3444.9742	0.0000	1.0000
F81I	15	-1707.7094	3445.4187	0.0000	1.0000
F81IG	16	-1706.7627	3445.5254	0.0000	1.0000
JC69G	12	-1710.7931	3445.5863	0.0000	1.0000
F81G	15	-1708.0963	3446.1927	0.0000	1.0000
K2P	12	-1715.4172	3454.8345	0.0000	1.0000
HKY	15	-1712.7987	3455.5975	0.0000	1.0000
SYM	16	-1713.7763	3459.5527	0.0000	1.0000
GTR	19	-1711.6544	3461.3088	0.0000	1.0000
JC69	11	-1734.6740	3491.3479	0.0000	1.0000
F81	14	-1732.0762	3492.1523	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
K2PI	13	-1690.9484	3468.3840	0.5723	0.5723
K2PG	13	-1691.3461	3469.1795	0.3845	0.9568
K2PIG	14	-1690.2245	3473.5891	0.0424	0.9992
HKYI	16	-1688.0775	3482.6008	0.0005	0.9997
HKYG	16	-1688.5164	3483.4786	0.0003	1.0000
HKYIG	17	-1687.4043	3487.9073	0.0000	1.0000
SYMI	17	-1688.8658	3490.8304	0.0000	1.0000
SYMG	17	-1689.3146	3491.7278	0.0000	1.0000
SYMIG	18	-1688.2950	3496.3416	0.0000	1.0000
JC69I	12	-1710.4236	3500.6816	0.0000	1.0000
JC69G	12	-1710.7931	3501.4206	0.0000	1.0000
JC69IG	13	-1709.4871	3505.4614	0.0000	1.0000
GTRI	20	-1686.4284	3505.9141	0.0000	1.0000
GTRG	20	-1686.9433	3506.9439	0.0000	1.0000
K2P	12	-1715.4172	3510.6688	0.0000	1.0000
GTRIG	21	-1685.9012	3511.5125	0.0000	1.0000
F81I	15	-1707.7094	3515.2117	0.0000	1.0000
F81G	15	-1708.0963	3515.9856	0.0000	1.0000
F81IG	16	-1706.7627	3519.9712	0.0000	1.0000
HKY	15	-1712.7987	3525.3904	0.0000	1.0000
SYM	16	-1713.7763	3533.9985	0.0000	1.0000
JC69	11	-1734.6740	3542.5294	0.0000	1.0000
GTR	19	-1711.6544	3549.7132	0.0000	1.0000
F81	14	-1732.0762	3557.2924	0.0000	1.0000



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

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

BEGIN MRBAYES;
 Charset 591-tomodel-phylip/chr5_11018.phylip = 1 - 775;
 Partition Dummy = 1:591-tomodel-phylip/chr5_11018.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;


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

BEGIN MRBAYES;
 Charset 591-tomodel-phylip/chr5_11018.phylip = 1 - 775;
 Partition Dummy = 1:591-tomodel-phylip/chr5_11018.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;


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

BEGIN MRBAYES;
 Charset 591-tomodel-phylip/chr5_11018.phylip = 1 - 775;
 Partition Dummy = 1:591-tomodel-phylip/chr5_11018.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

