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

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

Minimum AIC  model: HKYG
Minimum AICc model: HKYG
Minimum BIC  model: HKYG


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
HKYG	60	-1524.1540	3187.1739	0.7828	0.7828
HKYIG	61	-1524.6783	3190.9019	0.1214	0.9042
HKYI	60	-1526.5038	3191.8735	0.0747	0.9788
GTRG	64	-1522.5499	3194.7664	0.0176	0.9964
GTRIG	65	-1523.3316	3199.0652	0.0020	0.9985
GTRI	64	-1525.0204	3199.7075	0.0015	0.9999
K2PG	57	-1537.9326	3206.7756	0.0000	1.0000
K2PIG	58	-1538.2725	3210.0937	0.0000	1.0000
K2PI	57	-1539.9483	3210.8071	0.0000	1.0000
SYMG	61	-1536.8894	3215.3241	0.0000	1.0000
SYMIG	62	-1537.2553	3218.7489	0.0000	1.0000
SYMI	61	-1539.0799	3219.7049	0.0000	1.0000
HKY	59	-1550.6180	3237.4365	0.0000	1.0000
GTR	63	-1549.1789	3245.3032	0.0000	1.0000
K2P	56	-1563.6275	3255.5407	0.0000	1.0000
F81G	59	-1559.8200	3255.8406	0.0000	1.0000
F81IG	60	-1559.5260	3257.9180	0.0000	1.0000
F81I	59	-1561.3837	3258.9678	0.0000	1.0000
SYM	60	-1562.6711	3264.2082	0.0000	1.0000
JC69G	56	-1570.8410	3269.9676	0.0000	1.0000
JC69IG	57	-1570.5368	3271.9841	0.0000	1.0000
JC69I	56	-1572.4885	3273.2626	0.0000	1.0000
F81	58	-1585.2320	3304.0126	0.0000	1.0000
JC69	55	-1596.0712	3317.8168	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
HKYG	60	-1524.1540	3168.3080	0.7006	0.7006
HKYIG	61	-1524.6783	3171.3566	0.1526	0.8532
HKYI	60	-1526.5038	3173.0076	0.0668	0.9200
GTRG	64	-1522.5499	3173.0997	0.0638	0.9838
GTRIG	65	-1523.3316	3176.6631	0.0107	0.9946
GTRI	64	-1525.0204	3178.0408	0.0054	1.0000
K2PG	57	-1537.9326	3189.8651	0.0000	1.0000
K2PIG	58	-1538.2725	3192.5450	0.0000	1.0000
K2PI	57	-1539.9483	3193.8966	0.0000	1.0000
SYMG	61	-1536.8894	3195.7789	0.0000	1.0000
SYMIG	62	-1537.2553	3198.5105	0.0000	1.0000
SYMI	61	-1539.0799	3200.1597	0.0000	1.0000
HKY	59	-1550.6180	3219.2360	0.0000	1.0000
GTR	63	-1549.1789	3224.3578	0.0000	1.0000
F81G	59	-1559.8200	3237.6400	0.0000	1.0000
F81IG	60	-1559.5260	3239.0520	0.0000	1.0000
K2P	56	-1563.6275	3239.2549	0.0000	1.0000
F81I	59	-1561.3837	3240.7673	0.0000	1.0000
SYM	60	-1562.6711	3245.3422	0.0000	1.0000
JC69G	56	-1570.8410	3253.6819	0.0000	1.0000
JC69IG	57	-1570.5368	3255.0736	0.0000	1.0000
JC69I	56	-1572.4885	3256.9769	0.0000	1.0000
F81	58	-1585.2320	3286.4639	0.0000	1.0000
JC69	55	-1596.0712	3302.1425	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	60	-1524.1540	3414.7293	0.8812	0.8812
HKYI	60	-1526.5038	3419.4289	0.0841	0.9652
HKYIG	61	-1524.6783	3421.8850	0.0246	0.9898
K2PG	57	-1537.9326	3423.9654	0.0087	0.9985
K2PI	57	-1539.9483	3427.9969	0.0012	0.9997
K2PIG	58	-1538.2725	3430.7523	0.0003	1.0000
GTRG	64	-1522.5499	3435.9492	0.0000	1.0000
GTRI	64	-1525.0204	3440.8903	0.0000	1.0000
GTRIG	65	-1523.3316	3443.6196	0.0000	1.0000
SYMG	61	-1536.8894	3446.3073	0.0000	1.0000
SYMI	61	-1539.0799	3450.6881	0.0000	1.0000
SYMIG	62	-1537.2553	3453.1459	0.0000	1.0000
HKY	59	-1550.6180	3461.5504	0.0000	1.0000
K2P	56	-1563.6275	3469.2482	0.0000	1.0000
F81G	59	-1559.8200	3479.9544	0.0000	1.0000
F81I	59	-1561.3837	3483.0817	0.0000	1.0000
GTR	63	-1549.1789	3483.1002	0.0000	1.0000
JC69G	56	-1570.8410	3483.6752	0.0000	1.0000
F81IG	60	-1559.5260	3485.4734	0.0000	1.0000
JC69I	56	-1572.4885	3486.9702	0.0000	1.0000
JC69IG	57	-1570.5368	3489.1739	0.0000	1.0000
SYM	60	-1562.6711	3491.7636	0.0000	1.0000
F81	58	-1585.2320	3524.6712	0.0000	1.0000
JC69	55	-1596.0712	3528.0288	0.0000	1.0000



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

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

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


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

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


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

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


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

