Jesse001 Jesse001 - 1 year ago 85
R Question

How do I change colours of confidence interval lines when using `matlines` for prediction plot?

I'm plotting a logarithmic regression's line of best fit as well as the confidence intervals around that line. The code I'm using works well enough, except I'd rather that the confidence intervals both be "gray" (rather than the default "red" and "green"). Unfortunately, I'm not seeing a way to isolate them when specifying colour changes. I'd like

  • for the regression line:
    lty = 1, col = "black"

  • for confidence intervals to have:
    lty=2, col = "gray"

How can I achieve this? my code is of the sort:

R6cl <- lm(log(R6$y) ~ R6$x)
pR6cl <- predict(R6cl, interval="confidence")
plot(R6$x, log(R6$y), type = "p")
matlines(x = R6$x, y = log(R6$y), lwd = 2, lty = 1, col = "black")

which produces:

enter image description here

Answer Source

col, lty and lwd are vectorized. You can use

R6cl <- lm(log(y) ~ x, data = R6)  ## don't use $ in formula
pR6cl <- predict(R6cl, interval = "confidence")
plot(log(y) ~ x, data = R6)  ## Read `?plot.formula`
matlines(R6$x, pR6cl, lwd = 2, lty = c(1, 2, 2), col = c(1, 2, 2))

You can check the last figure in Piecewise regression with a quadratic polynomial and a straight line joining smoothly at a break point for what this code would produce.

If you are unclear why I advise against the use of $ in model formula, read Predict() - Maybe I'm not understanding it.

A side notice for other readers

OP has a dataset where x is sorted. If your x is not sorted, make sure you sort it first. See Messy plot when plotting predictions of a polynomial regression using lm() in R for more.

Recommended from our users: Dynamic Network Monitoring from WhatsUp Gold from IPSwitch. Free Download