Я не уверен, как элегантно добавить некоторые доверительные интервалы загрузочного ремешка, полученные из начальной загрузки, в мой график, используя R.
Hello. Я пытаюсь добавить доверительные интервалы к моему графику в R.
У меня было мало данных, поэтому я использовал bootstrap. Мой исходный код выглядит так:
data.df <- read.csv("E:\\data\\data.csv")
library(boot)
matplot(data.df[, 1], data.df[, -1], type="l",
main="Data",xlab="Year", ylab=" ", yaxt="n",
xlim=c(1970, 2015), ylim=c(15, 95),
lwd=c(1, 4, 4, 4),
col=c("black", "red", "burlywood", "cadetblue"),
lty=c("solid", "dotted", "dotted", "dotted" ))
lmodel6=lm(data.df$avg~data.df$X)
abline(lmodel6$coefficients[1],lmodel6$coefficients[2],col=1,lwd=4)
Я использовал matplot, а не просто обычный сюжет. В прошлом я привык к использованию matplot, поэтому решил продолжить их использование.
Далее я сделал бутстрап:
results <- boot(data=data.df, statistic=bs, R=1000, formula=avg~X)
results
bs <- function(formula, data, indices) {
d <- data[indices,]
fit <- lm(formula, data=d)
return(coef(fit))
}
boot.ci(results, type="bca", index=1)
boot.ci(results, type="bca", index=2)
Я надеюсь получить какой-нибудь результат, похожий на:
Но в большинстве примеров добавления доверительных интервалов к графикам онлайн используется обычный график или ggplot. Или они не имеют дело с начальной загрузкой. В других местах, таких как: https://thomasleeper.com/Rcourse/Tutorials/olsresultplots.html, создаются хорошие графики, но они используют другую процедуру и другой код.
Любое руководство или помощь будет принята с благодарностью!