matlab使用樣條插值重采樣估計INR資料研究

原文連結：http://tecdat.cn/?p=7019

易患血液凝固的人用華法林治療，血液稀釋劑。國際标準化比率（INR）衡量藥物的效果。較大劑量會增加INR，較小劑量會降低INR。患者由護士定期監測，當他們的INR超出目标範圍時，他們的劑量和測試頻率會發生變化。

該檔案

INR.mat

包含在五年内對患者進行的INR測量。該檔案包括一個

datetime

數組，其中包含每次測量的日期和時間，以及一個帶有相應INR讀數的矢量。加載資料。 

plot(Date,INR,'o','DatetimeTickFormat','MM/dd/yy')
 
plot([xlim;xlim]',[2 3;2 3],'k:')           

重新采樣資料以使INR讀數均勻分布。 使用

resample

當時在以後每星期五估計病人的INR。指定每周一次讀數的采樣率，或等效地，每秒讀數1 / （7 × 8 6 4 0 0 ）。使用樣條插值進行重采樣。

plot(tee,rum,'.-','DatetimeTickFormat','MM/dd/yy')

title('INR')
xlim([Date(1) Date(end)])
hold off           

每次INR讀數确定何時必須對患者進行測試。使用

diff

建構測量之間的時間間隔的向量。以周為機關表示間隔，并使用與以前相同的x軸繪制它們。 

plot(Date,diff(datenum([Date;nxt]))/7,'o-', ...
    'DatetimeTickFormat','MM/dd/yy')

title('Time Until Next Reading')
xlim([Date(1) Date(end)])
ylabel('Weeks')           

當INR超出範圍時，INR讀數之間的時間仍然很短。當INR過低時，患者會更頻繁地獲得讀數，因為血栓形成的風險會升高。當患者的INR在範圍内時，讀數之間的時間穩定增加，直到比率變得太小或太大。

重采樣的大幅波動可能是過沖的迹象。然而，華法林對身體有很大的影響。華法林劑量的微小變化可以大大改變INR，飲食，飛機上花費的時間或其他因素也會發生變化。此外，當比例非常低時（如2010年末，波動最大），華法林通過緊急注射依諾肝素來補充，其效果甚至更大。

參考文獻

1.matlab使用經驗模式分解emd 對信号進行去噪

2.Matlab使用Hampel濾波去除異常值

3.matlab偏最小二乘回歸(PLSR)和主成分回歸(PCR)

4.matlab預測ARMA-GARCH 條件均值和方差模型

5.matlab中使用VMD(變分模态分解) 

6.matlab使用貝葉斯優化的深度學習