# NOT RUN {
## Read PDB structure
pdb <- read.pdb( system.file("examples/1hel.pdb", package="bio3d") )
## Perform NMA
modes.a <- nma(pdb, ff="calpha")
modes.b <- nma(pdb, ff="anm")
## Calculate and plot RMSIP
r <- rmsip(modes.a, modes.b)
plot(r)
# }
Run the code above in your browser using DataLab