Hello,
I am fitting data using different methods e.g. Local Polynomial and Smoothing
splines. The data is generated out of a true function model with added normally
distributed noise.
I would like to know "how often the confidence band for all points
simultaneously contain all true values". I can answer the question for one
point in the following way:
e.g.
#
========================================================================================#
How many times the pointwise confidence interval at x=0.5 contains the true
value at 0.5
# i.e. what is the so called "coverage rate"?
#
========================================================================================pos
= which(x==0.5)
sum(abs(estlp[pos,] - m(x[pos])) <= 1.96*selp[pos,]) # equidistant x
outputs 946
# non-equidistant x
outputs 938
sum(abs(estss[pos,] - m(x[pos])) <= 1.96*sess[pos,]) # equidistant x
outputs 895
# non-equidistant x
outputs 936
This basically tells me that out of 1000 simulation runs with different random
noise, 946 times the true value i.e. m(x) for x=0.5 is contained within the 95%
confidence interval. The estlp Local Polynomial performs better than Smoothing
Splines under this criteria ...
Now is there any specific way to answer "how often the confidence band for
all points simultaneously contain all true values" other than this below?
#
========================================================================================#
How often does the confidence band for all points simultaneously contain all
true values?
#
========================================================================================sum(abs(estlp[,]
- m(x[])) <= 1.96*selp[,]) # equidistant x outputs 92560
# non-equidistant x
outputs 92109
sum(abs(estss[,] - m(x[])) <= 1.96*sess[,]) # equidistant x
outputs 90804
# non-equidistant x
outputs 94641
Is there a dedicated function in R for this purpose i.e. to build confidence
bands around a given fit ... maybe a way to plot it nicely too given that the
Estimated SE are calculated.
Many thanks in advance,
Best regards,
Giovanni
[[alternative HTML version deleted]]
