Можно ли раскрасить точки данных в воронкообразном графике на основе категории в пакете Metafor / R?

Question

Можно ли покрасить точки данных в графике воронки?Я использую пакет метафор.

 # Load package
 library(metafor)

 # Load data
 data(dat.bcg)

 # Code 
 ex <- rma(ai=tpos, bi=tneg, ci=cpos, 
 di=cneg, data=dat.bcg, measure="OR",
 slab=paste(author, year, sep=", "), method="FE")

# Funnel plot
funnel(ex, transf = exp, ylab="Test group", xlab="Control")

Можно ли раскрасить точки на основе переменной?

Fx синего цвета точки данных, которые представляют:

dat.bcg$alloc==random

Спасибо, С.

Answer 1

Я изучил исходный код метода funnel и раздраженно, нет, вы не можете передать вектор цветов в аргумент col.

Простой обходной путь - назначить вывод funnel() переменной, которая даст вам координаты x и y точек на воронкообразном графике, а затем нанесите поверх исходных точек на воронкообразном графике цветные точки.

# Create vector of colors
my_colors <- c('red','blue')[(dat.bcg$alloc == 'random') + 1]

# Create funnel plot, catching output
funnelplotdata <- funnel(ex, transf = exp, ylab="Test group", xlab="Control")

# Plot over points
with(funnelplotdata, points(x, y, col = my_colors, pch = 19))