Я в настоящее время самостоятельно изучаю многомерное шкалирование.Среди прочего я изучаю Borg & Groenen (2005): Современное многомерное масштабирование: теория и приложения.
На странице 10 они представляют реальный набор данных, о котором сообщает Wish (1971).Wish (1971) попросил 18 студентов оценить глобальное сходство разных пар наций, таких как Франция и Китай, по 9-балльной шкале от 1 = очень разные до 9 = очень похожие.Поскольку набор данных общедоступен, я хотел повторить результат в R для практических целей.В качестве первого шага я хотел воспроизвести следующую конфигурацию, также представленную в Borg & Groenen (2005, с. 10).
Я действовал следующим образом:
library(smacof) ### this package contains the data set
data(wish) ### that is the data set
Поскольку наборы данных содержат рейтинги сходства, я применил неметрическое многомерное масштабирование, используякоманда isoMDS пакета MASS.Хотя авторы учебника говорят о «двумерной конфигурации MDS», я также попробовал многомерные решения.Поэтому я закодировал цикл, который выполняет многомерное масштабирование для конфигураций, содержащих измерения от 2 до 9.
X <- c()
for (i in 2:9) {
MDS <- isoMDS(wish, k = i)
X <- c(X, MDS$stress)
plot(MDS$points[,c(1,2)])
text(MDS$points[, 1], MDS$points[, 2], colnames(as.matrix(wish)), cex=.6,
pos = 1)
}
plot(X, type = "b") ### this allowed me to plot the stress levels associated with each configuration
Ни один из полученных графиков не был похож на тот, что представлен в Borg & Groenen (2005, стр. 10).Например, карта для двух измерений выглядит следующим образом:
Я проверил, что набор данных идентичен тому, о котором сообщили Borg & Groenen (2005, стр. 10).Я также попробовал метрическое масштабирование следующим образом:
for (i in 2:9) {
plot(smacofSym(wish, ndim=i))
}
Опять же, я не смог повторить результаты, сообщенные Borg & Groenen (2005, с. 10).Однако я не уверен, что допустил ошибку при попытке повторить результаты.