Clik here to view.
大森・宇津公式を強度関数とした非定常ポアソン過程の重ね合わせによる熊本地震の余震の分析
スクレイピング使用するツールは R です。データは 気象庁｜地震情報から取得しました。rvest パッケージを使うとかんたんです。また XMLパッケージの readHTMLTable 関数を使う方法もあります。R から HTML の表を読み込む - 廿TT速報は一定時間たつと流れていってしまうため、数日おきにデータをダウンロードしてきて dplyr の full_join...
View ArticleClik here to view.
Clik here to view.
［dlm］状態空間モデルでトレンドと広告の効果を分離して推定する
はじめにStanで統計モデリングを学ぶ(7): 時系列の「トレンド」を目視ではなくきちんと統計的に推定する - 東京で働くデータサイエンティストのブログをみてください。上記事では Stan で状態空間モデルを推定しているので、ここでは R の dlm パッケージを使ってやってみます。モデル化コンバージョン（CV）がベースラインとなるトレンドと三種類の広告の効果の和で表わされるというモデルです。ここで...
View ArticleClik here to view.
岩波データサイエンス Vol.1 の年輪の例題を dlm でやる
DLMノイズが正規分布し、変数間の関係が線形の状態空間モデル動的線形モデル（Dynamic Linear Model; DLM）と呼ばれる。動的線形モデルは以下の式で表現できる。推定問題現時刻を kとするとき、推定問題はつぎの三つに分類できる。予測（prediction）：時刻（k-n）までのデータに基づいて現在の値 を推定する。フィルタリング（filtering）：時刻...
View ArticleClik here to view.
気温・湿度と熱中症の危険度（暑さ指数）の表
厚生労働省：職場における熱中症の予防についてより暑さ指数（WBGT）の表をプロットしてみた。WBGT値警戒水準 25℃未満注意25℃～28℃警戒28℃～31℃厳重警戒31℃以上危 険そろそろ熱中症に気をつけましょう。library(pdftools)library(pipeR)library(tidyr)library(dplyr)library(ggplot2) path1...
View ArticleClik here to view.
［rstan］横浜駅SFの移入と離脱のモデル
経緯この減衰曲線はどういう関数で説明するのが適切か考えている（指数関数ではない）https://t.co/PYdROur8XJpic.twitter.com/nFkikWjH1x— イスカリオテの湯葉 (@yubais)...
View ArticleClik here to view.
Alluvial diagram で Google アナリティクスのデータをプロットして遊ぼう
Alluvial diagram（アルビアルダイアグラム）をかんたんにプロットできる ggalluvial パッケージが、Rで解析：Alluvial diagramsをプロットしませんか「ggalluvial」パッケージで紹介されていた。これを使ってみようと思う。入口ページと出口ページまずは RGA パッケージで Googleアナリティクスのデータを R...
View ArticleClik here to view.
自殺者数の推移（原因・年代別）
自殺白書に関する報道は東京新聞と産経新聞で見出しの付け方が対照的だった。東京新聞:若者と高齢者の自殺深刻 政府、１６年版白書:社会(TOKYO Web)自殺白書 「経済」理由の自殺半減 厚労省「法改正や法律相談の充実が奏功」 （産経新聞） -...
View ArticleClik here to view.
失業と自殺は関係あるのか
男性に関して言えば、失業率と自殺率は強い相関を持つことが指摘されている。例えば データえっせい: 失業と自殺の時系列的相関や『この経済政策が民主主義を救う』などを参照。この経済政策が民主主義を救う: 安倍政権に勝てる対案作者:松尾匡出版社/メーカー:大月書店発売日: 2016/01/20メディア:単行本この商品を含むブログ (8件)...
View ArticleClik here to view.
ggplot2 で接線場や方向場を描く
接線場微分方程式の解の振る舞いをみるために接線場を描いてみる。以下で与えられるロジスティック方程式を考える。接線場は t-x平面の各点に傾き の小さな線分を描いたものである。logis <-function(t,x,a=1){ a*x*(1-x)}#点(t1,x1)を通り, 傾きが slope の直線 TheLine <-function(t1,x1,slope,t){...
View ArticleClik here to view.
ロトカ・ヴォルテラの方程式のパラメータ推定
ロトカ・ヴォルテラの方程式については、abrahamcow.hatenablog.comを見てください。hare–lynxデータここに一組のデータ・セットがあります。Hare=c(30,47.2,70.2,77.4,36.3,20.6,18.1,21.4,22,25.4,27.1,40.3,57,76.6,52.3,19.5,11.2,7.6,14.6,16.2,24.7)...
View ArticleClik here to view.
［KFAS］0-1データの状態空間モデル（打者の調子の波のモデル化 3）
とりあえずプロットAlbert (2008) 打者の調子の波のモデル化 - 廿TT打者の調子の波のモデル化（幾何分布編） - 廿TTと同じく、カルロス・ギーエンの打撃成績のデータを使います。y <-...
View ArticleClik here to view.
ぼくの考えたグラフ三原則
目次Request for Comments（ご意見求む）の精神です。良いグラフはフラットデザインである良いグラフはものさしで測れる良いグラフは色に頼りすぎない1....
View ArticleClik here to view.
原点がゼロでないグラフ、２軸のグラフ
前置きここに書くことはすでに言い尽くされているような気がしますが、どうもなかなかご理解頂けていないようなので、改めてしつこく述べます。経緯下記のツイートを参照。なんだこりゃ。TV屋は相変わらずヘボなグラフ作ってんな・・・・原点がゼロじゃないグラフを作るのはいいかげんにやめろよと ＞ 原油価格はこんなに下がっているのにガソリン価格なぜ下がらない ......
View ArticleClik here to view.
ロングテールSEOのためのパレート図の紹介と RGoogleAnalytics による実践
モチベーションロングテールSEOとは検索エンジンからのサイトへの訪問数を増やしたい場合、みんながいっぱい検索してる人気のあるキーワード（ビッグワード）からのアクセス数を稼ぐのは競合とか多くて大変。それよりも、ニッチなキーワード（複合キーワード、かけあわせワード）を積み重ねてあくせくアクセスしてもらおうみたいな話。例えば「家」とか「猫」とかのキーワードはいっぱい検索されているだろうけど、これらのワード...
View ArticleClik here to view.
3D円グラフとともに滅するべき棒グラフの省略表現
本稿の主張棒グラフをかくときは、縦軸をつけ、原点を 0 にしましょう。あと、大学の先生の教育面での取り組みがもっと評価される世の中になって欲しいです。はじめに：ちゃんと叱ってくれる先生は評価しよう3D円グラフふと思い立って「3D円グラフ」で検索してみたところ、なんとトップに奥村先生のページが！3D円グラフを使うのはやめよう | Okumura's...
View ArticleClik here to view.
縦軸に目盛のないグラフ（インターネットもデータ分析も好きなぼくが、入社一年未満でWeb解析コンサルタントを辞めた理由）
主張縦軸に目盛のない棒グラフや折れ線グラフを描くのはやめよう。棒グラフ以前の職場でよく見かけた棒グラフがこのようなものでした。これは万能棒グラフです。2007-10-22 - 定義すれば存在する目盛のベースラインが 0...
View ArticleClik here to view.
Clik here to view.
（R + Google アナリティクス）モザイクプロットでユーザーの構成比を把握
モザイクプロットの紹介モザイクプロットはクロス集計表を図示するための統計グラフです。例えば2015年1月にこのブログにアクセスしたユーザーを、年代、性別ごとにクロス集計すると下表のようになります。 18-24 25-34 35-44 45-54 55-64 65+ female 292 487 253 93 49 28 male 1155 4053 2166 590 217...
View ArticleClik here to view.