R/p.WR.r

#' @export

p.WR <- function(N, m, pk){
p <- rep(0,N)
I <- nk(N,m)
N <- dim(I)[1]
for(i in 1:N){
ni <- c(I[i,])
p[i] <- dmultinom(ni, prob=pk)
}
p
}
damarals/TeachingSampling documentation built on June 2, 2019, 9:06 p.m.