controls.py 文件源码

def build_control_IPD_dict( self, motifs, bi_motifs ):
        """

        """
        control_ipds_fn   = glob.glob( "control_ipds.tmp" )
        control_ipds_N_fn = glob.glob( "control_ipdsN.tmp")
        control_kmers_fn  = glob.glob( "control_ipdskmers.tmp")

        if (len(control_ipds_fn)>1 or len(control_ipds_N_fn)>1 or len(control_kmers_fn)>1):
            raise Exception("*** Double check the control files. There should not be multiples for a file type.")

        control_means,not_found = self.chunk_control_matrices(control_ipds_fn[0], control_ipds_N_fn[0], control_kmers_fn[0])

        if not_found > 0:
            logging.info("")
            logging.warning("WARNING: could not find sufficient instances (>=%s) for %s motifs (out of %s total) in control data!" % (self.opts.min_motif_count, not_found, (len(motifs)+len(bi_motifs))))
            logging.warning("   * If this is alarming, try reducing --min_motif_count or increasing --N_reads, although you just might not have those motifs in your reference sequence.")

        logging.info("")
        logging.info("Writing control data to a pickled file: %s" % self.opts.control_pkl_name)
        pickle.dump( control_means, open( self.opts.control_pkl_name, "wb" ) )

        return control_means