mslowess
用非参数法平滑峰值信号
语法
你
= mslowess (X
,强度
)
mslowess(…“秩序”,OrderValue
,……)
mslowess(…“跨越”,SpanValue
,……)
mslowess(…“内核”,KernelValue
,……)
mslowess(…,'RobustIterations',RobustIterationsValue
,……)
mslowess(…,'ShowPlot',ShowPlotValue
,……)
参数
X |
一组具有峰值的信号的分离单位值的向量。向量中的元素数等于矩阵中的行数强度 .分离单元可以量化波长、频率、距离、时间或m/z,这取决于产生信号数据的仪器。 |
强度 |
具有相同分离单位范围的一组峰值的强度值矩阵。每一行对应一个分离单元值,每一列对应一组有峰值的信号或保留时间。行数等于向量中的元素数X . |
描述
提示
使用以下语法处理来自产生信号数据的任何分离技术的数据,如光谱学、核磁共振、电泳、色谱或质谱。
平滑原始噪声信号数据,你
= mslowess (X
,强度
)强度
,使用默认跨度为的局部加权线性回归(Lowess)方法10
样本。
请注意
mslowess
假设输入向量,X
,可能没有均匀间距的分离单元。因此,用于平滑的滑动窗口的居中使用最接近的样本X
价值而非以之而论X
索引。
请注意
当输入向量,X
,没有重复值或NaN值,算法的速度大约是原来的两倍。
mslowess (
调用X
,强度
,……”PropertyName
',PropertyValue
,……)mslowess
使用属性名/属性值对的可选属性。可以以任意顺序指定一个或多个属性。每一个PropertyName
必须用单引号括起来,不区分大小写。这些属性名/属性值对如下:
mslowess(…“秩序”,
指定顺序(OrderValue
,……)
)。输入OrderValue
1
(线性多项式拟合或Lowess),2
(二次多项式拟合或黄土),或0
(相当于加权局部均值估计器,而且可能更快,因为只执行平均值计算,而不是最小二乘回归)。默认值为1
.
请注意
曲线拟合工具箱™软件也指低阶平滑2
黄土平滑。
mslowess(…“跨越”,
指定平滑内核的窗口大小。如果SpanValue
,……)
大于SpanValue
1
,窗口等于SpanValue
独立于分离单位向量的样本数,X
.默认值为10
样本。较高的值将使信号更加平滑,但要牺牲计算时间。如果SpanValue
小于1
,窗口大小被认为是数据中点数的一个分数。例如,当SpanValue
是0.005
,窗口大小等于0.50
中点数的%X
.
mslowess(…“内核”,
所指定的函数KernelValue
,……)KernelValue
用于对观测到的强度进行加权。靠近被平滑的分离单元位置的样本在确定估计值时具有最大的权重。KernelValue
可以是以下字符向量(或字符串)中的任何一个:
“tricubic”
(默认)(1 - (dist /距离)。^ 3)。^ 3
“高斯”
- - - - - -exp (- (2 * dist /距离)^ 2)。
“线性”
- - - - - -1-dist /距离
mslowess(…,'RobustIterations',
指定迭代次数(RobustIterationsValue
,……)
)适合强健的身体。如果RobustValue
RobustIterationsValue
是0
(默认),不执行健壮拟合。对于鲁棒平滑,在每个跨度上的小残差值被忽略,以改进新的估计。1
或2
健壮的迭代通常就足够了,而较大的值在计算上可能会比较昂贵。
请注意
对于一个X
具有均匀间距分离单元的向量,具有非鲁棒平滑OrderValue
等于0
相当于用核向量对信号进行滤波。
mslowess(…,'ShowPlot',
在原始信号上绘制平滑信号。当你打电话ShowPlotValue
,……)mslowess
没有输出参数,信号被绘制,除非ShowPlotValue
是假
.当ShowPlotValue
是真正的
,只有第一个信号强度
是绘制。ShowPlotValue
也可以包含到?中的一个信号的索引强度
.
例子
加载一个mat文件,包括生物信息学工具箱™软件,其中包含一些示例数据。
负载sample_lo_res
平滑频谱,画出原始信号和平滑信号的第一个频谱图。
y = mslowess (MZ_lo_res Y_lo_res,“Showplot”,真正的);
放大图中的一个区域,看看原始信号和平滑信号的区别。
轴([7350 7550 0.1 1.0])
另请参阅
mspalign
|msbackadj
|msdotplot
|msalign
|msheatmap
|mslowess
|msnorm
|mspeaks
|msresample
|msppresample
|mssgolay
|msviewer