在R

时间:2017-06-15 13:06:30

标签: r plot legend

我绘制了五个图表和一个图例。图表工作正常,但是legens消失没有错误。 My preview in RStudio looks like this

当我放大时,图例应该是的区域是空白的。 我使用以下代码:

opar <- par (no.readonly = TRUE)
par (mfrow = c(3, 2))

library(deSolve)

# Plot A
LotVmod <- function (Time, State, Pars) {
    with(as.list(c(State, Pars)), {
        dx = (b*x) - (b*x*x/K) - (y*(x^k/(x^k+C^k)*(l*x/(1+l*h*x))))
        dy = (y*e*(x^k/(x^k+C^k)*(l*x/(1+l*h*x)))) - (m*y)
        return(list(c(dx, dy)))
    })
}

Pars <- c(b = 1.080, e = 2.200, K = 130.000, k = 20.000, l = 2.000, 
          h = 0.030, C = 2.900, m = 0.050)

State <- c(x = 0.25, y = 2.75)  
Time <- seq(1, 9, by = 1)

out <- as.data.frame(ode(func = LotVmod, y = State, parms = Pars, times = Time))

matplot(out[,-1], type = "l", xlim = c(1, 9), ylim = c(0, 45),  
        xlab = "time", 
        ylab = "population",
        main = "Compartment A")
mtext ( "Coefficient of Variance 4.96", cex = 0.8 )

x <- c(# Validation data)
y <- c(# Validation data)

lines (Time, x,  type="l", lty=1, lwd=2.5, col="black") 
lines (Time, y, type="l", lty=1, lwd=2.5, col="red")

# Legend
plot.new()
legend("center", c(expression (italic ("F. occidentalis")*" observed"), 
                   expression (italic ("M. pygmaeus")*" observed"), 
                   expression (italic ("F. occidentalis")*" simulated"),
                   expression (italic ("M. pygmaeus")*" simulated")),
       lty = c(1, 1, 1, 2), 
       col = c(1, 2, 1, 2), 
       lwd = c(2.5, 2.5, 1, 1),
       box.lwd = 0, bty = "n")

# Plot C to F = same as A

par(opar)

我的输出没有出错。我之前使用过完全相同的代码没有任何问题,因此我重新启动了R,删除了所有对象,清除了所有的图并重新启动了RStudio和我的计算机。

2 个答案:

答案 0 :(得分:2)

尝试在您的图例语句中添加xpd=TRUE。即。

legend("center", c(expression (italic ("F. occidentalis")*" observed"), 
                   expression (italic ("M. pygmaeus")*" observed"), 
                   expression (italic ("F. occidentalis")*" simulated"),
                   expression (italic ("M. pygmaeus")*" simulated")),
       lty = c(1, 1, 1, 2), 
       col = c(1, 2, 1, 2), 
       lwd = c(2.5, 2.5, 1, 1),
       box.lwd = 0, bty = "n", xpd=TRUE)

默认情况下,图例会被绘图区域截断。此xpd参数可以在绘图区域外绘图。参见例如关于?par的更多信息,请xpd

答案 1 :(得分:2)

这是由于绘图画布的设置方式以及该设备如何重新缩放。按照您的方式,您可以在右上图的绘图区域中添加图例。然而,绘图区域不是完整的设备,而是仅由轴形成的空间内的部分。如果您重新缩放,该绘图区域也将重新缩放。绘图区域周围的边距不会改变大小,因此放大会使您的绘图区域变得如此之小,以至于它不再适合图例。它隐藏在绘图区周围的边缘。

因此,AEBilgrau非常正确,您需要添加xpd = TRUE。这允许图例延伸到绘图区域之外,因此在调整绘图设备大小时它不会消失在边距后面。