Question: Seqaln Function Error While Using Bio3D Package In R
0
gravatar for Peter
7.2 years ago by
Peter0
Peter0 wrote:

Hello everyone,

I am using a bio3d package from R for sequence alignment. But I am getting some errors while I was running seqaln function, which calls for the installed MUSCLE program. I got the following example from Bio3D tutorial.

pdb <- read.pdb("1bg2")  

s1 <- aa321(pdb$seqres)
s2 <- seq.pdb(pdb)
aln <- seqaln( seqbind(s1, s2), id=c("seqres","atom") )

muscle -in C:\Users\Peter \AppData\Local\Temp\RtmpaGTZAD\filee116a23 -out aln.fa  -
seqtype proteinError in file(file, "r") : cannot open the connection
In addition: Warning message:
In file(file, "r") : cannot open file 'aln.fa': No such file or directory

I don't understand about this error. your help will be much appreciated.

R sequence alignment • 2.1k views
ADD COMMENTlink modified 7.2 years ago by Neilfws48k • written 7.2 years ago by Peter0
2
gravatar for Neilfws
7.2 years ago by
Neilfws48k
Sydney, Australia
Neilfws48k wrote:

R cannot read the alignment file aln.fa which should have been generated by muscle. This is because the file was not written to disk. There are 2 possible reasons:

  1. muscle did not run (R cannot find it in your path)
  2. muscle did run, but you do not have permission to write file aln.fa in your current working directory

I was able to reproduce the error on a Linux machine by running R from a directory where I did not have write permission. However, before the error, I saw:

file7c082f8b0189 2 seqs, max length 325, avg  length 324
00:00:00    12 MB(-4%)  Iter   1  100.00%  K-mer dist pass 1
00:00:00    12 MB(-4%)  Iter   1  100.00%  K-mer dist pass 2
00:00:00    14 MB(-4%)  Iter   1  100.00%  Align node       
00:00:00    14 MB(-4%)  Iter   1  100.00%  Root alignment

Since you did not mention seeing these lines, I suspect that muscle did not run at all.

I note that you are working in Windows; unfortunately I have no idea how to add muscle to the path for this OS. Perhaps someone else can help with that.

ADD COMMENTlink written 7.2 years ago by Neilfws48k
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 2.3.0
Traffic: 2104 users visited in the last hour