主要内容

msppresample

与同时保留峰峰值重采样信号

句法

[X强度] = msppresample(PeakList.N
msppresample(PeakList.N, ...'范围',RangeValue.,......)
msppresample(PeakList.N,... 'FWHH',fwhhvalue.,......)
msppresample(PeakList.N,... 'ShowPlot',ShowPlotValue,......)

输入参数

PeakList.

以下任一操作:

  • 双列矩阵,第一列包含分离单元值,第二列包含强度值。分离单元可以根据产生信号数据的仪器量化波长,频率,距离,时间或m / z。

  • 单元阵列峰列表,其中每个元素是分离单元的值和强度值,并且每个元素对应的两列的矩阵的信号或保留时间。

提示

你可以使用mzxml2peaks.功能或mspeaks.创建的功能PeakList.矩阵或单元阵列。

N 整数,指定在二次采样信号相等间隔分(分离单元值)的数目。
RangeValue. 1-by-2向量,指定输出矩阵的最小和最大分离单元值强度RangeValue.必须在[分钟(inputSU) 最大限度(inputSU, 在哪里inputSU是来自输入的连接分离单元值PeakList.。默认是全系列[分钟(inputSU) 最大限度(inputSU
fwhhvalue. 指定分离单元中半高(FWHH)的全宽的值。FWHH用于将每个峰值转换为高斯形曲线。默认为中位数(DIFF(inputSU))/ 2, 在哪里inputSU是来自输入的连接分离单元值PeakList.。默认值是输入数据中观察到的分辨率的粗略近似,PeakList.

提示

为了确保峰的分辨率被保留,集fwhhvalue.感兴趣的两个峰是最接近彼此之间距离的一半。

ShowPlotValue 控制的原始和重新采样的信号的曲线图的显示。选择是真的错误的, 或者一世的整数,指定的信号的索引中的强度。如果你设置为真的,第一个信号强度被绘制了。默认为:
  • 错误的- 当指定的返回值。

  • 真的- 未指定返回值时。

输出参数

X 矢量与峰值一组信号的同等间隔,常见的分离单元值。矢量中的元素数量等于N,在矩阵中的行数或强度
强度 用于共享相同分离单元范围的一组峰的重建强度值的矩阵。每行对应于分离单元值,并且每列对应于具有峰值或保留时间的一组信号。行数等于N,或矢量中的元素数量X

描述

提示

从任何分离技术,其产生的信号数据,如光谱学,NMR,电泳,层析,或质谱法使用带有数据如下语法。

[X强度] = msppresample(PeakList.N重建PeakList.,一个峰值列表,通过将蜈蚣峰转换为半连续的原始信号,以保留峰值信息。重采样信号具有N同样间隔的点。输出X是的向量N元素指定具有峰值的一组信号集的同等间隔的常见分离单元值。输出强度是一组共享相同分离单元范围的一组峰的重建强度值的矩阵。每行对应于分离单元值,并且每列对应于具有峰值或保留时间的一组信号。行数等于N

msppresample使用高斯内核重建信号。任何给定分离单元值的强度取自任何贡献(重叠)峰的最大强度。

提示

msppresample可用于制备一组信号用于成像功能,如msheatmap和预处理功能,如msbackadjmsnorm.

msppresample(PeakList.N,......propertyname.',适当的价值,......)呼叫msppresample部分可选的属性是使用属性名称/属性值对。您可以按任意顺序指定一个或多个属性。每个propertyname.必须用单引号括起来,而是不论不敏感。这些属性名称/属性值对如下:

msppresample(PeakList.N, ...'范围',RangeValue.,......)指定输出矩阵的分离单元范围强度使用1-by-2向量中指定的最小和最大分离值RangeValue.RangeValue.必须在[分钟(inputSU) 最大限度(inputSU, 在哪里inputSU是来自输入的连接分离单元值PeakList.。默认是全系列[分钟(inputSU) 最大限度(inputSU

msppresample(PeakList.N,... 'FWHH',fwhhvalue.,......)在分离单元中设置半高(FWHH)的全宽度。FWHH用于将每个峰值转换为高斯形曲线。默认为中位数(DIFF(inputSU))/ 2, 在哪里inputSU是来自输入的连接分离单元值PeakList.。默认值是输入数据中观察到的分辨率的粗略近似,PeakList.

提示

为了确保峰的分辨率被保留,集fwhhvalue.感兴趣的两个峰是最接近彼此之间距离的一半。

msppresample(PeakList.N,... 'ShowPlot',ShowPlotValue,......)控制原始和重采样信号的曲线的显示。选择是真的错误的, 或者一世的整数,指定的信号的索引中的强度。如果你设置为真的,第一个信号强度被绘制了。默认为:

  • 错误的- 当指定的返回值。

  • 真的- 未指定返回值时。

例子

  1. 加载MAT文件,包含在生物信息工具箱™软件中,包含液相色谱/质谱(LC / MS)数据变量。这包括,峰列表,其中每个元素是m / z值和离子强度值,并且每个元素对应的两列的矩阵的光谱或保留时间的单元阵列。

    加载lcmsdata.
  2. 重新取样的数据,在重新采样的信号指定5000米/ z值。然后创建LC / MS数据的热图。

    [MZ,Y] = msppresample(ms_peaks,5000);msheatmap(MZ,ret_time,日志(Y))

  3. 画出两个保持时间之间重建的轮廓光谱。

    图T1 = 3370;T2 = 3390;H =查找(RET_TIME> T1&RET_TIME '质量/充电(m / z)')ylabel('保留时间')Zlabel(“相对强度”

  4. 重新确定数据以绘制总离子色谱图(TIC)。

    图[MZ,Y] = MSPPRESAMPLE(MS_PEAKS,5000);plot(ret_time,sum(y))标题('总离子色谱图(TIC)')xlabel('保留时间')ylabel(“相对强度”

  5. 重新确定数据以在450至500 m / z范围内绘制提取的离子色谱图(XIC)。

    图[MZ,Y] = MSPPRESAMPLE(MS_PEAKS,5000,'范围',[450 500]);plot(ret_time,sum(y))标题('提取离子色谱(XIC)从450到500 M / Z')xlabel('保留时间')ylabel(“相对强度”

在R2007A介绍