Exponenta Banner

msppresample

Передискретизируйте сигнал с peaks при сохранении peaks

Синтаксис

[X, Intensities] = msppresample(Peaklist, N)
msppresample(Peaklist, N, ...'Range', RangeValue, ...)
msppresample(Peaklist, N, ...'FWHH', FWHHValue, ...)
msppresample(Peaklist, N, ...'ShowPlot', ShowPlotValue, ...)

Входные параметры

Peaklist

Любое из следующего:

  • Матрица 2D столбца, где первый столбец содержит разделительные стоимости единицы и второй столбец, содержит значения интенсивности. Разделительный модуль может определить количество длины волны, частоты, расстояния, время или m/z в зависимости от инструмента, который генерирует данные сигнала.

  • Массив ячеек пиковых списков, где каждый элемент является матрицей 2D столбца разделительных стоимостей единицы и значений интенсивности и каждого элемента, соответствует время задержания или сигналу.

Совет

Можно использовать функцию mzxml2peaks или функцию mspeaks, чтобы создать матрицу Peaklist или массив ячеек.

NЦелое число, задающее количество равномерно распределенных точек (разделительные стоимости единицы) в передискретизируемом сигнале.
RangeValue1 2 вектор, задающий минимальные и максимальные разделительные стоимости единицы для выходной матрицы Intensities. RangeValue должен быть в [min(inputSU) max(inputSU)], где inputSU является конкатенированными разделительными стоимостями единицы от входа Peaklist. Значением по умолчанию является полный спектр [min(inputSU) max(inputSU)].
FWHHValueЗначение, которое задает полную ширину на половине высоты (FWHH) в разделительных модулях. FWHH используется, чтобы преобразовать каждый пик в кривую гауссовой формы. Значением по умолчанию является median(diff(inputSU))/2, где inputSU является конкатенированными разделительными стоимостями единицы от входа Peaklist. Значением по умолчанию является грубое приближение разрешения, наблюдаемого во входных данных, Peaklist.

Совет

Чтобы гарантировать, что разрешение peaks сохраняется, устанавливает FWHHValue на половину расстояния между двумя peaks интереса, которые являются самыми близкими друг к другу.

ShowPlotValueУправляет отображением графика исходного и передискретизируемого сигнала. Выбором является true, false или I, целое число, задающее индекс сигнала в Intensities. Если вы устанавливаете на true, первый сигнал в Intensities построен. Значение по умолчанию:

  • ложь Когда возвращаемые значения заданы.

  • tRUE Когда возвращаемые значения не заданы.

Выходные аргументы

X Вектор равномерно распределенных, общих разделительных стоимостей единицы для набора сигналов с peaks. Число элементов в векторе равняется N или количеству строк в матричном Intensities.
Intensities Матрица восстановленных значений интенсивности для набора peaks, который совместно использует ту же область значений разделительного модуля. Каждая строка соответствует разделительной стоимости единицы, и каждый столбец соответствует или набору сигналов с peaks или время задержания. Количество строк равняется N или числу элементов в векторном X.

Описание

Совет

Используйте следующие синтаксисы с данными из любого разделительного метода, который производит данные сигнала, такие как спектроскопия, NMR, электрофорез, хроматография или масс-спектрометрия.

[X, Intensities] = msppresample(Peaklist, N) передискретизирует Peaklist, пиковый список, путем преобразования центроидного peaks в полунепрерывный, необработанный сигнал, что консервы достигают максимума информация. Передискретизируемый сигнал имеет N равномерно распределенные точки. Вывод X является вектором элементов N, задающих равномерно распределенные, общие разделительные стоимости единицы для набора сигналов с peaks. Вывод Intensities является матрицей восстановленных значений интенсивности для набора peaks, который совместно использует ту же область значений разделительного модуля. Каждая строка соответствует разделительной стоимости единицы, и каждый столбец соответствует или набору сигналов с peaks или время задержания. Количество строк равняется N.

msppresample использует Гауссово ядро, чтобы восстановить сигнал. Интенсивность по любой данной разделительной стоимости единицы взята из максимальной интенсивности любого содействия (наложение) peaks.

