data(paracou16)
plot(paracou16)
# Calculate D
r <- 0:30
ActualValues.X <- D.r(paracou16, r, "V. Americana", "Q. Rosea", Intertype = TRUE)
# Plot
PlotResults(r, DivideByPiR2(ActualValues.X, r), , ylab="D / (pi R^2)", ReferenceValue=0)
Run the code above in your browser using DataLab