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

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

Minimum AIC  model: HKYIG
Minimum AICc model: HKYIG
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
HKYIG	61	-1534.6280	3210.1187	0.9385	0.9385
GTRIG	65	-1532.0447	3215.7015	0.0576	0.9961
K2PIG	58	-1544.2669	3221.4744	0.0032	0.9993
SYMIG	62	-1541.0198	3225.5697	0.0004	0.9997
HKYG	60	-1544.4123	3227.0336	0.0002	0.9999
HKYI	60	-1545.7748	3229.7585	0.0001	1.0000
GTRG	64	-1541.1083	3231.1211	0.0000	1.0000
GTRI	64	-1541.9826	3232.8696	0.0000	1.0000
K2PG	57	-1553.5954	3237.5167	0.0000	1.0000
SYMG	61	-1548.9373	3238.7374	0.0000	1.0000
K2PI	57	-1554.4423	3239.2105	0.0000	1.0000
SYMI	61	-1552.2032	3245.2693	0.0000	1.0000
F81IG	60	-1557.1111	3252.4312	0.0000	1.0000
JC69IG	57	-1566.0214	3262.3687	0.0000	1.0000
F81I	59	-1567.0505	3269.6692	0.0000	1.0000
F81G	59	-1568.3588	3272.2859	0.0000	1.0000
JC69I	56	-1575.6121	3278.9484	0.0000	1.0000
JC69G	56	-1577.2391	3282.2024	0.0000	1.0000
HKY	59	-1598.6477	3332.8637	0.0000	1.0000
GTR	63	-1596.8190	3339.8484	0.0000	1.0000
K2P	56	-1609.3657	3346.4556	0.0000	1.0000
SYM	60	-1605.9041	3350.0172	0.0000	1.0000
F81	58	-1624.3332	3381.6070	0.0000	1.0000
JC69	55	-1634.0588	3393.2527	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
HKYIG	61	-1534.6280	3191.2559	0.8034	0.8034
GTRIG	65	-1532.0447	3194.0894	0.1948	0.9982
K2PIG	58	-1544.2669	3204.5338	0.0011	0.9993
SYMIG	62	-1541.0198	3206.0397	0.0005	0.9998
HKYG	60	-1544.4123	3208.8246	0.0001	0.9999
GTRG	64	-1541.1083	3210.2165	0.0001	0.9999
HKYI	60	-1545.7748	3211.5496	0.0000	1.0000
GTRI	64	-1541.9826	3211.9651	0.0000	1.0000
SYMG	61	-1548.9373	3219.8746	0.0000	1.0000
K2PG	57	-1553.5954	3221.1908	0.0000	1.0000
K2PI	57	-1554.4423	3222.8846	0.0000	1.0000
SYMI	61	-1552.2032	3226.4064	0.0000	1.0000
F81IG	60	-1557.1111	3234.2222	0.0000	1.0000
JC69IG	57	-1566.0214	3246.0428	0.0000	1.0000
F81I	59	-1567.0505	3252.1010	0.0000	1.0000
F81G	59	-1568.3588	3254.7177	0.0000	1.0000
JC69I	56	-1575.6121	3263.2243	0.0000	1.0000
JC69G	56	-1577.2391	3266.4782	0.0000	1.0000
HKY	59	-1598.6477	3315.2954	0.0000	1.0000
GTR	63	-1596.8190	3319.6379	0.0000	1.0000
K2P	56	-1609.3657	3330.7314	0.0000	1.0000
SYM	60	-1605.9041	3331.8082	0.0000	1.0000
F81	58	-1624.3332	3364.6664	0.0000	1.0000
JC69	55	-1634.0588	3378.1175	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYIG	61	-1534.6280	3443.6573	0.6052	0.6052
K2PIG	58	-1544.2669	3444.5219	0.3927	0.9979
K2PG	57	-1553.5954	3457.0412	0.0008	0.9987
HKYG	60	-1544.4123	3457.0882	0.0007	0.9994
K2PI	57	-1554.4423	3458.7350	0.0003	0.9997
HKYI	60	-1545.7748	3459.8132	0.0002	0.9999
SYMIG	62	-1541.0198	3462.5787	0.0000	1.0000
GTRIG	65	-1532.0447	3463.0416	0.0000	1.0000
SYMG	61	-1548.9373	3472.2759	0.0000	1.0000
GTRG	64	-1541.1083	3475.0311	0.0000	1.0000
GTRI	64	-1541.9826	3476.7797	0.0000	1.0000
SYMI	61	-1552.2032	3478.8078	0.0000	1.0000
JC69IG	57	-1566.0214	3481.8933	0.0000	1.0000
F81IG	60	-1557.1111	3482.4859	0.0000	1.0000
JC69I	56	-1575.6121	3494.9370	0.0000	1.0000
F81I	59	-1567.0505	3496.2269	0.0000	1.0000
JC69G	56	-1577.2391	3498.1909	0.0000	1.0000
F81G	59	-1568.3588	3498.8436	0.0000	1.0000
HKY	59	-1598.6477	3559.4213	0.0000	1.0000
K2P	56	-1609.3657	3562.4442	0.0000	1.0000
SYM	60	-1605.9041	3580.0718	0.0000	1.0000
GTR	63	-1596.8190	3580.3147	0.0000	1.0000
F81	58	-1624.3332	3604.6545	0.0000	1.0000
JC69	55	-1634.0588	3605.6925	0.0000	1.0000



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

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

BEGIN MRBAYES;
 Charset 183-phylip/chr8_4241.phylip = 1 - 463;
 Partition Dummy = 1:183-phylip/chr8_4241.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;


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

BEGIN MRBAYES;
 Charset 183-phylip/chr8_4241.phylip = 1 - 463;
 Partition Dummy = 1:183-phylip/chr8_4241.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;


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

BEGIN MRBAYES;
 Charset 183-phylip/chr8_4241.phylip = 1 - 463;
 Partition Dummy = 1:183-phylip/chr8_4241.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

