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

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

Minimum AIC  model: GTRI
Minimum AICc model: GTRI
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
GTRI	20	-1847.9154	3737.2691	0.4175	0.4175
GTRG	20	-1847.9879	3737.4142	0.3883	0.8059
GTRIG	21	-1847.9031	3739.3912	0.1445	0.9504
HKYG	16	-1854.8832	3742.6916	0.0277	0.9781
HKYI	16	-1855.7368	3744.3988	0.0118	0.9899
HKYIG	17	-1854.8398	3744.7222	0.0101	1.0000
SYMG	17	-1877.4671	3789.9767	0.0000	1.0000
GTR	19	-1875.4265	3790.1521	0.0000	1.0000
SYMI	17	-1877.7277	3790.4979	0.0000	1.0000
SYMIG	18	-1877.5033	3792.1738	0.0000	1.0000
HKY	15	-1882.5993	3796.0136	0.0000	1.0000
K2PG	13	-1886.9747	3800.5652	0.0000	1.0000
K2PI	13	-1887.7571	3802.1301	0.0000	1.0000
K2PIG	14	-1886.9804	3802.6727	0.0000	1.0000
SYM	16	-1901.7321	3836.3894	0.0000	1.0000
F81G	15	-1905.8685	3842.5519	0.0000	1.0000
F81I	15	-1906.6442	3844.1034	0.0000	1.0000
F81IG	16	-1905.8350	3844.5952	0.0000	1.0000
K2P	12	-1911.7659	3848.0588	0.0000	1.0000
JC69G	12	-1932.3100	3889.1470	0.0000	1.0000
F81	14	-1930.6824	3890.0766	0.0000	1.0000
JC69I	12	-1933.0268	3890.5805	0.0000	1.0000
JC69IG	13	-1932.3164	3891.2487	0.0000	1.0000
JC69	11	-1955.1249	3932.6950	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRI	20	-1847.9154	3735.8307	0.4174	0.4174
GTRG	20	-1847.9879	3735.9758	0.3882	0.8057
GTRIG	21	-1847.9031	3737.8063	0.1555	0.9611
HKYG	16	-1854.8832	3741.7664	0.0215	0.9826
HKYI	16	-1855.7368	3743.4736	0.0091	0.9918
HKYIG	17	-1854.8398	3743.6796	0.0082	1.0000
GTR	19	-1875.4265	3788.8529	0.0000	1.0000
SYMG	17	-1877.4671	3788.9341	0.0000	1.0000
SYMI	17	-1877.7277	3789.4553	0.0000	1.0000
SYMIG	18	-1877.5033	3791.0065	0.0000	1.0000
HKY	15	-1882.5993	3795.1987	0.0000	1.0000
K2PG	13	-1886.9747	3799.9493	0.0000	1.0000
K2PI	13	-1887.7571	3801.5142	0.0000	1.0000
K2PIG	14	-1886.9804	3801.9609	0.0000	1.0000
SYM	16	-1901.7321	3835.4642	0.0000	1.0000
F81G	15	-1905.8685	3841.7369	0.0000	1.0000
F81I	15	-1906.6442	3843.2885	0.0000	1.0000
F81IG	16	-1905.8350	3843.6701	0.0000	1.0000
K2P	12	-1911.7659	3847.5318	0.0000	1.0000
JC69G	12	-1932.3100	3888.6200	0.0000	1.0000
F81	14	-1930.6824	3889.3648	0.0000	1.0000
JC69I	12	-1933.0268	3890.0535	0.0000	1.0000
JC69IG	13	-1932.3164	3890.6328	0.0000	1.0000
JC69	11	-1955.1249	3932.2498	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	16	-1854.8832	3812.2500	0.6784	0.6784
HKYI	16	-1855.7368	3813.9573	0.2889	0.9673
HKYIG	17	-1854.8398	3818.5685	0.0288	0.9961
GTRI	20	-1847.9154	3823.9353	0.0020	0.9981
GTRG	20	-1847.9879	3824.0804	0.0018	0.9999
GTRIG	21	-1847.9031	3830.3161	0.0001	1.0000
K2PG	13	-1886.9747	3857.2173	0.0000	1.0000
K2PI	13	-1887.7571	3858.7821	0.0000	1.0000
HKY	15	-1882.5993	3861.2771	0.0000	1.0000
K2PIG	14	-1886.9804	3863.6341	0.0000	1.0000
SYMG	17	-1877.4671	3863.8230	0.0000	1.0000
SYMI	17	-1877.7277	3864.3442	0.0000	1.0000
SYMIG	18	-1877.5033	3870.3006	0.0000	1.0000
GTR	19	-1875.4265	3872.5523	0.0000	1.0000
K2P	12	-1911.7659	3900.3945	0.0000	1.0000
SYM	16	-1901.7321	3905.9479	0.0000	1.0000
F81G	15	-1905.8685	3907.8154	0.0000	1.0000
F81I	15	-1906.6442	3909.3669	0.0000	1.0000
F81IG	16	-1905.8350	3914.1537	0.0000	1.0000
JC69G	12	-1932.3100	3941.4828	0.0000	1.0000
JC69I	12	-1933.0268	3942.9163	0.0000	1.0000
JC69IG	13	-1932.3164	3947.9008	0.0000	1.0000
F81	14	-1930.6824	3951.0380	0.0000	1.0000
JC69	11	-1955.1249	3980.7073	0.0000	1.0000



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

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

BEGIN MRBAYES;
 Charset 591-tomodel-phylip/chr7_10530.phylip = 1 - 605;
 Partition Dummy = 1:591-tomodel-phylip/chr7_10530.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=6 rates=propinv;
 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) pinvarpr=Uniform(0.0,1.0);
END;


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

BEGIN MRBAYES;
 Charset 591-tomodel-phylip/chr7_10530.phylip = 1 - 605;
 Partition Dummy = 1:591-tomodel-phylip/chr7_10530.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=6 rates=propinv;
 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) pinvarpr=Uniform(0.0,1.0);
END;


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

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

