Another useful view of our data is to compute for each strain the average log ratio between the ChIP signal and the NUC signal, so that any acetylated nucleosome should exhibit positive values while non acetylated nucleosomes should exhibit negative values. To do so, we will use the option --substract
of NMtdb
as the intensities are already log transformed during the normalization procedure:
rm
$NMtdb -i Data/BY_S288c/Array/NucH3K14_norm.txt \ -o Data/BY_S288c/Array/NucH3K14_log.db \ -c Data/BY_S288c/Array/NucH3K14_log.dbconf \ --substract