msnorm

Нормируйте набор сигналов с пиками

Описание

пример

yOut = msnorm(XIntensities) нормирует группу сигналов с peaks путем стандартизации области под кривой (AUC) к медиане группы и возвращает нормированные данные yOut.

пример

[yOut,normParams] = msnorm(XIntensities) также возвращает параметры нормализации normParams то, что можно использовать, чтобы нормировать другую группу сигналов.

пример

[___] = msnorm(XIntensities,Name,Value) дополнительные опции использования, заданные одним или несколькими аргументами пары "имя-значение" и, возвращают любой выходной аргумент в предыдущих синтаксисах. Например, out = msnorm(X,Y,'Quantile',[0.9 1]) устанавливает ниже (0.9) и верхний (1) предел квантиля использовать только самые большие 10% интенсивности в каждом сигнале вычислить AUC.

Примеры

свернуть все

В этом примере показано, как нормировать область под кривой каждого массового спектра из массовых данных о спецификации.

Загрузите MAT-файл, включенный с программным обеспечением Bioinformatics Toolbox™, которое содержит демонстрационные массовые данные о спецификации, включая MZ_lo_res, вектор m/z значений, и Y_lo_res, матрицу значений интенсивности.

load sample_lo_res

Создайте подмножество (четыре сигнала) данных.

MZ = MZ_lo_res;
Y = Y_lo_res(:,[1 2 5 6]);

Постройте эти четыре спектра.

plot(MZ, Y)
axis([-1000 20000 -20 105])
xlabel('Mass-charge Ratio')
ylabel('Relative Ion Intensities')
title('Original Spectra')

Нормируйте область под кривой (AUC) каждого спектра к медиане, устранив малую массу (m/z <1,000) шум, и постповторно масштабировавшись таким образом, что максимальная интенсивность равняется 100. Постройте эти четыре спектра.

Y1 = msnorm(MZ,Y,'Limits',[1000 inf],'Max',100);
plot(MZ, Y1)
axis([-1000 20000 -20 105])
xlabel('Mass-charge Ratio')
ylabel('Relative Ion Intensities')
title('AUC Normalized Spectra')

В этом примере показано, как нормировать ионную интенсивность каждого спектра из массовых данных о спецификации.

Загрузите MAT-файл, включенный с программным обеспечением Bioinformatics Toolbox™, которое содержит демонстрационные массовые данные о спецификации, включая MZ_lo_res, вектор m/z значений, и Y_lo_res, матрицу значений интенсивности.

load sample_lo_res

Создайте подмножество (четыре сигнала) данных.

MZ = MZ_lo_res;
Y = Y_lo_res(:,[1 2 5 6]);

Нормируйте ионную интенсивность каждого спектра к максимальной интенсивности одного самого высокого пика от любого из спектров в области значений выше 1000 m/z. Постройте эти четыре спектра.

Y2 = msnorm(MZ,Y,'QUANTILE', [1 1],'LIMITS',[1000 inf]);
plot(MZ, Y2)
axis([-1000 20000 -20 105])
xlabel('Mass-charge Ratio')
ylabel('Relative Ion Intensities')
title('Maximum-Intensity Normalized Spectra')

В этом примере показано, как выполнить нормализацию квантиля для массовых данных о спецификации.

Загрузите MAT-файл, включенный с программным обеспечением Bioinformatics Toolbox™, которое содержит демонстрационные массовые данные о спецификации, включая MZ_lo_res, вектор m/z значений, и Y_lo_res, матрицу значений интенсивности.

load sample_lo_res

Создайте подмножество (четыре сигнала) данных.

MZ = MZ_lo_res;
Y = Y_lo_res(:,[1 2 5 6]);

Нормируйте использование данных в m/z областях, где интенсивность в четвертом квартиле по крайней мере в 90% спектрограмм. Обратите внимание на то, что можно использовать параметры нормализации во втором выходе, чтобы нормировать другой набор данных в тех же m/z областях. Постройте эти четыре спектра.

[Y3,S] = msnorm(MZ,Y,'Quantile',[0.75 1],'Consensus',0.9);
area(MZ,S.Xh.*1000,'LineStyle','None','FaceColor',[.8 .8 .8])
hold on
plot(MZ, Y3)
hold off
axis([-1000 20000 -20 105])
xlabel('Mass-charge Ratio')
ylabel('Relative Ion Intensities')
title('Fourth-quartile Normalized Spectra')

