[ DATT-A9D > データで遊ぶ > 新逢坂山田前 ]
PCゲーム「A列車で行こう9」で遊びながら学ぶ、Rのハローワールド。(最終更新:2025年2月19日)
パッケージ「nnclust」をソースからコンパイルしてインストール(管理者として実行) | # ダウンロード url <- "https://cran.r-project.org/src/contrib/Archive/nnclust/nnclust_2.2.tar.gz" pkgFile <- "nnclust_2.2.tar.gz" download.file(url = url, destfile = pkgFile) # インストール install.packages(pkgs=pkgFile, type="source", repos=NULL) |
---|
「適応的新逢坂山田前」のRスクリプト(「ベクトル円海山」を実行したあと続けて実行してください) |
---|
myakari3dc <- subset(mynewakari3d,V3 < 30 & V3 > -20) # -10m~20mの地点を使う plot(myakari3dc) # 所望のデータになったか確認 library(nnclust) # 起動するたびに必要 myakari2d0 <- myakari3d0[ ,c(1, 2)] mycmpakari0elv <- prcomp(myakari3dc, scale=TRUE) mykmakari90elvc3 <- kmeans(mycmpakari0elv$x, as.integer(nrow(myakari3dc) / 41625 * 90), nstart=50) original <- mykmakari90elvc3$centers %*% t(mycmpakari0elv$rotation) original <- scale(original, center = FALSE, scale = 1 / mycmpakari0elv$scale) original <- scale(original, center = -mycmpakari0elv$center, scale = FALSE) mykmakari90elvc3o <- as.table(original) mykmakari90to30 <- kmeans(original, as.integer(nrow(myakari3dc) / 41625 * 30), nstart=50) #write.table(mykmakari90elvc3o, file="clipboard", sep="\t") #write.table(mykmakari90to30$centers, file="clipboard", sep="\t") # いくつ? nrow(mykmakari90elvc3o) nrow(mykmakari90to30$centers) |
鉄道駅(およそ30地点)の最小全域木(MST)を描く |
par(cex=0.8) # 文字を小さく Y <- matrix( mykmakari90to30$centers, ncol=3 ) mst <- mst(Y) plot(myakari2d0, type="n") segments( Y[mst$from,1], Y[mst$from,2], Y[mst$to,1], Y[mst$to,2], col="green") text(mykmakari90to30$centers) rect(0, 0, 256, -256, border = "red") # プロットされた図を右クリックしてコピーや保存をします |
路線バス(およそ90地点)の最小全域木(MST)を描く |
par(cex=0.6) # 文字を小さく X <- matrix( mykmakari90elvc3o, ncol=3 ) mst <- mst(X) plot(myakari2d0, type="n") text(mykmakari90elvc3o) segments( X[mst$from,1], X[mst$from,2], X[mst$to,1], X[mst$to,2], col="black") rect(0, 0, 256, -256, border = "red") # プロットされた図を右クリックしてコピーや保存をします |
※諸般の事情に鑑み、マップによらず変数名を「akari」にしています。Windowsでの実行を念頭にしていますので入出力に「クリップボード」を使用しています。
※プログラミングの学習との混同がないよう、for文などの使用を避けています。1つずつ確かめながら対話的に実行してください。
※クラスターの重心の座標をコピーして使いたいときは「write.table」の行のコメントアウトを外して実行してください。
水面に揺れる工場の灯 (k=81・k=27) | ![]() |
---|---|
文明開化の薫る街 (k=82・k=27) | ![]() |
宵明けの大都 (k=86・k=28) | ![]() |
空港連絡線は続くよ (k=91・k=30) | ![]() |
古都再興 (k=103・k=34) | ![]() |
夕日町計画 (k=105・k=35) | ![]() |
広域都市計画 (k=106・k=35) | ![]() |
おどる広告都市 (k=110・k=36) | ![]() |
城下町の復興 (k=110・k=36) | ![]() |
大都市構想 (k=121・k=40) | ![]() |
統計数理研究所のミラーサイトから「R」を入手しよう 統計数理研究所とは |
ARX小粒あずきケット(41439) All Rights Reserved. ©2018-2025, tht.