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

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

Minimum AIC  model: GTRG
Minimum AICc model: GTRG
Minimum BIC  model: K2PG


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
GTRG	64	-1436.1144	3051.2717	0.4113	0.4113
K2PIG	58	-1447.7004	3051.8978	0.3007	0.7120
K2PG	57	-1449.8712	3052.6364	0.2079	0.9198
HKYG	60	-1446.1706	3056.1735	0.0355	0.9553
GTRIG	65	-1437.4672	3057.8974	0.0150	0.9703
SYMIG	62	-1443.3639	3058.0734	0.0137	0.9840
SYMG	61	-1445.2849	3058.1360	0.0133	0.9973
HKYIG	61	-1446.8711	3061.3084	0.0027	1.0000
JC69IG	57	-1467.0603	3087.0147	0.0000	1.0000
JC69G	56	-1468.9698	3087.2729	0.0000	1.0000
K2PI	57	-1468.7258	3090.3457	0.0000	1.0000
F81G	59	-1466.1274	3092.3977	0.0000	1.0000
SYMI	61	-1462.4448	3092.4558	0.0000	1.0000
GTRI	64	-1457.8399	3094.7227	0.0000	1.0000
F81IG	60	-1465.6939	3095.2200	0.0000	1.0000
HKYI	60	-1465.8988	3095.6299	0.0000	1.0000
GTR	63	-1465.6575	3106.4858	0.0000	1.0000
K2P	56	-1479.3624	3108.0581	0.0000	1.0000
HKY	59	-1474.8947	3109.9323	0.0000	1.0000
SYM	60	-1473.6983	3111.2289	0.0000	1.0000
JC69I	56	-1486.4689	3122.2711	0.0000	1.0000
F81I	59	-1484.4570	3129.0568	0.0000	1.0000
JC69	55	-1497.8909	3141.5958	0.0000	1.0000
F81	58	-1495.0083	3146.5136	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRG	64	-1436.1144	3000.2288	0.9008	0.9008
GTRIG	65	-1437.4672	3004.9345	0.0857	0.9865
SYMIG	62	-1443.3639	3010.7279	0.0047	0.9912
K2PIG	58	-1447.7004	3011.4007	0.0034	0.9946
HKYG	60	-1446.1706	3012.3412	0.0021	0.9967
SYMG	61	-1445.2849	3012.5698	0.0019	0.9986
K2PG	57	-1449.8712	3013.7423	0.0010	0.9996
HKYIG	61	-1446.8711	3015.7421	0.0004	1.0000
GTRI	64	-1457.8399	3043.6798	0.0000	1.0000
SYMI	61	-1462.4448	3046.8895	0.0000	1.0000
JC69IG	57	-1467.0603	3048.1206	0.0000	1.0000
JC69G	56	-1468.9698	3049.9395	0.0000	1.0000
F81G	59	-1466.1274	3050.2548	0.0000	1.0000
F81IG	60	-1465.6939	3051.3877	0.0000	1.0000
K2PI	57	-1468.7258	3051.4516	0.0000	1.0000
HKYI	60	-1465.8988	3051.7975	0.0000	1.0000
GTR	63	-1465.6575	3057.3150	0.0000	1.0000
SYM	60	-1473.6983	3067.3965	0.0000	1.0000
HKY	59	-1474.8947	3067.7894	0.0000	1.0000
K2P	56	-1479.3624	3070.7247	0.0000	1.0000
JC69I	56	-1486.4689	3084.9378	0.0000	1.0000
F81I	59	-1484.4570	3086.9139	0.0000	1.0000
JC69	55	-1497.8909	3105.7819	0.0000	1.0000
F81	58	-1495.0083	3106.0165	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
K2PG	57	-1449.8712	3209.2150	0.6247	0.6247
K2PIG	58	-1447.7004	3210.3028	0.3626	0.9874
HKYG	60	-1446.1706	3218.1019	0.0073	0.9947
GTRG	64	-1436.1144	3219.7069	0.0033	0.9980
SYMG	61	-1445.2849	3221.7598	0.0012	0.9992
SYMIG	62	-1443.3639	3223.3473	0.0005	0.9997
HKYIG	61	-1446.8711	3224.9322	0.0002	0.9999
GTRIG	65	-1437.4672	3227.8419	0.0001	1.0000
JC69G	56	-1468.9698	3241.9829	0.0000	1.0000
JC69IG	57	-1467.0603	3243.5933	0.0000	1.0000
K2PI	57	-1468.7258	3246.9243	0.0000	1.0000
F81G	59	-1466.1274	3252.5862	0.0000	1.0000
SYMI	61	-1462.4448	3256.0796	0.0000	1.0000
F81IG	60	-1465.6939	3257.1484	0.0000	1.0000
HKYI	60	-1465.8988	3257.5583	0.0000	1.0000
K2P	56	-1479.3624	3262.7681	0.0000	1.0000
GTRI	64	-1457.8399	3263.1579	0.0000	1.0000
HKY	59	-1474.8947	3270.1208	0.0000	1.0000
SYM	60	-1473.6983	3273.1573	0.0000	1.0000
GTR	63	-1465.6575	3273.3638	0.0000	1.0000
JC69I	56	-1486.4689	3276.9811	0.0000	1.0000
F81I	59	-1484.4570	3289.2453	0.0000	1.0000
JC69	55	-1497.8909	3294.3959	0.0000	1.0000
F81	58	-1495.0083	3304.9186	0.0000	1.0000



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

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

BEGIN MRBAYES;
 Charset 183-phylip/chr12_5828.phylip = 1 - 228;
 Partition Dummy = 1:183-phylip/chr12_5828.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=6 rates=gamma;
 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);
END;


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

BEGIN MRBAYES;
 Charset 183-phylip/chr12_5828.phylip = 1 - 228;
 Partition Dummy = 1:183-phylip/chr12_5828.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=6 rates=gamma;
 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);
END;


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

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


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

