Как нарисовать сюжет скрипки с цветом, показывающим выражение генной ценности? - PullRequest
0 голосов
/ 23 декабря 2018

Я пытаюсь представить генную экспрессию «гена А» среди нескольких групп.Я использую ggplot2 для рисования, но у меня ничего не получается

p <- ggplot(MAPK_plot, aes(x = group, y = gene_A)) + geom_violin(trim = FALSE , aes( colour = gene_A))  + theme_classic()

И я хочу получить такую ​​фигуру, как this из https://www.researchgate.net/publication/313728883_Neuropilin-1_Is_Expressed_on_Lymphoid_Tissue_Residing_LTi-like_Group_3_Innate_Lymphoid_Cells_and_Associated_with_Ectopic_Lymphoid_Aggregates

1 Ответ

0 голосов
/ 23 декабря 2018

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

Хорошая работа по выяснению того, что выдолжен использовать geom_violin.Далее вам понадобится некоторая форма огранки / мультипанели.Наконец, чтобы сделать полную аннотацию, как в приведенной ссылке, вам необходимо использовать функциональность пакета grid (которую я здесь не использую).

Я не знаком с наборами данных генной экспрессии, но я использую набор данных рейтинга фильмов IMDB для этого примера (хранится в пакете ggplot2movies).

library(ggplot2)

library(ggplot2movies)
library(data.table)
mv <- copy(movies)
setDT(mv)
# make some variables for our plotting example
mv[, year_10 := cut_width(year, 10)]
mv[, rating_10yr_avg := mean(rating), by = year_10]
mv[, length_3gr := cut_number(length, 3)]

ggplot(mv, 
       aes(x = year_10, 
           y = rating)) + 
  geom_violin(aes(fill = rating_10yr_avg),
              scale = "width") +
  facet_grid(rows = vars(length_3gr))

enter image description here

Пожалуйста, не принимайте этоответьте как форму, поощряющую не публиковать данные, относящиеся к вашей проблеме.

...