当前位置: 首页>>代码示例>>Python>>正文


Python Helper.median方法代码示例

本文整理汇总了Python中Helper.median方法的典型用法代码示例。如果您正苦于以下问题:Python Helper.median方法的具体用法?Python Helper.median怎么用?Python Helper.median使用的例子?那么恭喜您, 这里精选的方法代码示例或许可以为您提供帮助。您也可以进一步了解该方法所在Helper的用法示例。


在下文中一共展示了Helper.median方法的3个代码示例,这些例子默认根据受欢迎程度排序。您可以为喜欢或者感觉有用的代码点赞,您的评价将有助于系统推荐出更棒的Python代码示例。

示例1: doStatistics

# 需要导入模块: import Helper [as 别名]
# 或者: from Helper import median [as 别名]
def doStatistics(idlist, outputcsv):
        seqlengths = Helper.getSequenceLengthsForAccessionsIds(idlist)
        tsvs = Helper.getTsvEntriesForIdList(idlist)
        
        domainCount = 0
        domainCountCutoff = 0
        domainCover = {}
        clanCount = {}
        
        cutoff = rcp.getint("Options", "domainlengthcutoff")
        
        lengths = []
        
        for seq in seqlengths:
            lengths.append(seqlengths[seq])
            domainlength = 0
            for domain in tsvs[seq]:
                domainlength = (domain.alignmentEnd - domain.alignmentStart)
                domainCount += 1
                if domainlength >= cutoff:
                    domainCountCutoff += 1
                if domain.clan not in clanCount:
                    clanCount[domain.clan] = 0
                clanCount[domain.clan] += 1
            domainCover[seq] = float(domainlength) / float(seqlengths[seq])

        if not os.path.exists(outputcsv):
            out = open(outputcsv, "w")
            out.write(";".join(["taxid", "lengthsum", "meanlength", "medianlength", "domaincount", "domaincountcutoff", "meanclans", "medianclans", "meancover", "mediancover"]))
            out.write("\n")
        
        meanClan = Helper.mean(list(clanCount.values()))
        medianClan = Helper.median(list(clanCount.values()))
        
        meanCover = Helper.mean(list(domainCover.values()))
        medianCover = Helper.median(list(domainCover.values()))
        
        lengthSum = str(sum(lengths))
        meanLength = str(Helper.mean(lengths))
        medianLength = str(Helper.median(lengths)) 
        
        out = open(outputcsv, "a")
        out.write(";".join([str(tsvs[tsvs.keys()[0]][0].tax), lengthSum, meanLength, medianLength, str(domainCount), str(domainCountCutoff), str(meanClan), str(medianClan), str(meanCover), str(medianCover)]))
        out.write("\n")
        
        print "lengthSum", lengthSum
        print "mean length", meanLength
        print "median length", medianLength
        print "domain count", domainCount
        print "domain count considering cutoff (", cutoff, ")", domainCountCutoff
        print "mean clans", meanClan
        print "median clans", medianClan
        print "mean coverage", meanCover
        print "median coverage", medianCover
        
        return clanCount
开发者ID:expectopatronum,项目名称:orth-scripts,代码行数:58,代码来源:statsAboutFullSequenceList.py

示例2: open

# 需要导入模块: import Helper [as 别名]
# 或者: from Helper import median [as 别名]
        
        handle.close()
        
        sumA = 0
        for key in inparalogsPerGroupA:
            sumA += inparalogsPerGroupA[key]
            
        sumB = 0
        for key in inparalogsPerGroupB:
            sumB += inparalogsPerGroupB[key]
        
        out = open(outputfile, 'w')
        out.write("Score summary\n")
        out.write(output)
        out.write("Mean " + str(Helper.mean(scores)) + '\n')
        out.write("Median " + str(Helper.median(scores)) + '\n')
        out.write("Max " + str(max(scores)) + '\n')
        out.write("Len " + str(len(scores)) + '\n\n')
        
        out.write("Organism A\n")
        out.write("Inparalgos " + str(sumA) + '\n')
        out.write("Inparalogs per group " + str(float(sumA)/float(len(inparalogsPerGroupA))) + '\n')
        out.write("Organism B\n")
        out.write("Inparalogs " + str(sumB) + '\n')
        out.write("Inparalogs per group " + str(float(sumB)/float(len(inparalogsPerGroupB))) + '\n')
        out.close()

    else:
        print __doc__
        sys.exit(0)
开发者ID:expectopatronum,项目名称:orth-scripts,代码行数:32,代码来源:analyse.py

示例3: addMinimumDomainCountToFullOrthologs

# 需要导入模块: import Helper [as 别名]
# 或者: from Helper import median [as 别名]
def addMinimumDomainCountToFullOrthologs(pairwise):    
    pairwiseCount = {}
    notTheSameA = []
    notTheSameB = []
    theSame = []
    for pair in pairwise:
        accPair = OrthologPair(first = pair.first.accession, second = pair.second.accession)
        if len(pair.first.domains) == len(pair.second.domains):
            theSame.append(len(pair.first.domains))
        else:
            notTheSameA.append(len(pair.first.domains))
            notTheSameB.append(len(pair.second.domains))
        pairwiseCount[accPair] = min(len(pair.first.domains), len(pair.second.domains))
    print len(theSame), "sequence pairs have the same amount of domains with a mean of", Helper.mean(theSame), "and a median of", Helper.median(theSame)
    print len(notTheSameA), "sequence pairs do not have the same amount of domains with a mean of", Helper.mean(notTheSameA), "/", Helper.mean(notTheSameB), "and a median of", Helper.median(notTheSameA), "/" , Helper.median(notTheSameB)
    return pairwiseCount
开发者ID:expectopatronum,项目名称:orth-scripts,代码行数:18,代码来源:compare.py


注:本文中的Helper.median方法示例由纯净天空整理自Github/MSDocs等开源代码及文档管理平台,相关代码片段筛选自各路编程大神贡献的开源项目,源码版权归原作者所有,传播和使用请参考对应项目的License;未经允许,请勿转载。