Как рассчитать индекс рассеяния в статистике R

Question

Я пытаюсь выяснить, как рассчитать индекс рассеяния на графике рассеяния в R. Хотя у меня есть формула, я все еще очень плохо знаком с использованием R и кодирования и не уверен, как написать его на языке программирования. , Формула:

Мои данные:

structure(list(S = c(0.619, 0.582, 0.653, 0.611, 0.649, 0.558 ), 
               O = c(0.58, 0.55, 0.52, 0.56, 0.49, 0.54)), 
          row.names = c(NA, 6L), class = "data.frame")

Любая помощь будет принята с благодарностью.

Answer 1

Вот удар по вашей формуле:

ScatterIndex <- function(s, o) sqrt( sum((s-mean(s)) - (o-mean(o)))^2 / sum(o^2) )

Указанные векторы должны быть одинаковой длины.

---

Ваши данные и вывод:

x <- structure(list(S = c(0.619, 0.582, 0.653, 0.611, 0.649, 0.558 ), O = c(0.58, 0.55, 0.52, 0.56, 0.49, 0.54)), row.names = c(NA, 6L), class = "data.frame")
ScatterIndex(x$S, x$O)
# [1] 1.676298e-16