フォーラム - neorail.jp R16

【自由研究】の話題

 → 茨城県内で必要な「主論文・野帳・掲示物」とは? 【詳細】(2017/10/19) NEW!
 → 【自由研究】の記事一覧(新着順)


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

研究を楽しく「追体験」! 真っ白のキャンバスに虹色の未来を描く方法、教えます。
R with Excel prcomp | kmeans・hclust | rpart | ggvoronoi | spantree | lm NEW!


2021年5月の話題
更新:2021/11/3

[4547]

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

研究ホワイトボックス(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

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

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

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

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://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級対策講座
 > リテールマーケティング(販売士)検定試験3級対策講座
 > 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

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


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

clusterwherePC3PC6prefcity
春日291.67591.52322238
大井281.30990.82672028
川内180.4088-0.87271416
日向1120.0340.79861316
千代田180.0252-0.55441320
愛宕180.618-0.03191117
鹿島1120.7911-0.51681114
若葉180.642-0.27351016
瑞穂58-1.38330.4906915
熊野510-0.6723-1.7972912
青葉180.5211-0.4237814
千歳512-0.6468-0.6699710
明石181.31650.174679
有明58-1.30790.0975612
石垣512-0.08430.37566
佐渡580.57990.81755
秋津580.08770.006355
秋月590.52541.292755
福江58-0.75060.40155
平戸580.0481-0.123745
長良58-0.3339-0.793745
由良58-0.0718-0.054445
羽黒580.26770.462235


 「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


関連する記事


[4969]

研究ホワイトボックス(49) とにかく「主成分分析つき回帰木」するには tht - 2023/1/1


[3574]

研究ホワイトボックス(28) 多変量データの全域を考察するには〜駅弁の商品構成を例に(中編) tht - 2017/12/3


[3639]

【自由研究】ふわコレ(6) tht - 2018/4/30


[3669]

【自由研究】ゆるシミュ(9) tht - 2018/8/16


[3868]

まだ見ぬ「ドリームかまち」を探して(インサイト編) tht - 2019/12/30


[3636]

Re:[3583] 或るナポリタンと停電のハフソク(則) tht - 2018/4/10


[4075]

「照焼大橋」不滅なれ(モスクワ編) tht - 2020/5/21


[4235]

きょうは天下茶屋で有楽町層。 tht - 2020/10/1






neorail.jp/は、個人が運営する非営利のウェブサイトです。広告ではありません。 All Rights Reserved. ©1999-2024, tht.