R/defineSeqRef.R

Defines functions defineSeqRef

defineSeqRef<- function(type, void){
  
  if (void!=1) print("Provide values (%) per each residue!")
  
  if(type == 1)
    referenceNames<-"ACDEFGHIKLMNPQRSTVWY"
  else if(type == 2)
    referenceNames<-"acgtu"
  else
    referenceNames<-"O"
  
  referenceVec<-vector(mode = "double")
  
  name<-strsplit(referenceNames, split = "")[[1]]
  
  for(i in name){
    if (void == 1)
      referenceVec[i]<-0
    else
      referenceVec[i]<-as.numeric(readline(paste(i, ": ")))
  }
    
  
  reference<-as.data.frame(referenceVec)
  
  reference
  
}

Try the PdPDB package in your browser

Any scripts or data that you put into this service are public.

PdPDB documentation built on Nov. 17, 2017, 5:09 a.m.