Hi Abhilash,

  The ChimeraX "volume subtract" command has an option "minrms true" (minimize root mean square) that will find the scale factor for a second map such that when the scaled second map is subtracted from the first map the sum of the root mean square residuals is minimum.  The scale factor is reported in the log.

    https://plato.cgl.ucsf.edu/chimerax/docs/user/commands/volume.html#subtract

Here's an example.

open 43675 from emdb
open 43676 from emdb
open 8vyq
open 8vyr
matchmaker #4 to #3
view position #2 sameas #4
        volume #1 level 0.5
volume #2 level 0.1
volume #1 region 0,0,0,255,255,115
volume subtract #2,1 minrms true
> Minimum RMS scale factor for "emdb 43675 #1" above level 0.5 is 0.14031

Several important points to note.  The maps need to be aligned with each other.  I do that in the example by opening fit atomic models, moving the second atomic model to align with the first, and then moving the second map to the same position as the moved atomic model.  For the subtraction to give a useful scale factor you want to limit it to a region where the two maps mostly overlap.  I cropped the first map to such a region in the example.  Then I subtracted the first map from the second -- the RMS is minimized within the contour level of the map that is being subtracted (#1 in this example), so you also might want to choose the contour level.  The result in the example is that map #1 needs to be multiplied by 0.13696 to minimize the RMS difference.  I think this is prone to giving a bad scaling factor if there is significant non-overlapping density for the region within the map contour level.  You can judge how well the scaling worked by looking at how much density is left in the difference map.

Tom

minrms.png

On May 30, 2025, at 10:58 AM, Abhilash Padavannil via ChimeraX-users <chimerax-users@cgl.ucsf.edu> wrote:

Hello,

I have five different cryoEM maps coming from similar samples. The maps are at different resolutions. I want to compare the density at the region of interest across these maps. How do I normalize the maps say to the first map to make accurate analysis.

I know I could scale them using vop scale, but I am wondering which of the parameters I should use.

Thanks in advance for your suggestions.

Abhilash.


_______________________________________________
ChimeraX-users mailing list -- chimerax-users@cgl.ucsf.edu
To unsubscribe send an email to chimerax-users-leave@cgl.ucsf.edu
Archives: https://mail.cgl.ucsf.edu/mailman/archives/list/chimerax-users@cgl.ucsf.edu/