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

Input data from file "183-phylip/chr4_9731.phylip" (29 taxa, 357 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	-1526.2491	3197.2279	0.7430	0.7430
HKYIG	61	-1526.0911	3199.8230	0.2030	0.9460
HKYI	60	-1529.1910	3203.1118	0.0392	0.9852
GTRG	64	-1524.5433	3205.5797	0.0114	0.9966
GTRIG	65	-1524.4097	3208.3040	0.0029	0.9995
GTRI	64	-1527.7409	3211.9750	0.0005	1.0000
K2PG	57	-1541.1400	3218.3938	0.0000	1.0000
K2PIG	58	-1540.8637	3220.6937	0.0000	1.0000
K2PI	57	-1543.8627	3223.8391	0.0000	1.0000
SYMG	61	-1540.1259	3227.8925	0.0000	1.0000
SYMIG	62	-1539.8378	3230.2471	0.0000	1.0000
SYMI	61	-1542.9101	3233.4609	0.0000	1.0000
HKY	59	-1567.3940	3276.6265	0.0000	1.0000
F81G	59	-1568.7510	3279.3405	0.0000	1.0000
F81IG	60	-1568.5611	3281.8519	0.0000	1.0000
F81I	59	-1571.1169	3284.0721	0.0000	1.0000
GTR	63	-1565.8326	3285.1874	0.0000	1.0000
JC69G	56	-1579.2228	3291.7256	0.0000	1.0000
JC69IG	57	-1578.9874	3294.0886	0.0000	1.0000
JC69I	56	-1581.5824	3296.4448	0.0000	1.0000
K2P	56	-1581.9044	3297.0889	0.0000	1.0000
SYM	60	-1580.9166	3306.5630	0.0000	1.0000
F81	58	-1607.8138	3354.5940	0.0000	1.0000
JC69	55	-1618.3372	3367.1396	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
HKYG	60	-1526.2491	3172.4982	0.6131	0.6131
HKYIG	61	-1526.0911	3174.1823	0.2642	0.8773
GTRG	64	-1524.5433	3177.0865	0.0618	0.9391
HKYI	60	-1529.1910	3178.3820	0.0324	0.9715
GTRIG	65	-1524.4097	3178.8194	0.0260	0.9975
GTRI	64	-1527.7409	3183.4818	0.0025	1.0000
K2PG	57	-1541.1400	3196.2801	0.0000	1.0000
K2PIG	58	-1540.8637	3197.7273	0.0000	1.0000
K2PI	57	-1543.8627	3201.7254	0.0000	1.0000
SYMG	61	-1540.1259	3202.2518	0.0000	1.0000
SYMIG	62	-1539.8378	3203.6756	0.0000	1.0000
SYMI	61	-1542.9101	3207.8203	0.0000	1.0000
HKY	59	-1567.3940	3252.7881	0.0000	1.0000
F81G	59	-1568.7510	3255.5021	0.0000	1.0000
F81IG	60	-1568.5611	3257.1221	0.0000	1.0000
GTR	63	-1565.8326	3257.6652	0.0000	1.0000
F81I	59	-1571.1169	3260.2338	0.0000	1.0000
JC69G	56	-1579.2228	3270.4456	0.0000	1.0000
JC69IG	57	-1578.9874	3271.9748	0.0000	1.0000
JC69I	56	-1581.5824	3275.1648	0.0000	1.0000
K2P	56	-1581.9044	3275.8089	0.0000	1.0000
SYM	60	-1580.9166	3281.8333	0.0000	1.0000
F81	58	-1607.8138	3331.6276	0.0000	1.0000
JC69	55	-1618.3372	3346.6745	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	60	-1526.2491	3405.1623	0.8949	0.8949
HKYIG	61	-1526.0911	3410.7242	0.0555	0.9504
HKYI	60	-1529.1910	3411.0462	0.0472	0.9976
K2PG	57	-1541.1400	3417.3110	0.0021	0.9997
K2PIG	58	-1540.8637	3422.6360	0.0001	0.9998
K2PI	57	-1543.8627	3422.7563	0.0001	1.0000
GTRG	64	-1524.5433	3425.2616	0.0000	1.0000
GTRIG	65	-1524.4097	3430.8723	0.0000	1.0000
GTRI	64	-1527.7409	3431.6569	0.0000	1.0000
SYMG	61	-1540.1259	3438.7937	0.0000	1.0000
SYMIG	62	-1539.8378	3444.0953	0.0000	1.0000
SYMI	61	-1542.9101	3444.3621	0.0000	1.0000
HKY	59	-1567.3940	3481.5745	0.0000	1.0000
F81G	59	-1568.7510	3484.2885	0.0000	1.0000
JC69G	56	-1579.2228	3487.5988	0.0000	1.0000
F81I	59	-1571.1169	3489.0202	0.0000	1.0000
F81IG	60	-1568.5611	3489.7863	0.0000	1.0000
JC69I	56	-1581.5824	3492.3180	0.0000	1.0000
K2P	56	-1581.9044	3492.9621	0.0000	1.0000
JC69IG	57	-1578.9874	3493.0058	0.0000	1.0000
GTR	63	-1565.8326	3501.9626	0.0000	1.0000
SYM	60	-1580.9166	3514.4974	0.0000	1.0000
F81	58	-1607.8138	3556.5363	0.0000	1.0000
JC69	55	-1618.3372	3559.9499	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/chr4_9731.phylip = 1 - 357;
 Partition Dummy = 1:183-phylip/chr4_9731.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/chr4_9731.phylip = 1 - 357;
 Partition Dummy = 1:183-phylip/chr4_9731.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/chr4_9731.phylip = 1 - 357;
 Partition Dummy = 1:183-phylip/chr4_9731.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

