- html - 出于某种原因,IE8 对我的 Sass 文件中继承的 html5 CSS 不友好?
- JMeter 在响应断言中使用 span 标签的问题
- html - 在 :hover and :active? 上具有不同效果的 CSS 动画
- html - 相对于居中的 html 内容固定的 CSS 重复背景?
我需要一些帮助来理解我想使用的函数,但我不完全确定它的某些部分是做什么的。我知道该函数是通过读取 Fasta 文件来创建字典。据我了解,这应该为最终扩展重叠群(重叠的 dna 序列)生成前和后缀词典。代码:
def makeSuffixDict(reads, lenSuffix = 20, verbose = True):
lenKeys = len(reads[0]) - lenSuffix
dict = {}
multipleKeys = []
i = 1
for read in reads:
if read[0:lenKeys] in dict:
multipleKeys.append(read[0:lenKeys])
else:
dict[read[0:lenKeys]] = read[lenKeys:]
if verbose:
print("\rChecking suffix", i, "of", len(reads), end = "", flush = True)
i += 1
for key in set(multipleKeys):
del(dict[key])
if verbose:
print("\nCreated", len(dict), "suffixes with length", lenSuffix, \
"from", len(reads), "Reads. (", len(reads) - len(dict), \
"unambigous)")
return(dict)
附加信息:reads = readFasta("smallReads.fna", verbose = True)
函数是这样调用的:
if __name__ == "__main__":
reads = readFasta("smallReads.fna", verbose = True)
suffixDicts = makeSuffixDicts(reads, 10)
smallReads.fna 文件包含碱基串 (Dna):
"> read 1
TTATGAATATTACGCAATGGACGTCCAAGGTACAGCGTATTTGTACGCTA
"> read 2
AACTGCTATCTTTCTTGTCCACTCGAAAATCCATAACGTAGCCCATAACG
"> read 3
TCAGTTATCCTATATACTGGATCCCGACTTTAATCGGCGTCGGAATTACT
以下是我不明白的部分:
lenKeys = len(reads[0]) - lenSuffix
值 [0] 是什么意思?据我了解,“len”返回列表中的元素数。为什么“阅读”会自动成为一个列表?编辑:似乎可以将 Fasta 文件声明为列表。有人可以证实吗?
if read[0:lenKeys] in dict:
这是否意味着“从 0 到‘lenKeys’”?仍然对值(value)感到困惑。在另一个函数中有类似的行:if read[-lenKeys:] in dict:
“-”是做什么的?
def makeSuffixDict(reads, lenSuffix = 20, verbose = True):
这里的参数我不明白:reads
怎么能是参数呢?除了从 len(reads[0])
中减去的值之外,此函数上下文中的 lenSuffix = 20
是什么?什么是冗长?我读过有关忽略空格的“详细模式”的信息,但我从未见过它被用作参数,后来又被用作变量。
最佳答案
你问题的语气让我觉得你混淆了程序特性(len
、函数等)和原始程序员定义的东西( 的类型)阅读
,冗长
等)。
def some_function(these, are, arbitrary, parameters):
pass
这个函数定义了一堆参数。除了我暗中赋予它们的值(value)外,它们没有任何意义。例如,如果我这样做:
def reverse_string(s):
pass
s
应该是字符串吧?在您的示例中,我们有:
def makeSuffixDict(reads, lenSuffix = 20, verbose = True):
lenKeys = len(reads[0]) - lenSuffix
...
从这两行我们可以推断出一些事情:
lenSuffix
是一个 int
,而 verbose
是一个 bool
(来自它们的默认参数)reads
可以被索引(字符串?列表?元组?)reads
中的项目有长度(字符串?列表?元组?)由于 Python 是动态类型的,所以这是目前我们所能知道的关于函数的信息。其余的将通过其文档或调用方式进行解释。
就是说:让我按顺序回答您所有的问题:
some_object[0]
is grabbing the first item in a container.[1,2,3][0] == 1
,"Hello, World!"[0] == "H"
. This is called indexing, and is governed by the__getitem__
magic method
len
is a built-in function that returns the length of an object. It is governed by the__len__
magic method.len('abc') == 3
, alsolen([1, 2, 3]) == 3
. Note thatlen(['abc']) == 1
, since it is measuring the length of the list, not the string inside it.
reads
is a parameter. It is whatever the calling scope passes to it. It does appear that it expects a list, but that's not a hard and fast rule!
Slicing is doing
some_container[start_idx : end_idx [ : step_size]]
. It does pretty much what you'd expect:"0123456"[0:3] == "012"
. Slice indexes are considered to be zero-indexed and lay between the elements, so[0:1]
is identical to[0]
, except that slices return lists, not individual objects (so'abc'[0] == 'a'
but'abc'[0:1] == ['a']
). If you omit either start or end index, it is treated as the beginning or end of the string respectively. I won't go into step size here.Negative indexes count from the back, so
'0123456'[-3:] == '456'. Note that
[-0]is not the last value,
[-1]is. This is contrasted with
[0]` being the first value.
Because the function is defined as
makeSuffixDict(reads, ...)
. That's what a parameter is.
Looks like it's the length of the expected suffix!
verbose
?
verbose
has no meaning on its own. It's just another parameter. Looks like the author included theverbose
flag so you could get output while the function ran. Notice all theif verbose
blocks seem to do nothing, just provide feedback to the user.
关于python - 了解 Python 函数,我们在Stack Overflow上找到一个类似的问题: https://stackoverflow.com/questions/30584108/
我正在处理一组标记为 160 个组的 173k 点。我想通过合并最接近的(到 9 或 10 个组)来减少组/集群的数量。我搜索过 sklearn 或类似的库,但没有成功。 我猜它只是通过 knn 聚类
我有一个扁平数字列表,这些数字逻辑上以 3 为一组,其中每个三元组是 (number, __ignored, flag[0 or 1]),例如: [7,56,1, 8,0,0, 2,0,0, 6,1,
我正在使用 pipenv 来管理我的包。我想编写一个 python 脚本来调用另一个使用不同虚拟环境(VE)的 python 脚本。 如何运行使用 VE1 的 python 脚本 1 并调用另一个 p
假设我有一个文件 script.py 位于 path = "foo/bar/script.py"。我正在寻找一种在 Python 中通过函数 execute_script() 从我的主要 Python
这听起来像是谜语或笑话,但实际上我还没有找到这个问题的答案。 问题到底是什么? 我想运行 2 个脚本。在第一个脚本中,我调用另一个脚本,但我希望它们继续并行,而不是在两个单独的线程中。主要是我不希望第
我有一个带有 python 2.5.5 的软件。我想发送一个命令,该命令将在 python 2.7.5 中启动一个脚本,然后继续执行该脚本。 我试过用 #!python2.7.5 和http://re
我在 python 命令行(使用 python 2.7)中,并尝试运行 Python 脚本。我的操作系统是 Windows 7。我已将我的目录设置为包含我所有脚本的文件夹,使用: os.chdir("
剧透:部分解决(见最后)。 以下是使用 Python 嵌入的代码示例: #include int main(int argc, char** argv) { Py_SetPythonHome
假设我有以下列表,对应于及时的股票价格: prices = [1, 3, 7, 10, 9, 8, 5, 3, 6, 8, 12, 9, 6, 10, 13, 8, 4, 11] 我想确定以下总体上最
所以我试图在选择某个单选按钮时更改此框架的背景。 我的框架位于一个类中,并且单选按钮的功能位于该类之外。 (这样我就可以在所有其他框架上调用它们。) 问题是每当我选择单选按钮时都会出现以下错误: co
我正在尝试将字符串与 python 中的正则表达式进行比较,如下所示, #!/usr/bin/env python3 import re str1 = "Expecting property name
考虑以下原型(prototype) Boost.Python 模块,该模块从单独的 C++ 头文件中引入类“D”。 /* file: a/b.cpp */ BOOST_PYTHON_MODULE(c)
如何编写一个程序来“识别函数调用的行号?” python 检查模块提供了定位行号的选项,但是, def di(): return inspect.currentframe().f_back.f_l
我已经使用 macports 安装了 Python 2.7,并且由于我的 $PATH 变量,这就是我输入 $ python 时得到的变量。然而,virtualenv 默认使用 Python 2.6,除
我只想问如何加快 python 上的 re.search 速度。 我有一个很长的字符串行,长度为 176861(即带有一些符号的字母数字字符),我使用此函数测试了该行以进行研究: def getExe
list1= [u'%app%%General%%Council%', u'%people%', u'%people%%Regional%%Council%%Mandate%', u'%ppp%%Ge
这个问题在这里已经有了答案: Is it Pythonic to use list comprehensions for just side effects? (7 个答案) 关闭 4 个月前。 告
我想用 Python 将两个列表组合成一个列表,方法如下: a = [1,1,1,2,2,2,3,3,3,3] b= ["Sun", "is", "bright", "June","and" ,"Ju
我正在运行带有最新 Boost 发行版 (1.55.0) 的 Mac OS X 10.8.4 (Darwin 12.4.0)。我正在按照说明 here构建包含在我的发行版中的教程 Boost-Pyth
学习 Python,我正在尝试制作一个没有任何第 3 方库的网络抓取工具,这样过程对我来说并没有简化,而且我知道我在做什么。我浏览了一些在线资源,但所有这些都让我对某些事情感到困惑。 html 看起来
我是一名优秀的程序员,十分优秀!