cryoET 3D classification

[Esther Bullitt]

Hi,

Is there somewhere I can read about how to choose the 3D classification scheme for a cryoET dataset? — EMAN2 includes ortho-proj that uses MSA, PCA-based, and multi-reference 3D classification.

Do I try them all, or is there a way to make an informed decision?

thanks very much,

Esther

[Muyuan Chen]

Unfortunately I don’t think there is one... A main reason is that I keep making new programs and changing existing one. 

I often start from e2spt_sgd_new.py. It has the —ncls and —classify options that does classification, and can also load from existing aligned particles. More recently it also supports —realspace and —mask that let you focus classification on local regions. Notably it only gives structures of multiple classes but not classify all particles. For that it requires e2spt_refinemulti_new.py. This one also works with or without initial alignment, with or without mask, and it split the particles to corresponding classes. It also taks subtilt refinement results for higher resolution classification. However, its convergence rate is much slower so it is better to run it with references. There is also now e2gmm_spt_heter_refine.py, which is probably better at continuous motion of local regions. 

Muyuan

On Apr 7, 2026, at 7:29 AM, 'Esther Bullitt' via EMAN2 <em...@googlegroups.com> wrote:

Hi,

--
--
----------------------------------------------------------------------------------------------
You received this message because you are subscribed to the Google
Groups "EMAN2" group.
To post to this group, send email to em...@googlegroups.com
To unsubscribe from this group, send email to eman2+un...@googlegroups.com
For more options, visit this group at
http://groups.google.com/group/eman2

---
You received this message because you are subscribed to the Google Groups "EMAN2" group.
To unsubscribe from this group and stop receiving emails from it, send an email to eman2+un...@googlegroups.com.
To view this discussion visit https://groups.google.com/d/msgid/eman2/310239d2-6a55-4d89-93fb-373aeebaf10an%40googlegroups.com.

[Steve Ludtke]

You left out the GMM-based dynamics classification. That's actually probably the one I'd use as a first attempt nowadays, unless you are confident that you know your data should divide into N<5 discrete classes (and you know N).  It's use on CryoET data is still a bit experimental, but it worked pretty well in some testing I did.

[Esther Bullitt]

The GMM-based dynamics classification looks promising.

For those interested, here is the documentation website:  https://blake.bcm.edu/emanwiki/doku.php?id=eman2:e2gmm

Thanks!

[Steve Ludtke]

I'd actually recommend starting with the video tutorial to explain some of the interface, though it's a few years old now, and probably needs an update, it will at least have the basics:

https://youtu.be/1OM-XIzkSWw

The earlier videos in that series also have good tutorials on tomography and single particle work (as of 2022). Most of it is still relevant.

To view this discussion visit https://groups.google.com/d/msgid/eman2/ee58ad05-3ec0-4a89-8101-cda2b6d83b82n%40googlegroups.com.

[Bullitt, Esther]

1.  $ e2version.py

EMAN 2.99.72 ( GITHUB: 2026-03-06 12:57 - commit: 1bb4904d8 )

Your EMAN2 is running on: Linux-6.5.0-28-generic-x86_64-with-glibc2.39 6.5.0-28-generic

Your Python version is: 3.12.13

2.  e2spt_refine_new.py --ptcls=sets/extracted_tomobox.lst --ref=clip72_CS20_mask.mrc --startres=30.0 --goldstandard --sym=H1:3:111.8:8.95 --iters=p3,t2,p,t,r,d --keep=0.95 --localrefine --maxres=0.0 --minres=0.0 --parallel=thread:4 --threads=6

3.  tomography, ~500K particles, 72 pix^3

We have been having problems with the computer hanging, and needing to be rebooted by pressing the power button (it works on smaller datasets).

The first problem was a full SSD drive, from old cached data (NOT eman2), which we deleted.

Our IT person wrote me to say:

Something is wrong with your eman2 installation or how you are running. You were running out of RAM just for test I created 200gb swap file and already over 100gb used. Eman2 suppose to create it is own cache but it does not.

He created a swap file on the SSD, but I don’t know how to get EMAN2 to use it — is there a configuration file somewhere that I can edit? 

Thanks!

Sincerely,

Esther

-- 

Esther Bullitt, PhD

Professor of Pharmacology, Physiology & Biophysics

Director, Cryogenic Electron Microscopy Core

Boston University Chobanian & Avedisian School of Medicine

700 Albany St, W-302

Boston, MA. 02118