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

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

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


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	-3041.3003	6224.6006	0.5911	0.5911
HKYIG	61	-3046.4089	6225.3380	0.4088	0.9998
GTRG	64	-3050.7638	6241.1477	0.0002	1.0000
HKYG	60	-3060.6817	6251.5300	0.0000	1.0000
SYMIG	62	-3061.8991	6258.6785	0.0000	1.0000
K2PIG	58	-3072.1150	6269.7092	0.0000	1.0000
GTRI	64	-3070.0604	6279.7410	0.0000	1.0000
SYMG	61	-3076.8971	6286.3144	0.0000	1.0000
HKYI	60	-3079.4232	6289.0130	0.0000	1.0000
K2PG	57	-3088.1892	6299.5236	0.0000	1.0000
SYMI	61	-3094.1412	6320.8025	0.0000	1.0000
K2PI	57	-3105.5633	6334.2718	0.0000	1.0000
F81IG	60	-3184.1240	6498.4148	0.0000	1.0000
F81G	59	-3198.1068	6524.0333	0.0000	1.0000
JC69IG	57	-3209.9765	6543.0983	0.0000	1.0000
F81I	59	-3213.5545	6554.9287	0.0000	1.0000
JC69G	56	-3223.7894	6568.3964	0.0000	1.0000
JC69I	56	-3239.0908	6598.9992	0.0000	1.0000
GTR	63	-3250.0628	6637.3725	0.0000	1.0000
HKY	59	-3265.7656	6659.3510	0.0000	1.0000
SYM	60	-3278.4303	6687.0273	0.0000	1.0000
K2P	56	-3290.7310	6702.2798	0.0000	1.0000
F81	58	-3390.1320	6905.7432	0.0000	1.0000
JC69	55	-3417.8374	6954.1713	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRIG	65	-3041.3003	6212.6006	0.7518	0.7518
HKYIG	61	-3046.4089	6214.8179	0.2481	0.9998
GTRG	64	-3050.7638	6229.5276	0.0002	1.0000
HKYG	60	-3060.6817	6241.3634	0.0000	1.0000
SYMIG	62	-3061.8991	6247.7982	0.0000	1.0000
K2PIG	58	-3072.1150	6260.2300	0.0000	1.0000
GTRI	64	-3070.0604	6268.1209	0.0000	1.0000
SYMG	61	-3076.8971	6275.7942	0.0000	1.0000
HKYI	60	-3079.4232	6278.8463	0.0000	1.0000
K2PG	57	-3088.1892	6290.3784	0.0000	1.0000
SYMI	61	-3094.1412	6310.2824	0.0000	1.0000
K2PI	57	-3105.5633	6325.1266	0.0000	1.0000
F81IG	60	-3184.1240	6488.2481	0.0000	1.0000
F81G	59	-3198.1068	6514.2136	0.0000	1.0000
JC69IG	57	-3209.9765	6533.9531	0.0000	1.0000
F81I	59	-3213.5545	6545.1090	0.0000	1.0000
JC69G	56	-3223.7894	6559.5788	0.0000	1.0000
JC69I	56	-3239.0908	6590.1816	0.0000	1.0000
GTR	63	-3250.0628	6626.1256	0.0000	1.0000
HKY	59	-3265.7656	6649.5313	0.0000	1.0000
SYM	60	-3278.4303	6676.8606	0.0000	1.0000
K2P	56	-3290.7310	6693.4621	0.0000	1.0000
F81	58	-3390.1320	6896.2639	0.0000	1.0000
JC69	55	-3417.8374	6945.6747	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYIG	61	-3046.4089	6499.1129	0.9997	0.9997
GTRIG	65	-3041.3003	6515.5380	0.0003	1.0000
HKYG	60	-3060.6817	6520.9979	0.0000	1.0000
GTRG	64	-3050.7638	6527.8044	0.0000	1.0000
K2PIG	58	-3072.1150	6530.5434	0.0000	1.0000
SYMIG	62	-3061.8991	6536.7539	0.0000	1.0000
K2PG	57	-3088.1892	6556.0311	0.0000	1.0000
HKYI	60	-3079.4232	6558.4808	0.0000	1.0000
SYMG	61	-3076.8971	6560.0893	0.0000	1.0000
GTRI	64	-3070.0604	6566.3977	0.0000	1.0000
K2PI	57	-3105.5633	6590.7794	0.0000	1.0000
SYMI	61	-3094.1412	6594.5775	0.0000	1.0000
F81IG	60	-3184.1240	6767.8826	0.0000	1.0000
F81G	59	-3198.1068	6789.1876	0.0000	1.0000
JC69IG	57	-3209.9765	6799.6058	0.0000	1.0000
F81I	59	-3213.5545	6820.0829	0.0000	1.0000
JC69G	56	-3223.7894	6820.5710	0.0000	1.0000
JC69I	56	-3239.0908	6851.1738	0.0000	1.0000
GTR	63	-3250.0628	6919.7419	0.0000	1.0000
HKY	59	-3265.7656	6924.5052	0.0000	1.0000
K2P	56	-3290.7310	6954.4543	0.0000	1.0000
SYM	60	-3278.4303	6956.4951	0.0000	1.0000
F81	58	-3390.1320	7166.5773	0.0000	1.0000
JC69	55	-3417.8374	7202.0064	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/chr3_6135.phylip = 1 - 781;
 Partition Dummy = 1:183-phylip/chr3_6135.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/chr3_6135.phylip = 1 - 781;
 Partition Dummy = 1:183-phylip/chr3_6135.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 (HKYIG)]

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


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

