# NOT RUN {
## Fetch PDB files and split to chain A only PDB files
ids <- c("1a70_A", "1czp_A", "1frd_A", "1fxi_A", "1iue_A", "1pfd_A")
raw.files <- get.pdb(ids, path = "raw_pdbs")
files <- pdbsplit(raw.files, ids, path = "raw_pdbs/split_chain")
## Sequence Alignement, and connectivity check
pdbs <- pdbaln(files)
cons <- inspect.connectivity(pdbs)
## omit files with missing residues
trim.pdbs(pdbs, row.inds=which(cons))
# }
Run the code above in your browser using DataLab