set.seed(123)
group1 <- rexp(18, rate = 1/10) # Exponential distribution
group2 <- runif(22, min = 10, max = 18) # Uniform distribution
group3 <- rchisq(20, df = 5) + 7 # Right-skewed Chi-square + shift
RWPNPBT(list(group1, group2, group3), B = 1000)
Run the code above in your browser using DataLab