powered by
bgpd
bgpdSetSeed(x)
gpd
method = "simulate"
# data <- rnorm(1000) # mod <- gpd(data, qu=.7, method="simulate") # bgpdSetSeed(mod) # mod1 <- gpd(data, qu=.7, method="simulate") # this produces the same MCMC output as mod
Run the code above in your browser using DataLab