随便看了一下,也不一定对;发现错了再修改。
话说基因表达是随时间变化的,比如有n个基因,在时间t表达情况A(t)就是一个n维向量,也可以看成是1*n的矩阵
基因之间是相互影响(调控)的,每个基因都可能对其他所有的基因调控。调控的影响力可以看成是一个n*n的矩阵X。
那么A(t+1)=A(t)*X
如果已知A(0),A(1),...,A(t),尽量推断出一个靠谱的、错误率较低的X,就是基因调控网络的构建。
当然,这个问题可以复杂化,就是A(t+1)不等于A(t)*X,而是A(t)*X后再通过sigmoid函数变换才是A(t+1);
这个问题也可以简单化,就是把A(t)和X都变成0,1;那么这个网络就称为布尔网络。X矩阵中1就代表调控,可以画出一条有向边,0代表不调控,没有边,这就画出了调控网络图。