# Using the forceR::df.all.200.tax dataset:
# sumarize by measurement and specimen
df.summary.specimen <- summarize_measurements(df = df.all.200.tax,
var1 = "measurement",
var2 = "specimen")
# plot results
if (FALSE) {
require(ggplot2)
ggplot(data = df.summary.specimen, mapping = aes(x=specimen,y=max.F.measurement)) +
geom_jitter(aes(color='blue'),alpha=0.7) +
geom_boxplot(fill="bisque",color="black",alpha=0.3) +
# scale_y_log10() +
labs(y="max(F)/specimen") +
guides(color="none") +
theme_minimal()
}
Run the code above in your browser using DataLab