Преобразуйте необработанные пиковые данные в пиковый список (централизованные данные)
Peaklist = mspeaks(X, Intensities)
[Peaklist, PFWHH]
= mspeaks(X, Intensities)
[Peaklist, PFWHH, PExt]
= mspeaks(X, Intensities)
mspeaks(X, Intensities,
...'Base', BaseValue, ...)
mspeaks(X, Intensities,
...'Levels', LevelsValue, ...)
mspeaks(X, Intensities,
...'NoiseEstimator', NoiseEstimatorValue,
...)
mspeaks(X, Intensities,
...'Multiplier', MultiplierValue, ...)
mspeaks(X, Intensities,
...'Denoising', DenoisingValue, ...)
mspeaks(X, Intensities,
...'PeakLocation', PeakLocationValue, ...)
mspeaks(X, Intensities,
...'FWHHFilter', FWHHFilterValue, ...)
mspeaks(X, Intensities,
...'OverSegmentationFilter', OverSegmentationFilterValue,
...)
mspeaks(X, Intensities,
...'HeightFilter', HeightFilterValue, ...)
mspeaks(X, Intensities,
...'ShowPlot', ShowPlotValue, ...)
mspeaks(X, Intensities,
...'Style', StyleValue, ...)
Peaklist = mspeaks(X, Intensities)Peaklist, двухколоночная матрица, содержащая значение оси разделения и интенсивность для каждого пика. X является вектором значений единиц разделения для набора сигналов с peaks. Intensities является матрицей значений интенсивности для набора пиков, которые имеют одну и ту же область значений единиц разделения.
[ возвращает Peaklist, PFWHH]
= mspeaks(X, Intensities)PFWHH, двухколоночная матрица, указывающая левое и правое положения полной ширины на маркерах половинной высоты (FWHH) для каждого пика. Для любого пика, не разрешенного в FWHH, mspeaks возвращает границы формы пика. Когда Intensities включает несколько сигналов, затем PFWHH - массив ячеек из матриц.
[ возвращает Peaklist, PFWHH, PExt]
= mspeaks(X, Intensities)PExt, матрицу из двух столбцов, указывающую левое и правое местоположения границ формы пика, определенных после шумоподавления вейвлета. Когда Intensities включает несколько сигналов, затем PExt - массив ячеек из матриц.
mspeaks вызывает (X, Intensities... 'PropertyName', PropertyValue, ...)mspeaks с необязательными свойствами, которые используют пары имя/значение свойства. Можно задать одно или несколько свойств в любом порядке. Заключайте каждую PropertyName в одинарных кавычках. Каждый PropertyName является нечувствительным к регистру. Эти имена свойства/пары значения свойств следующие:
mspeaks( задает основу вейвлета. X, Intensities,
...'Base', BaseValue, ...)
mspeaks( задает количество уровней для разложения вейвлета. X, Intensities,
...'Levels', LevelsValue, ...)
mspeaks( задает метод оценки порога, X, Intensities,
...'NoiseEstimator', NoiseEstimatorValue,
...)T, для фильтрации шумных компонентов при первом высокополосном разложении (y_h).
mspeaks( задает пороговую константу умножения. X, Intensities,
...'Multiplier', MultiplierValue, ...)
mspeaks( управляет использованием шумоподавления вейвлет, чтобы сглаживать сигнал. Варианты X, Intensities,
...'Denoising', DenoisingValue, ...)true (по умолчанию) или false.
mspeaks( задает долю высоты пика, используемую для выбора точек, используемых для вычисления значения оси центроида для соответствующего пика. X, Intensities,
...'PeakLocation', PeakLocationValue, ...)PeakLocationValue должно быть значением ≥ 0 и ≤ 1. По умолчанию это 1.0.
mspeaks( задает минимальную полную ширину на половине высоты (FWHH), в единицах разделения, для сообщаемого peaks. Peaks с FWHH ниже этого значения исключаются из списка выхода X, Intensities,
...'FWHHFilter', FWHHFilterValue, ...)Peaklist.
mspeaks( задает минимальное расстояние между соседним peaks в единицах разделения. Когда сигнал не сглаживается должным образом, может оказаться, что несколько максимумов представляют один и тот же пик. Увеличьте это значение фильтра, чтобы объединить перегранный peaks в один пик.X, Intensities,
...'OverSegmentationFilter', OverSegmentationFilterValue,
...)
mspeaks( задает минимальную высоту для сообщаемого peaks. Peaks с высотами ниже этого значения исключаются из списка выхода X, Intensities,
...'HeightFilter', HeightFilterValue, ...)Peaklist.
mspeaks( управляет отображением графика исходного и сглаженного сигналов с peaks, включенными в выходную матрицу X, Intensities,
...'ShowPlot', ShowPlotValue, ...)Peaklist отмечен.
mspeaks( задает стиль маркировки peaks на графике. X, Intensities,
...'Style', StyleValue, ...)
mspeaks находит peaks в данных любого метода разделения, который формирует сигнальные данные, такие как спектроскопия, ядерный магнитный резонанс (ЯМР), электрофорез, хроматография или масс-спектрометрия.
|
Вектор значений единиц разделения для набора сигналов с peaks. Количество элементов в векторе равняется количеству строк в матрице |
|
Матрица значений интенсивности для набора пиков, которые имеют одну и ту же область значений единиц разделения. Каждая строка соответствует значению единицы разделения, и каждый столбец соответствует либо набору сигналов с peaks, либо времени удержания. Количество строк равняется количеству элементов в векторе |
|
Целое число от По умолчанию: |
|
Целое число от По умолчанию: |
|
Вектор символов, строка или скаляр, который задает метод для оценки порога,
|
|
Положительное действительное значение, которое задает пороговую константу умножения. По умолчанию: |
|
Управляет использованием шумоподавления вейвлета, чтобы сглаживать сигнал. Варианты Совет Если ваши данные были ранее сглажены, например, с
|
|
Значение, которое задает долю высоты пика, используемую для выбора точек, для вычисления значения оси разделения центроида соответствующего пика. Значение должно быть Примечание Когда По умолчанию: |
|
Положительное действительное значение, которое задает минимальную полную ширину на половине высоты (FWHH), в единицах разделения, для сообщаемого peaks. Peaks с FWHH ниже этого значения исключаются из списка выхода По умолчанию: |
|
Положительное действительное значение, которое задает минимальное расстояние в единицах разделения между соседним peaks. Когда сигнал не сглаживается должным образом, может оказаться, что несколько максимумов представляют один и тот же пик. Увеличьте это значение фильтра, чтобы объединить перегранный peaks в один пик. По умолчанию: |
|
Положительное действительное значение, которое задает минимальную высоту для сообщаемого peaks. По умолчанию: |
|
Управляет отображением графика исходного сигнала и сглаженного сигнала с peaks, включенными в выходную матрицу
|
|
Вектор символов или строка, задающая стиль для маркировки peaks на графике. Варианты:
|
|
Двухколоночная матрица, где каждая строка соответствует пику. Первый столбец содержит значения единиц разделения (указывающие на расположение peaks вдоль оси разделения). Второй столбец содержит значения интенсивности. Когда |
|
Двухколоночная матрица, указывающая левое и правое положения полной ширины на маркерах половинной высоты (FWHH) для каждого пика. Для любого пика, не разрешенного в FWHH, |
|
Двухколоночная матрица, указывающая левое и правое положения границ формы пика, определенных после шумоподавления вейвлета. Когда |
Загрузите MAT-файл, включенный в программное обеспечение Bioinformatics Toolbox™, который содержит две переменные данных масс-спектрометрии, MZ_lo_res и Y_lo_res. MZ_lo_res является вектором значений m/z для множества спектров. Y_lo_res является матрицей значений интенсивности для множества массовых спектров, которые имеют одинаковые области значений m/z.
load sample_lo_resОтрегулируйте базовую линию восьми спектров, сохраненных в Y_lo_res.
YB = msbackadj(MZ_lo_res,Y_lo_res);
Преобразовать необработанные данные масс-спектрометрии в список пиков путем нахождения соответствующего peaks в каждом спектре.
P = mspeaks(MZ_lo_res,YB);
Постройте график третьего спектра в YB, матрица скорректированных по базовому уровню значений интенсивности с отмеченным обнаруженным peaks.
P = mspeaks(MZ_lo_res,YB,'SHOWPLOT',3);
Сглаживайте сигнал, используя mslowess функция. Затем преобразуйте сглаженные данные в список пиков путем нахождения соответствующего peaks и построения третьего спектра.
YS = mslowess(MZ_lo_res,YB,'SHOWPLOT',3);
P = mspeaks(MZ_lo_res,YS,'DENOISING',false,'SHOWPLOT',3);
Используйте cellfun функция для удаления всего peaks со значениями m/z менее 2000 из восьми пиков, перечисленных в выходных P. Затем постройте график peaks третьего спектра (в красном) по его сглаженному сигналу (в синем).
Q = cellfun(@(p) p(p(:,1)>2000,:),P,'UniformOutput',false); figure plot(MZ_lo_res,YS(:,3),'b',Q{3}(:,1),Q{3}(:,2),'rx') xlabel('Mass/Charge (M/Z)') ylabel('Relative Intensity') axis([0 20000 -5 95])
mspeaks преобразует необработанные пиковые данные в пиковый список (центроидальные данные) путем:
Сглаживание сигнала с помощью неопределенного преобразования вейвлета с коэффициентами Daubechies
Назначение пиковых местоположений
Оценка шума
Исключение peaks, которые не удовлетворяют заданным критериям
[1] Morris, J.S., Coombes, K.R., Koomen, J., Baggerly, K.A., and Kobayash, R. (2005) Редукция данных и количественная оценка для масс-спектрометрии в биомедицинских применениях с использованием среднего спектра. Биоинфоматика 21:9, 1764-1775.
[2] Ясуй, Я., Пепе, М., Томпсон, М.Л., Адам, Б.Л., Райт, Г.Л., Цюй, Я., Поттер, Ж.Д., Winget, M., Thornquist, M. and Feng, Z. (2003) Data-analytic strategy for belek biomarker discovery: profiling of high-dimensional proteomic data for cancer detection. Биостатистика 4:3, 449-463.
[3] Donoho, D.L., and Johnstone, I.M. (1995) Адаптация к неизвестной гладкости через усадку вейвлета. J. Am. Statist. Asso. 90, 1200–1224.
[4] Strang, G. and Nguyen, T. (1996) Wavelets and Filter Banks (Wellesley: Cambridge Press).
[5] Coombes, K.R., Tsavachidis, S., Моррис, J.S., Baggerly, K.A., Хун, Член конгресса, и Kuerer, H.M. (2005) Улучшенное пиковое обнаружение и определение количества данных масс-спектрометрии, приобретенных от увеличенной поверхностью лазерной десорбции и ионизации спектрами шумоподавления с неподкошенным дискретным вейвлетом, преобразовывают. Протеомика 5 (16), 4107-4117 .
msalign | msbackadj | msdotplot | msheatmap | mslowess | msnorm | mspalign | mspeaks | msppresample | msresample | mssgolay | msviewer