程序如下:
- library(kohonen)
- cities<-read.table("clipboard",sep="\t",header=T); cities
- names(cities)=c("city","SO2","NO2","PM10","CO","O3","PM2.5","day")
- city=cities[,-1]
- rownames(city)=cities[,1];city
- city=as.matrix(cities)
- city
- SOM=som(city,grid=somgrid(xdim=1,ydim=9,topo="rectangular"))
- summary(SOM)
- table(SOM$unit.classif)
- par(mfrow=c(2,2))
- plot(city,main="模拟数据观测点的分布",xlab="x1",ylab="x2")
- plot(SOM,type="mapping",main="SOM网络输出层示意图",pch=SOM$unit.classif)
- plot(SOM,type="changes",main="SOM网络聚类评价图")
- plot(SOM$data,pch=SOM$unit.classif,main="模拟数据SOM网络聚类结果",xlab="x1",ylab="x2")#运行不出来
- points(SOM$code,col=2,pch=10)