基于希爾伯特空間核方法的生物調(diào)控關(guān)系研究.pdf

1、生物醫(yī)學(xué)領(lǐng)域中經(jīng)常需要研究變量間的影響關(guān)系，例如，基因關(guān)聯(lián)的因果關(guān)系問題。本文基于再生核希爾伯特空間獨立性、條件獨立性理論，研究基因變量的調(diào)控關(guān)系。與一些傳統(tǒng)的方法相比，獨立性和條件獨立性更適合于因果性關(guān)系的描述，比如用來描述基因的本質(zhì)調(diào)控關(guān)系。通過在再生核希爾伯特空間（Reproducing Kernel Hilbert Space--RKHS）中定義類歐式空間的統(tǒng)計量，把對原空間統(tǒng)計特性描述的維數(shù)擴(kuò)展到無窮維，從而可以對分布做更精確

2、的定量描述，例如在RKHS中定義的協(xié)方差算子可以用來描述獨立性和條件獨立性。它是一種基于分布的統(tǒng)計獨立性度量方法，適合于關(guān)系復(fù)雜，數(shù)據(jù)維數(shù)高而樣本點少和高噪聲的生物醫(yī)學(xué)數(shù)據(jù)；它也是一種非參數(shù)的方法，不要求數(shù)據(jù)符合某種特殊的分布，有比較好的適應(yīng)性；另外，核方法能帶來計算上的高效性。對于協(xié)方差算子的度量問題，本文著重介紹了一種標(biāo)準(zhǔn)的度量方法，即Hilbert-Schmidt范數(shù)度量方法，該方法是一種基于RKHS上算子特征譜的度量方法。根據(jù)此

3、度量方法可得到獨立性和條件獨立性的判別準(zhǔn)則。特別地，在使用基于Hilbert-Schmidt范數(shù)的統(tǒng)計獨立性度量準(zhǔn)則（HSIC）進(jìn)行獨立性檢驗時，可以利用獨立性假設(shè)下的近似分布--Gamma分布，與普通的Monte Carlo重采樣法相比能極大地提高計算效率。在對該問題進(jìn)行仿真實驗表明，HSIC方法能很好地檢測出獨立性或非獨立性的關(guān)系。使用Gamma分布進(jìn)行檢驗的結(jié)果與permutation檢驗的結(jié)果很吻合，說明在獨立性假設(shè)下Gamma

4、分布對原分布的近似效果很好，并且它的參數(shù)可以直接從樣本數(shù)據(jù)中估計出來，無需像permutation檢驗?zāi)菢幼龃罅侩S機(jī)實驗來獲取分布數(shù)據(jù)。
　　對于時間序列數(shù)據(jù)的因果性推斷問題，本文在獨立性和條件獨立性理論下對Granger因果性方法進(jìn)行了擴(kuò)展，使其能適用于非線性的情況，給出了時間序列數(shù)據(jù)因果性推斷方法。對于大樣本的問題，給出了基于Subsampling的檢驗方法，該方法能顯著地降低整體的計算量，并且能提高檢驗的可靠性。另外，還給出

5、了更一般意義下的系統(tǒng)辨識問題的算法步驟，并對算法進(jìn)行了實驗測試。
　　最后，使用HSIC方法對Dialogue for Reverse Engineering Assessments and Methods第3次競賽項目（DREAM3）中的基因調(diào)控網(wǎng)絡(luò)問題的數(shù)據(jù)進(jìn)行了計算，結(jié)果顯示 HSIC方法在計算準(zhǔn)確性上和計算效率上都有較好的表現(xiàn)。DREAM3中的大腸桿菌（E.coli）基因調(diào)控網(wǎng)絡(luò)結(jié)構(gòu)辨識問題具有關(guān)系復(fù)雜（基因連接成網(wǎng)絡(luò)化）