Совет

msppresample полезен, чтобы подготовить набор сигналов для обработки изображений функций, таких как msheatmap и предварительная обработка функций, таких как msbackadj и msnorm.

msppresample(Peaklist, N, ... 'PropertyName', PropertyValue, ...)   вызывает msppresample с дополнительными свойствами, которые используют имя свойства / пары значения свойства. Можно задать одно или несколько свойств в любом порядке. Каждый PropertyName должен быть заключен в одинарные кавычки и нечувствительный к регистру. Это имя свойства / пары значения свойства следующие:

msppresample(Peaklist, N, ...'Range', RangeValue, ...) задает область значений разделительного модуля для выходной матрицы Intensities с помощью минимальных и максимальных разделительных значений, заданных в 1 2 векторном RangeValue. RangeValue должен быть в [min(inputSU) max(inputSU)], где inputSU является конкатенированными разделительными стоимостями единицы от входа Peaklist. Значением по умолчанию является полный спектр [min(inputSU) max(inputSU)]

msppresample(Peaklist, N, ...'FWHH', FWHHValue, ...) устанавливает полную ширину на половине высоты (FWHH) в разделительных модулях. FWHH используется, чтобы преобразовать каждый пик в кривую гауссовой формы. Значением по умолчанию является median(diff(inputSU))/2, где inputSU является конкатенированными разделительными стоимостями единицы от входа Peaklist. Значением по умолчанию является грубое приближение разрешения, наблюдаемого во входных данных, Peaklist.

Совет

Чтобы гарантировать, что разрешение peaks сохраняется, устанавливает FWHHValue на половину расстояния между двумя peaks интереса, которые являются самыми близкими друг к другу.

msppresample(Peaklist, N, ...'ShowPlot', ShowPlotValue, ...) управляет отображением графика исходного и передискретизируемого сигнала. Выбором является true, false или I, целое число, задающее индекс сигнала в Intensities. Если вы устанавливаете на true, первый сигнал в Intensities построен. Значение по умолчанию:

  • ложь Когда возвращаемые значения заданы.

  • tRUE Когда возвращаемые значения не заданы.

Примеры

  1. Загрузите MAT-файл, включенный с программным обеспечением Bioinformatics Toolbox™, которое содержит жидкостную хроматографию / переменные данных (LC/MS) масс-спектрометрии. Это включает peaks, массив ячеек пиковых списков, где каждый элемент является матрицей 2D столбца m/z значений и ионных значений интенсивности, и каждый элемент соответствует время задержания или спектру. 

    load lcmsdata

  2. Передискретизируйте данные, задавая 5 000 m/z значений в передискретизируемом сигнале. Затем создайте карту тепла данных LC/MS.

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

  3. Постройте восстановленные спектры профиля между двумя разами задержания.

    figure
t1 = 3370;
t2 = 3390;
h = find(ret_time>t1 & ret_time<t2);
[MZ,Y] = msppresample(ms_peaks(h),10000);
plot3(repmat(MZ,1,numel(h)),repmat(ret_time(h)',10000,1),Y)
xlabel('Mass/Charge (M/Z)')
ylabel('Retention Time')
zlabel('Relative Intensity')

  4. Передискретизируйте данные, чтобы построить Общую ионную хроматограмму (TIC).

    figure
[MZ,Y] = msppresample(ms_peaks,5000);
plot(ret_time,sum(Y))
title('Total Ion Chromatogram (TIC)')
xlabel('Retention Time')
ylabel('Relative Intensity')

  5. Передискретизируйте данные, чтобы построить Извлеченную Ионную Хроматограмму (XIC) в 450 - 500 областях значений m/z.

    figure
[MZ,Y] = msppresample(ms_peaks,5000,'Range',[450 500]);
plot(ret_time,sum(Y))
title('Extracted Ion Chromatogram (XIC) from 450 to 500 M/Z')
xlabel('Retention Time')
ylabel('Relative Intensity')

Смотрите также

| | | | | | |

Темы

Представленный в R2007a

Документация Bioinformatics Toolbox
Поддержка