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

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

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


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	-2499.5730	5040.5097	0.3742	0.3742
HKYG	16	-2504.1228	5041.1230	0.2753	0.6495
GTRIG	21	-2499.3739	5042.2503	0.1567	0.8062
HKYIG	17	-2503.9714	5042.9314	0.1115	0.9177
K2PG	13	-2509.7334	5046.0511	0.0234	0.9411
GTRI	20	-2502.4555	5046.2746	0.0210	0.9621
HKYI	16	-2507.0491	5046.9756	0.0148	0.9768
K2PIG	14	-2509.5926	5047.8605	0.0095	0.9863
SYMG	17	-2506.5623	5048.1134	0.0084	0.9947
SYMIG	18	-2506.3571	5049.8211	0.0036	0.9982
K2PI	13	-2512.7528	5052.0899	0.0011	0.9994
SYMI	17	-2509.1430	5053.2748	0.0006	1.0000
GTR	19	-2547.0759	5133.3835	0.0000	1.0000
HKY	15	-2551.8453	5134.4635	0.0000	1.0000
K2P	12	-2558.3149	5141.1299	0.0000	1.0000
SYM	16	-2554.4135	5141.7044	0.0000	1.0000
F81G	15	-2566.9209	5164.6147	0.0000	1.0000
F81I	15	-2567.6136	5166.0001	0.0000	1.0000
F81IG	16	-2566.7323	5166.3419	0.0000	1.0000
JC69G	12	-2571.6596	5167.8191	0.0000	1.0000
JC69I	12	-2572.3617	5169.2234	0.0000	1.0000
JC69IG	13	-2571.4637	5169.5116	0.0000	1.0000
F81	14	-2609.8365	5248.3482	0.0000	1.0000
JC69	11	-2614.8174	5252.0572	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRG	20	-2499.5730	5039.1460	0.4081	0.4081
HKYG	16	-2504.1228	5040.2455	0.2355	0.6435
GTRIG	21	-2499.3739	5040.7479	0.1832	0.8267
HKYIG	17	-2503.9714	5041.9427	0.1008	0.9275
GTRI	20	-2502.4555	5044.9110	0.0228	0.9504
K2PG	13	-2509.7334	5045.4668	0.0173	0.9677
HKYI	16	-2507.0491	5046.0982	0.0126	0.9803
SYMG	17	-2506.5623	5047.1247	0.0076	0.9878
K2PIG	14	-2509.5926	5047.1853	0.0073	0.9952
SYMIG	18	-2506.3571	5048.7143	0.0034	0.9986
K2PI	13	-2512.7528	5051.5057	0.0008	0.9994
SYMI	17	-2509.1430	5052.2861	0.0006	1.0000
GTR	19	-2547.0759	5132.1517	0.0000	1.0000
HKY	15	-2551.8453	5133.6905	0.0000	1.0000
K2P	12	-2558.3149	5140.6299	0.0000	1.0000
SYM	16	-2554.4135	5140.8270	0.0000	1.0000
F81G	15	-2566.9209	5163.8418	0.0000	1.0000
F81I	15	-2567.6136	5165.2271	0.0000	1.0000
F81IG	16	-2566.7323	5165.4645	0.0000	1.0000
JC69G	12	-2571.6596	5167.3191	0.0000	1.0000
JC69I	12	-2572.3617	5168.7234	0.0000	1.0000
JC69IG	13	-2571.4637	5168.9274	0.0000	1.0000
F81	14	-2609.8365	5247.6729	0.0000	1.0000
JC69	11	-2614.8174	5251.6348	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
K2PG	13	-2509.7334	5103.4048	0.8983	0.8983
K2PI	13	-2512.7528	5109.4437	0.0439	0.9422
K2PIG	14	-2509.5926	5109.5801	0.0410	0.9831
HKYG	16	-2504.1228	5111.5539	0.0153	0.9984
HKYI	16	-2507.0491	5117.4065	0.0008	0.9992
HKYIG	17	-2503.9714	5117.7078	0.0007	0.9999
SYMG	17	-2506.5623	5122.8898	0.0001	1.0000
SYMI	17	-2509.1430	5128.0511	0.0000	1.0000
GTRG	20	-2499.5730	5128.2814	0.0000	1.0000
SYMIG	18	-2506.3571	5128.9361	0.0000	1.0000
GTRI	20	-2502.4555	5134.0464	0.0000	1.0000
GTRIG	21	-2499.3739	5134.3400	0.0000	1.0000
K2P	12	-2558.3149	5194.1111	0.0000	1.0000
HKY	15	-2551.8453	5200.5421	0.0000	1.0000
SYM	16	-2554.4135	5212.1353	0.0000	1.0000
GTR	19	-2547.0759	5216.8303	0.0000	1.0000
JC69G	12	-2571.6596	5220.8004	0.0000	1.0000
JC69I	12	-2572.3617	5222.2046	0.0000	1.0000
JC69IG	13	-2571.4637	5226.8654	0.0000	1.0000
F81G	15	-2566.9209	5230.6933	0.0000	1.0000
F81I	15	-2567.6136	5232.0787	0.0000	1.0000
F81IG	16	-2566.7323	5236.7728	0.0000	1.0000
JC69	11	-2614.8174	5300.6592	0.0000	1.0000
F81	14	-2609.8365	5310.0677	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 591-tomodel-phylip/chr3_6253.phylip = 1 - 637;
 Partition Dummy = 1:591-tomodel-phylip/chr3_6253.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 591-tomodel-phylip/chr3_6253.phylip = 1 - 637;
 Partition Dummy = 1:591-tomodel-phylip/chr3_6253.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 (K2PG)]

BEGIN MRBAYES;
 Charset 591-tomodel-phylip/chr3_6253.phylip = 1 - 637;
 Partition Dummy = 1:591-tomodel-phylip/chr3_6253.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=2 rates=gamma;
 Prset applyto=(1) statefreqpr=Fixed(Equal) tratiopr=Beta(1.0,1.0) shapepr=Uniform(0.1,50.0);
END;


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

