Bioequivalence and Bioavailability Forum

Trapezoidal rule - Hmisc [🇷 for BE/BA]

posted by d_labes  – Berlin, Germany, 2010-04-19 15:41 (5901 d 11:12 ago) – Posting: # 5168
Views: 12,067

Dear All,

just discovered in Frank Harrels library (package) Hmisc:
Area under curve defined by arbitrary x and y vectors, using trapezoidal rule and assuming x sorted appropriately

trap.rule <- function(x, y) sum(diff(x) * (y[-1] + y[-length(y)]))/2

Wow ! The cream of the crop of programming art!

Complete thread:

• Trapezoidal rule - Hmiscd_labes 2010-04-19 13:41 [🇷 for BE/BA] 
  • Trapezoidal rule - PK Helmut 2010-04-20 12:53
    • package PK d_labes 2010-04-20 14:35
      • package PK Jack 2010-04-20 15:21
      • Syntax... Helmut 2010-04-20 15:31
      • package PK Jack 2010-04-23 17:00
        • package PK d_labes 2010-04-26 08:04