У меня есть набор данных микрочипа, и я хотел бы представить его в соответствии с прикрепленным рисунком (сделано с помощью PowerPoint). Я пробовал различные пакеты, доступные для R (VennDiagram, venneuler, limma). Однако мне не удалось нарисовать данные пропорционально или показать значения.
Кто-нибудь знает, как это сделать?
Код для игры можно найти ниже!
Спасибо!
# you might need these:
source("http://www.bioconductor.org/biocLite.R")
biocLite("limma")
library(VennDiagram)
library(limma)
library(venneuler)
# Trial 1: Kind of proportional but difficult to customize for labels and other stuff
su <- venneuler(c(A=162, B=104, C=86, "A&B"=206, "A&C"=112, "B&C"=90 ,"A&B&C"=2433))
plot(su)
# Trial 2: looks nice but not proportional
hw <- c(F,F,F,F,T,T,T,T)
hm <- c(F,F,T,T,F,F,T,T)
hr <- c(F,T,F,T,F,T,F,T)
c4 <- cbind(hw,hm,hr)
e <- vennCounts(c4)
e[1:8,4] <- c(3193,86,104,90,162,112,206,2433)
vennDiagram(e)
venneuler
очень бедна, но вы можете изменить метки с помощьюsu$labels = c("Foo", "Bar", "Baz")
, цвета с помощьюsu$colors = c(0.1, 0.44, 0.3)
и т. д. Однако я не думаю, что возможно нарисовать пропорциональные диаграммы Венна для трех наборов с помощью кругов. Есть только шесть параметров (три радиуса и три стороны треугольника, соединяющие центры кругов), которые пытаются представить семь различных областей. Однако это может быть (?) Возможно с другими формами, например. эллипсы. Я нашел eulerAPE, но мне это не помогло. - person Sparhawk   schedule 06.10.2016