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

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

Minimum AIC  model: GTRG
Minimum AICc model: GTRG
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
GTRG	20	-2100.9745	4243.1753	0.3741	0.3741
GTRI	20	-2101.3994	4244.0250	0.2446	0.6187
GTRIG	21	-2100.9745	4245.2999	0.1293	0.7480
HKYG	16	-2106.2962	4245.3819	0.1241	0.8721
HKYI	16	-2106.6902	4246.1699	0.0837	0.9558
HKYIG	17	-2106.2974	4247.4843	0.0434	0.9992
GTR	19	-2108.6538	4256.4154	0.0005	0.9997
HKY	15	-2113.3487	4257.3930	0.0003	1.0000
SYMG	17	-2128.5195	4291.9285	0.0000	1.0000
SYMI	17	-2128.7193	4292.3282	0.0000	1.0000
K2PG	13	-2133.5568	4293.6397	0.0000	1.0000
SYMIG	18	-2128.5381	4294.0718	0.0000	1.0000
K2PI	13	-2133.7779	4294.0818	0.0000	1.0000
K2PIG	14	-2133.5779	4295.7637	0.0000	1.0000
SYM	16	-2135.2318	4303.2531	0.0000	1.0000
K2P	12	-2140.2881	4305.0263	0.0000	1.0000
F81G	15	-2162.4820	4355.6596	0.0000	1.0000
F81I	15	-2162.7635	4356.2227	0.0000	1.0000
F81IG	16	-2162.4892	4357.7680	0.0000	1.0000
F81	14	-2168.8996	4366.4071	0.0000	1.0000
JC69G	12	-2185.1955	4394.8413	0.0000	1.0000
JC69I	12	-2185.3994	4395.2491	0.0000	1.0000
JC69IG	13	-2185.2099	4396.9457	0.0000	1.0000
JC69	11	-2191.4067	4405.1937	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRG	20	-2100.9745	4241.9490	0.3894	0.3894
GTRI	20	-2101.3994	4242.7987	0.2546	0.6440
GTRIG	21	-2100.9745	4243.9491	0.1432	0.7872
HKYG	16	-2106.2962	4244.5924	0.1038	0.8911
HKYI	16	-2106.6902	4245.3803	0.0700	0.9611
HKYIG	17	-2106.2974	4246.5947	0.0382	0.9993
GTR	19	-2108.6538	4255.3075	0.0005	0.9998
HKY	15	-2113.3487	4256.6974	0.0002	1.0000
SYMG	17	-2128.5195	4291.0389	0.0000	1.0000
SYMI	17	-2128.7193	4291.4386	0.0000	1.0000
SYMIG	18	-2128.5381	4293.0761	0.0000	1.0000
K2PG	13	-2133.5568	4293.1137	0.0000	1.0000
K2PI	13	-2133.7779	4293.5558	0.0000	1.0000
K2PIG	14	-2133.5779	4295.1558	0.0000	1.0000
SYM	16	-2135.2318	4302.4635	0.0000	1.0000
K2P	12	-2140.2881	4304.5761	0.0000	1.0000
F81G	15	-2162.4820	4354.9639	0.0000	1.0000
F81I	15	-2162.7635	4355.5270	0.0000	1.0000
F81IG	16	-2162.4892	4356.9785	0.0000	1.0000
F81	14	-2168.8996	4365.7992	0.0000	1.0000
JC69G	12	-2185.1955	4394.3911	0.0000	1.0000
JC69I	12	-2185.3994	4394.7988	0.0000	1.0000
JC69IG	13	-2185.2099	4396.4197	0.0000	1.0000
JC69	11	-2191.4067	4404.8133	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	16	-2106.2962	4317.5462	0.5762	0.5762
HKYI	16	-2106.6902	4318.3342	0.3885	0.9647
HKYIG	17	-2106.2974	4324.1082	0.0217	0.9864
HKY	15	-2113.3487	4325.0916	0.0132	0.9996
GTRG	20	-2100.9745	4333.1413	0.0002	0.9998
GTRI	20	-2101.3994	4333.9910	0.0002	1.0000
GTRIG	21	-2100.9745	4339.7010	0.0000	1.0000
GTR	19	-2108.6538	4341.9402	0.0000	1.0000
K2PG	13	-2133.5568	4352.3887	0.0000	1.0000
K2PI	13	-2133.7779	4352.8308	0.0000	1.0000
K2PIG	14	-2133.5779	4358.9905	0.0000	1.0000
K2P	12	-2140.2881	4359.2915	0.0000	1.0000
SYMG	17	-2128.5195	4368.5524	0.0000	1.0000
SYMI	17	-2128.7193	4368.9521	0.0000	1.0000
SYMIG	18	-2128.5381	4375.1492	0.0000	1.0000
SYM	16	-2135.2318	4375.4174	0.0000	1.0000
F81G	15	-2162.4820	4423.3581	0.0000	1.0000
F81I	15	-2162.7635	4423.9212	0.0000	1.0000
F81	14	-2168.8996	4429.6339	0.0000	1.0000
F81IG	16	-2162.4892	4429.9323	0.0000	1.0000
JC69G	12	-2185.1955	4449.1065	0.0000	1.0000
JC69I	12	-2185.3994	4449.5142	0.0000	1.0000
JC69	11	-2191.4067	4454.9691	0.0000	1.0000
JC69IG	13	-2185.2099	4455.6947	0.0000	1.0000



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

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

BEGIN MRBAYES;
 Charset 591-tomodel-phylip/chr7_10800.phylip = 1 - 706;
 Partition Dummy = 1:591-tomodel-phylip/chr7_10800.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=Dirichlet(1.0,1.0,1.0,1.0) shapepr=Uniform(0.1,50.0);
END;


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

BEGIN MRBAYES;
 Charset 591-tomodel-phylip/chr7_10800.phylip = 1 - 706;
 Partition Dummy = 1:591-tomodel-phylip/chr7_10800.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=Dirichlet(1.0,1.0,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/chr7_10800.phylip = 1 - 706;
 Partition Dummy = 1:591-tomodel-phylip/chr7_10800.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

