Clik here to view.
竜・龍・辰が使われている地名地図
はじめに 「竜」・「龍」・「辰」を含む地名を日本地図上にプロットしたもの。 2024年は辰年だ。ところで、長野県には「辰野町」、茨城県には「龍ケ崎市」と、「辰」や「龍」が用いられている地名は日本各地にある。私の調べでは、2024年1月現在、「竜」・「龍」・「辰」という漢字が用いられている市区町村名は8つある 。このほか、大字（おおあざ）、小字（こあざ）レベルだともっとたくさんある。...
View ArticleClik here to view.
ggplot2 パッケージでの地図出力でポイントとラベルの重なりを防ぐ方法
はじめに ggplot2 パッケージを使って sf オブジェクトから地図を出力する際、ポイントで表される地物とそのラベルの両方を出力したいことがある。このとき何も考えないと、ポイントと点が重なってしまう。 例えば、以下の地図では横浜市の市役所・区役所のポイント（赤色）がそのラベルと重なってしまっている。...
View ArticleClik here to view.
R言語のことを「統計のリンガフランカ」と表現する例
R言語はリンガフランカ 統計解析でよく使われるR言語のことを「統計のリンガフランカ」あるいは「データサイエンスのリンガフランカ」と呼ぶことがある。「リンガフランカ」(lingua franca) とは共通言語のこと。統計解析に関わる人が統計解析について議論するときの共通言語として、Rが使われているということだ。...
View ArticleR で1つのセルに複数の値が入っている場合に、複数の列にわたるダミー変数にする
今回扱いたいデータ 表型のデータで、1つのセルに複数の値が入っている場合がある。例えば、以下の例を見てみよう。太郎の好きな飲み物は「牛乳,紅茶」というコンマ区切りで1つのセルに値が入ってしまっている。 名前 性別 好きな飲み物 太郎 男 牛乳,紅茶 次郎 男 紅茶 三郎 男 コーヒー,緑茶 春子 女 緑茶,紅茶,コーヒー 夏子 女 コーヒー 秋子 女 冬子 女 緑茶,紅茶...
View ArticleR における加算代入演算子など
加算代入演算子とは プログラミングの際に、ある変数に数を足したうえでもとの変数に代入するということがよくある。例えば、以下の JavaScript では、もともと x に5が入っていたのが、それに 2 を足したうえで x に代入しなおしている。 var x = 5; x = x + 2; console.log(x); // 出力は「7」 そして、JavaScript では、x = x + 2...
View ArticleRによるツイート分析(1)：指定したTwitterアカウントのツイートをリツイート数が多い順に並べる
今回したいこと 今回したいことは、ある Twitter アカウントのツイートを、リツイート数が多い順に並べることだ。これは、Rの rtweet パッケージを使えばできる。 例えば、NHKニュースの Twitter アカウント (@nhk_news) の直近の1000個のツイートの中で、リツイートされた回数が多いものから順に並べたいとしよう。このことは、次のスクリプトで実行できる。...
View ArticleClik here to view.
Rのバージョン3.6.0がリリース
バージョン3.6.0のリリース 2019年4月26日、R のバージョン3.6.0がリリースされた。コードネームは Planting of a Tree（木を植えること）である。バージョン3.6.0での変更点の詳細は、Rアナウンスメーリングリストでのバージョン3.6.0リリースの通知（R 3.3.0 is released、英語）を参照されたい。...
View ArticleClik here to view.
Rの geosphere パッケージ(3)：ある地点から一定の方位・距離にある地点を求める
一定の方位・距離にある地点の計算 このパッケージの今までの紹介記事として以下のものがある。 Rの geosphere パッケージ(1)：2点間の距離と方位角を計算するRの geosphere パッケージ(2)：大圏航路を描く R の geoshpere パッケージでは、destPoint()...
View ArticleClik here to view.
Rの geosphere パッケージ(2)：大圏航路を描く
はじめに この記事では、R言語の geosphere パッケージを使って大圏航路を求め、それを地図上に描いてみたいと思う。 このパッケージの紹介記事の第1弾は「Rの geosphere パッケージ(1)：2点間の距離と方位角を計算する」を参照のこと。 準備 まずは、必要となるパッケージを読みこんでおく。 library("geosphere") library("tidyverse") また、東京...
View ArticleRの geosphere パッケージ(1)：2点間の距離と方位角を計算する
はじめに これから何回かに分けて、R言語の geosphere パッケージについて説明していきたいと思う。このパッケージは、球面三角法 (spherical trigonometry) に関する計算を行うパッケージであり、地球上の2点間の距離を計算したり、大圏航路を求めたりすることができる。 今回はまずパッケージの導入と、2点間の距離と方位角の計算方法について説明していきたいと思う。...
View ArticleRで複数のオブジェクトからリストを作るときに、オブジェクト名をリストの要素の名前にする
やりたいこと R言語では、あらかじめ複数個のオブジェクトを用意して、それらを1個のリストにまとめるということがよくある。以下の例では、animal, number, eraという3つのオブジェクトを作り、それらを unified という1個のリストにまとめている。 animal <- c("dog", "cat", "bear", "pig") number <- c(24, 35,...
View ArticleRで重複している箇所を取り出す方法
やりたいこと 1つのIDに対して1つのレコードしか存在しない設定であるにもかかわらず、何らかの理由で1つのIDに対して2つ以上のレコードが存在してしまっている場合がある。つまり、何らかの理由で重複が発生しているということだ。そういったときに原因を究明するには、重複しているところだけを取り出すことが有用である。 例えば、次の表では、ID が 2435...
View ArticleClik here to view.
Rによるデータクリーニング実践――政府統計からのグラフ作成を例として
はじめに データクリーニングは、データ分析の際に非常に重要なプロセスの1つであるが、データ分析の教科書では必ずしも十分に扱われていない。そこで、現実のデータクリーニングがどのように行われるかについて、一事例を紹介したいと思う。具体的には、統計処理に適したプログラミング言語のRを用いて、粗悪なデータから簡単な折れ線グラフが作成できる程度のきれいなデータにするまでのデータクリーニングを実施していく。...
View ArticleClik here to view.
R言語の上級者になりたければ『R言語徹底解説』を読もう
はじめに 『R言語徹底解説』という本を紹介したいと思う。この本は、R言語の中級者から上級者になるための助けになる本だ。この本を通じて、上級者になるために不可欠なR言語の仕組みを学べるのだ。 Wickham, H.〔著〕、石田基広・市川太祐・高柳慎一・福島真太朗〔訳〕．(2016).『R言語徹底解説』東京：共立出版． 原著：Wickham, H. (2015). Advanced R. Boca...
View ArticleClik here to view.
2017年の統計関係のアドベントカレンダー
アドベントカレンダーとは アドベントカレンダーでは1日に1つ記事を書いていく。 。 本来、アドベントカレンダー (advent calendar)...
View ArticleClik here to view.
fortunes パッケージでRに関する名言を楽しもう
fortunes パッケージとは R には fortunes というパッケージがある。これは、Rに関するさまざまな名言を収録したパッケージで、CRAN から手に入れることができる。 このパッケージに載っている名言の例をいくつか紹介しよう。 I think that it’s generally a good idea not to resist the most natural way of...
View ArticleR言語で連続する日付のベクトルを作る
開始日と終了日とを指定する R言語で、連続する日付のベクトルを作りたければ、seq() 函数を使えばよい 。例えば、2017年12月25日から2018年1月8日までの連続する15日間の日付のベクトルは以下のようにして作ることができる。 seq(as.Date("2017-12-25"), as.Date("2018-01-08"), by = "day") # [1] "2017-12-25"...
View ArticleRで文字列の先頭や末尾から決まった文字数だけ削る方法
やりたいこと データ分析の際に、文字列の先頭から決まった文字数だけ削りたいということがある。例えば、"AS101", "BD2372A", "7G66B" といった文字列の先頭から2文字分だけ削って、"101", "2372A", "66B"としたいことがある。 あるいは、文字列の末尾から決まった文字数だけ削りたいということがある。例えば、"福岡県知事", "神奈川県知事", "東京都知事",...
View Article米国のトランプ大統領の謎の covfefe ツイートをRで再現する
はじめに 米国のドナルド・トランプ大統領は Twitter を結構よく使っている 。何しろ、大統領に就任してからも日に4回ツイートしている 。 そのトランプ大統領が、2017年5月31日に以下のようなツイートをした 。このツイートでは、covfefe という謎の単語が使われていた。 Despite the constant negative press covfefe...
View ArticleCRAN にある R パッケージへのリンクとして選ぶべきURL
1つのパッケージに複数のURL CRAN では複数のURLが同じ R パッケージの紹介ページを示すことがある。例えば、2017年3月30日現在、以下の3つのURLのどれを選んでも Nippon というパッケージの紹介ページに行き着く。 https://CRAN.R-project.org/package=Nippon...
View Article
