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使用贝叶斯优化的深度学习

7.matlab贝叶斯隐马尔可夫hmm模型

8.matlab中的隐马尔可夫模型(HMM)实现

9.matlab实现MCMC的马尔可夫切换ARMA – GARCH模型

标签：读数,xlim,样条,华法林,INR,matlab,Date
来源： https://blog.51cto.com/u_15198753/2770834