Nothing
nearest.neighbour.envelope<-function(X,Y,Z,X2,Y2,Z2,psz=25,n=100,range=.5)
{
do<-function(i,X,Y,Z,X2,Y2,Z2,psz,range)
{
X2<-X2+runif(1,-range,range)
Y2<-Y2+runif(1,-range,range)
Z2<-Z2+runif(1,-range,range)
redgreen<-nearest.neighbour.distribution(X,Y,Z,X2,Y2,Z2,same=TRUE,psz=psz)
hist<-hist(redgreen,n=100,plot=FALSE)
return(list("x"=hist$mids,"y"=cumsum(hist$density)/100))
}
runs<-lapply(1:n,do,X,Y,Z,X2,Y2,Z2,psz=psz,range=range)
#run.max<-apply(runs,1,max)
#run.min<-apply(runs,1,min)
true<-nearest.neighbour.distribution(X,Y,Z,X2,Y2,Z2,same=TRUE,psz=psz)
truehist<-hist(true,n=100,plot=FALSE)
true<-cumsum(hist$density)/100
return(list("true"=true,"runs"=runs,"x"=truehist$mids))
}
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.