Skip to content

[R-es] Area bajo la curva

3 messages · WILMER CONTRERAS SEPULVEDA, Jorge I Velez, Javier Marcuzzi

#
Buenas noches.

Quisiera saber si existe una libreria o función que me calcule el área bajo
la curva de dos vectores de números. Ademas esta debe tener en cuenta si
existen datos faltantes o NaN.

Por ejemplo calcular el área bajo la curva del vector en el eje x:

x<- c(1,2,3,4,5,6,7,8,9,10)

y el vector en el eje y:

y<-c(0,2,15,30,50,NaN,NaN,10,2,1)


Muchas gracias.
#
Buenas noches, Wilmer.

Puedes intentar la aproximación sugerida en
https://stackoverflow.com/questions/4954507/calculate-the-area-under-a-curve


## preparación
x <- c(1,2,3,4,5,6,7,8,9,10)
y <- c(0,2,15,30,50,NaN,NaN,10,2,1)
d <- data.frame(x, y)
d <- d[complete.cases(d),]
id <- order(d$x)

## cálculo
require(zoo)
with(d, sum(diff(x[id])*rollmean(y[id],2)))
## [1] 169.5

Espero sea de utilidad.
Jorge.-


2017-08-31 19:17 GMT-05:00 WILMER CONTRERAS SEPULVEDA <wilmercs en ufps.edu.co>
:

  
  
#
Estimado Wilmer Contreras Sepulveda

Me parece buena la sugerencia de Jorge, hay muchas formas, yo aparte de la sugerencia que le hicieron le agrego otras dos, en cualquier buscador las palabras claves son R AUC, y el artículo https://www.r-bloggers.com/calculating-auc-the-area-under-a-roc-curve/ , lógicamente aparecerán tantas alternativas que puede llegar a confundirse respecto a cuál es conveniente para su caso.

Javier Rubén Marcuzzi

De: Jorge I Velez
Enviado: jueves, 31 de agosto de 2017 22:43
Para: WILMER CONTRERAS SEPULVEDA
CC: R-help-es
Asunto: Re: [R-es] Area bajo la curva

Buenas noches, Wilmer.

Puedes intentar la aproximación sugerida en
https://stackoverflow.com/questions/4954507/calculate-the-area-under-a-curve


## preparación
x <- c(1,2,3,4,5,6,7,8,9,10)
y <- c(0,2,15,30,50,NaN,NaN,10,2,1)
d <- data.frame(x, y)
d <- d[complete.cases(d),]
id <- order(d$x)

## cálculo
require(zoo)
with(d, sum(diff(x[id])*rollmean(y[id],2)))
## [1] 169.5

Espero sea de utilidad.
Jorge.-


2017-08-31 19:17 GMT-05:00 WILMER CONTRERAS SEPULVEDA <wilmercs en ufps.edu.co>
:
_______________________________________________
R-help-es mailing list
R-help-es en r-project.org
https://stat.ethz.ch/mailman/listinfo/r-help-es