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

 
Немає коментарів:
Дописати коментар