## 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