R (по-видимому) производит ошибки в расчетах

Это почти наверняка человеческие ошибки, но я не могу понять, что я делаю не так. Буду признателен за любые мысли.

Воспроизводимый пример:

data <- structure(list(pre_IB1 = c(-3L, -2L, 2L, 2L), pre_IB2 = c(3L, 
-1L, 1L, -3L), pre_IB3 = c(3L, 0L, 0L, -2L), pre_IB4 = c(-3L, 
-1L, -1L, -3L), pre_IB5 = c(3L, 0L, -3L, 2L), pre_IH1 = c(3L, 
1L, 2L, -3L), pre_IH2 = c(3L, 2L, -3L, -3L), pre_IH3 = c(3L, 
1L, 2L, 2L), pre_IH4 = c(3L, 2L, 1L, -3L), foot = c(3L, 3L, 3L, 
3L), switch = c(3L, 3L, 3L, 3L), condition = c("ff", "ff", "ff", 
"ff"), post_IB1 = c(3L, -2L, 0L, 0L), post_IB2 = c(3L, 0L, -1L, 
-1L), post_IB3 = c(3L, 1L, -1L, -1L), post_IB4 = c(3L, -1L, 0L, 
-2L), post_IB5 = c(3L, -3L, -3L, 2L), post_IH1 = c(3L, 1L, 1L, 
-1L), post_IH2 = c(3L, 1L, -3L, -1L), post_IH3 = c(3L, 1L, 2L, 
1L), post_IH4 = c(3L, 2L, 2L, -1L), pol_gen = c(0L, 0L, -1L, 
-1L), age = c(76L, 34L, 66L, 35L), gender = c(1L, 1L, 2L, 2L), 
    trolley = c(2L, 1L, 1L, 1L), participant = c("5c58244f62290e00012fa64c", 
    "5f50db296d9ce726e639d020", "5d48673efef71c0016d8bf40", "61017941304cf7ac38c2e110"
    )), row.names = c(NA, 4L), class = "data.frame")

Я пытаюсь выполнить два расчета, первый простой, второй более сложный. В первом случае я пересчитываю значения ответов, а затем вычисляю средний балл:

data[,1:9] <- apply(data[,1:9], 2, likert_convert,
            top.x = 3, bot.x = -3, top.y = 1, bot.y = -1)
data[,13:21] <- apply(data[,13:21], 2, likert_convert,
            top.x = 3, bot.x = -3, top.y = 1, bot.y = -1)

data <- mutate(data, post_consTotal=(post_IB1 + post_IB2 + post_IB3 + post_IB4 + post_IB5 + post_IH1 + post_IH2 + post_IH3 + post_IH4)/9)

Этот код создает значения post_consTotal, которые неверны. Например, во второй строке должно быть (-.66 + 0 + .33 + -.33 + -1 + .33 + .33 + .33 + .66)/9= -0.1, но R выдает -4.934325e-17.

Интересно, что эта проблема не возникает, если я пропускаю шаг пересчета.

Второй расчет, который я выполняю, включает в себя получение чего-то похожего (но не совсем) на среднее абсолютное отклонение между набором переменных. Сначала я пересчитываю, как указано выше, чтобы переменные находились между 1 и -1:

data[,1:9] <- apply(data[,1:9], 2, likert_convert,
            top.x = 3, bot.x = -3, top.y = 1, bot.y = -1)
data[,13:21] <- apply(data[,13:21], 2, likert_convert,
            top.x = 3, bot.x = -3, top.y = 1, bot.y = -1)

arrayData = array(data = NA, dim = c(5, 5, nrow(data)))
finalMeanScore = c()
for( numberOfParticipants in 1:nrow(data) ){
  tempData = unlist(as.vector(data[ numberOfParticipants, c("pre_IB1", "pre_IB2", "pre_IB3", "pre_IB4", "pre_IB5") ]))
  for ( i in 1:length(tempData)){
    tempCoh = 1 - abs(tempData[i] - tempData)
    arrayData[i,, numberOfParticipants] = tempCoh
  }
  diag(arrayData[,, numberOfParticipants]) = 0
  finalMeanScore[numberOfParticipants] = (sum(arrayData[,, numberOfParticipants]))/20
}
data <- data %>% 
  mutate(pre_cohIB=finalMeanScore)

Учитывая, что все входные значения находятся между 1 и -1, максимальное абсолютное отклонение составляет 2. Таким образом, строка

tempCoh = 1 - abs(tempData[i] - tempData)

не должна давать отрицательных значений, однако такая ситуация происходит довольно часто.

Я предположил, что это может быть связано с точностью представления чисел с плавающей запятой, что может объяснить первую проблему. Однако расхождение между -0.01 и -0.00000000000000004934325 кажется слишком большим, чтобы его можно было объяснить ошибкой округления.

Извините за длинный и насыщенный пост, любые предложения будут весьма полезны.