Calculating characteristics such as median, mean,… of a subset of data is quite straightforward in R:
For a data set containing results from several “models”, a subset for the model “base” is created by
base_only <- subset(all, all$model == "base")
Then, the median of the variable “data” can be computed via
median(base_only$data)
Here are a few more plotting options for boxplots:
Let’s start plotting the full set
plot(b$mod, b$x)
Plot labels for a subset in full set plot (label all points x < -1)
text(subset(b$mod, b$x < -1), subset(b$x, b$x < -1), subset(b$site, b$x < -1), cex=0.6, pos=4, col="red")
Plot subset with x > -1
plot(subset(b$mod, b$x > -1), subset(b$x, b$x > -1))
Plot horizontal gridlines
grid(nx = NA, ny = NULL)