【研究ホワイトボックス】

研究ホワイトボックス（48）複雑なデータと「対話」するには～郵便番号データを例に（回帰木編）

・「R with Excel」（続き）
・「R with Excel」（続き）

（約6000字）

　[4544],[4545],[4546]の続きです。

・[4544]
　> 本日「このフォーラム」、▼[3707],[3724]と▼[3708]、それに▼[4417],[4418],[4421]と▼[4419],[4431]などからの続きです。

　https://blog-imgs-128-origin.fc2.com/o/p/e/opentp2/blog_import_5d5dc9d02741f.jpeg

　やってきました京浜東北線ホーム。（違）

・[3697]
　> 京浜東北線の停車駅をそらんじてみせる
　> …なぜ京浜東北線なのかと問えば答えに詰まるっ。

　https://livedoor.blogimg.jp/tamudouraku/imgs/e/9/e965ab37.jpg
　https://images.tetsudo.com/report/278/11.jpg

　地上時代の京浜東北線。（違）

・[4140]
　> 「京浜東北線の快速電車」は1単語です。これはもう、ゼッタイに1単語です。

・[3990]
　> そういうふうに思っている。

　https://stat.ameba.jp/user_images/20200920/19/hamidashi103/81/1c/j/o1819181914822403532.jpg

　> のびーる
　> ひろがーる
　> たのしいー！

　そこは「たーのしー！」じゃないのか。えー…（てんてんてん）。

■「R with Excel」（続き）

ライブラリを読み込みます	library(rpart) library(rpart.plot)
データフレームを用意します	myindexzipall <- data.frame(cbind(myzipall[,c(2)], mycmpzipall$x))
rpartを実行する	myrpartzipall <- rpart(V1~., myindexzipall , minsplit = 5)
rpartの実行結果を rpart.plotでプロットする	rpart.plot(myrpartzipall , type = 2, uniform = TRUE, under = 1) # プロットされた図を右クリックしてコピーや保存をします

　2列目「city」をcbindしたら「V1」になるけれど、▼そのV1を、▼PC1からPC49までを使って説明する回帰木をつくれー。…お、おぅ。

　https://neorail.jp/forum/uploads/r_myzipall_rpart.png

　件数非公表なので「％」にしてみました。…ほぅ？

データフレームを用意します	myindexzipall8 <- data.frame(cbind(myzipall[,c(2)], mycmpzipall$x[, c(3,4,5,6,7,8)]))
rpartを実行する	myrpartzipall8 <- rpart(V1~., myindexzipall8 , minsplit = 20)
rpartの実行結果を rpart.plotでプロットする	rpart.plot(myrpartzipall8 , type = 2, uniform = TRUE, under = 1) # プロットされた図を右クリックしてコピーや保存をします

　https://neorail.jp/forum/uploads/r_myzipall_rpart_city_by_pc3to8.png

　うーん。PC3からPC8までを使って「city」を説明せよ。…それなんて回帰？（ぐぇ）しつれいしました。この限りにおいてはPC7が、PC1とPC2で外れ値的なものを見分けるのと似たことができるとおっしゃる。PC3を何度もぐちゃぐちゃと細切れに使う。すごーく“ご都合主義”のにおいがする。（※個人の感想です。）

　https://neorail.jp/forum/uploads/r_myzipall_rpart_pref_by_pc3to8.png

　スクリプトは省略しますが「city」ではなく「pref」をPC3からPC8で説明せよという回帰木。えー…（てんてんてん）。

　https://neorail.jp/forum/uploads/r_myzipall_rpart_pref_by_pc3to8_x8.png

　そんなに細かく分けないでいいよと言ったから記念日。（※詠嘆）PC8って、そんなに“優秀な”（＊）主成分だったのかしら。本当かしら。（※ジト目）

＊真っ先にそれを見るとざっくりした分けかたができる、の意。

　『PC3－PC6平面』で見ようかなと思った、その気持ちを説明するなら、たぶんこんなかんじ（※ひらがな）ということが表現できた回帰木ではありませんか。えー…（てんてんてん）。この意味では、PC8は外れ値的なものを分けられればなんでもよかった的なことになってきます。PC3とPC6だけではよく見えない中心部を細かく見たければPC7とPC5を見るといいよ。そんなふうに見えてきませんか。えー…（てんてんてん）。

