TokyoR#107_R_GeoData - Speaker Deck

Slide 1

Slide 1 text

#107 @kilometer00 2023.03.15 BeginneR Session R & Geostatistics

Slide 2

Slide 2 text

Who！？ Who?

Slide 3

Slide 3 text

Who！？名前：三村 @kilometer 職業：特任教員 (こうがくはくし) 専⾨：⾏動神経科学(霊⻑類) 脳イメージング医療システム⼯学 R歴： ~ 10年ぐらい近況: 3 & 0歳児, あたらしい職場

Slide 4

Slide 4 text

宣伝!!（書籍の翻訳に参加しました。）絶賛販売中！増刷出来!!

Slide 5

Slide 5 text

→ スライド資料は全てSpeakerDeckで公開しています。 Tokyo.Rでいっぱい話してきました。

Slide 6

Slide 6 text

BeginneR Session

Slide 7

Slide 7 text

BeginneR

Slide 8

Slide 8 text

Before After BeginneR Session BeginneR BeginneR

Slide 9

Slide 9 text

BeginneR Advanced Hoxo_m If I have seen further it is by standing on the shoulders of Giants. -- Sir Isaac Newton, 1676

Slide 10

Slide 10 text

@kilometer00 BeginneR Session R & Geostatistics

Slide 11

Slide 11 text

□ Rの基本 □ Geostatisticsの基本

Slide 12

Slide 12 text

No content

Slide 13

Slide 13 text

IEEE Spectrum’s Top Programming Languages 2022 https://spectrum.ieee.org/top-programming-languages-2022

Slide 14

Slide 14 text

https://spectrum.ieee.org/top-programming-languages-methods

Slide 15

Slide 15 text

https://spectrum.ieee.org/top-programming-languages-methods

Slide 16

Slide 16 text

No content

Slide 17

Slide 17 text

TIOBE Index for July 2023 https://www.tiobe.com/tiobe-index/

Slide 18

Slide 18 text

https://www.tiobe.com/tiobe-index/

Slide 19

Slide 19 text

https://twitter.com/hadleywickham/status/1279405379449913344?s=20

Slide 20

Slide 20 text

https://gist.github.com/daroczig/ Number of R packages ever published on CRAN

Slide 21

Slide 21 text

https://cran.r-project.org/web/packages/ The Comprehensive R Archive Network (CRAN)

Slide 22

Slide 22 text

h"ps://www.datasciencemeta.com/rpackages

Slide 23

Slide 23 text

① 興味を持つ ③ 周りを⾒渡す。 ④ 隣⼈がみつかる。 ② 始めてみる。 ⑤ 隣の芝が⻘く⾒える。 ⑥ 始めてみる。無駄にはならない（賢く⽴ち回ろうとして、外から恐る恐る⾒てるより、⼿を動かして始めてみよう）⼈⽣は有限で⻘春は瞬きかもしれないけど、⼤丈夫。プログラミング⾔語の2つや3つを学ぶ時間ぐらいある。

Slide 24

Slide 24 text

の始め⽅ 1. Rをインストール 2. RStudioをインストール https://www.r-project.org/ https://posit.co/download/rstudio-desktop/

Slide 25

Slide 25 text

の始め⽅ 3. RStudioを⽴ち上げる 4. 新規プロジェクトを作る 6. tidyverseをインストール install.packages(“tidyverse”) 5. 新規スクリプトを開く File -> New Project... ⌘ ⇧ N

Slide 26

Slide 26 text

① ここにRのコードを書く ② 選択して実⾏(⌘↩) ③ 実⾏結果が表⽰スクリプト → 名前をつけて保存コンソール RStudioの操作

Slide 27

Slide 27 text

① ここにRのコードを書く ② 選択して実⾏(⌘↩) ③ 実⾏結果が表⽰スクリプト → 名前をつけて保存コンソール環境変数などが表⽰ファイル/プロット/ヘルプなど⾃動で反映 RStudioの操作

Slide 28

Slide 28 text

の始め⽅おすすめパッケージ達・ tidyverse (データ科学総合ツール) ・ data.table (大きいデータ取り扱うなら) ・ cmdstanr (ベイズ統計やるなら) ・ patchwork (データ可視化やるなら) 最初から全部いれる必要はないよ！

Slide 29

Slide 29 text

の基本

Slide 30

Slide 30 text

例えばExcelでは、セルの内部に関数を書いて結果だけを表⽰する

Slide 31

Slide 31 text

例えばExcelでは、セルの内部に関数を書いて結果だけを表⽰するどこまでが共通の属性なのか直感的な配置から読み取る

Slide 32

Slide 32 text

の基本 a <- 5 x <- 1:5 y <- a * x 代⼊演算⼦オブジェクトコード

Slide 33

Slide 33 text

選択して実⾏(⌘↩) ⾃動で反映実⾏内容書く

Slide 34

Slide 34 text

の基本 Q. オブジェクトの中⾝が⾒たい時はどうするの？ A. オブジェクト名だけを選択して実⾏します。

Slide 35

Slide 35 text

実⾏(⌘↩)

Slide 36

Slide 36 text

の基本 Q. オブジェクトの中⾝が⾒たい時はどうするの？ A. オブジェクト名だけを選択して実⾏します。 Q. 既にある表データをRに読み込みたいんですが？ A. read_csv()関数を使います。

Slide 37

Slide 37 text

の基本 path <- “data/Book1.csv” read_csv(path)

Slide 38

Slide 38 text

実⾏(⌘↩) エラー：そんな関数ないよ

Slide 39

Slide 39 text

Don’t worry!!!

Slide 40

