Joe Cocker
2012-Nov-05 16:33 UTC
[R] Customly low standard deviation in fracdiff.var function
Hi,I have a question about the fracdiff.var function (package fracdiff) which
goal is to recompute more precise confidence intervals for the parameters
estimated by fracdiff (or arfima). More precisely, it deals with the standard
error of the "d" coefficient :
Is it normal that the standard error of the "d" coefficient can be
brought customly close to zero by decreasing the "h" argument ?
The example of the help file shows this fact :
{Quote}
## Generate a fractionally-differenced ARIMA(1,d,1) model :ts.test <-
fracdiff.sim(10000, ar = .2, ma = .4, d = .3)## estimate the parameters in an
ARIMA(1,d,1) model for the simulated seriesfd.out <- fracdiff(ts.test$ser,
nar= 1, nma = 1)
## Modify the covariance estimate by changing the finite-difference
interval(fd.o2 <- fracdiff.var(ts.test$series, fd.out, h = .0001))## looks
identical as print(fd.out),## however these (e.g.) differ
:vcov(fd.out)vcov(fd.o2)
{end of quote}
Looking at vcov(fd.o2) gives that sd's are indeed lower than in
vcov(fd.out). However, by testing :
vcov(fracdiff.var(ts.test$series, fd.o2, h = 0.00000000001))or
evenvcov(fracdiff.var(ts.test$series, fd.o2, h = 1e-90))
the sd of "d" doesn't stop decreasing, while the others have
reached their "computational" limit.
Is that normal ? If not, which value of sd(d) can I use ?
Hopping I have been clear enough,Many thanks !
[[alternative HTML version deleted]]
