rmabackadj

Выполните фоновую корректировку на данных тестового уровня Affymetrix микромассивов с помощью процедуры Устойчивого среднего значения мультимассивов (RMA)

Синтаксис

BackAdjustedMatrix = rmabackadj(PMData)
BackAdjustedMatrix = rmabackadj(..., 'Method', MethodValue, ...)
BackAdjustedMatrix = rmabackadj(..., 'Truncate', TruncateValue, ...)
BackAdjustedMatrix = rmabackadj(..., 'Showplot', ShowplotValue, ...)

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

PMData

Матрица значений интенсивности, где каждая строка соответствует идеальной паре (PM), зонд и каждый столбец соответствуют файлу Affymetrix® CEL. (Каждый файл CEL сгенерирован от отдельного чипа. Все микросхемы должны иметь тот же тип.)

MethodValue

Задает метод оценки для фоновых параметров модели корректировки. Введите любой 'RMA' (чтобы использовать метод оценки, описанный Bolstad, 2005) или 'MLE' (чтобы оценить параметры с помощью наибольшего правдоподобия). Значением по умолчанию является 'RMA'.

TruncateValue

Задает модель фонового шума. Войдите любой true (используйте усеченное Распределение Гаусса), или false (используйте неусеченное Распределение Гаусса). Значением по умолчанию является true.

ShowplotValue

Управляет графическим выводом гистограммы, показывающей распределение (синих) значений интенсивности зонда PM и замысловатая (красная) функция распределения вероятностей, с предполагаемыми параметрами mu, сигмой и альфой. Введите или 'all' (постройте гистограмму для каждого столбца или чипа), или задайте подмножество столбцов (микросхемы) путем ввода номера столбца, списка чисел или области значений чисел.

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

BackAdjustedMatrixМатрица настроенных фоном тестовых значений интенсивности.

Описание

BackAdjustedMatrix = rmabackadj(PMData) возвращается фон настроил значения тестовых значений интенсивности в матрице, PMData. Обратите внимание на то, что каждая строка в PMData соответствует идеальной паре (PM), зонд и каждый столбец в PMData соответствуют файлу Affymetrix CEL. (Каждый файл CEL сгенерирован от отдельного чипа. Все микросхемы должны иметь тот же тип.) Детали о фоновой корректировке описаны Bolstad, 2005.

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

BackAdjustedMatrix = rmabackadj(..., 'Method', MethodValue, ...) задает метод оценки для фоновых параметров модели корректировки. Когда MethodValue является 'RMA', rmabackadj реализует метод оценки, описанный Bolstad, 2005. Когда MethodValue является 'MLE', rmabackadj оценивает параметры с помощью наибольшего правдоподобия. Значением по умолчанию является 'RMA'.

BackAdjustedMatrix = rmabackadj(..., 'Truncate', TruncateValue, ...) задает используемую модель фонового шума. Когда TruncateValue является false, rmabackadj использует неусеченный Гауссов в качестве модели фонового шума. Значением по умолчанию является true.

BackAdjustedMatrix = rmabackadj(..., 'Showplot', ShowplotValue, ...) позволяет вам построить гистограмму, показывающую распределение (синих) значений интенсивности зонда PM и замысловатая (красная) функция распределения вероятностей, с предполагаемыми параметрами mu, сигмой и альфой. Когда ShowplotValue является 'all', rmabackadj строит гистограмму для каждого столбца или чипа. Когда ShowplotValue является номером, списком чисел или областью значений чисел, rmabackadj строит гистограмму для обозначенного номера столбца (чип).

Например:

  • (..., 'Showplot', 3,...) строит значения интенсивности в столбце 3 PMData.

  • (..., 'Showplot', [3,5,7],...) строит значения интенсивности в столбцах 3, 5 и 7 из PMData.

  • (..., 'Showplot', 3:9,...) строит значения интенсивности в столбцах 3 - 9 PMData.

Примеры

  1. Загрузите MAT-файл, включенный с программным обеспечением Bioinformatics Toolbox™, которое содержит данные тестового уровня Affymetrix, включая pmMatrix, матрицу значений интенсивности зонда PM из нескольких файлов CEL.

    load prostatecancerrawdata
  2. Выполните фоновую корректировку на значениях интенсивности зонда премьер-министра в матрице, pmMatrix, создав новую матрицу, BackgroundAdjustedMatrix.

     BackgroundAdjustedMatrix = rmabackadj(pmMatrix);
  3. Выполните фоновую корректировку на значениях интенсивности зонда премьер-министра только в столбце 3 матрицы, pmMatrix, создав новую матрицу, BackgroundAdjustedChip3.

     BackgroundAdjustedChip3 = rmabackadj(pmMatrix(:,3));

Файл prostatecancerrawdata.mat, используемый в предыдущем примере, содержит данные из Лучшего и др., 2005.

Ссылки

[1] Irizarry, R.A., Хоббс, B., Коллин, F., Бизер-Барклай, Y.D., Antonellis, K.J., Scherf, U., скорость, T.P. (2003). Исследование, нормализация и сводные данные данных об уровне зонда олигонуклеотида высокой плотности массивов. Биостатистика 4, 249–264.

[2] Bolstad, B. (2005). “affy: встроенные методы обработки” https://www.bioconductor.org/packages/2.1/bioc/vignettes/affy/ inst/doc/builtinMethods.pdf

[3] Лучше всего, C.J.M., Гиллеспи, J.W., И, Y., Chandramouli, G.V.R., Perlmutter, M.A., Gathright, Y., Эриксон, H.S., Георгевич, L., Tangrea, M.A., Duray, P.H., Гонсалес, S., Веласко, A., Linehan, W.M., Matusik, R.J., Цена, D.K., Figg, W.D., Emmert-маркер, M.R., и Chuaqui, R.F. (2005). Молекулярные изменения при первичном раке простаты после терапии абляции андрогена. Клинические Исследования рака 11, 6823–6834.

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