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

Input data from file "591-tomodel-phylip/chr2_24728.phylip" (7 taxa, 655 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	-2804.8266	5650.9782	0.5282	0.5282
GTRIG	21	-2804.5839	5652.6275	0.2316	0.7598
GTRG	20	-2806.2043	5653.7334	0.1332	0.8930
HKYI	16	-2811.0118	5654.8764	0.0752	0.9682
HKYIG	17	-2811.0050	5656.9707	0.0264	0.9946
HKYG	16	-2813.6432	5660.1390	0.0054	1.0000
GTR	19	-2832.4883	5704.1734	0.0000	1.0000
HKY	15	-2839.5021	5709.7553	0.0000	1.0000
K2PI	13	-2852.7385	5732.0448	0.0000	1.0000
SYMI	17	-2848.7905	5732.5418	0.0000	1.0000
SYMG	17	-2849.0482	5733.0572	0.0000	1.0000
SYMIG	18	-2848.3268	5733.7291	0.0000	1.0000
K2PIG	14	-2852.6259	5733.9081	0.0000	1.0000
K2PG	13	-2854.0765	5734.7208	0.0000	1.0000
F81I	15	-2874.0919	5778.9349	0.0000	1.0000
F81IG	16	-2874.1133	5781.0793	0.0000	1.0000
SYM	16	-2874.5384	5781.9294	0.0000	1.0000
K2P	12	-2879.2385	5782.9629	0.0000	1.0000
F81G	15	-2877.0427	5784.8365	0.0000	1.0000
F81	14	-2899.5926	5827.8415	0.0000	1.0000
JC69I	12	-2905.1250	5834.7359	0.0000	1.0000
JC69IG	13	-2905.1320	5836.8318	0.0000	1.0000
JC69G	12	-2907.3830	5839.2520	0.0000	1.0000
JC69	11	-2929.7968	5882.0041	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRI	20	-2804.8266	5649.6532	0.5310	0.5310
GTRIG	21	-2804.5839	5651.1678	0.2490	0.7800
GTRG	20	-2806.2043	5652.4085	0.1339	0.9139
HKYI	16	-2811.0118	5654.0237	0.0597	0.9736
HKYIG	17	-2811.0050	5656.0099	0.0221	0.9957
HKYG	16	-2813.6432	5659.2864	0.0043	1.0000
GTR	19	-2832.4883	5702.9766	0.0000	1.0000
HKY	15	-2839.5021	5709.0042	0.0000	1.0000
K2PI	13	-2852.7385	5731.4770	0.0000	1.0000
SYMI	17	-2848.7905	5731.5810	0.0000	1.0000
SYMG	17	-2849.0482	5732.0965	0.0000	1.0000
SYMIG	18	-2848.3268	5732.6537	0.0000	1.0000
K2PIG	14	-2852.6259	5733.2518	0.0000	1.0000
K2PG	13	-2854.0765	5734.1529	0.0000	1.0000
F81I	15	-2874.0919	5778.1837	0.0000	1.0000
F81IG	16	-2874.1133	5780.2266	0.0000	1.0000
SYM	16	-2874.5384	5781.0767	0.0000	1.0000
K2P	12	-2879.2385	5782.4769	0.0000	1.0000
F81G	15	-2877.0427	5784.0853	0.0000	1.0000
F81	14	-2899.5926	5827.1853	0.0000	1.0000
JC69I	12	-2905.1250	5834.2499	0.0000	1.0000
JC69IG	13	-2905.1320	5836.2640	0.0000	1.0000
JC69G	12	-2907.3830	5838.7660	0.0000	1.0000
JC69	11	-2929.7968	5881.5935	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYI	16	-2811.0118	5725.7779	0.8986	0.8986
HKYG	16	-2813.6432	5731.0405	0.0647	0.9633
HKYIG	17	-2811.0050	5732.2487	0.0354	0.9987
GTRI	20	-2804.8266	5739.3459	0.0010	0.9997
GTRG	20	-2806.2043	5742.1012	0.0003	0.9999
GTRIG	21	-2804.5839	5745.3451	0.0001	1.0000
HKY	15	-2839.5021	5776.2737	0.0000	1.0000
GTR	19	-2832.4883	5788.1846	0.0000	1.0000
K2PI	13	-2852.7385	5789.7772	0.0000	1.0000
K2PG	13	-2854.0765	5792.4532	0.0000	1.0000
K2PIG	14	-2852.6259	5796.0367	0.0000	1.0000
SYMI	17	-2848.7905	5807.8198	0.0000	1.0000
SYMG	17	-2849.0482	5808.3353	0.0000	1.0000
SYMIG	18	-2848.3268	5813.3771	0.0000	1.0000
K2P	12	-2879.2385	5836.2926	0.0000	1.0000
F81I	15	-2874.0919	5845.4532	0.0000	1.0000
F81G	15	-2877.0427	5851.3548	0.0000	1.0000
F81IG	16	-2874.1133	5851.9808	0.0000	1.0000
SYM	16	-2874.5384	5852.8309	0.0000	1.0000
JC69I	12	-2905.1250	5888.0655	0.0000	1.0000
F81	14	-2899.5926	5889.9702	0.0000	1.0000
JC69G	12	-2907.3830	5892.5817	0.0000	1.0000
JC69IG	13	-2905.1320	5894.5642	0.0000	1.0000
JC69	11	-2929.7968	5930.9245	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/chr2_24728.phylip = 1 - 655;
 Partition Dummy = 1:591-tomodel-phylip/chr2_24728.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/chr2_24728.phylip = 1 - 655;
 Partition Dummy = 1:591-tomodel-phylip/chr2_24728.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/chr2_24728.phylip = 1 - 655;
 Partition Dummy = 1:591-tomodel-phylip/chr2_24728.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