Slide 40 text

の基本 CRAN The Comprehensive R Archive Network

Slide 41

Slide 41 text

install.packages(“tidyverse”) の基本 1. パッケージのインストール 2. インストールされたパッケージを使う library(tidyverse)

Slide 42

Slide 42 text

実⾏(⌘↩) 読み込めたけど... コレとコレは省きたい (データじゃない)

Slide 43

Slide 43 text

Don’t worry!!!

Slide 44

Slide 44 text

の基本 library(tidyverse) path <- “data/Book1.csv” read_csv(path)

Slide 45

Slide 45 text

の基本 library(tidyverse) path <- “data/Book1.csv” read_csv(path, skip = 2)

Slide 46

Slide 46 text

No content

Slide 47

Slide 47 text

いちいちcsvに出⼒するの⾯倒だからエクセルファイルのまま読み込みたいなぁ。

Slide 48

Slide 48 text

Don’t worry!!!

Slide 49

Slide 49 text

library(tidyverse) path <- “data/Book1.xls” read_xlsx(path, skip = 2) の基本

Slide 50

Slide 50 text

エラー：そんな関数ないよ

Slide 51

Slide 51 text

Don’t worry!!!

Slide 52

Slide 52 text

library(tidyverse) library(readxl) path <- “data/Book1.xls” read_xlsx(path, skip = 2) の基本

Slide 53

Slide 53 text

そんな名前のファイル無いよ

Slide 54

Slide 54 text

無いことないだろ！！！？？さっき作ったばっかりだぞ！！

Slide 55

Slide 55 text

無いことないだろ！！！？？さっき作ったばっかりだぞ！！ほら！！！！！！

Slide 56

Slide 56 text

無いことないだろ！！！？？さっき作ったばっかりだぞ！！ほら！！！！！！ 😇 ほら...ぁ...

Slide 57

Slide 57 text

library(tidyverse) library(readxl) path <- “data/Book1.xlsx” read_xlsx(path, skip = 2) の基本

Slide 58

Slide 58 text

😇 できた

Slide 59

Slide 59 text

プログラムは思った通りには動かない。書いた通りに動くのだ。誰が⾔ったか知らないが、

Slide 60

Slide 60 text

表データの加⼯と可視化の概観 Long Wide Nested plot Figures Data table read_csv write_csv pivot_longer pivot_wider group_nest unnest ggplot ggsave wrap_plots map

Slide 61

Slide 61 text

□ Rの基本 □ Geostatisticsの基本 ✔

Slide 62

Slide 62 text

データ視覚化の⼈類史グラフの発明から時間と空間の可視化まで Michael Friendly Howrad Wainer 飯嶋貴⼦ (訳) ⻘⼟社, 2021

Slide 63

Slide 63 text

データ視覚化の⼈類史グラフの発明から時間と空間の可視化まで Michael Friendly Howrad Wainer 飯嶋貴⼦ (訳) ⻘⼟社, 2021 何かを理解しようとするならば、その始まりと推移を観察しなさい -- アリストテレス『形⽽上学』本書 p.21

Slide 64

Slide 64 text

主題図 Thematic map テーマに対して、必要な主題データを利⽤して分析した結果や、複数の主題データを重ね合せて表現した地図 h"ps://www.esrij.com/gis-guide/maps/thema6c-map/ h"ps://www.ny6mes.com/interac6ve/2016/06/24/world/europe/how-britain-voted-brexit-referendum.html

Slide 65

Slide 65 text

最初の主題図「新ヨーロッパ地図」アウグスト・クローム (独) ⽇本最初の主題地図帳「⼤⽇本帝国地産要覧図」最古の世界図バビロニアの粘⼟版 B.C.700頃 1782 1889 1650頃近代的な「データ」の概念の確⽴ジョン・グラント(英) ウィリアム・ペティ(英) 1829 最初の影付き⽐較地図「教育状況と犯罪数の⽐較統計」アンドレ=ミシェル・ゲリー(仏) アントニオ・バルビ(仏) 「回帰」概念の確⽴フランシス・ゴルトン(英) 1886 「散布図」概念の誕⽣ウィリアム・ハーシェル(英) 1833

Slide 66

Slide 66 text

Slide 67

Slide 67 text

「形状の緯度について」ニコラス・オレーム(仏) 1482 “緯度” “経度” 物理量の関数関係の可視化法 Wikipedia-en

Slide 68

Slide 68 text

𝑋 𝑌 𝑦! 𝑥! 𝑦" 𝑥" 𝑋 𝑌 𝑥! 𝑥" 𝑦! 𝑦" σʔλՄࢹԽ ࣸ૾ mapping

Slide 69

Slide 69 text

࣮ଘ ࣸ૾ʢ؍࡯ʣ σʔλ ࣸ૾ʢσʔλՄࢹԽʣ άϥϑ 𝑋 𝑌 𝑦! 𝑥! 𝑦" 𝑥" 𝑋 𝑌 𝑥! 𝑥" 𝑦! 𝑦" EBUB mapping σʔλՄࢹԽ

Slide 70

Slide 70 text

写像 (mapping) 𝑓: 𝑋 → 𝑌 𝑋 𝑌 ある情報の集合の要素を、別の情報の集合のただ1つの要素に対応づけるプロセス

Slide 71

Slide 71 text

「ナポレオンの⾏進」シャルル・ジョセフ・ミナール, 1812 -- これは、かつて作成された統計グラフィックのなかで最良のものと⾔えるだろう by エドワード・タフテ, 1983 ・地理情報（経度, 緯度）・軍の規模・軍の移動⽅向・気温