2017-01-04 3 views
3

Я новичок в R и работать над созданием лесной участок с набором данных с использованием forest.default реализации в metafor пакета '(я только есть x, ci.ub и ci.ub в данных для лесной участок)Добавление пробелов в данных для лесного участка в METAFOR

Я почти закончен, и мой сюжет похож на this, за исключением того, что мои данные проходят черную горизонтальную границу сверху графика. Я работал с этим, но у меня есть один лесной участок всегда на том же ряду, что и на горизонтальной границе, и наложение белой линией также не помогает, как упомянуто here.

Мой вопрос заключается в следующем: Как добавить пробелы в набор данных, чтобы сам лесной участок сдвинул черную рамку вверх? This is the image of the plot currently with the border (currently superimposed with a white line) cutting through my data

Спасибо!

EDIT: Вот фрагмент кода

require(metafor) 
    metafor::forest(
    x = as.vector(t(data_set[, forest_plot$mean])), 
    xlim = c(min_x_value - 2 * (plot_span), max_x_value), 
    alim = c(min_x_value , max_x_value), 
    ci.lb = as.vector(t(data_set[, forest_plot$low])), 
    ci.ub = as.vector(t(data_set[, forest_plot$high])), 
    col = "darkgoldenrod4", 
    ilab = display_data_matrix, 
    ilab.pos = 4, 
    annotate = FALSE, 
    slab = NA, 
    ilab.xpos = c(min_x_value - 2 * plot_span, min_x_value - 1.5 * plot_span, min_x_value - 0.75 * plot_span, min_x_value - 0.25 * plot_span), 
    rows = row_groupings, 
    xlab = paste(measure, ep_type, ep) 
) 
    # print("e") 
    text(min_x_value - 2* plot_span, subgroup_rows, col = "brown3", pos = 4 , subgroup_titles, cex = 1.2) 
    header_line <- dim(data_set)[1] + gaps + 1 
    text(c(min_x_value - 2 * plot_span, 
      min_x_value - 1.5 * plot_span, 
      min_x_value - 0.75 * plot_span, 
      min_x_value - 0.25 * plot_span), 
     pos = 4, 
     cex = 1.4, 
     col = "darkorchid4", 
     header_line, c("Trial", "Randomized Treatment", "N Arm", "Time")) 

    abline(h=dim(data_set)[1]+1, lwd=2, col="white") 

ответ

1

Вам нужно настроить ylim аргумент. Выполните следующие действия:

sav <- forest(x = [...]) 
sav 

(где [...] все остальные части кода). Затем посмотрите, какие значения по умолчанию для ylim. Отрегулируйте это, добавив 2 к верхнему значению (поскольку вы добавили две строки для заголовков подгрупп), а затем используйте ylim=c(lower,upper) в своем обращении к forest().

+0

Спасибо, профессор! Вы действительно живое доказательство того, что статистики могут быть довольно дружелюбны и социально адаптированы! – rushsangs

Смежные вопросы