## symmetric smacof
data <- na.omit(PVQ40[,1:5])
diss <- dist(t(data)) ## Euclidean distances
fit <- mds(diss)
res.jk <- jackmds(fit)
plot(res.jk, col.p = "black", col.l = "gray")
plot(res.jk, hclpar = list(c = 80, l = 40))
plot(res.jk, hclpar = list(c = 80, l = 40), plot.lines = FALSE)
Run the code above in your browser using DataLab