Коды R: вероятность охвата интервалов прогноза для cdf
Я пытаюсь получить коды R для вероятности охвата 95% интервалов прогнозирования для cdf. Я построил доверительный интервал в R, но я не знаю, как получить интервалы прогнозирования и вероятность покрытия интервалов прогнозирования для cdf. Примечание: я знаю, что есть другой способ получить доверительный интервал для этих моделей, но позже я хотел бы изменить термин ошибки на двухрежимный. Подскажите, пожалуйста, как решить мои проблемы? Заранее спасибо!
Коды R:
set.seed(100)
n=1000
e=rnorm(n,0,1)
t=1:n
#x=ts(t/n)or
x=(t-0.5)/n
y <- 1 - 2*x + e # model
FF <- ecdf(y) #cdf
#now we are going to create the confidence intervals for y
alpha <- 0.05
eps <- sqrt(log(2/alpha)/(2*n))
xx <- sort(y)
ll <- pmax(FF(xx)-eps,0) # lower interval
uu <- pmin(FF(xx)+eps,1) # upper interval
plot(FF,cex=0.5)
lines(xx, ll, col="blue") # lines (also points) doesn't need the add=TRUE command
lines(xx, uu, col="blue")