## Not run:
# ##-- Generate a small kinesin alignment and read corresponding structures
# pdbfiles <- get.pdb(c("1bg2","2ncd","1i6i","1i5s"), URLonly=TRUE)
# pdbs <- pdbaln(pdbfiles)
#
# ##-- Find 'core' positions
# core <- core.find(pdbs)
# plot(core)
#
# ##-- Fit on these relatively invarient subset of positions
# core.inds <- print(core, vol=0.5)
#
# print(core, vol=0.7)
# print(core, vol=1.0)
#
# ## End(Not run)
Run the code above in your browser using DataLab