par(mfrow=c(2,2))
#=====> Example 1 <=====#
F1 <- seq(0.001, 0.999, 0.001)[-sample(999, 300)]
x1 <- quantile(rt(10000, 15), F1)
F2 <- seq(0.001, 0.999, 0.001)[-sample(999, 300)]
x2 <- qnorm(F2)
hold <- cdfDist(x1, F1, x2, F2)
plot(hold)
summary(hold)
#=====> Example 2 <=====#
F1 <- seq(0.001, 0.999, 0.001)[-sample(999, 300)]
x1 <- exp(quantile(rnorm(10000, 1, sd=1), F1))
F2 <- seq(0.001, 0.999, 0.001)[-sample(999, 300)]
x2 <- qchisq(F2, mean(x1))
hold <- cdfDist(x1, F1, x2, F2)
plot(hold)
summary(hold)
#=====> Example 3 <=====#
F1 <- seq(0.001, 0.999, 0.001)[-sample(999, 300)]
x1 <- exp(quantile(rnorm(10000, 0.5, sd=0.5), F1))
F2 <- seq(0.001, 0.999, 0.001)[-sample(999, 300)]
x2 <- qchisq(F2, mean(x1))
hold <- cdfDist(x1, F1, x2, F2)
plot(hold)
summary(hold)
#=====> Example 4 <=====#
F1 <- seq(0.001, 0.999, 0.001)[-sample(999, 300)]
x1 <- exp(quantile(rnorm(10000, 0.5, sd=0.5), F1))
F2 <- seq(0.001, 0.999, 0.001)[-sample(999, 300)]
x2 <- qchisq(F2, mean(x1)+1)
hold <- cdfDist(x1, F1, x2, F2)
plot(hold)
summary(hold)
Run the code above in your browser using DataLab