File: Rcode-3-4.R

package info (click to toggle)

urca 1.3-3-1

links: PTS, VCS
area: main
in suites: bookworm
size: 2,432 kB
sloc: fortran: 501; ansic: 15; makefile: 2

file content (9 lines) | stat: -rwxr-xr-x 285 bytes

parent folder | download | duplicates (7)

library(fracdiff)
set.seed(123456)
y <- fracdiff.sim(n=1000, ar=0.0, ma=0.0, d=0.3)
y.spec <- spectrum(y$series, plot=FALSE)
lhs <- log(y.spec$spec)
rhs <- log(4*(sin(y.spec$freq/2))^2)
gph.reg <- lm(lhs ~ rhs)
gph.sum <- summary(gph.reg)
sqrt(gph.sum$cov.unscaled*pi/6)[2,2]