Я получаю следующую ошибку при попытке вычислить размер эффекта с помощью функции cohen.d в psych. См. примеры данных ниже.
library(psych)
x <- c(1, 2, 3, 4, 5, 6, 7, 8)
y <- c(1, 1, 1, 1, 2, 2, 2, 2)
xy <- data.frame(x,y)
names(xy) <- c("x", "group")
CohensD <- cohen.d(xy, "group", alpha=.05)
Ошибка в .rowNamesDF<-(x, значение = значение): недопустимая длина 'row.names'
Вероятно, это ошибка, так как в 1.8.12 он работает нормально, а также выдает ту же ошибку с
cohen.d(sat.act[1:8, c('education','gender')], "gender")
но не с
cohen.d(sat.act[1:8, c('education','gender','age')], "gender")
Итак, должно нормально работать с фреймом данных, имеющим более двух столбцов.
x <- c(1, 2, 3, 4, 5, 6, 7, 8)
y <- c(1, 1, 1, 1, 2, 2, 2, 2)
xy <- data.frame(x=x, group=y, z=rnorm(8))
cohen.d(xy, "group")
Call: cohen.d(x = xy, group = "group")
Cohen d statistic of difference between two means
lower effect upper
x 0.53 3.58 6.60
z -0.71 0.88 2.37
Multivariate (Mahalanobis) distance between groups
[1] 3.6
r equivalent of difference between two means
x z
0.87 0.40
Это действительно ошибка, которую я допустил при добавлении расстояний Махалабониса в cohen.d. Исправление доступно в psych_1.9.12.21, который доступен на моем сервере. Вы можете установить его оттуда:
install.packages("psych",repos="http://personality-project.org/r",type="source")
Я опубликую эту исправленную версию в CRAN, как только она появится на следующей неделе. Тем временем, если вам не нужна исправленная версия, просто используйте набор данных с более чем двумя столбцами (как указал А. Сулиман).
Вы должны использовать пакет effsize вместо psych.
install.packages(effsize)
библиотека (эффсайз)
Затем используйте коэн.д effsize.
x <- c(1, 2, 3, 4, 5, 6, 7, 8)
y <- c(1, 1, 1, 1, 2, 2, 2, 2)
коэн.d(x,y)
