É possível reproduzir esse gráfico de Regressão no R ou Python?
Obs. Esse gráfico foi produzido usando o objeto QAQC do acQuire 4. Link
E esse é o modelo que que criei até o momento.
dados = read.csv("C:\\Users\\....\\Desktop\\data3.csv", header = T)
# 0. Build linear model
model <- lm(CK ~ OR, data = dados)
summary(model)
# 1. Add predictions
pred.int <- predict(model, interval = "prediction")
mydata <- cbind(dados, pred.int)
# 2. Regression line
library("ggplot2")
p <- ggplot(mydata, aes(OR, CK)) +
geom_point() +
stat_smooth(method = lm)
# 3. Add prediction intervals
p + geom_line(aes(y = lwr), color = "red", linetype = "dashed")+
geom_line(aes(y = upr), color = "red", linetype = "dashed")
Minha base de dados: Link