data(paracou16)
# Keep only 50% of points to run this example
X <- rthin(paracou16, 0.5)
plot(X)
# Calculate Kmm
r <- seq(0, 30, 2)
ActualValues.X <- Kmm.r(X, r)
# Plot
PlotResults(r, DivideByPiR2(ActualValues.X, r), ylab="Kmm / (pi R^2)", ReferenceValue=1)
Run the code above in your browser using DataLab