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

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

Minimum AIC  model: GTRIG
Minimum AICc model: GTRIG
Minimum BIC  model: GTRIG


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
GTRIG	65	-3100.1285	6343.4367	0.9998	0.9998
GTRG	64	-3110.3823	6361.5254	0.0001	1.0000
HKYIG	61	-3115.2379	6364.0240	0.0000	1.0000
HKYG	60	-3123.0433	6377.2452	0.0000	1.0000
GTRI	64	-3127.5575	6395.8757	0.0000	1.0000
K2PIG	58	-3140.9560	6408.3131	0.0000	1.0000
SYMIG	62	-3137.5923	6411.1296	0.0000	1.0000
HKYI	60	-3142.5219	6416.2022	0.0000	1.0000
K2PG	57	-3147.8071	6419.6477	0.0000	1.0000
SYMG	61	-3143.3696	6420.2874	0.0000	1.0000
K2PI	57	-3164.3751	6452.7835	0.0000	1.0000
SYMI	61	-3162.4296	6458.4074	0.0000	1.0000
F81IG	60	-3197.3262	6525.8109	0.0000	1.0000
F81G	59	-3203.8415	6536.4593	0.0000	1.0000
JC69IG	57	-3218.8853	6561.8040	0.0000	1.0000
F81I	59	-3220.5559	6569.8881	0.0000	1.0000
JC69G	56	-3224.2793	6570.2314	0.0000	1.0000
JC69I	56	-3238.5421	6598.7569	0.0000	1.0000
GTR	63	-3269.0692	6676.4875	0.0000	1.0000
HKY	59	-3283.9875	6696.7513	0.0000	1.0000
SYM	60	-3305.8417	6742.8419	0.0000	1.0000
K2P	56	-3312.1126	6745.8979	0.0000	1.0000
F81	58	-3361.1365	6848.6743	0.0000	1.0000
JC69	55	-3384.8677	6889.0546	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRIG	65	-3100.1285	6330.2570	0.9999	0.9999
GTRG	64	-3110.3823	6348.7647	0.0001	1.0000
HKYIG	61	-3115.2379	6352.4759	0.0000	1.0000
HKYG	60	-3123.0433	6366.0867	0.0000	1.0000
GTRI	64	-3127.5575	6383.1149	0.0000	1.0000
K2PIG	58	-3140.9560	6397.9119	0.0000	1.0000
SYMIG	62	-3137.5923	6399.1846	0.0000	1.0000
HKYI	60	-3142.5219	6405.0437	0.0000	1.0000
SYMG	61	-3143.3696	6408.7393	0.0000	1.0000
K2PG	57	-3147.8071	6409.6143	0.0000	1.0000
K2PI	57	-3164.3751	6442.7501	0.0000	1.0000
SYMI	61	-3162.4296	6446.8593	0.0000	1.0000
F81IG	60	-3197.3262	6514.6523	0.0000	1.0000
F81G	59	-3203.8415	6525.6831	0.0000	1.0000
JC69IG	57	-3218.8853	6551.7706	0.0000	1.0000
F81I	59	-3220.5559	6559.1118	0.0000	1.0000
JC69G	56	-3224.2793	6560.5587	0.0000	1.0000
JC69I	56	-3238.5421	6589.0842	0.0000	1.0000
GTR	63	-3269.0692	6664.1384	0.0000	1.0000
HKY	59	-3283.9875	6685.9750	0.0000	1.0000
SYM	60	-3305.8417	6731.6833	0.0000	1.0000
K2P	56	-3312.1126	6736.2252	0.0000	1.0000
F81	58	-3361.1365	6838.2730	0.0000	1.0000
JC69	55	-3384.8677	6879.7354	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
GTRIG	65	-3100.1285	6627.6369	0.8746	0.8746
HKYIG	61	-3115.2379	6631.5555	0.1233	0.9978
HKYG	60	-3123.0433	6640.5912	0.0013	0.9992
GTRG	64	-3110.3823	6641.5695	0.0008	1.0000
K2PIG	58	-3140.9560	6663.2663	0.0000	1.0000
K2PG	57	-3147.8071	6670.3936	0.0000	1.0000
GTRI	64	-3127.5575	6675.9198	0.0000	1.0000
HKYI	60	-3142.5219	6679.5483	0.0000	1.0000
SYMIG	62	-3137.5923	6682.8393	0.0000	1.0000
SYMG	61	-3143.3696	6687.8189	0.0000	1.0000
K2PI	57	-3164.3751	6703.5294	0.0000	1.0000
SYMI	61	-3162.4296	6725.9389	0.0000	1.0000
F81IG	60	-3197.3262	6789.1569	0.0000	1.0000
F81G	59	-3203.8415	6795.6125	0.0000	1.0000
JC69IG	57	-3218.8853	6812.5499	0.0000	1.0000
JC69G	56	-3224.2793	6816.7629	0.0000	1.0000
F81I	59	-3220.5559	6829.0413	0.0000	1.0000
JC69I	56	-3238.5421	6845.2885	0.0000	1.0000
GTR	63	-3269.0692	6952.3682	0.0000	1.0000
HKY	59	-3283.9875	6955.9045	0.0000	1.0000
K2P	56	-3312.1126	6992.4294	0.0000	1.0000
SYM	60	-3305.8417	7006.1879	0.0000	1.0000
F81	58	-3361.1365	7103.6274	0.0000	1.0000
JC69	55	-3384.8677	7131.3645	0.0000	1.0000



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

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

BEGIN MRBAYES;
 Charset 183-phylip/chr12_5999.phylip = 1 - 717;
 Partition Dummy = 1:183-phylip/chr12_5999.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=6 rates=invgamma;
 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) pinvarpr=Uniform(0.0,1.0);
END;


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

BEGIN MRBAYES;
 Charset 183-phylip/chr12_5999.phylip = 1 - 717;
 Partition Dummy = 1:183-phylip/chr12_5999.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=6 rates=invgamma;
 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) pinvarpr=Uniform(0.0,1.0);
END;


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

BEGIN MRBAYES;
 Charset 183-phylip/chr12_5999.phylip = 1 - 717;
 Partition Dummy = 1:183-phylip/chr12_5999.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=6 rates=invgamma;
 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) pinvarpr=Uniform(0.0,1.0);
END;


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

