01數學模型到底是什麼?
非線性回歸的關鍵在於將模型擬合至數據中,這就提出了一個問題:什麼是模型?
數學模型是對物理、化學或生物狀態或過程的描述。使用模型可幫助思考化學和生理過程或機制,因此你可設計更好的實驗並理解結果。將一個模型擬合到你的數據,將獲得可在模型背景中進行解讀的最佳擬合值。
數學模型既不是假設,也不是理論。與科學假說不同,模型不能通過實驗直接驗證。因為所有模型既真又假……模型的驗證並非說它是「真的」,而是它產生了與重要問題相關的良好的可檢驗假設。——R. Levins, Am.,《科學家》54:421-31,1966
使用模型的目的不一定是完美描述系統。一個完美模型可能有太多有用的參數。相反,目的是找到儘可能簡單的模型來描述自己的實驗體系。我們希望一個模型足夠簡單,這樣就能使模型擬合到數據,但又得足夠複雜,以與的數據很好地擬合併給到一些參數,以幫助自己理解體系,得出有效的科學結論,並設計新的實驗。
我們一起來了解數學模型的作用概念,下文簡要描述了三個常用模型:
作為濃度函數的光密度背景:比色化學分析基於一個簡單的原理。向樣本中加入適當的反應物,引發化學反應,反應產物是有色物質。終止反應時,有色產物的濃度與待檢測的物質的初始濃度成正比。
模型:由於光密度與有色物質的濃度成正比,光密度也與所檢測的物質的濃度成正比。
現實檢查:從數學上講,該方程適用於任何x值。但結果僅當某些值下才有意義。
1)負的X值沒有意義,因為濃度不能是負值;
2)反應不再受物質濃度限制時,該模型可能在高物質濃度下失效;
3)如果溶液變得非常暗(光密度非常高),以至於幾乎沒有光到達檢測器,則模型在高濃度下也可能失敗。此時,儀器的噪聲可能超過信號;
一個模型只適用於一定範圍的值是很正常的。只需要意識到這些限制,不要試圖使用超出其有用範圍的模型。
指數衰減每當某事發生的速度與剩餘的數量成正比時,使用指數方程。例如,配體與受體分離、放射性同位素的衰變和藥物的代謝。表示為微分方程:
將微分方程轉換成在可定義任何時間的Y值的模型需要一些演算。只有一個函數的導數與指數函數Y成正比。同時對微分方程的兩側進行積分,得到一個新的指數方程,該方程將Y定義為X(時間)、速率常數k和Y在時間為0時的值(Y0)的函數。
平衡結合配體與受體相互作用時,或基質與酶相互作用時,結合遵循質量作用定律。
測量結合的量,也就是RL複合物的濃度,將其繪製在Y軸上。改變添加配體的量,我們可假設它與游離配體的濃度(L)相同,因此形成了X軸。一些簡單(但乏味)的代數可得出該方程:
02不會選擇模型怎麼辦?借力軟體程序
非線性回歸目的是將模型擬合到數據中。程序會找到模型中參數的最佳擬合值(可能是速率常數、親和力、受體數量等),從而方便科學地解讀這些最佳擬合值。模型的選擇是一個基於科學的決定。應根據自己對化學或生理學(或遺傳學等)的理解作出判斷,而非僅僅基於圖表的形狀。
一些程序是會自動將數據擬合成數千個方程,然後向你推薦與數據擬合最佳的方程。這樣的程序確實很具有吸引力。但有一個問題,這些程序並不了解你的實驗科學背景。與數據擬合最佳的方程可能不符合有科學意義的模型,這些方程無法解讀參數的最佳擬合值,相應的結果就不大可靠。
如果你的目的是簡單地為模擬或插值創建一條平滑的曲線,那可以讓程序幫你推薦選擇模型。這樣,我們可以不關注參數的值或模型的意義,只需要關注曲線是否很好地擬合數據,且不會產生過大擺動,就OK了。
但,曲線擬合的目的是為了將數據擬合到基於化學、物理或生物原理的模型時,不建議直接使用程序推薦的模型。科學決策還是需要基於實驗體系來決定,而非計算機的算法。
這裡推薦大家了解一下GraphPad Prism提供的理解模型含義的提示,或許讓你的選擇更科學。這裡用描述酶活性的米氏方程作為底物濃度的函數舉例:
Y=Vmax*X/(Km +X)
提示01確信你知道X和Y的含義和單位
在此示例中,Y為酶活性,根據不同的酶,可用不同單位表示。X為底物濃度,單位為摩爾或微摩爾或其他濃度單位。
提示02計算參數的單位
在示例方程中,參數Km加到了X上。只需添加用相同單位表示的東西就有意義了,因此Km必須用與X相同的濃度單位表示。這意味著單位在X/(Km+X)中被抵消,因此Vmax必須用與Y相同的酶活性單位表示。
提示03算出在X極限值下Y的值
因為X為濃度,不能為負。但其可為零。將X=0代入方程,你會發現Y也是零。讓我們也來看看X變得很大時會發生什麼。X相比於Km變大時,分母(X+Km)的值與X非常接近。因此,X/(X+Km)接近1.0,Y接近Vmax。因此,X變得非常大時,模型的圖表必須在Y=Vmax處趨於平穩。
提示04算出在X特定值下Y的值
因為Km是以與X相同的單位表示的,可想像下如果X等於Km會發生什麼?在此情況下,X/(Km+X)等於0.5,因此Y等於Vmax的一半。這意味著Km是導致速度等於最大速度Vmax一半的底物濃度。
提示05用不同的參數值繪製模型
繪製一系列具有不同參數值的曲線可幫助我們直觀了解參數的含義。如需使用Prism實現這一點,可以使用分析「創建一系列理論曲線」