view visualize_pc.R @ 4:4354d92dc1aa draft

planemo upload for repository https://github.com/galaxycomputationalchemistry/galaxy-tools-compchem/tree/master/tools/bio3d commit 3ff06e3182c3a1546ea0a3b29e0d4383e12169e1
author chemteam
date Wed, 03 Apr 2019 15:45:27 -0400
parents 50a2b87784b5
children
line wrap: on
line source

#!/usr/bin/env Rscript

options(stringAsfactors = FALSE)
args <- commandArgs(trailingOnly = TRUE)

library(bio3d)

dcdfile <- args[1]
pdbfile <- args[2]

dcd <- read.dcd(dcdfile)
pdb <- read.pdb(pdbfile)

method <- args[3]
selection <- args[4]
domain <- args[5]
id <- args[6] 
pcid <- as.integer(id)

pdbout <- args[7]


if (selection == "string") {
    inds <- atom.select(pdb, string = domain)
}
if (selection == "elety") {
    inds <- atom.select(pdb, elety = domain)
}
if (selection == "resid") {
    inds <- atom.select(pdb, resid = domain)
}
if (selection == "segid") {
    inds <- atom.select(pdb, segid = domain)
}
xyz <- fit.xyz(fixed=pdb$xyz, mobile=dcd, fixed.inds=inds$xyz, mobile.inds=inds$xyz)

if (method == "FALSE") {
    pc <- pca.xyz(xyz[,inds$xyz], use.svd=FALSE)
}
if (method == "TRUE") {
    pc <- pca.xyz(xyz[,inds$xyz], use.svd=TRUE)
}

mktrj.pca(pc, pc=pcid, b=pc$au[,pcid], file=pdbout)