R/nearestneighbourenvelope.R

In bioimaginggroup/bioimagetools: Tools for Microscopy Imaging

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))
}