使用 seqRFLP

编程问答 2022-05-27

问题描述

我对 seqRFLP 上的 frag.dat 函数感兴趣，但该函数仅适用于每次运行的单个酶（酶 = X）。

frag.dat(fil = fas,enznames = X,enzdata = enzdata)

这是一个简单的例子

directory <- system.file("extdata",package = "seqRFLP")
path <- file.path(directory,"try.fasta")
fas <- read.fasta(path)
frag.dat(fil = fas,enznames = "AluI",enzdata = enzdata)

结果

enznames recogSite          cutting_Site
>try     AluI     AG'CT 55,94,244,256,343,553

但是如果我想要一次运行多种酶，那么

frag.dat(fil = fas,"TaqI",enzdata = enzdata)
Error in frag.dat(fil = fas,enzdata = enzdata) : 
  unused argument ("TaqI")

是否可以同时运行多种酶？提前致谢。

解决方法

暂无找到可以解决该程序问题的有效方法，小编努力寻找整理中！

如果你已经找到好的解决方法，欢迎将解决方案带上本链接一起发送给小编。

小编邮箱:dio#foxmail.com (将#修改为@）

dataframeenzymerestrictionsequence