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年末，波动最大），华法林通过紧急注射依诺肝素来补充，其效果甚至更大。