全部版块 我的主页
论坛 计量经济学与统计论坛 五区 计量经济学与统计软件
1043 1
2025-10-13

本文将讲解如何对曲线进行平滑处理、将曲线转换为其导数或积分,以及如何计算曲线下面积。

曲线的平滑、求导与积分

一次Prism分析可对曲线进行平滑处理,和/或将曲线转换为其导数或积分。


1.gif


求曲线的导数或积分

一阶导数是曲线在每个X值处的陡峭程度。当曲线向上倾斜时,导数为正;向下倾斜时,导数为负。曲线的峰和谷处导数为零。计算出数值导数后,若您选择,Prism可对结果进行平滑处理。

二阶导数是导数曲线的导数。曲线的拐点处二阶导数为零。

积分是曲线与Y=0线(或您输入的其他值)之间的累积面积。


注意事项:


  • Prism不支持符号代数或符号微积分运算。若您向Prism提供一系列定义曲线的XY点,它可以通过这些点计算数值导数(或积分)。但如果您向Prism提供一个方程,它无法计算出定义该导数或积分的新方程。
  • 该分析对曲线进行积分,会得到另一条显示累积面积的曲线。请不要将其与Prism中另一个单独计算曲线下面积单一数值的分析混淆。



曲线平滑

若您从仪器导入一条曲线,您可能希望对数据进行平滑处理以改善图表的外观。由于平滑曲线时会丢失数据,因此在进行非线性回归或其他分析之前,不应对曲线进行平滑处理。平滑处理并非一种数据分析方法,而纯粹是一种创建更美观图表的方式。Prism提供两种调整曲线平滑度的方式:您可以选择用于平均的相邻点数量,以及平滑多项式的“阶数”。由于平滑的唯一目的是让曲线看起来更美观,您可以简单尝试几种设置,直到满意结果的外观为止。如果设置过高,一些峰将会被平滑掉;如果设置过低,曲线则不够平滑。这种平衡是主观的,可以通过反复尝试来确定。


结果表的行数比原始数据少。


勿分析已平滑的数据


对曲线进行平滑处理可能具有误导性。其本意是减少“杂点”,以便您能看清实际趋势。但问题在于,您可能会看到并不真实存在的“趋势”。下方上排的三个图表是模拟数据,每个数值均来自均值为50、标准差为10的高斯分布,且每个数值都是独立抽取的,与之前的数值无关。当您观察这三个图表时,会看到围绕一条水平线的随机散点,这与数据的生成方式完全一致。


2.png


上方下排的三个图表展示了相同数据经平滑处理(每侧平均10个数值,并使用二阶平滑多项式)后的结果。当您观察这些图表时,会看到趋势;第一个图表呈现下降趋势,第二个图表似乎又规律地波动,第三个图表呈现上升趋势。这些趋势都是平滑处理的产物,每个趋势都与上方未处理的图表显示的内容一致。


对数据进行平滑处理会通过放大随机出现的大幅高低值波动、同时削弱点与点之间的变异性,从而营造出趋势的假象。相关分析、线性回归喝非线性回归的一个关键假设是数据点相互独立。对于已平滑的数据,这一假设不成立。如果一个数值出现极高或极低的情况,平滑后相邻的点也会如此。由于随机趋势被放大,随机散点被削弱,对已平滑数据的任何分析(若未考虑平滑处理)都将无效。


数学细节


  • 一阶导数的计算方式如下(x和Y是数据数组,x'和y'是包含结果的数组):

x' = (x[i+1] + x) / 2


y' at x' = (y[i+1] - y) / (x[i+1] - x)


  • 二阶导数的计算是通过两次运行该算法,本质上是计算一阶导数的一阶导数。
  • Prism使用梯形法则对曲线进行积分。结果的X值与您分析的数据的X值相同。结果的第一个Y值等于您指定的数值(通常为0.0)。对于其他行,结果的Y值等于前一个结果加上通过该点新增的面积。该面积等于X值的差值乘以之前和当前Y值的平均值
  • 平滑处理采用Savitsky和Golay的方法
  • 若您要求Prism同时进行平滑处理并转换为导数(一阶或二阶)或积分,Prism会按顺序执行这些步骤:首先生成导数或积分,然后进行平滑处理



曲线下面积

操作方法:曲线下面积

曲线下面积是对可测量效应或现象的综合测量。在药代动力学种,它被用作药物效应的累积测量;在色谱分析中,它被用作比较峰的手段。

请注意,Prism还会再单独的受试者工作特征(ROC)分析中计算ROC曲线下的面积。


从代表曲线的数据表或结果表开始。点击“分析”,并从XY分析列表中选择“曲线下面积”。


3.gif


曲线下面积结果解读

若您的数据来自色谱或光谱分析,Prism可将数据划分为不同区域,并确定每个区域的最高点(峰)。不过,Prism仅在区域定义清晰时才能执行此操作:信号(或效应/现象的图形表示)在区域之间必须低于基线,且峰之间不能重叠。

对于每个区域,Prism会以X轴单位乘以Y轴单位的形式显示面积。同时,Prism还会将每个区域显示为所有区域总面积的占比。面积的计算采用梯形法则,即通过在定义曲线的每一组相邻点之间连接直线,然后对这些直线下方的面积求和。

接下来,Prism会识别每个区域的峰,报告该区域最高点的X和Y坐标,以及代表区域起始和结束的两个X坐标。


