- 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/
我开始在 Ethereum blockchain 上了解如何开发智能合约以及如何写 web-script用于与智能合约交互(购买、销售、统计......)我得出了该怎么做的结论。我想知道我是否正确理解
我正在 UIView 中使用 CATransform3DMakeRotation,并且我正在尝试进行 45º,变换就像向后放置一样: 这是我拥有的“代码”,但显然没有这样做。 CATransform3
我目前正在测试 WebRTC 的功能,但我有一些脑逻辑问题。 WebRTC 究竟是什么? 我只读了“STUN”、“P2P”和其他...但是在技术方面什么是正确的 WebRTC(见下一个) 我需要什么
我在看 DelayedInit在 Scala in Depth ... 注释是我对代码的理解。 下面的 trait 接受一个非严格计算的参数(由于 => ),并返回 Unit .它的行为类似于构造函数
谁能给我指出一个用图片和简单的代码片段解释 WCF 的资源。我厌倦了谷歌搜索并在所有搜索结果中找到相同的“ABC”文章。 最佳答案 WCF 是一项非常复杂的技术,在我看来,它的文档记录非常少。启动和运
我期待以下 GetArgs.hs打印出传递给它的参数。 import System.Environment main = do args main 3 4 3 :39:1: Coul
private int vbo; private int ibo; vbo = glGenBuffers(); ibo = glGenBuffers(); glBindBuffer(GL_ARRAY_
我正在尝试一个 for 循环。我添加了一个 if 语句以在循环达到 30 时停止循环。 我见过i <= 10将运行 11 次,因为循环在达到 10 次时仍会运行。 如果有设置 i 的 if 语句,为什
我正在尝试了解 WSGI 的功能并需要一些帮助。 到目前为止,我知道它是一种服务器和应用程序之间的中间件,用于将不同的应用程序框架(位于服务器端)与应用程序连接,前提是相关框架具有 WSGI 适配器。
我是 Javascript 的新手,我正在尝试绕过 while 循环。我了解它们的目的,我想我了解它们的工作原理,但我在使用它们时遇到了麻烦。 我希望 while 值自身重复,直到两个随机数相互匹配。
我刚刚偶然发现Fabric并且文档并没有真正说明它是如何工作的。 我有根据的猜测是您需要在客户端和服务器端都安装它。 Python 代码存储在客户端,并在命令运行时通过 Fabric 的有线协议(pr
我想了解 ConditionalWeakTable .和有什么区别 class ClassA { static readonly ConditionalWeakTable OtherClass
关闭。这个问题需要更多focused .它目前不接受答案。 想改善这个问题吗?更新问题,使其仅关注一个问题 editing this post . 5年前关闭。 Improve this questi
我还没有成功找到任何可以引导我理解 UIPickerView 和 UIPickerView 模型的好例子。有什么建议吗? 最佳答案 为什么不使用默认的 Apple 文档示例?这是来自苹果文档的名为 U
我在看foldM为了获得关于如何使用它的直觉。 foldM :: Monad m => (a -> b -> m a) -> a -> [b] -> m a 在这个简单的例子中,我只返回 [Just
答案What are _mm_prefetch() locality hints?详细说明提示的含义。 我的问题是:我想要哪一个? 我正在处理一个被重复调用数十亿次的函数,其中包含一些 int 参数。
我一直在读这个article了解 gcroot 模板。我明白 gcroot provides handles into the garbage collected heap 然后 the handle
提供了一个用例: 流处理架构;事件进入 Kafka,然后由带有 MongoDB 接收器的作业进行处理。 数据库名称:myWebsite集合:用户 并且作业接收 users 集合中的 user 记录。
你好 我想更详细地了解 NFS 文件系统。我偶然发现了《NFS 图解》这本书,不幸的是它只能作为谷歌图书提供,所以有些页面丢失了。有人可能有另一个很好的资源,这将是在较低级别上了解 NFS 的良好开始
我无法理解这个问题,哪个更随机? rand() 或: rand() * rand() 我发现这是一个真正的脑筋急转弯,你能帮我吗? 编辑: 凭直觉,我知道数学答案是它们同样随机,但我忍不住认为,如果您
我是一名优秀的程序员,十分优秀!