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

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

Minimum AIC  model: GTRI
Minimum AICc model: HKYI
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
HKYI	16	-2293.9689	4620.7041	0.2458	0.2458
GTRI	20	-2289.9188	4621.0274	0.2091	0.4548
GTRG	20	-2289.9927	4621.1752	0.1942	0.6490
HKYG	16	-2294.3068	4621.3798	0.1753	0.8243
HKYIG	17	-2293.9245	4622.7121	0.0900	0.9144
GTRIG	21	-2289.7510	4622.8126	0.0856	1.0000
HKY	15	-2324.9437	4680.5624	0.0000	1.0000
GTR	19	-2321.1440	4681.3630	0.0000	1.0000
SYMI	17	-2325.4488	4685.7608	0.0000	1.0000
SYMG	17	-2325.5274	4685.9180	0.0000	1.0000
SYMIG	18	-2325.6628	4688.2917	0.0000	1.0000
K2PI	13	-2337.1972	4700.9050	0.0000	1.0000
K2PG	13	-2337.5783	4701.6670	0.0000	1.0000
K2PIG	14	-2337.1931	4702.9762	0.0000	1.0000
F81I	15	-2350.9409	4732.5570	0.0000	1.0000
F81G	15	-2351.2605	4733.1961	0.0000	1.0000
F81IG	16	-2350.8934	4734.5530	0.0000	1.0000
SYM	16	-2355.7756	4744.3174	0.0000	1.0000
K2P	12	-2367.3421	4759.1212	0.0000	1.0000
F81	14	-2383.2679	4795.1257	0.0000	1.0000
JC69I	12	-2394.1651	4812.7672	0.0000	1.0000
JC69G	12	-2394.6375	4813.7120	0.0000	1.0000
JC69IG	13	-2394.1667	4814.8440	0.0000	1.0000
JC69	11	-2425.2516	4872.8724	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRI	20	-2289.9188	4619.8376	0.2294	0.2294
HKYI	16	-2293.9689	4619.9379	0.2182	0.4476
GTRG	20	-2289.9927	4619.9854	0.2131	0.6606
HKYG	16	-2294.3068	4620.6136	0.1556	0.8163
GTRIG	21	-2289.7510	4621.5020	0.0998	0.9161
HKYIG	17	-2293.9245	4621.8490	0.0839	1.0000
HKY	15	-2324.9437	4679.8873	0.0000	1.0000
GTR	19	-2321.1440	4680.2880	0.0000	1.0000
SYMI	17	-2325.4488	4684.8976	0.0000	1.0000
SYMG	17	-2325.5274	4685.0548	0.0000	1.0000
SYMIG	18	-2325.6628	4687.3256	0.0000	1.0000
K2PI	13	-2337.1972	4700.3944	0.0000	1.0000
K2PG	13	-2337.5783	4701.1565	0.0000	1.0000
K2PIG	14	-2337.1931	4702.3863	0.0000	1.0000
F81I	15	-2350.9409	4731.8819	0.0000	1.0000
F81G	15	-2351.2605	4732.5210	0.0000	1.0000
F81IG	16	-2350.8934	4733.7868	0.0000	1.0000
SYM	16	-2355.7756	4743.5512	0.0000	1.0000
K2P	12	-2367.3421	4758.6842	0.0000	1.0000
F81	14	-2383.2679	4794.5358	0.0000	1.0000
JC69I	12	-2394.1651	4812.3302	0.0000	1.0000
JC69G	12	-2394.6375	4813.2750	0.0000	1.0000
JC69IG	13	-2394.1667	4814.3335	0.0000	1.0000
JC69	11	-2425.2516	4872.5032	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYI	16	-2293.9689	4693.3607	0.5707	0.5707
HKYG	16	-2294.3068	4694.0364	0.4071	0.9777
HKYIG	17	-2293.9245	4699.8607	0.0221	0.9999
GTRI	20	-2289.9188	4711.6161	0.0001	0.9999
GTRG	20	-2289.9927	4711.7639	0.0001	1.0000
GTRIG	21	-2289.7510	4717.8694	0.0000	1.0000
HKY	15	-2324.9437	4748.7212	0.0000	1.0000
K2PI	13	-2337.1972	4760.0505	0.0000	1.0000
K2PG	13	-2337.5783	4760.8126	0.0000	1.0000
SYMI	17	-2325.4488	4762.9094	0.0000	1.0000
SYMG	17	-2325.5274	4763.0666	0.0000	1.0000
K2PIG	14	-2337.1931	4766.6313	0.0000	1.0000
GTR	19	-2321.1440	4767.4776	0.0000	1.0000
SYMIG	18	-2325.6628	4769.9263	0.0000	1.0000
F81I	15	-2350.9409	4800.7158	0.0000	1.0000
F81G	15	-2351.2605	4801.3549	0.0000	1.0000
F81IG	16	-2350.8934	4807.2096	0.0000	1.0000
K2P	12	-2367.3421	4813.7514	0.0000	1.0000
SYM	16	-2355.7756	4816.9740	0.0000	1.0000
F81	14	-2383.2679	4858.7808	0.0000	1.0000
JC69I	12	-2394.1651	4867.3973	0.0000	1.0000
JC69G	12	-2394.6375	4868.3421	0.0000	1.0000
JC69IG	13	-2394.1667	4873.9895	0.0000	1.0000
JC69	11	-2425.2516	4922.9814	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/chr7_4965.phylip = 1 - 727;
 Partition Dummy = 1:591-tomodel-phylip/chr7_4965.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 (HKYI)]

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


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

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

