Внедрить критерий MAP Величера

Question

Я смотрю на пакет psych и учебник VSS . Могу ли я просто заменить VSS на MAP? Как это:

MAP(x, n = 8, rotate = "varimax", diagonal = FALSE, fm = "pa", n.obs=NULL,plot=TRUE,title="Very Simple Structure",...)

или есть другой способ сделать это?

Я сейчас занимаюсь факторным анализом и использую метод локтя на графике осыпей. Я пытаюсь проверить, могу ли я использовать критерий MAP Величера

Answer 1

@ Райан Росарио

Я наконец понял:

install.packages("psych")
library("psych")

Я выполнил шаги образца

my.VSS <- VSS(test.data,title="VSS of 24 mental tests")

и

VSS(sim.circ(nvar=24),fm="mle", title="VSS of 24 circumplex variables")

и

VSS(sim.item(nvar=24),fm="mle", title="VSS of 24 circumplex variables")

вы получите что-то вроде этого в качестве вывода (для последнего ввода):

Очень простая структура VSS 24 циркумплексных переменных Вызов: VSS (x = sim.item (nvar = 24), fm = "mle", title = "VSS из 24 переменных циркумплекса") VSS сложность 1 достигает максимума 0,84 с 3 факторами VSS сложность 2 достигает максимума 0,87 с 8 коэффициентами

Критерий Velicer MAP достигает как минимум 0,05 с 2 факторами

Velicer MAP 1 0,05 0,01 0,01 0,01 0,01 0,02 0,02 0,02

Очень простая структура сложности 1 1 0,44 0,84 0,84 0,80 0,75 0,76 0,80 0,80

Очень простая структура сложности 2 1 0,00 0,85 0,85 0,85 0,86 0,86 0,86 0,87

как сказано в документации (выделение выше - мое):

"Критерий MAP Уэйна Великера имеет был добавлен в качестве дополнительного теста для оптимальное количество компонентов для извлечь. Обратите внимание, что VSS и MAP будут не всегда согласны относительно оптимального номер. "

В этом случае VSS со сложностью 1 и 2 дает ответ 3 и 8 факторов соответственно, в то время как критерий MAP Велицера дает 2.

Answer 2

Поздний ответ, но я подумал, что последую за комментарием подателя страницы. Допустим, вы создаете объект с помощью функции vss:

my.vss <- vss(test.data)

Функция сводки будет предоставлять результаты критерия VSS и MAP, например,

summary(my.vss)

но вы также можете легко извлечь результаты MAP из объекта (как мне было нужно при выполнении критериальных тестов VSS и MAP для многих наборов данных), например:

#returns the number of factors recommended by MAP
which(my.vss$map == min(my.vss$map))   

#returns the number of factors recommended by VSS for complexity 1
which(my.vss$cfit.1 == max(my.vss$cfit.1)

#returns the number of factors recommended by VSS for complexity 2
which(my.vss$cfit.2 == max(my.vss$cfit.2))

В объекте vss.stats data.frame с объектом vss есть целая куча забавных статистических данных, т. Е.

class(my.vss$vss.stats)

Answer 3

Ознакомьтесь с документацией для этой функции, используя ?MAP.

Если там ничего нет, возможно, отправьте электронное письмо автору, указанному на странице CRAN пакета: CRAN Page для Psych

Этот пакет сделан людьми на этом сайте. У них, вероятно, есть список рассылки: PersonalityProject.org