若Prism识别的区域数量超出您的需求,可返回“参数”对话框,增大区域最小宽度和/或峰最小高度的数值。


该分析的局限性


请注意以下局限性:


  • 基线必须是水平的
  • 无平滑或曲线拟合操作
  • Prism无法分离重叠的峰。除非信号在两个峰之间完全降至基线,否则程序无法区分相邻的两个峰。同样,Prism也无法识别一个峰在另一个峰肩部的情况
  • 若信号在基线以上开始(或结束),第一个(或最后一个)峰将不完整,Prism会报告其“看到”的尾部下方的面积
  • 若您的第一个X值大于0,Prism不会外推回X=0
  • Prism不会外推超出数据集中的最大X值,因此也不会将曲线外推至基线以下
  • Prism不再要求X值等距。在对梯形面积求和时,即使部分梯形更“宽”,也不会有影响



Prism如何计算曲线下面积


Prism采用梯形法则计算曲线下面积,如下图所示。


4.gif


在Prism中,由非线性回归创建的曲线本质上是一系列相连的XY点,且X值等距分布。Prism也可对您输入的XY表格计算曲线下面积,不要求X值必须等距。上图左侧展示了其中两个点以及作为虚线的基线。该部分曲线下的梯形区域被阴影标注。图的中间部分展示了Prism计算面积的方式:中间面板中的两个三角形面积相同,因此左侧梯形的面积与右侧矩形的面积(更易计算)相等。因此,面积为ΔX*([(Y1+Y2)/2]-基线]。Prism会对定义曲线的每对相邻点重复使用该公式。


面积计算基于您指定的基线以及两个X值之间的曲线。具体取哪些X值呢?


  • 若所有数据点均大于基线,AUC(曲线下面积)计算从数据集中的最小X值开始,到最大X值结束。注意,Prism不会将曲线延伸至数据的X范围之外
  • 若最小X值对应的Y值低于您的基线:Prism会找到数据集中与大于基线的Y值对应的最小X值。它会在该点与数据集中下一个最小X值的点之间画一条直线,然后使用线性内插法找到该直线与基线的交点,并将该内插值作为计算AUC的第一个X值
  • 若最大X值对应的Y值低于您的基线:Prism会找到数据集中与大于基线的Y值对应的最大X值。它会在该点与数据集中下一个最大X值的点之间画一条直线,然后使用线性内插法找到该直线与基线的交点,并将该内插值作为计算AUC的最后一个X值



AUC的标准误与置信区间

若您输入带有重复Y值的数据,或以均值、标准差或标准误的形式输入数据,Prism会采用Gagnon所述的方法报告AUC的标准误和置信区间。若您在子列中输入重复的Y值,Prism会将其视为同一实验中的重复测量。若每个子列实际是不同的重复实验,Prism不会按子列分别计算AUC再取平均值。95%置信区间等于AUC加减1.96倍的标准误。它使用z分布(因此始终为1.96)而非t分布(其数值取决于样本量),因为参考文献1-3采用了这种方式。当定义曲线的点超过几十个时,t方法和z方法的结果几乎无差异。


Prism不会通过比较峰来提供差异的置信区间或对应的P值。但您可以通过以下步骤让Prism实现这一功能:


  • 创建一个新的分组表,格式设置为输入均值、SEM和样本量(n)。您只需在该表的第一行输入数值
  • 将AUC值作为均数输入
  • 将AUC值的标准误作为“SEM”输入
  • 为每个组定义自由度(df):该组的数据点数减去浓度数
  • 对于n,输入比自由度大1的数值。当Prism进行t检验时,会从输入的n中减1以获得正确的自由度
  • 点击“分析”,若要比较两个AUC则选择t检验,若要比较三个或更多则选择单因素方差分析



何为“峰”?

默认情况下,Prism仅将基线上方的点视为峰的一部分,因此仅报告高于基线的峰。您可以选择将低于基线的峰也纳入考虑。

默认情况下,Prism会忽略高度小于Y值最小值到最大值距离的10%的峰,但您可以在“曲线下面积参数”对话框中修改这一定义。您也可以设置让Prism忽略非常窄的峰。

总峰面积vs.总面积vs.净面积


Prism以两种或三种方式报告峰下面积:


  • 总面积。它是阳性峰、阴性峰、峰高不足以计数的峰以及峰宽过窄不足以计数的峰的面积总和。在分析对话框中,唯一会影响总面积定义的选择是基线的定义
  • 总峰面积。您要求Prism纳入考虑的峰的面积总和。该数值会受到分析对话框中若干选择的影响:基线的定义、关于纳入或忽略阴性峰的选择,以及关于峰小到不足以计算的定义
  • 净面积。只有当您要求Prism将基线以下的峰定义为峰时,才会看到该数值。它是通过将基线上方峰的面积减去基线下方峰的面积计算得出的差值

附件列表
3.gif

原图尺寸 11.78 KB

3.gif

二维码

扫码加我 拉你入群

请注明:姓名-公司-职位

以便审核进群资格,未注明则拒绝

全部回复
2025-10-14 20:38:27
二维码

扫码加我 拉你入群

请注明:姓名-公司-职位

以便审核进群资格,未注明则拒绝

相关推荐
栏目导航
热门文章
推荐文章

分享

扫码加好友,拉您进群