R/IkWR.r

In damarals/TeachingSampling: Selection of Samples and Parameter Estimation in Finite Population

#' @export

IkWR <- function(N, m)
{
Q <- SupportWR(N, m, ID = FALSE)
I <- matrix(0, choose(N+m-1, m), N)
for (i in 1:m) {
for (j in 1:choose(N+m-1, m)) {
for (k in 1:N) {
if (Q[j, i] == k)
I[j, k] <- 1
}
}
}
I
}