Используйте параметры нормализации во втором выходе предыдущего шага, чтобы нормировать различное подмножество данных (четыре сигнала) использование данных в тех же m/z областях как предыдущий набор данных. Постройте эти четыре спектра.

Y4 = msnorm(MZ,Y_lo_res(:,[3 4 7 8]),S);
 
area(MZ,S.Xh.*1000,'LineStyle','None','FaceColor',[.8 .8 .8])
hold on
plot(MZ, Y4)
hold off
axis([-1000 20000 -20 105])
xlabel('Mass-charge Ratio')
ylabel('Relative Ion Intensities')
title('Fourth-quartile Normalized Spectra')

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

свернуть все

Вектор разделительных стоимостей единицы для набора сигналов с peaks, заданным как вектор.

Типы данных: double

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

Типы данных: double

Аргументы в виде пар имя-значение

Задайте дополнительные разделенные запятой пары Name,Value аргументы. Name имя аргумента и Value соответствующее значение. Name должен появиться в кавычках. Вы можете задать несколько аргументов в виде пар имен и значений в любом порядке, например: Name1, Value1, ..., NameN, ValueN.

Пример: out = msnorm(X,Y,'Quantile',[0.9 1])

Пределы квантиля, чтобы уменьшать набор разделительных стоимостей единицы в X, заданный как 1- 2 вектор или скаляр между 0 и 1.

Если вы задаете вектор, первый элемент является нижним пределом, и второй элемент является верхним пределом. Например, [0.9 1] средние значения, что функция использует только самые большие 10% интенсивности в каждом сигнале вычислить AUC. Значение по умолчанию [0 1] средние значения, что функция использует целый AUC, вместо того, чтобы ограничить интенсивность конкретным квантилем.

Если вы задаете скалярное значение, оно представляет более низкий предел квантиля. Верхний предел квантиля автоматически устанавливается к 1.

Пример: 'Quantile',[0.8 1]

Типы данных: double

Область значений разделительного модуля, чтобы выбрать точки нормализации, заданные как 1- 2 вектор. Значение по умолчанию [min(X) max(X)] выбирает все доступные точки из X. Если вы задаете нижний или верхний предел значению, которое не является в доступной области значений [min(X) max(X)], функция устанавливает нижний предел на min(X) и верхний предел max(X), соответственно.

Этот параметр полезен, чтобы устранить шум из вычисления AUC. Например, можно исключить матричный шум, который появляется в области малой массы (m/z <1000) из массового спектрометра SELDI путем установления предела к [1000 max(X)].

Пример: 'Limits',[900 max(X)]

Типы данных: double

Минимальный процент значений интенсивности в пределах квантиля, что положение разделительного модуля, должно быть, придется включить в вычисление AUC, задал как скаляр между 0 и 1. Те же положения разделительного модуля затем используются, чтобы нормировать все сигналы. Используйте этот параметр, чтобы устранить peaks низкой интенсивности и шум от нормализации.

Например, чтобы выбрать области MZ, интенсивность которых в третьем квантиле по крайней мере в 90% спектрограмм, устанавливает 'Quantile' и 'Consensus' к следующему: yOut = msnorm(MZ,Y,'Quantile',[0.5 0.75],'Consensus',0.9).

Пример: 'Consensus',0.8

Типы данных: double

Метод для нормализации AUC каждого сигнала, заданного как 'Median' или 'Mean'.

Пример: 'Method','Mean'

Типы данных: char | string

Общая максимальная интенсивность, чтобы масштабироваться к после нормализации каждого сигнала индивидуально, заданный как скаляр. Если вы не задаете этот параметр, никакое постмасштабирование не выполняется.

Примечание

Если вы задаете это значение и также устанавливаете 'Quantile' к [1 1], затем одна точка (пиковая высота самого высокого пика) нормирована к заданному максимальному значению.

Пример: 'Max'

Типы данных: double

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

свернуть все

Нормированные значения интенсивности, возвращенные как матрица.

Параметры нормализации, которые можно использовать, чтобы нормировать другую группу сигналов, возвратились как структура.

Представлено до R2006a

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