Рассчитать достоверность факторного анализа

Question

Краткое содержание моего вопроса: В настоящее время я делаю факторный анализ с 23 пунктами. Процедура предложена Энди Филдом в Обнаружение статистики с использованием R. Однако у меня есть проблема с вычислением достоверности факторного анализа.

Подробнее:

performance_problem <- .dataset2[, c(17,7,2,11,21)]
A <- dataset2 [, c(8,15,13,21)]
B <- dataset2 [, c(12,6,9,1)]
C <- dataset2[,c(20,14,4)]
D <- dataset2[,c(10,16,3)]
E <- dataset2 [,c(23,19,5)]

Я определил разные предметы под разными факторами. Например, первый фактор назвал проблемы с производительностью, включая 17-й, 7-й, 2-й, 11-й, 21-й элемент в наборе данных2. И я не получил ошибки при запуске кода R.

В моем факторном анализе нет отрицательных факторных нагрузок. Таким образом, я непосредственно рассчитал свою надежность с этим кодом R:

alpha(performance_problem)

И я тоже попробовал это:

psych:alpha(performance_problem)

Но ни один из них не работает. И сообщение об ошибке, которое я получил от R, было:

Error in grDevices::col2rgb(colour, TRUE) : 
  invalid color name 'c(2, 1, 3, 3, 1, 4, 4, 1, 4, 1, 1, 1, 5, 5, 5, 6, 3, 4, 2, 3, 1, 2, 4, 5, 4, 2, 4, 1, 5, 4, 1, 6, 3, 6, 5, 1, 2, 1, 1, 5, 6, 2, 2, 6, 3, 5, 4, 5, 4, 6, 1, 6, 1, 5, 6, 6, 1, 4, 1, 6, 1, 5, 2, 1, 2, 3, 4, 5, 1, 1, 5, 4, 1, 1, 5, 5, 2, 1, 3, 1, 1, 5, 6, 6, 3, 1, 6, 1, 3, 1, 6, 4, 3, 3, 1, 1, 5, 3, 5, 4, 1, 1, 2, 1, 5, 1, 3, 2, 3, 1, 2, 1, 3, 3, 6, 2, 6, 1, 6, 2, 2, 6, 1, 3, 6, 5, 1, 1, 6, 6, 3, 5, 4, 1, 4, 5, 4, 5, 1, 5, 2, 3, 2, 5, 2, 1, 5, 6, 1, 6, 2, 5, 2, 1, 1, 3, 6, 1, 4, 4, 1, 2, 4, 6, 1, 1, 1, 
6, 1, 5, 6, 1, 5, 3, 1, 3, 1, 4, 3, 5, 1, 1, 2, 2, 2, 1, 5, 5, 5, 6, 2, 3, 4, 6, 1, 1, 1, 5, 4, 2, 5, 6, 1, 1, 5, 5, 5, 1, 5, 4, 4, 1, 3, 2, 2, 6, 1, 6, 1, 4, 5, 5, 1, 1, 5, 4, 1, 1, 5, 1, 5, 1, 2, 5, 2, 1, 1, 1, 6, 2, 4, 5, 6, 6, 1, 1, 4, 4, 5, 6, 1, 4, 6, 3, 6, 2, 2, 1, 2, 1, 2, 3, 6, 1, 4, 6, 5, 1, 2, 6, 2, 5, 1, 5, 4, 5, 1, 3, 1, 4, 3, 3, 1, 5, 2, 4, 1, 2, 2, 4, 2, 5, 6, 6, 6, 5, 1, 2, 4, 1, 5, 3, 1, 3, 1, 2, 6, 4, 4, 1, 1, 4, 1, 2, 1, 3, 1, 2, 5, 1, 5, 1, 1

Но я не очень понимаю, почему у меня эта проблема в R? Любая помощь будет очень благодарна!

Answer 1

Попробуйте другой метод:

library(dplyr)
library(psych)

Создайте df только с переменными, которые необходимы для анализа надежности (называемый df1);df2 - это ваш df со всеми начальными переменными;и переменные, которые вы хотите выбрать, например, V1: V10

df1 <- dplyr::select(df2, V1:V10) 

psych::omega(df1, nfactors = 1) # you get both alpha and omega