bioinformatics 相关问题

仅将此标记用于与生物信息学相关的编程相关问题。其他问题不属于此处,但可能是https://bioinformatics.stackexchange.com/上的主题。有关更多信息,请参阅标记维基。

R:根据来自另一个data.table的条件“标记”一行

我有一个超过100,000行的data.table(A)。共有3列。 chrom start end 1:chr1 6484847 6484896 2:chr1 6484896 6484945 3:chr1 ...

回答 2 投票 1

在Django和Django REST Framework中使用保留字“class”作为字段名称

问题描述分类学是在共同特征的基础上定义和命名生物有机体组的科学。生物被组合成分类群(单数:...

回答 3 投票 9

如何修改Smith和Waterman perl脚本以加快速度?

我找到了这个perl脚本,但我有太多的序列需要分析。我想知道是否可以对其进行优化?我在它上面推出了NYTProf,看到了部分“计算匹配分数”,“计算差距...

回答 2 投票 0

在Python中替换核苷酸

问题是首先反转字符串。然后将“T”换成“A”,反之亦然,将“G”换成“C”,反之亦然。我没有问题扭转字符串。问题是将A替换为A和副......

回答 5 投票 -1

使用命名列表重命名列表列表

所以我正在使用包含其他列表的列表,其结构如下:库(图)库(RBGL)库(Rgraphviz)显示(树)$`SO:0001968` $`SO:0001968` $`SO: 0001622` $`SO:0001968` $`...

回答 2 投票 4

无法解析FASTA文件中的序列

如何从序列中删除'> gi | 2765658 | emb | Z78533.1 | CIZ78533 C.irapeanum 5.8S rRNA基因和ITS1和ITS2 DNA \ n'等ID?我有这个代码:open('sequence.fasta','r')为f:...

回答 1 投票 1

使用perl脚本计算核苷酸频率[关闭]

我在下面有这个perl脚本来计算序列长度及其频率以及核苷酸频率(A,T,G和C)。此脚本适用于具有大量序列的文件,但它...

回答 1 投票 -4

如何绘制R中每个染色体读数的密度直方图

我在这里有这个数据文件https://www.dropbox.com/s/hei2fa25r9ezvkx/myfile.csv?dl=0我想绘制chr1到chr16起始位置的密度图。所以我想把X轴分成chr1,chr2,......

回答 1 投票 -3

使用bash在每列中组合具有相同名称的行

我有一个如下文件(但有52列和4,000行):1NA2 1NB2 2RA2 2RB2 Vibrionaceae 0.22 0.25 0.36 1.02 Bacillaceae 2.0 1.76 0.55 0.23 ...

回答 1 投票 1

安装Bioconductor软件包自动化

我试图在R脚本的Ubuntu 16.04上自动安装Bioconductor及其不同的软件包。但是,每个安装都需要交互来说 - 是,是的我想安装它和-...

回答 1 投票 1

将一个数据框中的两列与另一个数据框中的两列匹配,与顺序无关

以两个数据帧为例:df1 = data.frame(V1 = c('JOHN','BRIAN','KATE','ERIC','CURT','ZACH'),V2 = c('ABIGAIL', 'ANDY','GEORGE','JOHN','MARY','FRANKLIN'),V3 = seq(1,6,1),V4 = seq(1,...

回答 3 投票 1

在Python中聚类相似字符串的算法?

我正在编写一个目前包含多个DNA序列列表的脚本(每个列表都有不同数量的DNA序列),我需要根据Hamming对每个列表中的序列进行聚类...

回答 1 投票 1

Python将RNA seq转换成单字母氨基酸序列

在编写将给定的RNA核苷酸序列转换成氨基酸序列的代码时,我需要一些帮助。我目前已经使用了2个字典:氨基酸密码子之一和......

回答 1 投票 -2

在python循环中打印来自交替文件的行

我试图使用python在两个单独的文件中找到感兴趣的四行块,然后按受控顺序打印出一些这些行。下面是两个输入文件和一个例子......

回答 3 投票 1

如何在Python中打印文件中的某些行

我想帮助弄清楚如何在.txt文件中只打印出给定数量的行。我创建了一个带有2个输入参数的函数文件(x,y),第一个'x'是文件,...

回答 2 投票 0

Python在上下文后打印行

如何在我对使用python感兴趣的上下文后打印两行。 Example.fastq @ read1 AAAGGCTGTACTTCGTTCCAGTTG +'(''%$'))%**)2 +'。(&&'/ 5- @ read2 CTGAGTTGAGTTAGTGTTGACTC +)(+ - 0 -...

回答 2 投票 1

如何使用tanimoto距离形成网络

我有这样的列数据-1列Column-2 Similarity_dist化学-1化学-2 0.6化学-1化学-3 0.1化学-1化学-4 0.7化学-1化学-5 0.3化学-2化学-...

回答 1 投票 3

perl从2个数组中提取常见元素(fastq文件中的常见序列)

我有2个文件(paired_end),其中包含来自同一生物体的Fastq格式的读取(File_R1.fastq和File_R2.fastq)。我想用bwa和samtools(生物信息学)确定深度覆盖率,因为......

回答 2 投票 1

如何使用lapply来计算r中列表中的唯一值

关于如何从数据帧中计算唯一值,我之前已经问了一个类似的问题,但是我需要使用“lapply”,因为我之前使用的方法不起作用或者我无法得到它...

回答 2 投票 1

如何使用python或R将三个字母的氨基酸代码转换为一个字母代码?

我有一个fasta文件,如下所示。我想将三个字母的代码转换为一个字母代码。我怎么能用python或R做到这一点? > 2ppo ARGHISLEULEULYS> 3oot METHISARGARGMET所需...

回答 11 投票 10

© www.soinside.com 2019 - 2024. All rights reserved.