### artificial 2-sample problem
df <- data.frame(y = rnorm(20), x = gl(2, 10))
### Ansari-Bradley test
at <- ansari_test(y ~ x, data = df, distribution = "exact")
### density of the exact distribution of the Ansari-Bradley statistic
dens <- sapply(support(at), dperm, object = at)
### plot density
plot(support(at), dens, type = "s")
### 95\% quantile
qperm(at, 0.95)
### one-sided p-value
pperm(at, statistic(at))
Run the code above in your browser using DataLab