　最下段で見て左から5番目の「11％」のところの山というかこぶのようなの（「65％」という「その他大勢」の集団がなかった場合の正規分布のピークみたいなの）が、実は「いい感じ？」なのではないか。もっと本当でしょうか。どんどんわからなくなってまいりました。（※遠い目）

　rpartの戻り値で色分けしたプロットはできないのかしら。な・・・なんだってー！！（※白目）

■「R with Excel」（続き）

library(mclust) # 起動するたびに必要
clPairs(mycmpzipall$x[,c(3,6)], cl=myrpartzipall8$where)
# プロットされた図を右クリックしてコピーや保存をします

　https://coub-anubis-a.akamaized.net/coub_storage/coub/simple/cw_timeline_pic/99a09ad0acd/abd5163aad2cfa8ec4579/1411640630_1401866495_image.jpg
　https://upload.wikimedia.org/wikipedia/en/b/bc/Aizawa_Inori.jpg
　https://e7.pngegg.com/pngimages/727/295/png-clipart-inori-aizawa-internet-explorer-web-browser-magical-girl-internet-explorer-cartoon-fictional-character.png
　https://neorail.jp/forum/uploads/r_myzipall_plot_pc3_pc6_rpart_pref_x8.png

　https://www.irishexaminer.com/cms_media/module_img/430/215159_1_googlediscover_bn-612409_71a17f0fb51e426ebc545ea92ae52e26.jpg

　うーん。他社ブラウザに侵食されながら微笑をたたえつつ宇宙のちりになりゆく某マイクロソフト社製ブラウザの図。…なにそれこわい。（違）

　http://www.sega-mj.com/arcade/img/bnr/bnr_symphogear-xv.jpg

　お好み焼き食べて早く帰ろ。…うん。（棒読み）▼多様な個性が響き合う130の専門店「ユアエルム八千代台」については[4543]を参照。

　https://r.gnavi.co.jp/eki/0002322/konamono/rs/

　> みいちゃん

　出ました「みいちゃん」。八千代台駅東口は徒歩4分であります。ぐるなびを見たとお伝えいただければ幸いです。八千代商工会議所パソコン教室すぐ先。（※ソコジャナイ。）

・[3944]
　> スライド1枚で1分ですよー。（棒読み）
　> 成毛さんが言ってるから絶対だ。（違）

　http://www.yachiyocci.jp/htmlyachiyo/PC.html

　> ワード・エクセルはもちろん、日商PC検定試験対策講座や日商簿記検定試験対策講座も充実しています。仕事の効率化を図りたい方にピッタリ。社員研修にもご利用いただけます。

　> 日商ビジネス英語検定3級対策講座
　> 日商ビジネス英語検定2級対策講座
　> リテールマーケティング(販売士)検定試験３級対策講座
　> CAD入門講座
　> 弥生会計講座
　> 弥生給与講座

　「八千代商工会議所パソコン教室」のご案内でした。

プロットする

myclwh <- cbind(mykmzipall5$cluster, myrpartzipall$where)
colnames(myclwh) <- c('cluster','where')
plot(myclwh)

　https://neorail.jp/forum/uploads/r_myzipall_plot_cluster_where.png

　そういうふうに、k-meansよりは柔軟に分けてくれている…かな、というふいんきであります。

	cluster	where	PC3	PC6	pref	city

