My Research Notes

This page is read only. You can view the source, but not change it. Ask your administrator if you think this is wrong.

===check the printing quota===
  lpquota

===convert tif to eps and pdf===

  convert -compress zip IMAGE.tif eps2:IMAGE.eps//this gives compressed version
  tiff2ps -e figS5d.tif > figS5d.eps//the output can be very large
  epstopdf figS5a.eps >figS5a.pdf
  
===convert pdf to tiff===

  gs -dNOPAUSE -q -r<resolution> -sDEVICE=tiff12nc(for color) -dBATCH -sOutputFile=Zuo-cssp.tiff Zuo-cssp.pdf
  
  === Merge Mappability Files ===

<code>
cd ...

for chr in 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 X Y;
do
 awk -v i=$chr '{print "chr"i"\t"$1"\t"$2}' "chr"$chr"_fragL200_bin200.txt" > "/scratch/czuo/power/example/chr"$chr"_tmp.txt"
done

cd /scratch/czuo/power/example/
rm -f hg19_map.txt
cat chr*_tmp.txt > hg19_map.txt
rm -f chr*_tmp.txt
</code>

=== Subsampling with Loop ===

<code>
for ((i=1;i<11;i=i+1))
do
    for ((j=0;j<=19;j++))
    do
	valc=$((i*2000+j))
#	valt1=$((valc + 100))
	valt2=$((valc + 200))
        prob=$(echo "scale=2;$i/50"|bc)
perl random_sampling.pl $valc $prob run62_WT-O2_B_Input_32bp_eland_Umatch.txt
#perl random_sampling.pl $valt1 $prob run62_Fnr_IP_32bp_eland_Umatch.txt
perl random_sampling.pl $valt2 $prob run101_Fnr_IP_32bp_eland_Umatch.txt
    done
done
</code>

My Research Notes

User Tools

Site Tools

Page Tools