data(transducin)
attach(transducin, warn.conflicts=FALSE)
# Ignore gap containing positions
gaps.pos <- gap.inspect(pdbs$xyz)
# PCA
pc.xray <- pca.xyz(pdbs$xyz[, gaps.pos$f.inds])
# Write PC trajectory of pc=1
outfile = tempfile()
a <- mktrj(pc.xray, file = outfile)
outfile
detach(transducin)
Run the code above in your browser using DataLab