if (FALSE) {
# Definition of an univariate hawkes model
provaHawkes2<-setHawkes()
str(provaHawkes2)
# Simulation
true.par <- list(nu1=0.5, c11=3.5, a11=4.5)
simprv1 <- simulate(object = provaHawkes2, true.parameter = true.par,
sampling = setSampling(Terminal =70, n=7000))
plot(simprv1)
# Computation of intensity
lambda1 <- Intensity.PPR(simprv1, param = true.par)
plot(lambda1)
# qmle
res1 <- qmle(simprv1, method="Nelder-Mead", start = true.par)
summary(res1)
}
Run the code above in your browser using DataLab