## Not run: 
# attach(transducin)
# 
# gaps.pos <- gap.inspect(pdbs$xyz)
# 
# #-- Do PCA without structures 2 and 7
# pc.xray <- pca.xyz(pdbs$xyz[-c(2,7), gaps.pos$f.inds])
# 
# #-- Project structures 2 and 7 onto the PC space
# d <- project.pca(pdbs$xyz[c(2,7), gaps.pos$f.inds], pc.xray)
# 
# plot(pc.xray$z[,1], pc.xray$z[,2],col="gray")
# points(d[,1],d[,2], col="red")
# 
# detach(transducin)
# ## End(Not run)
Run the code above in your browser using DataLab