Exponenta Banner

dimercount

Считайте димеры в последовательности нуклеотида

свернуть все на странице

Синтаксис

Dimers = dimercount(SeqNT)
[Dimers, Percent] = dimercount(SeqNT)
... = dimercount(SeqNT, 'Ambiguous', AmbiguousValue)
... = dimercount(SeqNT, 'Chart', ChartValue)

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

SeqNT

Одно из следующего:

Примеры: 'ACGT' или [1 2 3 4]

AmbiguousValue

Вектор символов или строка, задающая, как обработать димеры, содержащие неоднозначные символы нуклеотида (R, Y, K, M, S, W, B, D, H, V или N). Выбор:

  • 'ignore' (значение по умолчанию) — Димеры пропусков, содержащие неоднозначные символы

  • 'bundle' — Димеры количеств, содержащие неоднозначные символы и отчеты общее количество в поле Ambiguous Dimers структура вывода.

  • 'prorate' — Димеры количеств, содержащие неоднозначные символы и, распределяют их пропорционально в соответствующих полях димера, содержащих стандартные символы нуклеотида. Например, счета для димера AR распределяются равномерно между полями AA и AG.

  • 'warn' — Димеры пропусков, содержащие неоднозначные символы и, выводят предупреждение.

ChartValue Вектор символов или строка, задающая тип диаграммы. Выбором является 'pie' или 'bar'.

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

DimersСтруктура MATLAB, содержащая поля AA, AC, AG, AT, CA, CC, CG, CT, GA, GC, GG, GT, TA, TC, TG и TT, которые содержат количества димера в SeqNT.
PercentМатрица 4 на 4 с относительными пропорциями димеров в SeqNT. Строки соответствуют A, C, G и T в первом элементе димера, и столбцы соответствуют A, C, G и T во втором элементе димера.

Описание

Dimers = dimercount(SeqNT) считает димеры нуклеотида в SeqNT, последовательности нуклеотида, и возвращает количества димера в Dimers, структура MATLAB, содержащая поля AA, AC, AG, AT, CA, CC, CG, CT, GA, GC, GG, GT, TA, TC, TG и TT.

  • Для последовательностей, которые имеют димеры с символьным U, эти димеры добавляются к соответствующим димерам, содержащим T.

  • Если последовательность содержит разрывы, обозначенные дефисом (-), разрывы проигнорированы, и эти два символа по обе стороны от разрыва считаются как димер.

  • Если последовательность содержит нераспознанные символы, то димеры, содержащие эти символы, проигнорированы, и следующее предупреждающее сообщение появляется:

    Warning: Unknown symbols appear in the sequence. These will be ignored.

[Dimers, Percent] = dimercount(SeqNT) возвращает Percent, матрицу 4 на 4 с относительными пропорциями димеров в SeqNT. Строки соответствуют A, C, G и T в первом элементе димера, и столбцы соответствуют A, C, G и T во втором элементе димера.

... = dimercount(SeqNT, 'Ambiguous', AmbiguousValue) задает, как обработать димеры, содержащие неоднозначные символы нуклеотида. Выбор:

  • 'ignore' (значение по умолчанию)

  • 'bundle'

  • 'prorate'

  • 'warn'

... = dimercount(SeqNT, 'Chart', ChartValue) создает график, показывающий относительные пропорции димеров. ChartValue может быть 'pie' или 'bar'.

Примеры

свернуть все

Скрипт Open Live Script
seq = randseq(100)
seq = 
'TTATGACGTTATTCTACTTTGATTGTGCGAGACAATGCTACCTTACCGGTCGGAACTCGATCGGTTGAACTCTATCACGCCTGGTCTTCGAAGTTAGCAC'

[Dimers, Percent] = dimercount(seq)
Dimers = struct with fields:
    AA: 4
    AC: 9
    AG: 3
    AT: 6
    CA: 3
    CC: 3
    CG: 8
    CT: 9
    GA: 8
    GC: 4
    GG: 4
    GT: 6
    TA: 7
    TC: 8
    TG: 7
    TT: 10


Percent = 4×4

    0.0404    0.0909    0.0303    0.0606
    0.0303    0.0303    0.0808    0.0909
    0.0808    0.0404    0.0404    0.0606
    0.0707    0.0808    0.0707    0.1010

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

| | | | |

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

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