powered by
rmsf(xyz, average=FALSE)
read.dcd
fit.xyz
read.fasta.pdb
attach(transducin) # Ignore Gaps gaps <- gap.inspect(pdbs$ali) r <- rmsf(pdbs$xyz) plot(r[gaps$f.inds], typ="h", ylab="RMSF (A)") detach(transducin)
Run the code above in your browser using DataLab