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