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

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

Minimum AIC  model: SYMG
Minimum AICc model: SYMG
Minimum BIC  model: SYMG


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
SYMG	17	-2787.1553	5609.1595	0.5090	0.5090
SYMIG	18	-2787.1551	5611.2602	0.1781	0.6871
SYMI	17	-2788.2247	5611.2982	0.1747	0.8618
GTRG	20	-2785.7571	5612.6840	0.0874	0.9491
GTRIG	21	-2785.7570	5614.8027	0.0303	0.9794
GTRI	20	-2787.2039	5615.5777	0.0206	1.0000
K2PI	13	-2802.3795	5631.2611	0.0000	1.0000
K2PG	13	-2802.8384	5632.1788	0.0000	1.0000
HKYI	16	-2800.5303	5633.8141	0.0000	1.0000
K2PIG	14	-2802.8383	5634.2567	0.0000	1.0000
HKYG	16	-2801.1612	5635.0759	0.0000	1.0000
HKYIG	17	-2801.1614	5637.1716	0.0000	1.0000
JC69I	12	-2865.5054	5755.4405	0.0000	1.0000
F81I	15	-2862.5555	5755.7750	0.0000	1.0000
JC69G	12	-2866.0093	5756.4483	0.0000	1.0000
F81G	15	-2863.1949	5757.0538	0.0000	1.0000
JC69IG	13	-2866.0093	5758.5206	0.0000	1.0000
F81IG	16	-2863.1951	5759.1437	0.0000	1.0000
SYM	16	-2864.7511	5762.2556	0.0000	1.0000
GTR	19	-2862.9440	5764.9450	0.0000	1.0000
K2P	12	-2873.9910	5772.4118	0.0000	1.0000
HKY	15	-2870.8770	5772.4180	0.0000	1.0000
F81	14	-2928.5491	5885.6783	0.0000	1.0000
JC69	11	-2932.3902	5887.1436	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
SYMG	17	-2787.1553	5608.3106	0.4916	0.4916
SYMIG	18	-2787.1551	5610.3102	0.1809	0.6724
SYMI	17	-2788.2247	5610.4494	0.1687	0.8411
GTRG	20	-2785.7571	5611.5141	0.0991	0.9402
GTRIG	21	-2785.7570	5613.5140	0.0364	0.9767
GTRI	20	-2787.2039	5614.4077	0.0233	1.0000
K2PI	13	-2802.3795	5630.7591	0.0000	1.0000
K2PG	13	-2802.8384	5631.6767	0.0000	1.0000
HKYI	16	-2800.5303	5633.0606	0.0000	1.0000
K2PIG	14	-2802.8383	5633.6766	0.0000	1.0000
HKYG	16	-2801.1612	5634.3225	0.0000	1.0000
HKYIG	17	-2801.1614	5636.3228	0.0000	1.0000
JC69I	12	-2865.5054	5755.0107	0.0000	1.0000
F81I	15	-2862.5555	5755.1111	0.0000	1.0000
JC69G	12	-2866.0093	5756.0186	0.0000	1.0000
F81G	15	-2863.1949	5756.3899	0.0000	1.0000
JC69IG	13	-2866.0093	5758.0185	0.0000	1.0000
F81IG	16	-2863.1951	5758.3902	0.0000	1.0000
SYM	16	-2864.7511	5761.5022	0.0000	1.0000
GTR	19	-2862.9440	5763.8879	0.0000	1.0000
HKY	15	-2870.8770	5771.7541	0.0000	1.0000
K2P	12	-2873.9910	5771.9820	0.0000	1.0000
F81	14	-2928.5491	5885.0982	0.0000	1.0000
JC69	11	-2932.3902	5886.7805	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
SYMG	17	-2787.1553	5686.6007	0.6245	0.6245
SYMI	17	-2788.2247	5688.7394	0.2143	0.8388
K2PI	13	-2802.3795	5690.6280	0.0834	0.9222
K2PG	13	-2802.8384	5691.5456	0.0527	0.9749
SYMIG	18	-2787.1551	5693.2056	0.0230	0.9979
K2PIG	14	-2802.8383	5698.1508	0.0019	0.9998
GTRG	20	-2785.7571	5703.6201	0.0001	0.9999
GTRI	20	-2787.2039	5706.5137	0.0000	1.0000
HKYI	16	-2800.5303	5706.7454	0.0000	1.0000
HKYG	16	-2801.1612	5708.0072	0.0000	1.0000
GTRIG	21	-2785.7570	5710.2253	0.0000	1.0000
HKYIG	17	-2801.1614	5714.6128	0.0000	1.0000
JC69I	12	-2865.5054	5810.2743	0.0000	1.0000
JC69G	12	-2866.0093	5811.2821	0.0000	1.0000
JC69IG	13	-2866.0093	5817.8874	0.0000	1.0000
F81I	15	-2862.5555	5824.1905	0.0000	1.0000
F81G	15	-2863.1949	5825.4693	0.0000	1.0000
K2P	12	-2873.9910	5827.2456	0.0000	1.0000
F81IG	16	-2863.1951	5832.0750	0.0000	1.0000
SYM	16	-2864.7511	5835.1869	0.0000	1.0000
HKY	15	-2870.8770	5840.8335	0.0000	1.0000
GTR	19	-2862.9440	5851.3886	0.0000	1.0000
JC69	11	-2932.3902	5937.4387	0.0000	1.0000
F81	14	-2928.5491	5949.5724	0.0000	1.0000



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

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

BEGIN MRBAYES;
 Charset 591-tomodel-phylip/chr2_13064.phylip = 1 - 739;
 Partition Dummy = 1:591-tomodel-phylip/chr2_13064.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=Fixed(Equal);
END;


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

BEGIN MRBAYES;
 Charset 591-tomodel-phylip/chr2_13064.phylip = 1 - 739;
 Partition Dummy = 1:591-tomodel-phylip/chr2_13064.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=Fixed(Equal);
END;


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

BEGIN MRBAYES;
 Charset 591-tomodel-phylip/chr2_13064.phylip = 1 - 739;
 Partition Dummy = 1:591-tomodel-phylip/chr2_13064.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=Fixed(Equal);
END;


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

