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

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

Minimum AIC  model: GTRG
Minimum AICc model: GTRG
Minimum BIC  model: GTRG


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
GTRG	20	-2438.1168	4917.3550	0.6232	0.6232
GTRIG	21	-2438.1171	4919.4696	0.2165	0.8397
GTRI	20	-2439.4745	4920.0705	0.1603	1.0000
HKYG	16	-2453.8337	4940.3898	0.0000	1.0000
GTR	19	-2450.9114	4940.8362	0.0000	1.0000
HKYIG	17	-2453.8337	4942.4812	0.0000	1.0000
HKYI	16	-2455.0151	4942.7526	0.0000	1.0000
SYMG	17	-2459.3673	4953.5485	0.0000	1.0000
SYMI	17	-2460.0776	4954.9690	0.0000	1.0000
SYMIG	18	-2459.3695	4955.6498	0.0000	1.0000
HKY	15	-2466.0509	4962.7384	0.0000	1.0000
SYM	16	-2469.5075	4971.7374	0.0000	1.0000
K2PG	13	-2474.0073	4974.4962	0.0000	1.0000
K2PI	13	-2474.7504	4975.9823	0.0000	1.0000
K2PIG	14	-2474.0073	4976.5709	0.0000	1.0000
K2P	12	-2484.3243	4993.0608	0.0000	1.0000
F81G	15	-2512.1509	5054.9384	0.0000	1.0000
F81I	15	-2513.0158	5056.6681	0.0000	1.0000
F81IG	16	-2512.1565	5057.0355	0.0000	1.0000
F81	14	-2522.1181	5072.7924	0.0000	1.0000
JC69G	12	-2527.3644	5079.1409	0.0000	1.0000
JC69I	12	-2528.0513	5080.5147	0.0000	1.0000
JC69IG	13	-2527.3731	5081.2276	0.0000	1.0000
JC69	11	-2536.6764	5095.7010	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRG	20	-2438.1168	4916.2335	0.6154	0.6154
GTRIG	21	-2438.1171	4918.2343	0.2263	0.8417
GTRI	20	-2439.4745	4918.9490	0.1583	1.0000
HKYG	16	-2453.8337	4939.6674	0.0000	1.0000
GTR	19	-2450.9114	4939.8229	0.0000	1.0000
HKYIG	17	-2453.8337	4941.6674	0.0000	1.0000
HKYI	16	-2455.0151	4942.0301	0.0000	1.0000
SYMG	17	-2459.3673	4952.7346	0.0000	1.0000
SYMI	17	-2460.0776	4954.1552	0.0000	1.0000
SYMIG	18	-2459.3695	4954.7390	0.0000	1.0000
HKY	15	-2466.0509	4962.1018	0.0000	1.0000
SYM	16	-2469.5075	4971.0150	0.0000	1.0000
K2PG	13	-2474.0073	4974.0147	0.0000	1.0000
K2PI	13	-2474.7504	4975.5008	0.0000	1.0000
K2PIG	14	-2474.0073	4976.0146	0.0000	1.0000
K2P	12	-2484.3243	4992.6486	0.0000	1.0000
F81G	15	-2512.1509	5054.3018	0.0000	1.0000
F81I	15	-2513.0158	5056.0315	0.0000	1.0000
F81IG	16	-2512.1565	5056.3131	0.0000	1.0000
F81	14	-2522.1181	5072.2361	0.0000	1.0000
JC69G	12	-2527.3644	5078.7287	0.0000	1.0000
JC69I	12	-2528.0513	5080.1026	0.0000	1.0000
JC69IG	13	-2527.3731	5080.7461	0.0000	1.0000
JC69	11	-2536.6764	5095.3527	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
GTRG	20	-2438.1168	5009.1614	0.7081	0.7081
GTRI	20	-2439.4745	5011.8768	0.1821	0.8902
HKYG	16	-2453.8337	5014.0096	0.0627	0.9529
GTRIG	21	-2438.1171	5015.8085	0.0255	0.9784
HKYI	16	-2455.0151	5016.3724	0.0192	0.9977
HKYIG	17	-2453.8337	5020.6561	0.0023	0.9999
GTR	19	-2450.9114	5028.1043	0.0001	1.0000
SYMG	17	-2459.3673	5031.7233	0.0000	1.0000
HKY	15	-2466.0509	5031.7977	0.0000	1.0000
SYMI	17	-2460.0776	5033.1439	0.0000	1.0000
K2PG	13	-2474.0073	5034.4178	0.0000	1.0000
K2PI	13	-2474.7504	5035.9039	0.0000	1.0000
SYMIG	18	-2459.3695	5038.3740	0.0000	1.0000
K2PIG	14	-2474.0073	5041.0641	0.0000	1.0000
SYM	16	-2469.5075	5045.3572	0.0000	1.0000
K2P	12	-2484.3243	5048.4053	0.0000	1.0000
F81G	15	-2512.1509	5123.9976	0.0000	1.0000
F81I	15	-2513.0158	5125.7274	0.0000	1.0000
F81IG	16	-2512.1565	5130.6553	0.0000	1.0000
JC69G	12	-2527.3644	5134.4854	0.0000	1.0000
JC69I	12	-2528.0513	5135.8593	0.0000	1.0000
F81	14	-2522.1181	5137.2856	0.0000	1.0000
JC69IG	13	-2527.3731	5141.1492	0.0000	1.0000
JC69	11	-2536.6764	5146.4630	0.0000	1.0000



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

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

BEGIN MRBAYES;
 Charset 683-tomodel-phylip/chr2_24774.phylip = 1 - 770;
 Partition Dummy = 1:683-tomodel-phylip/chr2_24774.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=6 rates=gamma;
 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) shapepr=Uniform(0.1,50.0);
END;


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

BEGIN MRBAYES;
 Charset 683-tomodel-phylip/chr2_24774.phylip = 1 - 770;
 Partition Dummy = 1:683-tomodel-phylip/chr2_24774.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=6 rates=gamma;
 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) shapepr=Uniform(0.1,50.0);
END;


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

BEGIN MRBAYES;
 Charset 683-tomodel-phylip/chr2_24774.phylip = 1 - 770;
 Partition Dummy = 1:683-tomodel-phylip/chr2_24774.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=6 rates=gamma;
 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) shapepr=Uniform(0.1,50.0);
END;


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

