library(dplyr)
## 
## Attaching package: 'dplyr'
## The following objects are masked from 'package:stats':
## 
##     filter, lag
## The following objects are masked from 'package:base':
## 
##     intersect, setdiff, setequal, union
library(ggplot2)

#Dimicks Trial

DimicksTrial %>%
  group_by(Tgt_Rate_s) %>%
  summarize(meanYield = mean(Yld_Vol_Dr))
ggplot(DimicksTrial, aes(Tgt_Rate_s, Yld_Vol_Dr)) +
  geom_point() +
  geom_smooth(method = "lm", color = "purple") +
  ggtitle("Yield vs Applied Rate (Dimicks)") +
  scale_y_continuous(limits = c(0,200)) +
  ##scale_x_continuous(limits = c(120,140)) +
  xlab("Applied Rate") +
  ylab("Yield (bu/ac)")
## `geom_smooth()` using formula = 'y ~ x'
## Warning: Removed 2 rows containing non-finite values (`stat_smooth()`).
## Warning: Removed 2 rows containing missing values (`geom_point()`).

#Marv East Trial

MarvETrial %>%
  group_by(Tgt_Rate_s) %>%
  summarize(meanYield = mean(Yld_Vol_Dr))
ggplot(MarvETrial, aes(Tgt_Rate_s, Yld_Vol_Dr)) +
  geom_point() +
  geom_smooth(method = "lm", color = "orange") +
  ggtitle("Yield vs Applied Rate (Marv East)") +
  ##scale_y_continuous(limits = c(0,120)) +
  ##scale_x_continuous(limits = c(120,140)) +
  xlab("Applied Rate") +
  ylab("Yield (bu/ac)")
## `geom_smooth()` using formula = 'y ~ x'

ggplot(MarvETrial, aes(Tgt_Rate_s, Yld_Vol_Dr)) +
  geom_point() +
  geom_smooth(method = "lm", color = "blue") +
  ggtitle("Yield vs Applied Rate") +
  scale_y_continuous(limits = c(0,200)) +
  ##scale_x_continuous(limits = c(120,140)) +
  xlab("Applied Rate") +
  ylab("Yield (bu/ac)")
## `geom_smooth()` using formula = 'y ~ x'
## Warning: Removed 17 rows containing non-finite values (`stat_smooth()`).
## Warning: Removed 17 rows containing missing values (`geom_point()`).