
Output from mraic.pl version 1.4.4 by Johan Nylander
-------------------------------------------------------------
00:39 02/17/2011

Input data from file "183-phylip/chr11_4723.phylip" (29 taxa, 315 characters)

Minimum AIC  model: HKYG
Minimum AICc model: HKYG
Minimum BIC  model: HKYG


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
HKYG	60	-3167.2838	6483.3866	0.6871	0.6871
HKYIG	61	-3166.5699	6485.0370	0.3011	0.9882
GTRIG	65	-3164.1276	6492.7131	0.0065	0.9947
GTRG	64	-3165.9231	6493.1263	0.0053	0.9999
K2PG	57	-3181.9313	6503.5902	0.0000	1.0000
K2PIG	58	-3180.7334	6504.2011	0.0000	1.0000
SYMG	61	-3176.8790	6505.6553	0.0000	1.0000
SYMIG	62	-3176.7883	6508.5765	0.0000	1.0000
HKYI	60	-3183.7796	6516.3782	0.0000	1.0000
GTRI	64	-3179.7341	6520.7481	0.0000	1.0000
K2PI	57	-3195.8039	6531.3354	0.0000	1.0000
SYMI	61	-3190.7446	6533.3864	0.0000	1.0000
GTR	63	-3191.8862	6541.8999	0.0000	1.0000
HKY	59	-3198.9574	6543.6795	0.0000	1.0000
SYM	60	-3200.5725	6549.9638	0.0000	1.0000
K2P	56	-3210.3445	6557.4333	0.0000	1.0000
JC69G	56	-3282.3203	6701.3848	0.0000	1.0000
F81G	59	-3278.5101	6702.7850	0.0000	1.0000
JC69IG	57	-3282.1834	6704.0945	0.0000	1.0000
F81IG	60	-3278.3562	6705.5312	0.0000	1.0000
JC69I	56	-3293.8351	6724.4143	0.0000	1.0000
F81I	59	-3290.0373	6725.8393	0.0000	1.0000
JC69	55	-3300.3841	6734.5519	0.0000	1.0000
F81	58	-3296.9123	6736.5591	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
HKYG	60	-3167.2838	6454.5677	0.5048	0.5048
HKYIG	61	-3166.5699	6455.1397	0.3792	0.8841
GTRIG	65	-3164.1276	6458.2553	0.0799	0.9639
GTRG	64	-3165.9231	6459.8463	0.0360	1.0000
SYMG	61	-3176.8790	6475.7580	0.0000	1.0000
K2PIG	58	-3180.7334	6477.4667	0.0000	1.0000
SYMIG	62	-3176.7883	6477.5765	0.0000	1.0000
K2PG	57	-3181.9313	6477.8625	0.0000	1.0000
GTRI	64	-3179.7341	6487.4681	0.0000	1.0000
HKYI	60	-3183.7796	6487.5593	0.0000	1.0000
SYMI	61	-3190.7446	6503.4892	0.0000	1.0000
K2PI	57	-3195.8039	6505.6078	0.0000	1.0000
GTR	63	-3191.8862	6509.7724	0.0000	1.0000
HKY	59	-3198.9574	6515.9148	0.0000	1.0000
SYM	60	-3200.5725	6521.1449	0.0000	1.0000
K2P	56	-3210.3445	6532.6891	0.0000	1.0000
F81G	59	-3278.5101	6675.0203	0.0000	1.0000
JC69G	56	-3282.3203	6676.6406	0.0000	1.0000
F81IG	60	-3278.3562	6676.7123	0.0000	1.0000
JC69IG	57	-3282.1834	6678.3669	0.0000	1.0000
F81I	59	-3290.0373	6698.0746	0.0000	1.0000
JC69I	56	-3293.8351	6699.6701	0.0000	1.0000
F81	58	-3296.9123	6709.8247	0.0000	1.0000
JC69	55	-3300.3841	6710.7682	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	60	-3167.2838	6679.7221	0.8945	0.8945
HKYIG	61	-3166.5699	6684.0467	0.1029	0.9974
K2PG	57	-3181.9313	6691.7592	0.0022	0.9995
K2PIG	58	-3180.7334	6695.1160	0.0004	0.9999
GTRG	64	-3165.9231	6700.0109	0.0000	1.0000
GTRIG	65	-3164.1276	6702.1725	0.0000	1.0000
SYMG	61	-3176.8790	6704.6650	0.0000	1.0000
SYMIG	62	-3176.7883	6710.2360	0.0000	1.0000
HKYI	60	-3183.7796	6712.7137	0.0000	1.0000
K2PI	57	-3195.8039	6719.5045	0.0000	1.0000
GTRI	64	-3179.7341	6727.6328	0.0000	1.0000
SYMI	61	-3190.7446	6732.3961	0.0000	1.0000
HKY	59	-3198.9574	6737.3166	0.0000	1.0000
K2P	56	-3210.3445	6742.8332	0.0000	1.0000
GTR	63	-3191.8862	6746.1845	0.0000	1.0000
SYM	60	-3200.5725	6746.2993	0.0000	1.0000
JC69G	56	-3282.3203	6886.7846	0.0000	1.0000
JC69IG	57	-3282.1834	6892.2635	0.0000	1.0000
F81G	59	-3278.5101	6896.4221	0.0000	1.0000
F81IG	60	-3278.3562	6901.8667	0.0000	1.0000
JC69I	56	-3293.8351	6909.8142	0.0000	1.0000
JC69	55	-3300.3841	6917.1597	0.0000	1.0000
F81I	59	-3290.0373	6919.4764	0.0000	1.0000
F81	58	-3296.9123	6927.4739	0.0000	1.0000



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

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

BEGIN MRBAYES;
 Charset 183-phylip/chr11_4723.phylip = 1 - 315;
 Partition Dummy = 1:183-phylip/chr11_4723.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=2 rates=gamma;
 Prset applyto=(1) statefreqpr=Dirichlet(1.0,1.0,1.0,1.0) tratiopr=Beta(1.0,1.0) shapepr=Uniform(0.1,50.0);
END;


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

BEGIN MRBAYES;
 Charset 183-phylip/chr11_4723.phylip = 1 - 315;
 Partition Dummy = 1:183-phylip/chr11_4723.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=2 rates=gamma;
 Prset applyto=(1) statefreqpr=Dirichlet(1.0,1.0,1.0,1.0) tratiopr=Beta(1.0,1.0) shapepr=Uniform(0.1,50.0);
END;


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

BEGIN MRBAYES;
 Charset 183-phylip/chr11_4723.phylip = 1 - 315;
 Partition Dummy = 1:183-phylip/chr11_4723.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=2 rates=gamma;
 Prset applyto=(1) statefreqpr=Dirichlet(1.0,1.0,1.0,1.0) tratiopr=Beta(1.0,1.0) shapepr=Uniform(0.1,50.0);
END;


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

