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

Input data from file "591-tomodel-phylip/chr6_7961.phylip" (7 taxa, 671 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	-2375.4917	4792.2756	0.5854	0.5854
GTRIG	21	-2375.4799	4794.3835	0.2040	0.7894
GTRG	20	-2376.5574	4794.4071	0.2017	0.9911
HKYI	16	-2384.5872	4802.0062	0.0045	0.9956
HKYG	16	-2385.0880	4803.0078	0.0027	0.9983
HKYIG	17	-2384.5383	4804.0138	0.0017	1.0000
K2PI	13	-2410.9948	4848.5436	0.0000	1.0000
K2PG	13	-2411.3226	4849.1992	0.0000	1.0000
K2PIG	14	-2410.9216	4850.4834	0.0000	1.0000
SYMI	17	-2408.1625	4851.2622	0.0000	1.0000
SYMG	17	-2408.9541	4852.8455	0.0000	1.0000
SYMIG	18	-2408.1636	4853.3762	0.0000	1.0000
GTR	19	-2407.3717	4853.9108	0.0000	1.0000
HKY	15	-2416.8287	4864.3902	0.0000	1.0000
K2P	12	-2445.5672	4915.6085	0.0000	1.0000
SYM	16	-2441.8741	4916.5800	0.0000	1.0000
F81I	15	-2445.4871	4921.7071	0.0000	1.0000
F81IG	16	-2445.5028	4923.8375	0.0000	1.0000
F81G	15	-2446.8033	4924.3395	0.0000	1.0000
JC69I	12	-2466.6778	4957.8297	0.0000	1.0000
JC69G	12	-2467.7064	4959.8870	0.0000	1.0000
JC69IG	13	-2466.6873	4959.9286	0.0000	1.0000
F81	14	-2476.8734	4982.3870	0.0000	1.0000
JC69	11	-2499.2179	5020.8364	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRI	20	-2375.4917	4790.9833	0.5784	0.5784
GTRIG	21	-2375.4799	4792.9597	0.2153	0.7937
GTRG	20	-2376.5574	4793.1148	0.1992	0.9929
HKYI	16	-2384.5872	4801.1744	0.0035	0.9965
HKYG	16	-2385.0880	4802.1760	0.0021	0.9986
HKYIG	17	-2384.5383	4803.0766	0.0014	1.0000
K2PI	13	-2410.9948	4847.9895	0.0000	1.0000
K2PG	13	-2411.3226	4848.6452	0.0000	1.0000
K2PIG	14	-2410.9216	4849.8431	0.0000	1.0000
SYMI	17	-2408.1625	4850.3250	0.0000	1.0000
SYMG	17	-2408.9541	4851.9083	0.0000	1.0000
SYMIG	18	-2408.1636	4852.3272	0.0000	1.0000
GTR	19	-2407.3717	4852.7433	0.0000	1.0000
HKY	15	-2416.8287	4863.6574	0.0000	1.0000
K2P	12	-2445.5672	4915.1343	0.0000	1.0000
SYM	16	-2441.8741	4915.7482	0.0000	1.0000
F81I	15	-2445.4871	4920.9743	0.0000	1.0000
F81IG	16	-2445.5028	4923.0057	0.0000	1.0000
F81G	15	-2446.8033	4923.6067	0.0000	1.0000
JC69I	12	-2466.6778	4957.3556	0.0000	1.0000
JC69IG	13	-2466.6873	4959.3746	0.0000	1.0000
JC69G	12	-2467.7064	4959.4128	0.0000	1.0000
F81	14	-2476.8734	4981.7467	0.0000	1.0000
JC69	11	-2499.2179	5020.4358	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYI	16	-2384.5872	4873.3147	0.5974	0.5974
HKYG	16	-2385.0880	4874.3163	0.3620	0.9594
HKYIG	17	-2384.5383	4879.7256	0.0242	0.9836
GTRI	20	-2375.4917	4881.1587	0.0118	0.9955
GTRG	20	-2376.5574	4883.2902	0.0041	0.9995
GTRIG	21	-2375.4799	4887.6439	0.0005	1.0000
K2PI	13	-2410.9948	4906.6035	0.0000	1.0000
K2PG	13	-2411.3226	4907.2592	0.0000	1.0000
K2PIG	14	-2410.9216	4912.9659	0.0000	1.0000
SYMI	17	-2408.1625	4926.9741	0.0000	1.0000
SYMG	17	-2408.9541	4928.5573	0.0000	1.0000
HKY	15	-2416.8287	4931.2889	0.0000	1.0000
SYMIG	18	-2408.1636	4933.4850	0.0000	1.0000
GTR	19	-2407.3717	4938.4099	0.0000	1.0000
K2P	12	-2445.5672	4969.2395	0.0000	1.0000
SYM	16	-2441.8741	4987.8885	0.0000	1.0000
F81I	15	-2445.4871	4988.6058	0.0000	1.0000
F81G	15	-2446.8033	4991.2382	0.0000	1.0000
F81IG	16	-2445.5028	4995.1460	0.0000	1.0000
JC69I	12	-2466.6778	5011.4608	0.0000	1.0000
JC69G	12	-2467.7064	5013.5181	0.0000	1.0000
JC69IG	13	-2466.6873	5017.9886	0.0000	1.0000
F81	14	-2476.8734	5044.8695	0.0000	1.0000
JC69	11	-2499.2179	5070.0322	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/chr6_7961.phylip = 1 - 671;
 Partition Dummy = 1:591-tomodel-phylip/chr6_7961.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 591-tomodel-phylip/chr6_7961.phylip = 1 - 671;
 Partition Dummy = 1:591-tomodel-phylip/chr6_7961.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 591-tomodel-phylip/chr6_7961.phylip = 1 - 671;
 Partition Dummy = 1:591-tomodel-phylip/chr6_7961.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

