Как сделать так, чтобы линия тренда на графике рассеяния соответствовала границам оси x?

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

Проблема в том, что, похоже, не соблюдается xlim = c(0, 20) в области графика, поскольку он выходит за пределы оси x. Я пробовал поиграть с опцией xpd, но безуспешно. Затем я попытался поиграть с различными par()options, но не нашел ничего, что могло бы помочь с этой проблемой.

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

График с нежелательным выводом

MWE - ПРИМЕЧАНИЕ. Я привожу только 15 точек данных, чтобы проиллюстрировать проблему, поэтому график будет отличаться от предоставленного изображения.

df.data <- data.frame(id = 1:15,
                      ll = c(-9.53026, -6.50640,-6.50640, -7.68535, -11.80899, -8.42790,
                             -6.50640, -6.50640, -7.92405, -6.50640, -8.95522, -9.99228,
                             -10.02286, -8.95969, -6.07313),
                      aspm = c(4.582104, 0.490244, 0.737765, 0.256699, 1.575931, 1.062693,
                               1.006984, 0.590355, 1.014370, 0.924855, 0.735989, 0.831025,
                               1.197886, 1.143220, 0.928068))

str.col.light.blue <- c(rgb(r = 110/255, g = 155/255, b = 225/255))
str.col.dark.blue <- c(rgb(r = 50/255, g = 100/255, b = 185/255))

layout(matrix(c(2, 4, 1, 3), 2, 2, byrow = TRUE), widths = c(5, 2), heights = c(2, 5))
layout.show(4)

par(omi = c(0.1, 0.1, 0.1, 0.1))
par(mar = c(2, 2, 0, 0))
par(mai = c(1, 1, 0, 0))

plot(df.data[, "ll"] ~ df.data[, "aspm"], col  = str.col.light.blue,
    xlim = c(0, 20), ylim = c(-15, -5), axes = FALSE,
    xlab = "X1", ylab = "X2",
    cex.lab = 1.25)

abline(a = -8.156670, b = -0.000879, lty = 5, col = "black", lwd = 2, xpd = FALSE)

axis(1, at = seq(0, 20, by = 5), labels = seq(0, 20, by = 5), cex.axis = 1)
axis(2, at = seq(-15, -5, by = 3), labels = seq(-15, -5, by = 3), cex.axis = 1, las = 1)

rect(0, -15, 20, log(1/3)*8, density = 10, angle = 45, lwd = 0.5, col = "gray")

par(mar = c(0, 2, 0, 0))
par(mai = c(0, 1, 0.25, 0))

x.hist <- hist(df.data[, "aspm"], plot = FALSE, breaks = 20)
barplot(x.hist$density, axes = FALSE, horiz = FALSE, space = 0, col = str.col.dark.blue)

par(mar = c(2, 0, 0, 0))
par(mai = c(1, 0, 0, 0.25))

y.hist <- hist(df.data[, "ll"], plot = FALSE, breaks = 20) 
barplot(y.hist$density, axes = FALSE, horiz = TRUE, space = 0, col = str.col.dark.blue)

r plot trendline
person edsandorf    schedule 05.09.2017    source источник
comment
Команда par(mfg) выдает ошибку Error in par(mfg = c(0, 0, 0, 0)) : invalid value specified for graphical parameter "mfg". Что касается вопроса, может быть, вы могли бы использовать segments вместо abline?   -  person Rui Barradas    schedule 05.09.2017
comment
Вы правы - par(mfg) выдала ошибку. Этого не было, когда я создавал MWE, но когда я запускал его в чистом сеансе, он работал. Похоже, что это не имеет видимого эффекта, поэтому я отредактировал MWE и удалил эту строку. Насчет использования segments() - возможно. Однако это требует от меня определения начальной и конечной точек линии. Теперь легко, когда линия тренда плоская, но при изменении наклона мне потребовалось бы пересчитать начальную и конечную точки, что было бы немного громоздко, поскольку я создаю несколько таких графиков с разными наклонами для линий тренда.   -  person edsandorf    schedule 05.09.2017

Ответы (1)


arrow_upward
0
arrow_downward

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

linear <- function(x, a, b) a + b*x

Затем я использовал ваш код со следующими изменениями. abline был заменен на segments со всеми графическими параметрами, которые вы использовали в исходном вызове.

x0 <- 0
y0 <- linear(x0, a = -8.156670, b = -0.000879)
x1 <- 20
y1 <- linear(x1, a = -8.156670, b = -0.000879)

segments(x0, y0, x1, y1, lty = 5, col = "black", lwd = 2, xpd = FALSE)

Этот вызов segment был размещен там, где было abline.
На последнем графике я вижу сегмент с хорошим поведением.

person Rui Barradas    schedule 05.09.2017
comment
Идеально! Блестяще решает проблему. - person edsandorf; 05.09.2017