library(CoSMoS)
## here we target to a process that has the Pareto type II
## marginal distribution with scale parameter 1 and shape parameter 0.3
## (note that all parameters have to be named)
dist <- 'paretoII'
distarg <- list(scale = 1, shape = .3)
x <- actpnts(margdist = dist, margarg = distarg, p0 = 0)
x
## you can see the points by using
ggplot(x,
aes(x = rhox,
y = rhoz)) +
geom_point(colour = 'royalblue4', size = 2.5) +
geom_abline(lty = 5) +
labs(x = bquote(Autocorrelation ~ rho[x]),
y = bquote(Gaussian ~ rho[z])) +
scale_x_continuous(limits = c(0, 1)) +
scale_y_continuous(limits = c(0, 1)) +
theme_classic()
Run the code above in your browser using DataLab