数据是R语言的数据包carData中的Status数据集。
install.packages("carData")
library(carData)
用的函数sammon()则是MASS函数包里的。没有的就要自己导入。
install.packages("MASS")
library(MASS)
直接上代码:
d2=as.matrix(States)
D2=dist(States)
MDS2=sammon(D2);MDS2
x=MDS2$points[,1];y=MDS2$points[,2]
plot(x,y);abline(h=0,v=0,lty=3)
text(x,y,adj=0.5,labels=rownames(d2))
思路参考:
《多元统计分析及R语言建模》王斌会编著 p326-p327,兼容参考p324