R: geom_lo go для нескольких генов из tibble - PullRequest
Новая
       31

R: geom_lo go для нескольких генов из tibble

2 голосов
/ 27 февраля 2020

Я не очень знаком с пакетом ggseqlo go, поэтому буду признателен за любую помощь.

Я приготовил тибль, который выглядит так: 

test <- tibble( gene = c("A", "B", "A", "C", "B"),
                seq = c("AAAAAAAAAAAAAAAAAAAA",
                        "GGGGGGGGGGGGGGGG",
                        "AAAAAATAAAAATAAAAAAA",
                        "AGTCGTCATGCATCAATCCCAATGGTGCA",
                        "GGGGGGGCCGGGGGGG") )

Я хочу подготовить seqlo go для каждого гена, основываясь на названии гена. Для каждого гена последовательности имеют одинаковую длину.

Насколько я попробовал это: 

ggplot() + 
 geom_logo(data = test$gene) +
 facet_grid(rows = ~ gene)

Но пока это лучшее, что я получил:

...