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Целое число, задающее количество равномерно распределенных точек (разделительные стоимости единицы) в передискретизируемом сигнале.
RangeValueВектор 1 на 2, задающий минимальные и максимальные разделительные стоимости единицы для выходной матрицы Intensities. RangeValue должен быть в [min (inputSU) макс. (inputSU)], где inputSU конкатенированные разделительные стоимости единицы от входа Peaklist. Значением по умолчанию является полный спектр [min (inputSU) макс. (inputSU)].
FWHHValueЗначение, которое задает полную ширину на половине высоты (FWHH) в разделительных модулях. FWHH используется, чтобы преобразовать каждый пик в кривую гауссовой формы. Значение по умолчанию является средним (diff (inputSU))/2, где inputSU конкатенированные разделительные стоимости единицы от входа Peaklist. Значением по умолчанию является грубое приближение разрешения, наблюдаемого во входных данных, Peaklist.

Совет

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

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

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

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

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

Описание

Совет

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

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

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

Совет

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

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

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

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

Совет

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

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

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

  • 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

Для просмотра документации необходимо авторизоваться на сайте