The commonest error with averaging is getting the mask wrong.
Check that the CCs after application of the averaging start at a
reasonable value - 0.3 at least and increase with each cycle ( by the
way why do ncycle 1?)
But in the end the density will not be identical, the Fobs are not
perfectly symmetric so there will be differences. The best idea is to
average (with correct matrices - I always find that takes several pases
before I get them all right - then build molecule A and refit it over
the others before starting refinement.
EleanorHailiang Zhang wrote:
Hi,
I am using the following DM script to perform a NCS averaging. I have a
fundemental question: after NCS averaging, are the density distrubitions
of different NCS unit being averaged supposed to be the same? I found they
are different by checking FCDM/PHICDM, and maybe I am wrong somewhere...
dm NCSIN ${PDB}.msk HKLIN ${PDBALL}.mtz HKLOUT ${PDBALL}-dm.mtz \
<<dmtest
mode AVER
ncycle 1
combine PERT
scheme ALL
solc 0.6213
#Identical
AVER REFI
NCSMASK NMER 1
ROTA MATRIX 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0
TRAN 0.0 0.0 0.0
#A>C
AVER REFI
NCSMASK NMER 1
ROTA MATRIX -0.22748 0.97259 0.04813 -0.97372 -0.22662 -0.02273
-0.01120 -0.05203 0.99858
TRAN 101.46837 81.74413 2.89341
#A>D
AVER REFI
NCSMASK NMER 1
ROTA MATRIX -0.90337 0.42792 0.02831 -0.42883 -0.90066 -0.07011
-0.00451 -0.07547 0.99714
TRAN 158.03317 36.91842 3.25853
#A>F
AVER REFI
NCSMASK NMER 1
ROTA MATRIX -0.21272 -0.97702 0.01352 0.97675 -0.21300 -0.02424 0.02657
0.00805 0.99961
TRAN 100.69797 -81.01860 -1.71365
#A>G
AVER REFI
NCSMASK NMER 1
ROTA MATRIX 0.61704 -0.78687 -0.01005 0.78590 0.61553 0.05899 -0.04023
-0.04429 0.99821
TRAN 32.50667 -65.76570 7.05504
#A>H
AVER REFI
NCSMASK NMER 1
ROTA MATRIX -0.85814 -0.51116 -0.04813 -0.51290 0.85771 0.03558
0.02310 0.05522 -0.99821
TRAN 156.14981 42.23873 48.93406
#A>I
AVER REFI
NCSMASK NMER 1
ROTA MATRIX -0.13703 -0.98975 -0.04032 -0.99048 0.13635 0.01902
-0.01332 0.04254 -0.99901
TRAN 95.97630 82.30948 52.82510
#A>J
AVER REFI
NCSMASK NMER 1
ROTA MATRIX 0.69662 -0.71695 -0.02645 -0.71716 -0.69691 0.00230
-0.02008 0.01737 -0.99965
TRAN 25.84588 59.76286 53.68224
#A>K
AVER REFI
NCSMASK NMER 1
ROTA MATRIX 0.99467 0.10258 -0.01072 0.10259 -0.99472 0.00088 -0.01057
-0.00197 -0.99994
TRAN 0.47215 -8.86082 52.78315
#A>L
AVER REFI
NCSMASK NMER 1
ROTA MATRIX 0.55782 0.82946 0.02875 0.82987 -0.55793 -0.00466 0.01218
0.02646 -0.99958
TRAN 36.68436 -68.63833 49.21587
#A>M
AVER REFI
NCSMASK NMER 1
ROTA MATRIX -0.30892 0.95102 -0.01113 0.95109 0.30890 -0.00370
-0.00008 -0.01173 -0.99993
TRAN 109.48987 -79.05550 52.39334
#A>N
AVER REFI
NCSMASK NMER 1
ROTA MATRIX -0.93676 0.34855 -0.03147 0.34600 0.93589 0.06627 0.05255
0.05119 -0.99731
TRAN 162.32979 -29.26800 45.41564
LABIN FP = FWT PHIO = PHIC FOMO = WCMB
LABOUT FDM=FDM PHIDM=PHIDM FOMDM=FOMDM FCDM=FCDM PHICDM=PHICDM
END
dmtest