Matlab, ¿cómo calcular el AUC (área debajo de la curva)?

tengo el archivo data.txt con dos columnas y N filas, algo como esto:Matlab, ¿cómo calcular el AUC (área debajo de la curva)?

0.009943796 0.4667975 
0.009795735 0.46777886 
0.009623984 0.46897832 
0.009564759 0.46941447 
0.009546991 0.4703958 
0.009428543 0.47224948 
0.009375241 0.47475737 
0.009298249 0.4767201 
[...]

Cada par de valores en el archivo corresponde a uno coordenadas del punto (x, y). Si se traza, estos puntos generan una curva. Me gustaría calcular el área bajo la curva (AUC) de esta curva.

Así que cargar los datos:

data = load("data.txt"); 
X = data(:,1); 
Y = data(:,2);

Así, X contiene todas las coordenadas x de los puntos, y Y todas las coordenadas y.

¿Cómo puedo calcular el área bajo la curva (AUC)?

Fuente

2011-12-28 DavideChicco.it

Depende. ¿Es la [regla trapezoidal] (http://en.wikipedia.org/wiki/Trapezium_rule) lo suficientemente buena para ti? –

http://www.mathworks.com/matlabcentral/fileexchange/24597-area-under-a-curve –

@Robert: parece que es el área bajo la curva de una * función * (Matlab tiene un montón de ' funciones quad'xxxx()). OP está buscando la integración numérica de datos. –

La manera más fácil es la función de regla trapezoidal trapz.

Si se sabe que sus datos son fáciles, puede intentar usar la regla de Simpson, pero no hay nada incorporado en MATLAB para integrar datos numéricos a través de la regla de Simpson. (& no estoy seguro de cómo usarlo para x/Y donde X datos no aumenta de manera constante)

Fuente

2011-12-28 19:49:59

Gracias chicos, lo he intentado con trapz(), pero extrañamente siempre me da valores negativos. ¿Por qué esto? Si es un área, debería ser siempre positiva ... ¿Alguna idea? Gracias! –

Si la curva está por debajo de 0, el área realmente disminuirá. Esto es solo integral, recuerda. Para obtener el AUC positivo, es posible que deba cambiar la línea base. Por ejemplo, restar 'min (Y)' de 'Y'.O puede usar 'abs (Y)' para resumir las áreas positivas y negativas. – yuk

Técnicamente, si usa 'trapz (x, y)', el signo del resultado depende del signo de yy del signo del cambio en x. (recuerde: esto es integral de y dx) Entonces, si sus valores y son positivos pero x está disminuyendo, obtendría un número negativo. En realidad, es un poco más complicado que eso: para curvas cerradas, el signo debe ser positivo para rodear en el sentido de las agujas del reloj y negativo para rodearlo en sentido antihorario (ver http://en.wikipedia.org/wiki/Green%27s_theorem#Area_Calculation). –

Puede hacer algo así:

AUC = sum((Y(1:end-1)+Y(2:end))/2.*... 
    (X(2:end)-X(1:end-1)));

Fuente

2012-01-02 13:27:15 Oli

Oli, ¿podría por favor agregar una explicación o un enlace para la ecuación que utilizó? – Tin

sólo tiene que añadir AUC = trapz (X , y) a su programa y obtendrá el área bajo la curva

Fuente

2012-03-01 18:46:55 Simon

Fuente: Link

Un ejemplo en MATLAB para ayudarle a obtener su respuesta ...

x=[3 10 15 20 25 30]; 
y=[27 14.5 9.4 6.7 5.3 4.5]; 
trapz(x,y)

En caso de tener valores negativos en y, se puede hacer como,

y=max(y,0)

Fuente

2013-05-02 11:33:49 Faheem

[~, ~, ~, AUC] = perfcurve (etiquetas, puntuaciones, posclass);

% posclass podría ser 1

http://www.mathworks.com/matlabcentral/newsreader/view_thread/252131

Fuente

2013-07-16 19:06:36 Nucular

Hay algunas opciones a trapz para la persona dispuesta a hacer algo de código por sí mismos. Este enlace muestra la implementación de Simpson's rule, con el código python incluido. También hay un File Exchange en la regla de los simpsons.

Fuente

2014-06-19 05:56:17 patrik

Matlab, ¿cómo calcular el AUC (área debajo de la curva)?

Respuesta

Cuestiones relacionadas