Well - if the CCs are 0 then no averaging can take place.
You know you can let DM make the mask itself - are you using the GUI?
It shows you what to set..
Eleanr
[log in to unmask] wrote:
> Hi,
>
> Thanks for reminding me checking the mask. I think their might be
> something wrong with the mask, since when DM read in the mask, it says:
>
>
> Number of columns, rows, sections ............... 84 74 69
> Map mode ........................................ 0
> Start and stop points on columns, rows, sections -53 30
> 80 153 -4 64
> Grid sampling on x, y, z ........................ 136 260 150
> Cell dimensions ................................. 135.57100
> 260.11200 150.20000 90.00000 101.14000 90.00000
> Fast, medium, slow axes ......................... Z X Y
> Minimum density ................................. 0.00000
> Maximum density ................................. 0.00000
> Mean density .................................... 0.00000
> Rms deviation from mean density ................. 0.00000
> Space-group ..................................... 4
> Number of titles ................................ 1
>
>
> It seems the mask is just null. However, I converted it to a map file, and
> coot clearly showed the mask, so I am not sure why the null mask was found
> by DM. Moreover, the NCS CCs are just 0s for the mask.
>
> Anyway, following is my NCSMASK script I used to generate the above mask,
> where XYZIN is the reorganized pdb containing only a single fixed NCS unit
> (chain A). and all the operations were generated by LSQKAB with Chain A
> mapped to other chains. Not sure whether there is something wrong here or
> not...
>
> ncsmask xyzin ${PDB}_A.pdb mskout ${PDB}.msk << eof
> SYMM P1211
> EXPAND 1.0
> OVERLAP 3
> AVERAGE 12
> #Identical
> 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
> 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
> 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
> 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
> 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
> 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
> 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
> 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
> 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
> 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
> 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
> 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
> eof
>
>
>> 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
>>
>
>
|