| 1 | 1/1 | 返回列表 |
| 查看: 4653 | 回復(fù): 0 | ||
1289622木蟲 (正式寫手)
|
[求助]
關(guān)于房室模型、權(quán)重的選擇
|
|
以下是網(wǎng)上總結(jié)的 關(guān)于房室模型、權(quán)重的選擇要求,不知是否正確,請專家指點(diǎn)。 還有,如何用winnonlin非房室模型進(jìn)行數(shù)據(jù)處理? 模型選擇: 不同模型計算出的AIC值(相同數(shù)據(jù),相同的數(shù)據(jù)擬合方法,相同權(quán)重條件下)AIC越小模型越好;但是有時候不一定能選好最佳模型。AIC第一。 BIC第二。 關(guān)于房室數(shù)的選擇 可根據(jù) F 檢驗和 AIC 值選擇房室數(shù). 對于相同的權(quán), 當(dāng)F 檢驗有顯著意義(P<0.05 或 P<0.01)時, 應(yīng)取 AIC 值較小的房室數(shù), 當(dāng) F 檢驗無顯著意義(P>0.05)時, 則取房室數(shù)少者為宜。模型選擇最好參照AIC 、BIC或HQIC,R2通常不能表征自由度的減少。 一、最佳房室數(shù)確定原則 1.測定值能夠均勻而隨機(jī)地分布在擬合曲線的兩側(cè)。 2.適當(dāng)?shù)厥箽埐钇椒胶停⊿)或加權(quán)殘差平方和(SW)達(dá)到最小。擬合的房室模型中,S或Sw最小者即為所求的房室模型。 二、確定房室數(shù)的具體方法 1.C- t散點(diǎn)圖判斷法 (1)iv后血藥濃度(C)對時間(t)在半對數(shù)坐標(biāo)紙上繪出散點(diǎn)圖,由散點(diǎn)圖形估計房室數(shù)。 如各數(shù)據(jù)點(diǎn)可用一條直線擬合,可初估為一室模型;如圖形在一處或兩處出現(xiàn)轉(zhuǎn)折,血藥濃度呈現(xiàn)先快后慢的衰減曲線,可初估為二室或三室模型,擬合雙指數(shù)或三指數(shù)方程。2)血管外給藥后的房室模型可根據(jù) C-T曲線中吸收后相的曲線形狀加以估計: 如吸收后相曲線為一直線,則可估計屬于口服一室模型,擬合雙指數(shù)方程;如吸收后相的曲線形狀表現(xiàn)為先快后慢的衰減曲線,則可估計屬于口服二室模型,擬合三指數(shù)方程式。 應(yīng)注意的是,一個靜脈注射為二室模型的藥物,如果分布速率并不十分快,其口服吸收曲線可能表現(xiàn)為口服一室模型,這是因為緩慢的分布相為吸收相所掩蓋;如該藥分布十分迅速,則可表現(xiàn)為口服二室模型。散點(diǎn)圖法簡單,但比較粗糙,不夠準(zhǔn)確,需采用以下方法進(jìn)一步確證。 2. 殘差平方和或加權(quán)殘差平方和判斷法 在應(yīng)用最小二乘法擬合藥-時曲線時,需先對數(shù)據(jù)進(jìn)行權(quán)重。原因如下:線性動力學(xué)中,藥代動力學(xué)參數(shù)的計算常將血藥濃度轉(zhuǎn)換為對數(shù)濃度后對時間進(jìn)行直線擬合(log C – t ),在計算中常用最小二乘法,但經(jīng)轉(zhuǎn)換后的數(shù)據(jù)殘差平方和達(dá)最小,不等于原來數(shù)據(jù)的殘差平方和達(dá)最;在測得的一組血藥濃度數(shù)值中,高、低濃度相差較大(如相差兩個數(shù)量級)時,若實驗值與計算值殘差不經(jīng)權(quán)重,則在最小二乘法計算過程中,低濃度數(shù)據(jù)的作用將會被忽視,而實際上低濃度實驗數(shù)據(jù)對曲線的擬合都是十分重要的。 多數(shù)文獻(xiàn)采用1/C2(濃度倒數(shù)的平方)進(jìn)行權(quán)重。 3.擬合優(yōu)度r12值大的為最佳房室模型。 4.F測驗判斷法 Sw1及Sw2分別為第一種和第二種模型的加權(quán)殘差平方和; df為自由度,即各自的實驗數(shù)據(jù)點(diǎn)的數(shù)目減去參數(shù)的數(shù)目; iv一室模型的C-T方程式需確2個參數(shù) (C=C0e-kt ); iv二室: 4個參數(shù) (C=Ae-αt + Be-βt ) ; iv三室: 6個參數(shù) (C=Ae-αt + Be-βt + Ge-γt ); 如某實驗測得12個實驗數(shù)據(jù)點(diǎn),則上述三種模型的df分別為10、8、6; 如算得的F值比相應(yīng)自由度的F界值(Fα,α=0.05,即5%顯著水平)大,也就是說P<0.05,便可認(rèn)為將參數(shù)的數(shù)目從i增至i+1是有意義的,模型2優(yōu)于模型1。Fα是查表得出的Fα(df1-df2,df2)。 . AIC 判斷法 AIC=N㏑(Re)+2P N為實驗數(shù)據(jù)點(diǎn)的數(shù)目; P為擬合的房室模型的相應(yīng)數(shù)學(xué)方程式中參數(shù)數(shù)目(P=2n,n為房室數(shù)); Re為加權(quán)殘差平方和(與Sw含義相同)。 AIC和權(quán)重的選擇有關(guān)。 所以,更換權(quán)重,不能用AIC作為評價指標(biāo),以相關(guān)系數(shù)更好。 同樣,評價兩個不同模型對一組數(shù)據(jù)的擬合效果,應(yīng)當(dāng)選擇相同權(quán)重, 數(shù)據(jù)擬合方法也要一致,才可以以AIC判斷。 在擬合的房室模型中,AIC小為好,AIC最小的數(shù)學(xué)方程式被認(rèn)為是對血藥濃度時程的最佳表達(dá),這種統(tǒng)計學(xué)方法謂之最低AIC測定(minimum AIC estimation,MAICE)。對于MAICE來說,不要求進(jìn)行F測驗及顯著水平測定。如按Re及AIC判斷結(jié)果不一致,而Re相差不大時,以AIC為判定標(biāo)準(zhǔn)。 6.其他 多數(shù)藥物屬于二室模型,其速率常數(shù)具有下列特征: α>β,且相差較大; k12、k21、ke(藥物自中央室消除的一級速率常數(shù))均為正值; k12+k21<20ke。 如α>>β 或 α=β,可視為一室模型加以處理。 屬于前者的藥物,向組織的分布相當(dāng)迅速,分布相比消除相快得多,以致在藥動學(xué)處理中,分布相可以忽略;屬于后者的藥物,分布相和消除相速率幾乎相等,故可認(rèn)為藥物在機(jī)體內(nèi)瞬即達(dá)平衡,因而也可作一室處理。 一般來說,當(dāng)k12+k21≥20ke,二室模型可作一室處理。 應(yīng)說明的是:①房室模型的確定主要取決于該藥本身的化學(xué)結(jié)構(gòu)及藥動學(xué)性質(zhì),然而,由于實驗?zāi)康囊约皩嶒炘O(shè)計(如采血時間,給藥途徑,檢測方法的靈敏度)等的不同,同一藥物也可被判斷為不同的房室模型。Saunders建議,在采集血樣的安排方面,峰濃度前應(yīng)有3個采樣點(diǎn),達(dá)峰時到二倍達(dá)峰時之間應(yīng)不少于2點(diǎn),此后到10倍達(dá)峰時應(yīng)不少于4點(diǎn)。②在具體確定某一藥物的房室數(shù)時,有時用不同方法得出不同的房室數(shù),這時應(yīng)采用上述幾種方法綜合判斷。如石吊蘭素給猴靜注石吊蘭素后,分別擬合一室及二室模型,F(xiàn)測驗值不顯著,但r12值以二室為大,另以兩種模型的計算值與實測值直觀比較,以二室模型的符合程度較佳,故以二室開放模型來描述較為恰當(dāng)。 當(dāng)數(shù)據(jù),權(quán)重都相同時: 更換一下模型,可以以AIC為準(zhǔn); 更換一下數(shù)據(jù)擬合方法,可以以AIC為準(zhǔn); 相關(guān)性任何時間都可以用。 兩個原則:AIC還有模型的簡單程度。參數(shù)越少,AIC越小,參數(shù)越少,同時也說明模型越簡單,AIC包含了模型簡單這個因素。 三 權(quán)重的選擇: 低濃度數(shù)據(jù)對于半衰期影響很大。 一般而言,對于DMPK來說,高濃度的數(shù)據(jù)要比低濃度的數(shù)據(jù)重要。 所以,實際情況中,多數(shù)以1作為權(quán)重。 但是對于蓄積評價來說,會影響到蓄積系數(shù),所以,要以1/y 或者1/yy為權(quán)重。 假如僅僅為了數(shù)據(jù)擬合得漂亮。 我一般選擇1/y,這樣既照顧了高濃度,也照顧了低濃度。 首先比較擬合優(yōu)度值(Goodness of fit), 數(shù)值越小越好, 如數(shù)值相等或相近,再同時結(jié)合最低檢測濃度, 參考最大絕對誤差和相誤差. 一般說來, 權(quán)重為 1 時, 最大絕對誤差較小, 對高濃度數(shù)據(jù)部份擬合較好, 有利于計算表觀分布容積及生物利用度,權(quán)重為 1/C/C 時, 最大相對誤差較小, 對低濃度數(shù)據(jù)部份擬合較好, 有利于計算消除相半衰期, 權(quán)重為 1/C 時, 通常能兼顧高低濃度兩部份數(shù)據(jù), 但最大絕對誤差及相對誤差往往并非最小。 當(dāng)高濃度小于理論值、而低濃度大于理論值時,用1/C或1/C^。 |
找到一些相關(guān)的精華帖子,希望有用哦~
| 1 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 347求調(diào)劑 +11 | 山頂見α 2026-03-25 | 11/550 |
|
|---|---|---|---|---|
|
[考研] 本科211安全工程,初試290分,求調(diào)劑 +3 | 2719846834 2026-03-28 | 3/150 |
|
|
[考研] 張芳銘-中國農(nóng)業(yè)大學(xué)-環(huán)境工程專碩-298 +8 | 手機(jī)用戶 2026-03-26 | 8/400 |
|
|
[考研] 085602 307分 求調(diào)劑 +8 | 不知道叫什么! 2026-03-26 | 8/400 |
|
|
[考研] 320分,材料與化工專業(yè),求調(diào)劑 +10 | 一定上岸aaa 2026-03-27 | 14/700 |
|
|
[考研] 269求調(diào)劑 +4 | 我想讀研11 2026-03-31 | 4/200 |
|
|
[考研] 293分求調(diào)劑,外語為俄語 +5 | 加一一九 2026-03-31 | 5/250 |
|
|
[考研] 085602化工求調(diào)劑(331分) +8 | 111@127 2026-03-30 | 8/400 |
|
|
[考研] 化學(xué)0703 調(diào)劑 306分 一志愿211 +8 | 26要上岸 2026-03-28 | 8/400 |
|
|
[考研] 求化學(xué)調(diào)劑 +11 | wulanna 2026-03-28 | 11/550 |
|
|
[考研] 材料與化工304求B區(qū)調(diào)劑 +4 | 邱gl 2026-03-26 | 7/350 |
|
|
[考研] 311求調(diào)劑 +6 | 冬十三 2026-03-24 | 6/300 |
|
|
[考研] 總分293求調(diào)劑 +8 | 加一一九 2026-03-25 | 11/550 |
|
|
[考研] 330分求調(diào)劑 +5 | qzenlc 2026-03-29 | 5/250 |
|
|
[碩博家園] 招收生物學(xué)/細(xì)胞生物學(xué)調(diào)劑 +4 | IceGuo 2026-03-26 | 5/250 |
|
|
[考研] 081200-11408-276學(xué)碩求調(diào)劑 +6 | 崔wj 2026-03-26 | 6/300 |
|
|
[考研] 283求調(diào)劑 +3 | A child 2026-03-28 | 3/150 |
|
|
[考研] 308求調(diào)劑 +7 | 墨墨漠 2026-03-27 | 7/350 |
|
|
[考研] 315分求調(diào)劑 +7 | 26考研上岸版26 2026-03-26 | 7/350 |
|
|
[考研] 265求調(diào)劑 +8 | 小木蟲085600 2026-03-27 | 8/400 |
|