## Not run:
#################
# A realisation of spatio-temporal homogeneous Poisson point processes
hpp <- rpp(lambda = 100, replace = FALSE)$xyt
# R plot
plot(hpp)
# This function provides an kernel estimator of the spatial mark variogram function
out <- gsp(hpp)
# R plot - Spatial mark variogram function
par(mfrow=c(1,1))
xl <- c(0,0.25)
yl <- c(0,max(out$gsptheo,out$egsp))
plot(out$ds,out$egsp,type="l",xlab="r = distance",ylab=expression(gamma[sp](r)),
xlim=xl,ylim=yl,col=1,cex.lab=1.5,cex.axis=1.5)
lines(out$ds,rep(out$gsptheo,length(out$ds)),col=11)
## End(Not run)
Run the code above in your browser using DataLab