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

Input data from file "591-tomodel-phylip/chrZ_7924.phylip" (7 taxa, 825 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	-2102.9694	4238.6121	0.4587	0.4587
HKYI	16	-2103.4019	4239.4771	0.2977	0.7564
HKYIG	17	-2102.9463	4240.6510	0.1655	0.9219
GTRG	20	-2101.2770	4243.5988	0.0379	0.9598
GTRI	20	-2101.6243	4244.2933	0.0268	0.9866
GTRIG	21	-2101.2610	4245.6727	0.0134	1.0000
HKY	15	-2115.5195	4261.6324	0.0000	1.0000
GTR	19	-2114.0117	4266.9675	0.0000	1.0000
K2PG	13	-2127.2406	4280.9301	0.0000	1.0000
K2PI	13	-2127.7143	4281.8774	0.0000	1.0000
K2PIG	14	-2127.1978	4282.9142	0.0000	1.0000
SYMG	17	-2125.7554	4286.2692	0.0000	1.0000
SYMI	17	-2126.2795	4287.3173	0.0000	1.0000
SYMIG	18	-2125.6981	4288.2449	0.0000	1.0000
K2P	12	-2139.6748	4303.7338	0.0000	1.0000
SYM	16	-2138.2794	4309.2320	0.0000	1.0000
F81G	15	-2151.3110	4333.2152	0.0000	1.0000
F81I	15	-2151.5740	4333.7414	0.0000	1.0000
F81IG	16	-2151.3112	4335.2956	0.0000	1.0000
F81	14	-2162.6199	4353.7582	0.0000	1.0000
JC69G	12	-2174.0107	4372.4056	0.0000	1.0000
JC69I	12	-2174.3091	4373.0025	0.0000	1.0000
JC69IG	13	-2174.0008	4374.4503	0.0000	1.0000
JC69	11	-2185.3192	4392.9632	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
HKYG	16	-2102.9694	4237.9388	0.4480	0.4480
HKYI	16	-2103.4019	4238.8038	0.2907	0.7386
HKYIG	17	-2102.9463	4239.8926	0.1686	0.9073
GTRG	20	-2101.2770	4242.5540	0.0446	0.9518
GTRI	20	-2101.6243	4243.2486	0.0315	0.9833
GTRIG	21	-2101.2610	4244.5220	0.0167	1.0000
HKY	15	-2115.5195	4261.0391	0.0000	1.0000
GTR	19	-2114.0117	4266.0234	0.0000	1.0000
K2PG	13	-2127.2406	4280.4813	0.0000	1.0000
K2PI	13	-2127.7143	4281.4286	0.0000	1.0000
K2PIG	14	-2127.1978	4282.3957	0.0000	1.0000
SYMG	17	-2125.7554	4285.5108	0.0000	1.0000
SYMI	17	-2126.2795	4286.5590	0.0000	1.0000
SYMIG	18	-2125.6981	4287.3963	0.0000	1.0000
K2P	12	-2139.6748	4303.3496	0.0000	1.0000
SYM	16	-2138.2794	4308.5587	0.0000	1.0000
F81G	15	-2151.3110	4332.6219	0.0000	1.0000
F81I	15	-2151.5740	4333.1481	0.0000	1.0000
F81IG	16	-2151.3112	4334.6223	0.0000	1.0000
F81	14	-2162.6199	4353.2397	0.0000	1.0000
JC69G	12	-2174.0107	4372.0214	0.0000	1.0000
JC69I	12	-2174.3091	4372.6183	0.0000	1.0000
JC69IG	13	-2174.0008	4374.0015	0.0000	1.0000
JC69	11	-2185.3192	4392.6384	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	16	-2102.9694	4313.3850	0.5936	0.5936
HKYI	16	-2103.4019	4314.2499	0.3852	0.9788
HKYIG	17	-2102.9463	4320.0542	0.0211	0.9999
HKY	15	-2115.5195	4331.7699	0.0001	1.0000
GTRG	20	-2101.2770	4336.8616	0.0000	1.0000
GTRI	20	-2101.6243	4337.5562	0.0000	1.0000
K2PG	13	-2127.2406	4341.7813	0.0000	1.0000
K2PI	13	-2127.7143	4342.7286	0.0000	1.0000
GTRIG	21	-2101.2610	4343.5451	0.0000	1.0000
K2PIG	14	-2127.1978	4348.4110	0.0000	1.0000
GTR	19	-2114.0117	4355.6157	0.0000	1.0000
K2P	12	-2139.6748	4359.9342	0.0000	1.0000
SYMG	17	-2125.7554	4365.6724	0.0000	1.0000
SYMI	17	-2126.2795	4366.7205	0.0000	1.0000
SYMIG	18	-2125.6981	4372.2732	0.0000	1.0000
SYM	16	-2138.2794	4384.0049	0.0000	1.0000
F81G	15	-2151.3110	4403.3527	0.0000	1.0000
F81I	15	-2151.5740	4403.8788	0.0000	1.0000
F81IG	16	-2151.3112	4410.0685	0.0000	1.0000
F81	14	-2162.6199	4419.2551	0.0000	1.0000
JC69G	12	-2174.0107	4428.6060	0.0000	1.0000
JC69I	12	-2174.3091	4429.2029	0.0000	1.0000
JC69IG	13	-2174.0008	4435.3015	0.0000	1.0000
JC69	11	-2185.3192	4444.5077	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_7924.phylip = 1 - 825;
 Partition Dummy = 1:591-tomodel-phylip/chrZ_7924.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_7924.phylip = 1 - 825;
 Partition Dummy = 1:591-tomodel-phylip/chrZ_7924.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_7924.phylip = 1 - 825;
 Partition Dummy = 1:591-tomodel-phylip/chrZ_7924.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

