Nothing
#' @export
S.SI<-function(N,n,e=runif(N))
{
c<-matrix(0,N,1)
dec<-matrix(0,N,1)
sam<-matrix(0,N,1)
for(k in 1:N){
c[k]<-(n-dec[k])/(N-k+1)
if(e[k]<c[k]){
dec[k:N]<-dec[k]+1
sam[k]<-k}
}
sam
}
Any scripts or data that you put into this service are public.
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.