上週二之後覺得很沮喪
但沮喪完還是得爬起來繼續加油
先講跑步
第一次堪堪跑1.5k就覺得要暈死在路邊，頭痛得不行
之後固定每天 5-6 pm 都去跑
昨天已經可以跑到 3k, 而且前2k 都可以在 7分速以內!!! H
(06:33, 06:45 [HR 182!? 其實當下沒啥感覺] ; 3rd km 就只能 07:13)
超感動要哭了，在台灣時怎麼跑都跑不進七分鐘阿奇怪 QQ
果真是心裡不開心就去訓練自己的身體吧
關於我的 AMK model
Dr. Natalia 給了我一篇 Monolix workflow 的 paper
我跟著做，學了 syncomore 覺得真是超好用
拿來做 model comparison 輕鬆很多，眼睛不會再抽筋
只差 prediction vs. observation plot 還是得慢慢看
(果然得用雙螢幕阿)
可是... 不管我怎麼改變項，預測不準的地方就是不準
永遠都有幾個點 underprediction 落在離線好遙遠的地方，不動如山
回頭去看data, 應該就是抽血時間和給藥時間有記錄錯誤
試問給了 08:46 IV 750 mg, 08:53 concentration= 3 mcg/mL 是怎樣 !?!?!?
(我都要哭了 .....)
而且當 covariate 加越多，Vd beta 的 RSE 就越大
但 Vd beta 本身縮得很小，所以即使 SE 有縮小，RSE 還是可以膨脹到接近100%
(如果把 UF volume 當成 regressor, RSE 變成 800%.......)
為此我還寫了email問 Dr. Natalia
因為照著數據畫常態分布圖，曲線並沒有變胖啊 QQ
如果照常態分佈圖，綠線應該比橘線或紅線好吧
但 Natalia 說，還是要以 RSE 為判斷依據
因為那才能代表該數據預測的「變異程度」
我後來好像想通了，就像 10+/- 1 (RSE 10%) 和 1 +/- 0.5 (RSE 50%)，終究是 10 +/- 的預測比較準，因為「變異程度比較小」
所以我用 Priscila's model 怎樣都會遇到 Vd beta 擴張的問題，暫時還想不到好解法 orz
另外是我想用 AMIDIAL-ICU published data 做 model fitting
問了 Priscila, 她是說不需要寫信去問作者
用他們 paper 上的數據應該就可以寫出模型
但是我不會阿.............. QQ
看了半天 Monolix ODE, 勉強拼湊了一點東西出來但完全不知道對不對?
應該說我不太懂我自己在寫甚麼
有種依樣畫葫蘆的感覺
DESCRIPTION:
The administration is via a bolus.
The PK model has a central compartment (volume V1), a peripheral compartment (volume V2, intercompartmental clearance Q), and a linear elimination (clearance Cl).
[LONGITUDINAL]
input = {Cl, Cld, V1, Q, V2}
Equation:
t_0=0
Ac_0=0
Ap_0=0
ddt_Ac= Ac-Q/V1*Ac-Cl/V1*Ac-Cld/V1*Ac+Q/V2*Ap
ddt_Ap=Q/V1*Ac-Q/V2*Ap
PK:
; Parameter transformations
V = V1
k12 = Q/V1
k21 = Q/V2
; PK model definition
Cc = Ac/V1
OUTPUT:
output = {Cc}
(Syntax 檢查說沒錯，但... Cld 沒有定義啊? 好吧我是真的不懂 orz)
心一橫，還是用了 ufl.edu 的信箱 email 給 AMIDIAL-ICU 的 corresponding author
反正大不了不回而已 QvQ
下週 Natalia 不在，我也要去 Jacksonville 了
再下一週就是 last week...... 怎麼辦我的 AMK model QQ
而且最後還要搞清楚建好 model 之後怎麼預測new subject的濃度
(應該是用 simulx?)
(然後跟我槓上的 BPS exam, 12月中要再去試筆，已瘋 XD)
ps. 心的某一個角落還是好想小孩喔 QQ...... 好想回台灣.......
想好好的抱一下孩子們阿阿阿
小孩10/20 管樂市賽加油!!
沒有留言:
張貼留言