пʼятниця, 21 листопада 2014 р.

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=1Npilog2pi

Немає коментарів:

Дописати коментар