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

Input data from file "183-phylip/chr1_32246.phylip" (29 taxa, 547 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	-1561.8124	3271.4627	0.9049	0.9049
GTRG	64	-1565.9958	3277.2530	0.0500	0.9549
HKYIG	61	-1570.2787	3278.1534	0.0319	0.9868
GTRI	64	-1567.3355	3279.9324	0.0131	0.9999
HKYG	60	-1577.9777	3291.0172	0.0001	1.0000
HKYI	60	-1578.1689	3291.3995	0.0000	1.0000
SYMIG	62	-1580.8593	3301.8591	0.0000	1.0000
SYMG	61	-1585.5333	3308.6625	0.0000	1.0000
SYMI	61	-1586.3280	3310.2518	0.0000	1.0000
K2PIG	58	-1593.2139	3316.4524	0.0000	1.0000
K2PI	57	-1600.4922	3328.5060	0.0000	1.0000
K2PG	57	-1600.8359	3329.1933	0.0000	1.0000
F81IG	60	-1626.6797	3388.4212	0.0000	1.0000
GTR	63	-1625.8712	3394.4380	0.0000	1.0000
F81I	59	-1633.7658	3400.0695	0.0000	1.0000
F81G	59	-1635.0000	3402.5380	0.0000	1.0000
HKY	59	-1635.2635	3403.0650	0.0000	1.0000
SYM	60	-1644.2117	3423.4852	0.0000	1.0000
JC69IG	57	-1648.7302	3424.9819	0.0000	1.0000
JC69I	56	-1655.1957	3435.4200	0.0000	1.0000
JC69G	56	-1655.2499	3435.5283	0.0000	1.0000
K2P	56	-1657.6176	3440.2638	0.0000	1.0000
F81	58	-1690.7529	3511.5304	0.0000	1.0000
JC69	55	-1709.7325	3542.0108	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRIG	65	-1561.8124	3253.6249	0.9400	0.9400
GTRG	64	-1565.9958	3259.9916	0.0390	0.9790
HKYIG	61	-1570.2787	3262.5575	0.0108	0.9898
GTRI	64	-1567.3355	3262.6710	0.0102	1.0000
HKYG	60	-1577.9777	3275.9554	0.0000	1.0000
HKYI	60	-1578.1689	3276.3378	0.0000	1.0000
SYMIG	62	-1580.8593	3285.7186	0.0000	1.0000
SYMG	61	-1585.5333	3293.0666	0.0000	1.0000
SYMI	61	-1586.3280	3294.6560	0.0000	1.0000
K2PIG	58	-1593.2139	3302.4278	0.0000	1.0000
K2PI	57	-1600.4922	3314.9845	0.0000	1.0000
K2PG	57	-1600.8359	3315.6719	0.0000	1.0000
F81IG	60	-1626.6797	3373.3595	0.0000	1.0000
GTR	63	-1625.8712	3377.7423	0.0000	1.0000
F81I	59	-1633.7658	3385.5315	0.0000	1.0000
F81G	59	-1635.0000	3388.0000	0.0000	1.0000
HKY	59	-1635.2635	3388.5270	0.0000	1.0000
SYM	60	-1644.2117	3408.4235	0.0000	1.0000
JC69IG	57	-1648.7302	3411.4604	0.0000	1.0000
JC69I	56	-1655.1957	3422.3914	0.0000	1.0000
JC69G	56	-1655.2499	3422.4997	0.0000	1.0000
K2P	56	-1657.6176	3427.2353	0.0000	1.0000
F81	58	-1690.7529	3497.5058	0.0000	1.0000
JC69	55	-1709.7325	3529.4649	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYIG	61	-1570.2787	3525.1289	0.9593	0.9593
GTRIG	65	-1561.8124	3533.4141	0.0152	0.9746
HKYG	60	-1577.9777	3534.2224	0.0102	0.9847
HKYI	60	-1578.1689	3534.6047	0.0084	0.9931
GTRG	64	-1565.9958	3535.4763	0.0054	0.9986
GTRI	64	-1567.3355	3538.1557	0.0014	1.0000
K2PIG	58	-1593.2139	3552.0858	0.0000	1.0000
SYMIG	62	-1580.8593	3552.5945	0.0000	1.0000
SYMG	61	-1585.5333	3555.6380	0.0000	1.0000
SYMI	61	-1586.3280	3557.2273	0.0000	1.0000
K2PI	57	-1600.4922	3560.3381	0.0000	1.0000
K2PG	57	-1600.8359	3561.0254	0.0000	1.0000
F81IG	60	-1626.6797	3631.6264	0.0000	1.0000
F81I	59	-1633.7658	3639.4940	0.0000	1.0000
F81G	59	-1635.0000	3641.9625	0.0000	1.0000
HKY	59	-1635.2635	3642.4895	0.0000	1.0000
GTR	63	-1625.8712	3648.9226	0.0000	1.0000
JC69IG	57	-1648.7302	3656.8140	0.0000	1.0000
JC69I	56	-1655.1957	3663.4405	0.0000	1.0000
JC69G	56	-1655.2499	3663.5489	0.0000	1.0000
SYM	60	-1644.2117	3666.6904	0.0000	1.0000
K2P	56	-1657.6176	3668.2844	0.0000	1.0000
F81	58	-1690.7529	3747.1638	0.0000	1.0000
JC69	55	-1709.7325	3766.2096	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/chr1_32246.phylip = 1 - 547;
 Partition Dummy = 1:183-phylip/chr1_32246.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/chr1_32246.phylip = 1 - 547;
 Partition Dummy = 1:183-phylip/chr1_32246.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/chr1_32246.phylip = 1 - 547;
 Partition Dummy = 1:183-phylip/chr1_32246.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

