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

Input data from file "683-tomodel-phylip/chrZ_11394.phylip" (7 taxa, 629 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	-2178.2647	4397.9110	0.3362	0.3362
HKYI	16	-2182.6597	4398.2082	0.2898	0.6260
GTRIG	21	-2178.2996	4400.1215	0.1113	0.7374
HKYIG	17	-2182.6947	4400.3911	0.0973	0.8347
GTRG	20	-2179.6794	4400.7404	0.0817	0.9164
HKYG	16	-2183.9506	4400.7901	0.0797	0.9961
HKY	15	-2188.3962	4407.5755	0.0027	0.9987
GTR	19	-2184.9228	4409.0935	0.0013	1.0000
K2PI	13	-2237.1137	4500.8193	0.0000	1.0000
K2PIG	14	-2237.1476	4502.9792	0.0000	1.0000
K2PG	13	-2238.3189	4503.2296	0.0000	1.0000
SYMI	17	-2235.4218	4505.8453	0.0000	1.0000
SYMIG	18	-2235.4578	4508.0369	0.0000	1.0000
SYMG	17	-2236.6779	4508.3574	0.0000	1.0000
K2P	12	-2242.6838	4509.8740	0.0000	1.0000
F81I	15	-2239.6977	4510.1784	0.0000	1.0000
F81IG	16	-2239.7316	4512.3521	0.0000	1.0000
F81G	15	-2240.8092	4512.4015	0.0000	1.0000
SYM	16	-2240.9909	4514.8708	0.0000	1.0000
F81	14	-2244.3275	4517.3391	0.0000	1.0000
JC69I	12	-2281.9472	4588.4009	0.0000	1.0000
JC69IG	13	-2281.9824	4590.5566	0.0000	1.0000
JC69G	12	-2283.0658	4590.6380	0.0000	1.0000
JC69	11	-2286.9386	4596.3051	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRI	20	-2178.2647	4396.5294	0.3696	0.3696
HKYI	16	-2182.6597	4397.3193	0.2490	0.6185
GTRIG	21	-2178.2996	4398.5992	0.1313	0.7498
GTRG	20	-2179.6794	4399.3588	0.0898	0.8396
HKYIG	17	-2182.6947	4399.3894	0.0884	0.9281
HKYG	16	-2183.9506	4399.9012	0.0685	0.9965
HKY	15	-2188.3962	4406.7925	0.0022	0.9987
GTR	19	-2184.9228	4407.8456	0.0013	1.0000
K2PI	13	-2237.1137	4500.2275	0.0000	1.0000
K2PIG	14	-2237.1476	4502.2951	0.0000	1.0000
K2PG	13	-2238.3189	4502.6377	0.0000	1.0000
SYMI	17	-2235.4218	4504.8436	0.0000	1.0000
SYMIG	18	-2235.4578	4506.9155	0.0000	1.0000
SYMG	17	-2236.6779	4507.3557	0.0000	1.0000
K2P	12	-2242.6838	4509.3675	0.0000	1.0000
F81I	15	-2239.6977	4509.3954	0.0000	1.0000
F81IG	16	-2239.7316	4511.4632	0.0000	1.0000
F81G	15	-2240.8092	4511.6184	0.0000	1.0000
SYM	16	-2240.9909	4513.9819	0.0000	1.0000
F81	14	-2244.3275	4516.6551	0.0000	1.0000
JC69I	12	-2281.9472	4587.8944	0.0000	1.0000
JC69IG	13	-2281.9824	4589.9648	0.0000	1.0000
JC69G	12	-2283.0658	4590.1315	0.0000	1.0000
JC69	11	-2286.9386	4595.8772	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYI	16	-2182.6597	4468.4254	0.7170	0.7170
HKYG	16	-2183.9506	4471.0073	0.1972	0.9142
HKY	15	-2188.3962	4473.4544	0.0580	0.9722
HKYIG	17	-2182.6947	4474.9397	0.0276	0.9998
GTRI	20	-2178.2647	4485.4120	0.0001	1.0000
GTRG	20	-2179.6794	4488.2414	0.0000	1.0000
GTRIG	21	-2178.2996	4491.9260	0.0000	1.0000
GTR	19	-2184.9228	4492.2841	0.0000	1.0000
K2PI	13	-2237.1137	4558.0012	0.0000	1.0000
K2PG	13	-2238.3189	4560.4114	0.0000	1.0000
K2P	12	-2242.6838	4562.6971	0.0000	1.0000
K2PIG	14	-2237.1476	4564.5130	0.0000	1.0000
F81I	15	-2239.6977	4576.0574	0.0000	1.0000
F81G	15	-2240.8092	4578.2804	0.0000	1.0000
F81	14	-2244.3275	4578.8729	0.0000	1.0000
SYMI	17	-2235.4218	4580.3939	0.0000	1.0000
F81IG	16	-2239.7316	4582.5693	0.0000	1.0000
SYMG	17	-2236.6779	4582.9060	0.0000	1.0000
SYM	16	-2240.9909	4585.0880	0.0000	1.0000
SYMIG	18	-2235.4578	4586.9099	0.0000	1.0000
JC69I	12	-2281.9472	4641.2240	0.0000	1.0000
JC69G	12	-2283.0658	4643.4611	0.0000	1.0000
JC69	11	-2286.9386	4644.7627	0.0000	1.0000
JC69IG	13	-2281.9824	4647.7385	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_11394.phylip = 1 - 629;
 Partition Dummy = 1:683-tomodel-phylip/chrZ_11394.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_11394.phylip = 1 - 629;
 Partition Dummy = 1:683-tomodel-phylip/chrZ_11394.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_11394.phylip = 1 - 629;
 Partition Dummy = 1:683-tomodel-phylip/chrZ_11394.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

