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

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

Minimum AIC  model: HKYG
Minimum AICc model: HKYG
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
HKYG	16	-2840.4500	5713.5791	0.5010	0.5010
HKYI	16	-2841.3482	5715.3755	0.2041	0.7051
HKYIG	17	-2840.4128	5715.5907	0.1832	0.8883
GTRG	20	-2838.3560	5717.7659	0.0618	0.9501
GTRI	20	-2839.1927	5719.4394	0.0267	0.9768
GTRIG	21	-2838.2818	5719.7244	0.0232	1.0000
K2PG	13	-2856.5318	5739.5164	0.0000	1.0000
K2PI	13	-2857.2669	5740.9866	0.0000	1.0000
K2PIG	14	-2856.4960	5741.5151	0.0000	1.0000
SYMG	17	-2855.6953	5746.1557	0.0000	1.0000
SYMI	17	-2856.3968	5747.5586	0.0000	1.0000
SYMIG	18	-2855.6631	5748.1822	0.0000	1.0000
HKY	15	-2868.1350	5766.8686	0.0000	1.0000
GTR	19	-2865.3164	5769.5851	0.0000	1.0000
K2P	12	-2884.5076	5793.4027	0.0000	1.0000
SYM	16	-2882.9665	5798.6121	0.0000	1.0000
F81G	15	-2891.3380	5813.2745	0.0000	1.0000
F81I	15	-2891.6309	5813.8602	0.0000	1.0000
F81IG	16	-2891.1840	5815.0471	0.0000	1.0000
JC69G	12	-2904.0708	5832.5291	0.0000	1.0000
JC69I	12	-2904.3337	5833.0551	0.0000	1.0000
JC69IG	13	-2903.9177	5834.2881	0.0000	1.0000
F81	14	-2917.1825	5862.8881	0.0000	1.0000
JC69	11	-2930.3518	5883.0311	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
HKYG	16	-2840.4500	5712.9000	0.4852	0.4852
HKYI	16	-2841.3482	5714.6963	0.1976	0.6828
HKYIG	17	-2840.4128	5714.8257	0.1852	0.8680
GTRG	20	-2838.3560	5716.7119	0.0721	0.9402
GTRI	20	-2839.1927	5718.3855	0.0312	0.9714
GTRIG	21	-2838.2818	5718.5636	0.0286	1.0000
K2PG	13	-2856.5318	5739.0637	0.0000	1.0000
K2PI	13	-2857.2669	5740.5339	0.0000	1.0000
K2PIG	14	-2856.4960	5740.9921	0.0000	1.0000
SYMG	17	-2855.6953	5745.3907	0.0000	1.0000
SYMI	17	-2856.3968	5746.7936	0.0000	1.0000
SYMIG	18	-2855.6631	5747.3261	0.0000	1.0000
HKY	15	-2868.1350	5766.2701	0.0000	1.0000
GTR	19	-2865.3164	5768.6327	0.0000	1.0000
K2P	12	-2884.5076	5793.0151	0.0000	1.0000
SYM	16	-2882.9665	5797.9329	0.0000	1.0000
F81G	15	-2891.3380	5812.6760	0.0000	1.0000
F81I	15	-2891.6309	5813.2617	0.0000	1.0000
F81IG	16	-2891.1840	5814.3680	0.0000	1.0000
JC69G	12	-2904.0708	5832.1415	0.0000	1.0000
JC69I	12	-2904.3337	5832.6675	0.0000	1.0000
JC69IG	13	-2903.9177	5833.8353	0.0000	1.0000
F81	14	-2917.1825	5862.3650	0.0000	1.0000
JC69	11	-2930.3518	5882.7035	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	16	-2840.4500	5788.2098	0.6909	0.6909
HKYI	16	-2841.3482	5790.0061	0.2814	0.9724
HKYIG	17	-2840.4128	5794.8423	0.0251	0.9974
K2PG	13	-2856.5318	5800.2529	0.0017	0.9991
K2PI	13	-2857.2669	5801.7231	0.0008	0.9999
K2PIG	14	-2856.4960	5806.8881	0.0001	1.0000
GTRG	20	-2838.3560	5810.8492	0.0000	1.0000
GTRI	20	-2839.1927	5812.5227	0.0000	1.0000
GTRIG	21	-2838.2818	5817.4077	0.0000	1.0000
SYMG	17	-2855.6953	5825.4073	0.0000	1.0000
SYMI	17	-2856.3968	5826.8103	0.0000	1.0000
SYMIG	18	-2855.6631	5832.0496	0.0000	1.0000
HKY	15	-2868.1350	5836.8730	0.0000	1.0000
K2P	12	-2884.5076	5849.4975	0.0000	1.0000
GTR	19	-2865.3164	5858.0631	0.0000	1.0000
SYM	16	-2882.9665	5873.2427	0.0000	1.0000
F81G	15	-2891.3380	5883.2789	0.0000	1.0000
F81I	15	-2891.6309	5883.8646	0.0000	1.0000
JC69G	12	-2904.0708	5888.6239	0.0000	1.0000
JC69I	12	-2904.3337	5889.1498	0.0000	1.0000
F81IG	16	-2891.1840	5889.6778	0.0000	1.0000
JC69IG	13	-2903.9177	5895.0246	0.0000	1.0000
F81	14	-2917.1825	5928.2611	0.0000	1.0000
JC69	11	-2930.3518	5934.4790	0.0000	1.0000



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

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

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


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

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


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

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

