R/UPMEsfromq.R

"UPMEsfromq" <-
function(q)
{
n=ncol(q)
N=nrow(q)
s=rep(0,times=N) 
for(k in 1:N) 
if(n!=0) if(runif(1)<q[k,n]) {s[k]=1;n=n-1} 
s 
}

Try the sampling package in your browser

Any scripts or data that you put into this service are public.

sampling documentation built on Nov. 2, 2023, 6:26 p.m.