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

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

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


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	-1515.6800	3189.5049	0.8167	0.8167
GTRIG	65	-1515.6463	3192.4926	0.1833	1.0000
GTRI	64	-1527.7326	3213.6102	0.0000	1.0000
HKYG	60	-1536.3305	3218.8039	0.0000	1.0000
HKYIG	61	-1536.1212	3221.3536	0.0000	1.0000
SYMG	61	-1536.8609	3222.8329	0.0000	1.0000
SYMIG	62	-1536.6341	3225.3689	0.0000	1.0000
K2PG	57	-1547.3003	3231.9645	0.0000	1.0000
K2PIG	58	-1546.9928	3234.2551	0.0000	1.0000
HKYI	60	-1545.2124	3236.5677	0.0000	1.0000
SYMI	61	-1546.1553	3241.4218	0.0000	1.0000
K2PI	57	-1555.4816	3248.3272	0.0000	1.0000
GTR	63	-1554.0784	3263.2687	0.0000	1.0000
HKY	59	-1572.3639	3287.9235	0.0000	1.0000
SYM	60	-1574.9990	3296.1409	0.0000	1.0000
K2P	56	-1585.0209	3304.5207	0.0000	1.0000
F81G	59	-1586.6073	3316.4103	0.0000	1.0000
F81IG	60	-1586.3546	3318.8521	0.0000	1.0000
JC69G	56	-1594.2585	3322.9959	0.0000	1.0000
JC69IG	57	-1593.9647	3325.2933	0.0000	1.0000
F81I	59	-1594.9072	3333.0101	0.0000	1.0000
JC69I	56	-1602.0420	3338.5630	0.0000	1.0000
F81	58	-1619.5395	3379.3486	0.0000	1.0000
JC69	55	-1629.6938	3391.0017	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRG	64	-1515.6800	3159.3599	0.7244	0.7244
GTRIG	65	-1515.6463	3161.2926	0.2756	1.0000
GTRI	64	-1527.7326	3183.4653	0.0000	1.0000
HKYG	60	-1536.3305	3192.6611	0.0000	1.0000
HKYIG	61	-1536.1212	3194.2424	0.0000	1.0000
SYMG	61	-1536.8609	3195.7218	0.0000	1.0000
SYMIG	62	-1536.6341	3197.2682	0.0000	1.0000
K2PG	57	-1547.3003	3208.6006	0.0000	1.0000
K2PIG	58	-1546.9928	3209.9856	0.0000	1.0000
HKYI	60	-1545.2124	3210.4249	0.0000	1.0000
SYMI	61	-1546.1553	3214.3107	0.0000	1.0000
K2PI	57	-1555.4816	3224.9633	0.0000	1.0000
GTR	63	-1554.0784	3234.1568	0.0000	1.0000
HKY	59	-1572.3639	3262.7278	0.0000	1.0000
SYM	60	-1574.9990	3269.9980	0.0000	1.0000
K2P	56	-1585.0209	3282.0419	0.0000	1.0000
F81G	59	-1586.6073	3291.2146	0.0000	1.0000
F81IG	60	-1586.3546	3292.7093	0.0000	1.0000
JC69G	56	-1594.2585	3300.5171	0.0000	1.0000
JC69IG	57	-1593.9647	3301.9293	0.0000	1.0000
F81I	59	-1594.9072	3307.8144	0.0000	1.0000
JC69I	56	-1602.0420	3316.0841	0.0000	1.0000
F81	58	-1619.5395	3355.0791	0.0000	1.0000
JC69	55	-1629.6938	3369.3876	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
GTRG	64	-1515.6800	3404.6004	0.9468	0.9468
GTRIG	65	-1515.6463	3410.3650	0.0530	0.9999
HKYG	60	-1536.3305	3422.5740	0.0001	1.0000
K2PG	57	-1547.3003	3427.0179	0.0000	1.0000
HKYIG	61	-1536.1212	3427.9873	0.0000	1.0000
GTRI	64	-1527.7326	3428.7058	0.0000	1.0000
SYMG	61	-1536.8609	3429.4666	0.0000	1.0000
K2PIG	58	-1546.9928	3432.2347	0.0000	1.0000
SYMIG	62	-1536.6341	3434.8449	0.0000	1.0000
HKYI	60	-1545.2124	3440.3378	0.0000	1.0000
K2PI	57	-1555.4816	3443.3806	0.0000	1.0000
SYMI	61	-1546.1553	3448.0555	0.0000	1.0000
GTR	63	-1554.0784	3475.5654	0.0000	1.0000
HKY	59	-1572.3639	3488.8088	0.0000	1.0000
K2P	56	-1585.0209	3496.6273	0.0000	1.0000
SYM	60	-1574.9990	3499.9109	0.0000	1.0000
JC69G	56	-1594.2585	3515.1025	0.0000	1.0000
F81G	59	-1586.6073	3517.2956	0.0000	1.0000
JC69IG	57	-1593.9647	3520.3466	0.0000	1.0000
F81IG	60	-1586.3546	3522.6222	0.0000	1.0000
JC69I	56	-1602.0420	3530.6695	0.0000	1.0000
F81I	59	-1594.9072	3533.8954	0.0000	1.0000
F81	58	-1619.5395	3577.3283	0.0000	1.0000
JC69	55	-1629.6938	3580.1412	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/chrZ_495.phylip = 1 - 341;
 Partition Dummy = 1:183-phylip/chrZ_495.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/chrZ_495.phylip = 1 - 341;
 Partition Dummy = 1:183-phylip/chrZ_495.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 (GTRG)]

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


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

