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

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

Minimum AIC  model: HKYG
Minimum AICc model: HKYG
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
HKYG	16	-2139.7612	4312.2962	0.3476	0.3476
GTRG	20	-2136.0595	4313.3207	0.2082	0.5558
HKYI	16	-2140.5312	4313.8362	0.1609	0.7167
HKYIG	17	-2139.5494	4313.9706	0.1505	0.8672
GTRIG	21	-2135.9867	4315.2972	0.0775	0.9447
GTRI	20	-2137.3857	4315.9731	0.0553	1.0000
HKY	15	-2167.2838	4365.2494	0.0000	1.0000
GTR	19	-2163.9167	4366.9191	0.0000	1.0000
F81G	15	-2172.1837	4375.0493	0.0000	1.0000
F81IG	16	-2172.0468	4376.8674	0.0000	1.0000
F81I	15	-2174.0386	4378.7590	0.0000	1.0000
K2PG	13	-2185.2127	4396.9411	0.0000	1.0000
K2PIG	14	-2185.1708	4398.9374	0.0000	1.0000
K2PI	13	-2186.4681	4399.4518	0.0000	1.0000
SYMG	17	-2182.9543	4400.7804	0.0000	1.0000
SYMIG	18	-2182.8875	4402.7507	0.0000	1.0000
SYMI	17	-2184.0337	4402.9392	0.0000	1.0000
F81	14	-2199.5514	4427.6984	0.0000	1.0000
JC69G	12	-2210.7565	4445.9543	0.0000	1.0000
JC69IG	13	-2210.6088	4447.7332	0.0000	1.0000
JC69I	12	-2212.6603	4449.7619	0.0000	1.0000
K2P	12	-2213.7706	4451.9824	0.0000	1.0000
SYM	16	-2211.4073	4455.5885	0.0000	1.0000
JC69	11	-2238.4608	4499.2945	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
HKYG	16	-2139.7612	4311.5224	0.3174	0.3174
GTRG	20	-2136.0595	4312.1190	0.2356	0.5530
HKYI	16	-2140.5312	4313.0624	0.1470	0.6999
HKYIG	17	-2139.5494	4313.0988	0.1443	0.8443
GTRIG	21	-2135.9867	4313.9734	0.0932	0.9375
GTRI	20	-2137.3857	4314.7714	0.0625	1.0000
HKY	15	-2167.2838	4364.5676	0.0000	1.0000
GTR	19	-2163.9167	4365.8334	0.0000	1.0000
F81G	15	-2172.1837	4374.3675	0.0000	1.0000
F81IG	16	-2172.0468	4376.0936	0.0000	1.0000
F81I	15	-2174.0386	4378.0772	0.0000	1.0000
K2PG	13	-2185.2127	4396.4255	0.0000	1.0000
K2PIG	14	-2185.1708	4398.3417	0.0000	1.0000
K2PI	13	-2186.4681	4398.9362	0.0000	1.0000
SYMG	17	-2182.9543	4399.9086	0.0000	1.0000
SYMIG	18	-2182.8875	4401.7750	0.0000	1.0000
SYMI	17	-2184.0337	4402.0674	0.0000	1.0000
F81	14	-2199.5514	4427.1027	0.0000	1.0000
JC69G	12	-2210.7565	4445.5130	0.0000	1.0000
JC69IG	13	-2210.6088	4447.2176	0.0000	1.0000
JC69I	12	-2212.6603	4449.3206	0.0000	1.0000
K2P	12	-2213.7706	4451.5411	0.0000	1.0000
SYM	16	-2211.4073	4454.8147	0.0000	1.0000
JC69	11	-2238.4608	4498.9217	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	16	-2139.7612	4384.7904	0.6626	0.6626
HKYI	16	-2140.5312	4386.3304	0.3068	0.9694
HKYIG	17	-2139.5494	4390.9461	0.0305	0.9999
GTRG	20	-2136.0595	4403.7040	0.0001	1.0000
GTRI	20	-2137.3857	4406.3564	0.0000	1.0000
GTRIG	21	-2135.9867	4410.1377	0.0000	1.0000
HKY	15	-2167.2838	4433.2563	0.0000	1.0000
F81G	15	-2172.1837	4443.0562	0.0000	1.0000
F81I	15	-2174.0386	4446.7660	0.0000	1.0000
F81IG	16	-2172.0468	4449.3616	0.0000	1.0000
GTR	19	-2163.9167	4452.8392	0.0000	1.0000
K2PG	13	-2185.2127	4455.9557	0.0000	1.0000
K2PI	13	-2186.4681	4458.4665	0.0000	1.0000
K2PIG	14	-2185.1708	4462.4512	0.0000	1.0000
SYMG	17	-2182.9543	4477.7559	0.0000	1.0000
SYMI	17	-2184.0337	4479.9147	0.0000	1.0000
SYMIG	18	-2182.8875	4484.2015	0.0000	1.0000
F81	14	-2199.5514	4491.2122	0.0000	1.0000
JC69G	12	-2210.7565	4500.4641	0.0000	1.0000
JC69I	12	-2212.6603	4504.2717	0.0000	1.0000
K2P	12	-2213.7706	4506.4921	0.0000	1.0000
JC69IG	13	-2210.6088	4506.7479	0.0000	1.0000
SYM	16	-2211.4073	4528.0827	0.0000	1.0000
JC69	11	-2238.4608	4549.2934	0.0000	1.0000



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

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

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


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

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


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

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

