Overlay boxplot and scatter.smooth line
Same thing using ggplot... # your data x <- c(rep(1,100),rep(2,100),rep(3,100)) y <- c(rnorm(100,1),rnorm(100,2),rnorm(100,3)) df <- data.frame(x=x, y=y) library(ggplot2) ggp <- ggplot(df) + labs(x="X", y="Y") ggp <- ggp + geom_boxplot(aes(x=factor(x), y=y)) ggp <- ggp + geom_smooth(formula=y~x, aes(x=x, y=y), method="loess", se=F) ggp Setting method="lm" above would plot linear model instead of loess, setting se=T (the default) would plot 95% CL for the model. Adding the line: ggp <- ggp + stat_summary(fun.data="mean_cl_normal", aes(x=factor(x),y=y), size=1.2, color="red") ggp will overlay mean and 95% CL (assuming normal distribution of error) at each level of x. -----Original Message----- From: Johannes Radinger [mailto:johannesradinger at gmail.com] Sent: Tuesday, November 19, 2013 10:01 AM To: Adams, Jean Cc: R help Subject: Re: [R] Overlay boxplot and scatter.smooth line Thanks...that works great! /J
On Tue, Nov 19, 2013 at 2:39 PM, Adams, Jean <jvadams at usgs.gov> wrote:
Use lines() with loess.smooth() instead of scatter.smooth() to add to the already existing boxplot. For example, boxplot(y~x) lines(loess.smooth(x, y)) Jean On Tue, Nov 19, 2013 at 4:43 AM, Johannes Radinger < johannesradinger at gmail.com> wrote:
Hi,
I'd like to plot a boxplot and use a scatter.smooth line plot as an
overlay for this plot.
This works except for the problem that the x-axis ticks for both
plots are differently spaced:
The boxplot ticks are closer and the space between the outer most
boxplots and the plot region (box) is larger for boxplots than for
the scatter.smooth plot. Is there any plot-option that can be changed
to produce the desired plot.
Here an example to illustrate what problem I am facing:
x <- c(rep(1,100),rep(2,100),rep(3,100))
y <- c(rnorm(100,1),rnorm(100,2),rnorm(100,3))
boxplot(y~x)
par(new=TRUE)
scatter.smooth(x, y, xaxt="n", yaxt="n", ann=FALSE)
Moreover, I'd like to plot just the smoothed line only and not the
datapoints.
Any suggestions?
best regards,
J.
[[alternative HTML version deleted]]
______________________________________________ R-help at r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code.