He calculado la probabilidad de que, para cada par de aminoácidos, la cantidad de hidden stops que albergan se deba al azar. Así que es una especie de p-value.
Es para el genoma completo (las variables son las mismas que ayer)
>> prob_noefecto=gen2prob_noefecto(genes,codigo);
>> subplot(1,2,1)
>> imagesc(prob_noefecto(:,:,1))
>> subplot(1,2,2)
>> imagesc(prob_noefecto(:,:,2))
Sale significativo para todos los pares de aminoácidos. Teniendo en cuenta que ayer salía que para algunos pares de aminoácidos los stops eran incluso menos probables en los genes que en la teoría, hay un error en alguna parte.
EN EFECTO: Hay un error. Ver siguiente post.
Cuaderno de laboratorio de al. 