2021-07-14线性拟合

 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")