def test_simulation(self): raw_motif_path = '../Example/jaspar_vertebrata.txt' #raw_motif_path = '/home/jiawei/yijing_motif/motif/CArG_motif_matrix_t.txt' genome_db_path = '/home/jiawei/.MotifScan/genome/tair10' sample_number = 100000 background = pd.read_pickle('%s/background'%genome_db_path)['background'] chromosome_size = pd.read_pickle('%s/chromosome_size'%genome_db_path) motif_table = motif.load_motif_matrix(raw_motif_path) motif_table = motif.compute_max_score(motif_table,background) motif_table = motif.simulation(motif_table,sample_number,genome_db_path,chromosome_size,background) print motif_table["score_cutoff_3"]
def test_load_motif_matrix(self): matrix_file = '/home/jiawei/MAmotif/motifscan_pkg/MotifScan/Example/jaspar_all_motif.txt' motif.load_motif_matrix(matrix_file)