EXERCISE BREAK: Compute the entropy of the NF-κB motif matrix (reproduced below).
Motifs T C G G G G g T T T t t
c C G G t G A c T T a C
a C G G G G A T T T t C
T t G G G G A c T T t t
a a G G G G A c T T C C
T t G G G G A c T T C C
T C G G G G A T T c a t
T C G G G G A T T c C t
T a G G G G A a c T a C
T C G G G t A T a a C C
На этом задании я очень поленилась. Достала массив Profile из теории, вставила его в блокнот по контрол-ф удалила все нули (единицы удалила вручную), заменила точки на "0.", а двойные пробелы на пробел с запятой. Получила строку ненулевых слагаемых в сумме для рассчета энтропии и загнала это в питон в таком виде:
Получила в итоге 9.916290005356972 и это коррект.
Формула энтропии: H(p1,…,pN)=−∑i=1Npi⋅log2pi
Немає коментарів:
Дописати коментар