; Apply reference-free alignment to the particles series


X40=2 ;number of classes (with & without ligand)

DO LB11 x50=1,83

VM
mkdir groups/grp{***x50}

UD N,x91,x92 ; Number of particles
views/docview/view{***x50}


; Prepare SEQuential file
; Code number
X95=x50

CA SI
views/view{**x50}/sar*****
1-X91
mask/mask1{***x50}
SEQX95
Y
(0.0)


CA S ; Perform PCA, get 5 factors ; before 20 factors
SEQX95
5,X95
Y
(0)


CL KM
2
X95
X40
1-5
(0)
1457
groups/grp{***x50}/sel***
groups/grp{***x50}/MAP{***x95}

; Calculate group averages


DO LB81 I=1,X40
X31=x0
AS DC
views/view{**x50}/sar*****
groups/grp{***x50}/sel{***X31}
A
groups/grp{***x50}/avg{***X31}
groups/grp{***x50}/var{***X31}
LB81

; Calculate the resolution for each group

DO LB82 I=1,X40
X31=x0
AS DC
views/view{**x50}/sar*****
groups/grp{***x50}/sel{***X31}
O
_1
_2
_3
_4
RF
_1
_3
(0.5)
(0.2,2.0)
groups/grp{***x50}/dres{***X31}

LB82

LB11

en