Homec4science
Diffusion scATAC-seq (master) Browse

scATAC-seq/scripts/10xgenomics_PBMC_5k_motifs_classification_5/classification_sp1_motif.sh
411e7adaffa4master

Last Change
Raw File

classification_sp1_motif.sh
View Options

# some paths
## directories
results_dir='results/10xgenomics_PBMC_5k_motifs_classification_5'
data_dir='data/10xgenomics_PBMC_5k_motifs'
## input
file_mat_seq="$data_dir/sp1_motifs_10e-7_sequences.mat"
file_mat_seq_rmsk="$data_dir/sp1_motifs_10e-7_sequences_rmsk.mat"
## file with seeds
file_seed=$results_dir'/sp1_motifs_10e-7_seed.txt'
file_seed_rmsk=$results_dir'/sp1_motifs_10e-7_seed_rmsk.txt'
mkdir -p $results_dir
touch $file_seed
touch $file_seed_rmsk
# parameters
n_iter='20'
n_shift='771'
n_core=32
# sequences
for k in 1 2 3 4 5 6 7 8 9 10
do
# all sequences
seed=$(< /dev/urandom tr -dc _A-Z-a-z-0-9 | head -c${1:-15};echo)
file_prob=$results_dir/'sp1_motifs_10e-7_open_bin1bp_sequences_'$k'class_prob.mat4d'
file_mod=$results_dir/'sp1_motifs_10e-7_sequences_'$k'class_model.mat'
echo "$file_prob $seed" >> $file_seed
bin/EMSequence --seq $file_mat_seq --class $k --shift $n_shift --flip --iter $n_iter --seed $seed --thread $n_core --out $file_prob
bin/ProbToModel --seq $file_mat_seq --prob $file_prob --thread $n_core 1> $file_mod
# repeat masked sequences
seed=$(< /dev/urandom tr -dc _A-Z-a-z-0-9 | head -c${1:-15};echo)
file_prob=$results_dir/'sp1_motifs_10e-7_open_bin1bp_sequences_rmsk_'$k'class_prob.mat4d'
file_mod=$results_dir/'sp1_motifs_10e-7_sequences_rmsk_'$k'class_model.mat'
echo "$file_prob $seed" >> $file_seed
bin/EMSequence --seq $file_mat_seq_rmsk --class $k --shift $n_shift --flip --iter $n_iter --seed $seed --thread $n_core --out $file_prob
bin/ProbToModel --seq $file_mat_seq_rmsk --prob $file_prob --thread $n_core 1> $file_mod
done
c4science · Help