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

Input data from file "683-tomodel-phylip/chrZ_6110.phylip" (7 taxa, 1088 characters)

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


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

Model	df	lnL		AIC		wAIC	cwAIC
GTRI	20	-2825.6926	5691.3851	0.3630	0.3630
GTRIG	21	-2824.9894	5691.9787	0.2698	0.6327
GTRG	20	-2826.0835	5692.1669	0.2455	0.8782
HKYI	16	-2831.2530	5694.5059	0.0762	0.9545
HKYG	16	-2832.1064	5696.2128	0.0325	0.9870
HKYIG	17	-2832.0650	5698.1300	0.0125	0.9994
GTR	19	-2833.3071	5704.6142	0.0005	0.9999
HKY	15	-2838.9327	5707.8653	0.0001	1.0000
F81I	15	-2873.6513	5777.3027	0.0000	1.0000
F81G	15	-2873.9862	5777.9724	0.0000	1.0000
F81IG	16	-2873.6585	5779.3170	0.0000	1.0000
F81	14	-2879.5491	5787.0982	0.0000	1.0000
SYMI	17	-2930.0274	5894.0549	0.0000	1.0000
SYMG	17	-2930.3628	5894.7255	0.0000	1.0000
SYMIG	18	-2930.0357	5896.0714	0.0000	1.0000
K2PI	13	-2936.8490	5899.6981	0.0000	1.0000
K2PG	13	-2937.1204	5900.2407	0.0000	1.0000
K2PIG	14	-2936.8544	5901.7088	0.0000	1.0000
SYM	16	-2938.1791	5908.3583	0.0000	1.0000
K2P	12	-2942.8775	5909.7549	0.0000	1.0000
JC69I	12	-2973.9159	5971.8318	0.0000	1.0000
JC69G	12	-2974.1720	5972.3441	0.0000	1.0000
JC69IG	13	-2973.9218	5973.8437	0.0000	1.0000
JC69	11	-2982.1709	5986.3417	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYI	16	-2831.2530	5774.3794	0.6873	0.6873
HKYG	16	-2832.1064	5776.0863	0.2927	0.9800
HKY	15	-2838.9327	5782.7467	0.0105	0.9905
HKYIG	17	-2832.0650	5782.9957	0.0092	0.9997
GTRI	20	-2825.6926	5791.2270	0.0002	0.9999
GTRG	20	-2826.0835	5792.0089	0.0001	1.0000
GTRIG	21	-2824.9894	5796.8127	0.0000	1.0000
GTR	19	-2833.3071	5799.4641	0.0000	1.0000
F81I	15	-2873.6513	5852.1841	0.0000	1.0000
F81G	15	-2873.9862	5852.8539	0.0000	1.0000
F81	14	-2879.5491	5856.9875	0.0000	1.0000
F81IG	16	-2873.6585	5859.1905	0.0000	1.0000
K2PI	13	-2936.8490	5964.5953	0.0000	1.0000
K2PG	13	-2937.1204	5965.1380	0.0000	1.0000
K2P	12	-2942.8775	5969.6601	0.0000	1.0000
K2PIG	14	-2936.8544	5971.5981	0.0000	1.0000
SYMI	17	-2930.0274	5978.9205	0.0000	1.0000
SYMG	17	-2930.3628	5979.5912	0.0000	1.0000
SYMIG	18	-2930.0357	5985.9291	0.0000	1.0000
SYM	16	-2938.1791	5988.2318	0.0000	1.0000
JC69I	12	-2973.9159	6031.7369	0.0000	1.0000
JC69G	12	-2974.1720	6032.2492	0.0000	1.0000
JC69IG	13	-2973.9218	6038.7409	0.0000	1.0000
JC69	11	-2982.1709	6041.2548	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 683-tomodel-phylip/chrZ_6110.phylip = 1 - 1088;
 Partition Dummy = 1:683-tomodel-phylip/chrZ_6110.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 683-tomodel-phylip/chrZ_6110.phylip = 1 - 1088;
 Partition Dummy = 1:683-tomodel-phylip/chrZ_6110.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 683-tomodel-phylip/chrZ_6110.phylip = 1 - 1088;
 Partition Dummy = 1:683-tomodel-phylip/chrZ_6110.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

