Question: setting database.path for extracting PSSM Feature in BioSeqClass
0
gravatar for greensandag
3.6 years ago by
Taiwan
greensandag0 wrote:

How can i use Feature PSSM in BioSeqClass packages?. i tried run this program:

 

if(interactive()){
  file = file.path(path.package("BioSeqClass"), "example", "acetylation_K.fasta")
  tmp = readAAStringSet(file)
  proteinSeq = as.character(tmp)
  ## Need "blastpgp" program and a formated database. Database can be formated by "formatdb" program.
  PSSM1 = featurePSSM(proteinSeq[1:2], start.pos=rep(1,2), stop.pos=rep(10,2), psiblast.path="blastpgp", database.path="./result1.fasta")
}

 

but, i got  error messages:

Error in featurePSSM(proteinSeq[1:2], start.pos = rep(1, 2), stop.pos = rep(10,  : 
  Do not have formated database file: ./result1.fasta.pin

 

How can i setting database.path? if my directory for database : C:\BLAST\db

ADD COMMENTlink written 3.6 years ago by greensandag0
Please log in to add an answer.

Help
Access

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 16.09
Traffic: 85 users visited in the last hour