R: fa.parallel пакета psych, глобальное значение false

У меня возникла проблема при вызове функции fa.parallel() из пакета psych в R.

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

Учитывая следующий пример набора данных:

  SQ100 SQ101 SQ102 SQ103 SQ104 SQ105 SQ106 SQ107 SQ108 SQ109 SQ110 SQ111 SQ112
1     1     2     3     2     1     1     2     2    NA    NA    NA     1     1
2     3     2     2     2     1     1     3     2     3     1     3     1     1
3     1     1     1     1     1     1     1     1     1     1     5     1     1
4    NA     1    NA     1     1     1     2     1    NA    NA    NA     1     1
5     1     5     4     2     3     2     3     2     1    NA    NA     2     2
6    NA    NA    NA    NA    NA     2     1     1     1    NA    NA     1     1

У меня есть категории от 1 до 5 и NA для отсутствующих значений.

Если я выполню

fa.parallel.poly(dataset,global=TRUE)

Вывод консоли говорит:

The items do not have an equal number of response alternatives, global set to FALSE

Как вы можете видеть в таблице выше, существуют переменные, такие как «SQ112», которые имеют только значения 1 и 2, или «SQ109» имеет 1 и NA в качестве ответов.

Однако это только наблюдаемые данные, и варианты ответов от 1 до 5 или NA возможны для всех переменных, но в этом случае просто не наблюдаются.

Итак, я хочу, чтобы для глобального значения было установлено значение TRUE, но каждый раз, когда я вызываю функцию, оно устанавливает глобальное значение false, даже если я указываю значение True в параметре функции.


r package psych factor-analysis
person Jannic Beck    schedule 30.12.2014    source источник

Ответы (1)


arrow_upward
3
arrow_downward

Сообщение указывает на то, что в ваших наблюдаемых данных не все элементы имеют одинаковое количество категорий. Параллельный анализ должен работать в любом случае. Начиная с psych 1.5.1, fa.parallel теперь имеет возможность напрямую выполнять полихорические корреляции, поэтому вам не нужен вызов fa.parallel.poly. Просто скажите fa.parallel(dataset,cor="poly"). (Он по-прежнему будет жаловаться на то, что данные не имеют одинакового количества вариантов ответа.)

person William Revelle    schedule 05.02.2015
comment
Но обратите внимание: однако это всего лишь наблюдаемые данные, и варианты ответов от 1 до 5 или NA возможны для всех переменных, но в этом случае просто не наблюдаются. Да, это работает, но я не хочу, чтобы элементы с менее наблюдаемыми категориями, чем 5, обрабатывались так, как если бы они имели только f.e. 3 категории ответов... - person Jannic Beck; 10.02.2015
comment
Проблема в том, что fa.parallel случайным образом пересэмплирует ваши наблюдаемые данные. Если ваши данные не содержат всех пяти альтернатив, они не могут повторно выбирать альтернативы, которых не существует. Однако вы можете попробовать параметр sim=TRUE, который будет генерировать смоделированные данные с той же размерностью, что и ваши данные, но с другой структурой элемента. - person William Revelle; 11.02.2015