R/data-raw/breakfastdissimilarities.R

data <- read.csv(system.file("extdata", "Breakfast_Dissimilarities.csv", package = "flipExampleData"))
breakfasts = names(data[,-1])
breakfastDissimilarities = vector("list", 4)
n = 4
k = 15
firstRow = seq(1,nrow(data),k)
names(breakfastDissimilarities) = data[firstRow,1]
for (i in 1:4) {
  zmat = as.matrix(data[firstRow[i]:(firstRow[i] + k - 1),-1])
  dimnames(zmat) = list(breakfasts,breakfasts)
  breakfastDissimilarities[[i]] = as.dist(zmat)
}
breakfast.dissimilarities <- breakfastDissimilarities
usethis::use_data(breakfast.dissimilarities, internal = FALSE, overwrite = TRUE)

# MDS - square
library(smacof)
data(breakfast.dissimilarities)
mdsInterval <- smacofSym(breakfast.dissimilarities[[4]],
     type = "interval", eps = 1e-12, itmax = 100000)

flipPlots::LabeledScatterPlot(mdsInterval)
NumbersInternational/flipExampleData documentation built on May 8, 2024, 5:15 a.m.