cluster_ctcf_mnase_k562.sh
faf6f1d50a68master

cluster_ctcf_mnase_k562.sh
View Options


			results_dir='results/bulk_sequencing'
			data_dir='data/bulk_sequencing/'

			mkdir -p $results_dir

			file_mnase=$data_dir'/ctcf_mnase_k562.mat'
			file_seed=$results_dir'/ctcf_mnase_k562_seed.txt'
			n_iter='20'
			n_shift='21'
			seeding='random'
			n_core=5

			for k in 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15
			do
			seed=$(< /dev/urandom tr -dc _A-Z-a-z-0-9 \| head -c${1:-15};echo)
			file_prob=$results_dir/'ctcf_mnase_k562_'$k'class_prob.mat4d'
			file_ref=$results_dir/'ctcf_mnase_k562_'$k'class_ref.mat'
			file_aic=$results_dir/'ctcf_mnase_k562_'$k'class_aic.txt'
			echo "$file_prob $seed" >> $file_seed
			bin/ChIPPartitioning --data $file_mnase --class $k --shift $n_shift --flip --iter $n_iter --seeding $seeding --seed $seed --parallel $n_core > $file_prob
			bin/probToRef --data $file_mnase --prob $file_prob --parallel $n_core 1> $file_ref 2> $file_aic
			done