He medido dos parámetros (carbono orgánico disuelto DOC = y, y descarga = x). Cuando estas dos variables se trazan una contra la otra, obtenemos un ciclo de histéresis (vea el ejemplo de código y la imagen).
Ahora, para un análisis más detallado, me gustaría determinar el área de este circuito histérico. Descubrí que esto se puede hacer usando el método de lanzamiento de Monte Carlo. Este método dice que el área de un área desconocida es proporcional al área de un rectángulo conocido multiplicado por los golpes en el campo interior (el bucle).
Mi problema ahora es cómo resolver el problema interno / externo con R. ¿Cómo puedo dibujar un rectángulo con un área conocida y cómo puedo sobresalir los golpes aleatorios dentro y fuera del circuito histérico?
Tenga en cuenta que estoy abierto a cualquier otro método ...
Busqué en Google y busqué en varios sitios estadísticos, pero no pude encontrar una respuesta. Cualquier ayuda directa o enlace a otros sitios web / publicaciones es muy apreciada.
Data <- read.table("http://dl.dropbox.com/u/2108381/DOC_Q_hystersis.txt", sep = ";",
header = T)
head(Data)
plot(Data$Q, Data$DOC, type = "o", xlab = "Discharge (m3 s-1)", ylab = "DOC (mg C l-1)",
main = "Hystersis loop of the C/Q relationship")