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

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

Minimum AIC  model: GTRI
Minimum AICc model: GTRI
Minimum BIC  model: HKYI


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	-1869.1266	3779.3808	0.3804	0.3804
GTRG	20	-1869.7158	3780.5592	0.2110	0.5914
GTRIG	21	-1869.1447	3781.5314	0.1298	0.7212
HKYI	16	-1874.5675	3781.8612	0.1101	0.8312
HKYG	16	-1875.1449	3783.0160	0.0618	0.8930
GTR	19	-1872.0288	3783.0763	0.0599	0.9530
HKYIG	17	-1874.5841	3783.9863	0.0380	0.9910
HKY	15	-1878.1126	3786.8651	0.0090	1.0000
K2PI	13	-1902.3225	3831.1291	0.0000	1.0000
K2PG	13	-1902.8000	3832.0841	0.0000	1.0000
K2PIG	14	-1902.3361	3833.2315	0.0000	1.0000
SYMI	17	-1900.4348	3835.6878	0.0000	1.0000
SYMG	17	-1900.9037	3836.6256	0.0000	1.0000
SYMIG	18	-1900.4485	3837.8126	0.0000	1.0000
K2P	12	-1907.0336	3838.4816	0.0000	1.0000
F81I	15	-1905.1843	3841.0086	0.0000	1.0000
SYM	16	-1904.4191	3841.5644	0.0000	1.0000
F81G	15	-1905.7332	3842.1064	0.0000	1.0000
F81IG	16	-1905.2002	3843.1267	0.0000	1.0000
F81	14	-1908.7685	3846.0962	0.0000	1.0000
JC69I	12	-1930.4772	3885.3688	0.0000	1.0000
JC69G	12	-1930.9666	3886.3476	0.0000	1.0000
JC69IG	13	-1930.4918	3887.4677	0.0000	1.0000
JC69	11	-1934.7516	3891.8533	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRI	20	-1869.1266	3778.2533	0.3938	0.3938
GTRG	20	-1869.7158	3779.4316	0.2185	0.6123
GTRIG	21	-1869.1447	3780.2895	0.1423	0.7546
HKYI	16	-1874.5675	3781.1349	0.0932	0.8478
GTR	19	-1872.0288	3782.0575	0.0588	0.9066
HKYG	16	-1875.1449	3782.2897	0.0523	0.9590
HKYIG	17	-1874.5841	3783.1682	0.0337	0.9927
HKY	15	-1878.1126	3786.2251	0.0073	1.0000
K2PI	13	-1902.3225	3830.6450	0.0000	1.0000
K2PG	13	-1902.8000	3831.6000	0.0000	1.0000
K2PIG	14	-1902.3361	3832.6722	0.0000	1.0000
SYMI	17	-1900.4348	3834.8696	0.0000	1.0000
SYMG	17	-1900.9037	3835.8075	0.0000	1.0000
SYMIG	18	-1900.4485	3836.8970	0.0000	1.0000
K2P	12	-1907.0336	3838.0673	0.0000	1.0000
F81I	15	-1905.1843	3840.3686	0.0000	1.0000
SYM	16	-1904.4191	3840.8381	0.0000	1.0000
F81G	15	-1905.7332	3841.4664	0.0000	1.0000
F81IG	16	-1905.2002	3842.4004	0.0000	1.0000
F81	14	-1908.7685	3845.5370	0.0000	1.0000
JC69I	12	-1930.4772	3884.9545	0.0000	1.0000
JC69G	12	-1930.9666	3885.9333	0.0000	1.0000
JC69IG	13	-1930.4918	3886.9836	0.0000	1.0000
JC69	11	-1934.7516	3891.5032	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYI	16	-1874.5675	3855.3938	0.4172	0.4172
HKY	15	-1878.1126	3855.8429	0.3333	0.7505
HKYG	16	-1875.1449	3856.5486	0.2342	0.9847
HKYIG	17	-1874.5841	3862.0683	0.0148	0.9995
GTR	19	-1872.0288	3870.2400	0.0002	0.9997
GTRI	20	-1869.1266	3871.0769	0.0002	0.9999
GTRG	20	-1869.7158	3872.2553	0.0001	1.0000
GTRIG	21	-1869.1447	3877.7543	0.0000	1.0000
K2PI	13	-1902.3225	3890.9804	0.0000	1.0000
K2PG	13	-1902.8000	3891.9354	0.0000	1.0000
K2P	12	-1907.0336	3893.7615	0.0000	1.0000
K2PIG	14	-1902.3361	3897.6488	0.0000	1.0000
F81I	15	-1905.1843	3909.9863	0.0000	1.0000
F81	14	-1908.7685	3910.5135	0.0000	1.0000
F81G	15	-1905.7332	3911.0841	0.0000	1.0000
SYMI	17	-1900.4348	3913.7697	0.0000	1.0000
SYMG	17	-1900.9037	3914.7076	0.0000	1.0000
SYM	16	-1904.4191	3915.0970	0.0000	1.0000
F81IG	16	-1905.2002	3916.6593	0.0000	1.0000
SYMIG	18	-1900.4485	3920.4382	0.0000	1.0000
JC69I	12	-1930.4772	3940.6487	0.0000	1.0000
JC69G	12	-1930.9666	3941.6274	0.0000	1.0000
JC69	11	-1934.7516	3942.5562	0.0000	1.0000
JC69IG	13	-1930.4918	3947.3190	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/chr4_4009.phylip = 1 - 766;
 Partition Dummy = 1:591-tomodel-phylip/chr4_4009.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/chr4_4009.phylip = 1 - 766;
 Partition Dummy = 1:591-tomodel-phylip/chr4_4009.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 (HKYI)]

BEGIN MRBAYES;
 Charset 591-tomodel-phylip/chr4_4009.phylip = 1 - 766;
 Partition Dummy = 1:591-tomodel-phylip/chr4_4009.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=2 rates=propinv;
 Prset applyto=(1) statefreqpr=Dirichlet(1.0,1.0,1.0,1.0) tratiopr=Beta(1.0,1.0) pinvarpr=Uniform(0.0,1.0);
END;


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

