基因转录调控网络推导方法——伪逆矩阵模型
在基因调控网络推导中,使用到的
基因芯片数据
通常具有
样本个数少(通常小于10)而
基因数目大(通常大于1000)的局限性,也就是说,实验样本个数远远小于基因的个数。另外,调控矩阵具有较强的
稀疏性,即每个基因只被少量的转录因子调控,而每个转录因子只调控少量的基因。
伪逆矩阵法
就是为了解决基因调控网络线性模型中
基因个数远远大于样本个数的问题而发展的一种有效的基因调控网络推导方法。Andrecut和Kauffman给出了推导基因调控网络的伪逆矩阵模型,并用贪婪算法进行了计算,得到了较好的效果。
在建模过程中,考虑含有N个基因在M次干扰实验中的基因调控网络,每一次干扰实验可以用一个N维动力系统建模
dxnmdt=∑i=1Nwnixim−μnxnm+unm,n=1,2,...,Nm=1,2,...,M(1)\frac{\mathrm{d} x_{nm}}{\mathrm{d} t} = \sum_{i=1}^{N}w_{ni}x_{im}-\mu _nx_{nm+u_{nm}}, n=1,2,...,N m=1,2,...,M (1)