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

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

Minimum AIC  model: SYMG
Minimum AICc model: SYMG
Minimum BIC  model: SYMG


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
SYMG	17	-2079.0393	4192.9989	0.4361	0.4361
SYMI	17	-2079.5911	4194.1025	0.2511	0.6872
SYMIG	18	-2079.0114	4195.0530	0.1561	0.8434
GTRG	20	-2077.5273	4196.3235	0.0827	0.9261
GTRI	20	-2078.1813	4197.6314	0.0430	0.9691
GTRIG	21	-2077.4486	4198.2951	0.0309	1.0000
K2PG	13	-2097.5173	4221.5787	0.0000	1.0000
K2PI	13	-2097.7572	4222.0585	0.0000	1.0000
K2PIG	14	-2097.4994	4223.6276	0.0000	1.0000
HKYG	16	-2097.0986	4227.0141	0.0000	1.0000
HKYI	16	-2097.3577	4227.5322	0.0000	1.0000
HKYIG	17	-2097.0824	4229.0851	0.0000	1.0000
SYM	16	-2108.6375	4250.0919	0.0000	1.0000
GTR	19	-2107.2298	4253.6059	0.0000	1.0000
K2P	12	-2127.3735	4279.2126	0.0000	1.0000
HKY	15	-2126.9117	4284.5430	0.0000	1.0000
JC69G	12	-2139.8690	4304.2036	0.0000	1.0000
JC69I	12	-2140.0084	4304.4824	0.0000	1.0000
JC69IG	13	-2139.8779	4306.2999	0.0000	1.0000
F81G	15	-2139.4756	4309.6709	0.0000	1.0000
F81I	15	-2139.6208	4309.9612	0.0000	1.0000
F81IG	16	-2139.4764	4311.7696	0.0000	1.0000
JC69	11	-2166.5596	4355.5126	0.0000	1.0000
F81	14	-2166.1788	4360.9862	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
SYMG	17	-2079.0393	4192.0786	0.4189	0.4189
SYMI	17	-2079.5911	4193.1822	0.2412	0.6601
SYMIG	18	-2079.0114	4194.0228	0.1585	0.8186
GTRG	20	-2077.5273	4195.0547	0.0946	0.9132
GTRI	20	-2078.1813	4196.3625	0.0492	0.9624
GTRIG	21	-2077.4486	4196.8972	0.0376	1.0000
K2PG	13	-2097.5173	4221.0346	0.0000	1.0000
K2PI	13	-2097.7572	4221.5144	0.0000	1.0000
K2PIG	14	-2097.4994	4222.9989	0.0000	1.0000
HKYG	16	-2097.0986	4226.1973	0.0000	1.0000
HKYI	16	-2097.3577	4226.7153	0.0000	1.0000
HKYIG	17	-2097.0824	4228.1648	0.0000	1.0000
SYM	16	-2108.6375	4249.2751	0.0000	1.0000
GTR	19	-2107.2298	4252.4596	0.0000	1.0000
K2P	12	-2127.3735	4278.7469	0.0000	1.0000
HKY	15	-2126.9117	4283.8233	0.0000	1.0000
JC69G	12	-2139.8690	4303.7379	0.0000	1.0000
JC69I	12	-2140.0084	4304.0167	0.0000	1.0000
JC69IG	13	-2139.8779	4305.7558	0.0000	1.0000
F81G	15	-2139.4756	4308.9512	0.0000	1.0000
F81I	15	-2139.6208	4309.2415	0.0000	1.0000
F81IG	16	-2139.4764	4310.9528	0.0000	1.0000
JC69	11	-2166.5596	4355.1191	0.0000	1.0000
F81	14	-2166.1788	4360.3575	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
SYMG	17	-2079.0393	4269.0290	0.6159	0.6159
SYMI	17	-2079.5911	4270.1326	0.3547	0.9706
SYMIG	18	-2079.0114	4275.4997	0.0242	0.9948
K2PG	13	-2097.5173	4279.8791	0.0027	0.9975
K2PI	13	-2097.7572	4280.3588	0.0021	0.9996
GTRG	20	-2077.5273	4285.5846	0.0002	0.9998
K2PIG	14	-2097.4994	4286.3698	0.0001	0.9999
GTRI	20	-2078.1813	4286.8924	0.0001	1.0000
GTRIG	21	-2077.4486	4291.9536	0.0000	1.0000
HKYG	16	-2097.0986	4298.6212	0.0000	1.0000
HKYI	16	-2097.3577	4299.1393	0.0000	1.0000
HKYIG	17	-2097.0824	4305.1152	0.0000	1.0000
SYM	16	-2108.6375	4321.6990	0.0000	1.0000
K2P	12	-2127.3735	4333.0649	0.0000	1.0000
GTR	19	-2107.2298	4338.4630	0.0000	1.0000
HKY	15	-2126.9117	4351.7207	0.0000	1.0000
JC69G	12	-2139.8690	4358.0559	0.0000	1.0000
JC69I	12	-2140.0084	4358.3347	0.0000	1.0000
JC69IG	13	-2139.8779	4364.6003	0.0000	1.0000
F81G	15	-2139.4756	4376.8486	0.0000	1.0000
F81I	15	-2139.6208	4377.1390	0.0000	1.0000
F81IG	16	-2139.4764	4383.3767	0.0000	1.0000
JC69	11	-2166.5596	4404.9106	0.0000	1.0000
F81	14	-2166.1788	4423.7284	0.0000	1.0000



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

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

BEGIN MRBAYES;
 Charset 591-tomodel-phylip/chr2_17644.phylip = 1 - 683;
 Partition Dummy = 1:591-tomodel-phylip/chr2_17644.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=Fixed(Equal);
END;


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

BEGIN MRBAYES;
 Charset 591-tomodel-phylip/chr2_17644.phylip = 1 - 683;
 Partition Dummy = 1:591-tomodel-phylip/chr2_17644.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=Fixed(Equal);
END;


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

BEGIN MRBAYES;
 Charset 591-tomodel-phylip/chr2_17644.phylip = 1 - 683;
 Partition Dummy = 1:591-tomodel-phylip/chr2_17644.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=Fixed(Equal);
END;


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

