
Output from mraic.pl version 1.4.4 by Johan Nylander
-------------------------------------------------------------
21:54 04/04/2011

Input data from file "591-tomodel-phylip/chrZ_4747.phylip" (7 taxa, 810 characters)

Minimum AIC  model: GTRI
Minimum AICc model: GTRI
Minimum BIC  model: HKYI


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
GTRI	20	-3133.5424	6308.1495	0.7176	0.7176
GTRIG	21	-3133.5751	6310.3228	0.2421	0.9597
GTRG	20	-3136.4487	6313.9620	0.0392	0.9989
HKYI	16	-3144.6168	6321.9196	0.0007	0.9997
HKYIG	17	-3144.6388	6324.0503	0.0003	0.9999
HKYG	16	-3146.9582	6326.6024	0.0001	1.0000
SYMI	17	-3164.0531	6362.8789	0.0000	1.0000
SYMIG	18	-3164.0644	6364.9936	0.0000	1.0000
SYMG	17	-3165.7937	6366.3600	0.0000	1.0000
K2PI	13	-3170.1933	6366.8439	0.0000	1.0000
K2PIG	14	-3170.1975	6368.9233	0.0000	1.0000
GTR	19	-3165.4118	6369.7856	0.0000	1.0000
K2PG	13	-3171.7918	6370.0408	0.0000	1.0000
HKY	15	-3175.7887	6382.1819	0.0000	1.0000
SYM	16	-3194.4696	6421.6251	0.0000	1.0000
K2P	12	-3201.2360	6426.8635	0.0000	1.0000
F81I	15	-3227.6257	6485.8559	0.0000	1.0000
F81IG	16	-3227.6507	6487.9874	0.0000	1.0000
F81G	15	-3229.7560	6490.1165	0.0000	1.0000
JC69I	12	-3247.5193	6519.4301	0.0000	1.0000
JC69IG	13	-3247.5382	6521.5338	0.0000	1.0000
JC69G	12	-3249.3424	6523.0763	0.0000	1.0000
F81	14	-3254.2523	6537.0330	0.0000	1.0000
JC69	11	-3273.8842	6570.0992	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRI	20	-3133.5424	6307.0848	0.7082	0.7082
GTRIG	21	-3133.5751	6309.1502	0.2522	0.9604
GTRG	20	-3136.4487	6312.8974	0.0387	0.9991
HKYI	16	-3144.6168	6321.2336	0.0006	0.9997
HKYIG	17	-3144.6388	6323.2776	0.0002	0.9999
HKYG	16	-3146.9582	6325.9164	0.0001	1.0000
SYMI	17	-3164.0531	6362.1062	0.0000	1.0000
SYMIG	18	-3164.0644	6364.1289	0.0000	1.0000
SYMG	17	-3165.7937	6365.5873	0.0000	1.0000
K2PI	13	-3170.1933	6366.3867	0.0000	1.0000
K2PIG	14	-3170.1975	6368.3950	0.0000	1.0000
GTR	19	-3165.4118	6368.8236	0.0000	1.0000
K2PG	13	-3171.7918	6369.5836	0.0000	1.0000
HKY	15	-3175.7887	6381.5774	0.0000	1.0000
SYM	16	-3194.4696	6420.9391	0.0000	1.0000
K2P	12	-3201.2360	6426.4720	0.0000	1.0000
F81I	15	-3227.6257	6485.2514	0.0000	1.0000
F81IG	16	-3227.6507	6487.3014	0.0000	1.0000
F81G	15	-3229.7560	6489.5120	0.0000	1.0000
JC69I	12	-3247.5193	6519.0387	0.0000	1.0000
JC69IG	13	-3247.5382	6521.0765	0.0000	1.0000
JC69G	12	-3249.3424	6522.6849	0.0000	1.0000
F81	14	-3254.2523	6536.5047	0.0000	1.0000
JC69	11	-3273.8842	6569.7684	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYI	16	-3144.6168	6396.3861	0.8080	0.8080
GTRI	20	-3133.5424	6401.0255	0.0794	0.8875
HKYG	16	-3146.9582	6401.0690	0.0777	0.9652
HKYIG	17	-3144.6388	6403.1271	0.0278	0.9930
GTRG	20	-3136.4487	6406.8381	0.0043	0.9973
GTRIG	21	-3133.5751	6407.7879	0.0027	1.0000
K2PI	13	-3170.1933	6427.4481	0.0000	1.0000
K2PG	13	-3171.7918	6430.6450	0.0000	1.0000
K2PIG	14	-3170.1975	6434.1534	0.0000	1.0000
SYMI	17	-3164.0531	6441.9558	0.0000	1.0000
SYMG	17	-3165.7937	6445.4369	0.0000	1.0000
SYMIG	18	-3164.0644	6448.6755	0.0000	1.0000
HKY	15	-3175.7887	6452.0329	0.0000	1.0000
GTR	19	-3165.4118	6458.0673	0.0000	1.0000
K2P	12	-3201.2360	6482.8364	0.0000	1.0000
SYM	16	-3194.4696	6496.0917	0.0000	1.0000
F81I	15	-3227.6257	6555.7069	0.0000	1.0000
F81G	15	-3229.7560	6559.9675	0.0000	1.0000
F81IG	16	-3227.6507	6562.4539	0.0000	1.0000
JC69I	12	-3247.5193	6575.4031	0.0000	1.0000
JC69G	12	-3249.3424	6579.0493	0.0000	1.0000
JC69IG	13	-3247.5382	6582.1379	0.0000	1.0000
F81	14	-3254.2523	6602.2632	0.0000	1.0000
JC69	11	-3273.8842	6621.4358	0.0000	1.0000



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

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

BEGIN MRBAYES;
 Charset 591-tomodel-phylip/chrZ_4747.phylip = 1 - 810;
 Partition Dummy = 1:591-tomodel-phylip/chrZ_4747.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=6 rates=propinv;
 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) pinvarpr=Uniform(0.0,1.0);
END;


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

BEGIN MRBAYES;
 Charset 591-tomodel-phylip/chrZ_4747.phylip = 1 - 810;
 Partition Dummy = 1:591-tomodel-phylip/chrZ_4747.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=6 rates=propinv;
 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) pinvarpr=Uniform(0.0,1.0);
END;


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

BEGIN MRBAYES;
 Charset 591-tomodel-phylip/chrZ_4747.phylip = 1 - 810;
 Partition Dummy = 1:591-tomodel-phylip/chrZ_4747.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=2 rates=propinv;
 Prset applyto=(1) statefreqpr=Dirichlet(1.0,1.0,1.0,1.0) tratiopr=Beta(1.0,1.0) pinvarpr=Uniform(0.0,1.0);
END;


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