春日	2	9	1.6759	1.5232	22	38
大井	2	8	1.3099	0.8267	20	28
川内	1	8	0.4088	-0.8727	14	16
日向	1	12	0.034	0.7986	13	16
千代田	1	8	0.0252	-0.5544	13	20
愛宕	1	8	0.618	-0.0319	11	17
鹿島	1	12	0.7911	-0.5168	11	14
若葉	1	8	0.642	-0.2735	10	16
瑞穂	5	8	-1.3833	0.4906	9	15
熊野	5	10	-0.6723	-1.7972	9	12
青葉	1	8	0.5211	-0.4237	8	14
千歳	5	12	-0.6468	-0.6699	7	10
明石	1	8	1.3165	0.1746	7	9
有明	5	8	-1.3079	0.0975	6	12
石垣	5	12	-0.0843	0.375	6	6
佐渡	5	8	0.5799	0.817	5	5
秋津	5	8	0.0877	0.0063	5	5
秋月	5	9	0.5254	1.2927	5	5
福江	5	8	-0.7506	0.401	5	5
平戸	5	8	0.0481	-0.1237	4	5
長良	5	8	-0.3339	-0.7937	4	5
由良	5	8	-0.0718	-0.0544	4	5
羽黒	5	8	0.2677	0.4622	3	5

　「pref」と「city」を整数という名の粒々で数えていたら見分けがつかないとか、整数で数えて1でも違えばぜったい違うとか厳密すぎていたところ、都道府県に重みをつけて、小数で扱えてくるといった趣き。（※体言止め）

　https://compe.japandesign.ne.jp/wp/wp-content/uploads/2019/04/02-9.jpg
　https://pbs.twimg.com/media/DMxSTS3VwAASlpU.png
　http://tabitetu-gate.net/wp-content/uploads/2019/02/yokh_10.jpg

　数字がいっぱい。

・[4534]
　> 「グッドテイストタマゴサラダ半熟風」をどこで区切って読みますか。

　> 「ファッションデリカフーズ」「グッドテイストタマゴサラダ」の（カタカナの）文字
　> 「ファッションデリカフーズ」「グッドテイストタマゴサラダ」のカタカナで「グランプリ」と描きたくなりませんか。なりますよね。

　https://upload.wikimedia.org/wikipedia/commons/e/e1/The_Tokyo_Stock_Exchange_-_main_room_3.jpg

　> これからはコンピューターの時代だとかいってね。

　JTB時刻表の「グッたいむ」とかね。（※ソコジャナイ。）

・[4440]
　> オリジナルセンターかいか家
　> 何の店かわからない名前だけど、きっと定食屋で服も売ってるんだ。

　https://atari-kamafuna.com/2018/09/10/mrr_gr_delicafoodshamamoto/

　「デリカフーズハマモト」スゲェ！（※言いかた）

　> 「濱本社長」が焼く焼き鳥

　古い言いかたをすれば「町のお肉屋さん」です。精肉店です。だからといって、工場で働く人が多く住む船橋という町で、だれもかれもが生の肉を購入して冷蔵庫で保管して、おもむろに調理して食べるという悠長なことができるわけもない。焼き鳥やコロッケなど、すぐ食べられるかたちで売る商売のほうがメインになっていった、その結果が現在の「デリカフーズハマモト」なのでしょう。遠い時代の音がする。（※表現は演出です。）

　> 濱本社長の釣りの写真が沢山。

　そっちが先で、それなら精肉店も、ということなんでしょうか。船橋は工場の町である前に漁師の町。というか、社長は何代目でいらっしゃる（しばらくお待ちください）ありがとうございました。このほか、船橋にはチキンハンバーグとミートボールで有名な某社もあることを申し添えます。

　https://www.ishiifood.co.jp/corp-gaiyo.php

　なんと佃煮と煮豆からのスタートでした。遠い時代の音がした気がしました。（※表現は演出です。）

　[4548]に続きます。

この記事のURL　https://neorail.jp/forum/4547/

この記事を参照している記事

[4540]	きょうは関ケ原で京成スカイライナー。	2021/5/1
[4548]	いわゆるMDSとPCAの違い【t-SNEのBarnes-Hutアルゴリズムあり】	2021/5/13
[4552]	きょうは八つ橋庵かけはしから東映太秦映画村まで。	2021/5/21
[4560]	【ぼよよん】月刊ログイン「けっきょく究極コンストラクション」（1985年7月）をココロに（談）【わしゃわしゃ】	2021/5/31
[4591]	レアパープルビッグエッグ（仮）	2021/6/1
[4653]	ツールとは何か	2021/11/16

【自由研究】の話題