Calculate the Area under a Curve

by

The AUC is approximated pretty easily by looking at a lot of trapezium figures, each time bound between x_i, x_{i+1}, y{i+1} and y_i. Using the rollmean of the zoo package, you can do:

library(zoo)

x <- 1:10
y <- 3*x+25
id <- order(x)

AUC <- sum(diff(x[id])*rollmean(y[id],2))

Make sure you order the x values, or your outcome won’t make sense. If you have negative values somewhere along the y axis, you’d have to figure out how exactly you want to define the area under the curve, and adjust accordingly (e.g. using abs() )

Regarding your follow-up : if you don’t have a formal function, how would you plot it? So if you only have values, the only thing you can approximate is a definite integral. Even if you have the function in R, you can only calculate definite integrals using integrate(). Plotting the formal function is only possible if you can also define it.

More Related Contents:

  1. double integral in R
  2. How to replace a value by another within a variable in R?
  3. Relative frequencies / proportions with dplyr
  4. How to add texture to fill colors in ggplot2
  5. Cleaning up factor levels (collapsing multiple levels/labels)
  6. Split a large dataframe into a list of data frames based on common value in column
  7. Save a plot in an object
  8. How to cbind or rbind different lengths vectors without repeating the elements of the shorter vectors?
  9. lapply vs for loop – Performance R
  10. ggplot, facet, piechart: placing text in the middle of pie chart slices
  11. YAML current date in rmarkdown
  12. Extract regression coefficient values
  13. How to plot a function curve in R
  14. How to make execution pause, sleep, wait for X seconds in R?
  15. Concatenate strings and expressions in a plot’s title
  16. Counting number of instances of a condition per row R [duplicate]
  17. Create a ranking variable with dplyr?
  18. Find value corresponding to maximum in other column [duplicate]
  19. Split delimited single value character vector
  20. R split numeric vector at position
  21. Creating a Plot Window of a Particular Size
  22. R: invalid multibyte string [duplicate]
  23. Reordering factor gives different results, depending on which packages are loaded
  24. Join two data frames in R based on closest timestamp
  25. Drawing a barchart to compare two sets of data using ggplot2 package?
  26. Applying a function to each row of a data.table
  27. R: Calculate cosine distance from a term-document matrix with tm and proxy
  28. Programming-safe version of subset – to evaluate its condition while called from another function
  29. ggplot2 – using two different color scales for overlayed plots
  30. Extract digit from numeric in r

Leave a Comment