library(ggplot2)
library(readxl)
library(dplyr)
ss_1 <- read_excel("SS浊度.xlsx")
head(ss_1)
ss_2 <- ss_1 %>% select(2,3)
head(ss_2)
lm.ss_2<-lm(turbidity ~ 0+conc_mg_L, data= ss_2)
summary(lm.ss_2)
summary(lm.ss_2)$coefficients
fit_x <- seq(0,1000,length.out = 1000)
fit_y <- 0.6767869*fit_x
d_fit <- data.frame(conc_mg_L=fit_x, turbidity=fit_y)
head(d_fit)
ggplot(ss_2, aes(conc_mg_L, turbidity))+
theme_bw()+
stat_summary(geom ="point",fun="mean",size=1, color="gray2")+
stat_summary(geom = "errorbar",fun.min = function(x) mean(x)-sd(x),
fun.max = function(x) mean(x)+sd(x),
width=2)+
geom_line(data = d_fit,color="red2")+
labs(x=expression("SS concentration (mg"~ L^"-1"*")"),
y="Turbidity")+
annotate(geom="text", x=900, y=150,label=expression("y=0.68x"),size=3)
ggsave("ss_t.png", width=280/90, height=224/90, dpi=600, units="in")