Clik here to view.
精神疾患の軽症化
精神疾患の軽症化が進んでいる。統合失調症の新規発症が激減し、入院を要するレベルの気分障害がほとんどいなくなり、病棟経営の支えである慢性の統合失調症患者は急速に高齢化が進んでいる。地域移行を進めるまでもなく精神科病床数は今後10年で確実に激減する。— 斎藤環 (@pentaxxx) November 14,...
View ArticleClik here to view.
［googleAnalyticsR］調整済みオッズ比の四半期比較プロット
とりあえず図とコードだけ貼ります。縦の棒はオッズ比の95%信頼区間です。library(googleAnalyticsR)library(tidyverse)library(lubridate)##### ga_auth() account_list <- ga_account_list() ga_id <- account_list$viewId[3] gadata <-...
View ArticleClik here to view.
Rcpp: カーネル密度推定のバンド幅を一個抜き交差検証法で決める
バンド幅 h をいろいろ変えて一個抜き交差検証法で評価した対数尤度が結構なめらかな形になったのでブレント法で最適なバンド幅を選んでみた。対象としたデータはこれ。h をいろいろ変えて一個抜き交差検証法で評価した対数尤度のプロットはこちら。選ばれた最適なバンド幅で推定した密度関数がこんな感じ。Rcpp のコード：#include <Rcpp.h>usingnamespace Rcpp; //...
View ArticleClik here to view.
目で見る尤度関数（『ベイズ統計の理論と方法』より）
『ベイズ統計の理論と方法』1.4節の例を R でやってみます。ベイズ統計の理論と方法作者:渡辺澄夫出版社/メーカー:コロナ社発売日: 2012/03/01メディア:単行本購入: 1人 クリック: 4回この商品を含むブログ (8件) を見る尤度関数が正規分布で近似できるとき、いろいろ良い性質がなりたちます。関数 を標準正規分布の密度関数とし、確率モデル (ただし )...
View ArticleClik here to view.
TikzのBayesNetライブラリでプレート表現を書いてみる
プレート表現については、たしか『実践 ベイズモデリング』などに解説がある。実践 ベイズモデリング -解析技法と認知モデル-作者:豊田秀樹出版社/メーカー:朝倉書店発売日: 2017/01/25メディア:単行本（ソフトカバー）この商品を含むブログ (2件) を見るプレート表現に慣れたいので、TikzのBayesNetライブラリ（GitHub - jluttine/tikz-bayesnet: TikZ...
View ArticleClik here to view.
ggplot2 でミニ棒グラフ行列を描画する geom_barmatrix を書きました
library(tidyverse) geom_barmatrix <-function(mapping =NULL, data =NULL, stat ="identity",position ="identity", ..., width =NULL, binwidth =NULL, na.rm =FALSE, show.legend =NA, inherit.aes =TRUE){...
View ArticleClik here to view.
ggplot2 のためのいくつかの関数をパッケージ化しました
まだ説明とかぜんぜん書いてない。GitHub - abikoushi/ggsomestat: Some stat in ggplot2インストールはdevtools::install_github("abikoushi/ggsomestat")で、たぶんいけます。練習用にsearchConsoleというsearchConsoleからとってきたデータが入ってます。stat_ma は移動平均を描画する...
View ArticleClik here to view.
打ち切りデータを扱うために前に書いた R の関数群をパッケージ化しました
R の自分用の関数がごちゃごちゃしてきたのでパッケージ化して整理しようとしているのです。tidysurv とかいう恐れ多い名前をつけてしまった……。GitHub - abikoushi/tidysurv: R package for censored data handlingまだ作りかけです。dplyr、ggplot2 と survival...
View ArticleClik here to view.
ggplot2とpatchworkでヒートマップの周辺度数をプロット
パッチワークパッケージ（GitHub - thomasp85/patchwork: The Composer of ggplots）は便利ですね。library(tidyverse)library(patchwork) data("author",package ="ca") author_t <- as_data_frame(author)%>%...
View ArticleClik here to view.
ggplot2 で内訳を左から右に詳細化していく棒グラフ (2)
以前にこんなのを書いたんだが、ggplot2 で内訳を左から右に詳細化していく棒グラフ - 廿TT新しいgeomを作るまでもないことに気がついた。library(tidyverse) make_drilldown <-function(data,dimensions,metrics){ len <- length(dimensions) out...
View ArticleClik here to view.
ggplot2で検索クエリのワードクラウド
ただパッケージ動かしてみたってだけなんですけど。はい。ぼくは原則的にはワードクラウドを使う必要ない（棒グラフのほうが優れている）と思っているんですけど、最近軟弱になってきてる。library(searchConsoleR)library(ggwordcloud) scr_auth() sc_websites <- list_websites() scdata <-...
View ArticleClik here to view.
殺人事件の最近の傾向（平成19年から28年）
「殺人事件が年々増えているというのはメディアが作った嘘で実際は減少している」みたいな言説を最近よく聞く気がします。ぼくはテレビ見ないし新聞も読まないので（見たほうがいいのはわかってるけど）メディアがどういう報道をしてるのかは知らない。そこで、警察庁Webサイト（http://www.npa.go.jp/toukei/seianki/h28zuhyousakuin.htm）より最近の殺人事件の被害者数...
View ArticleClik here to view.
［ggplot2］ヒストグラムを箱ひげ図風に並べるプロット
こんな感じです。ggplot(data = iris)+ geom_grid(aes(y = Sepal.Length, x=Species),binwidth =0.1)思い通りに動かないことも多いけど公開します。皆様の暖かいアドバイスをお待ちしております。coord_flip に対応したい四角が大きいとき隣の四角とかぶってしまうのをなんとかしたいbins...
View ArticleClik here to view.
［ggplot2］ヒストグラムを箱ひげ図風に並べるプロット（2）
［ggplot2］ヒストグラムを箱ひげ図風に並べるプロット - 廿TTをちょっと修正しました。こんな感じです。source("https://gist.githubusercontent.com/abikoushi/d7e443774155de3934b754cf1a66625a/raw/74d451393ec9e1a6537122d2c39d5ee2d7559965/geom_grid.R")...
View ArticleClik here to view.
［ggplot2］ヒストグラムを箱ひげ図風に並べるプロット（3）
追記：パッケージ化したのでこちらを参照してください→ggplot2でヒストグラムを箱ひげ図風に並べてプロットするパッケージggbrickを書いた - 廿TT進捗です。［ggplot2］ヒストグラムを箱ひげ図風に並べるプロット - 廿TT［ggplot2］ヒストグラムを箱ひげ図風に並べるプロット（2） -...
View ArticleClik here to view.
ggplot2でヒストグラムを箱ひげ図風に並べてプロットするパッケージggbrickを書いた
github.comggtetrisって名前にしようかと思ったんだけど、もうある（GitHub - EmilHvitfeldt/ggtetris: Create Tetris Chart Visualizations in...
View ArticleClik here to view.
ggbrick: histogram like dot plot based on ggplot2
Hello everyone. I introduce an R package ggbrick. github.comMy English is poor. If you don't understand my writing, please use comment field (コメントを書く).ggbrick provides the function geom_brick which is...
View ArticleClik here to view.
腸内細菌叢のデータで遊ぼう（Kostic et al, 2015）
Kostic et al, 2015 (https://www.cell.com/cell-host-microbe/fulltext/S1931-3128(15)00021-9)...
View ArticleClik here to view.
可視化で理解するKPIツリー：コンバージョンの寄与度分解
背景昔Web系コンサルをやっていたときKPIツリーとかを習ったんだけどぼくはこれが苦手だった。例えばあるサイトのコンバージョン（CV）はサイトへの訪問数（SS）×コンバージョンレート（CVR）だからコンバージョンというゴールを増やしたかったら、SSを増やすかCVRを増やすしかない。KPIツリーこうやって順々に要素を順々に分解していけばコンバージョンを増やす方策をロジカルに考えられますね、というわけ。...
View ArticleClik here to view.
棒グラフをとにかくいっぱい並べて書きたい
set.seed(2) y <- lapply(1:7, rexp, n=26) names(y)<- sapply(1:7,function(x)paste0(sample(LETTERS,10),collapse ="")) tab <-list(y=y,x=LETTERS) class(tab)<-"minibarTable" plot.minibarTable...
View Article