Площадь R под кривой ZOO пакета
Я хотел бы рассчитать AUC для кадра данных, используя пакет ZOO, как описано в предыдущем посте AUC:
library(zoo)
x <- hdt_br$time
y <- hdt_br$OD
id <- order(x)
AUC <- sum(diff(x[id])*rollmean(y[id],2))
Использование вышеуказанного метода, когда данные упорядочены (id) с точки зрения x
Я получаю относительно большой AUC 44197.59
, Тем не менее, при ручном расчете я получаю AUC 11019.2
, Вероятно, из-за сочетания + ive и -ive точек данных и тренда, являющегося неправильным интегралом.
Я пробовал несколько разных решений, но я действительно застрял и был бы признателен за любые идеи для решения этой проблемы.