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)