—2019/10/17 更新

前回の06.Google Colaboratoryで時系列データで折れ線グラフを描画では、GoogleColaboratory上で、自作データをアップロードして、ggplot2を用いて折れ線グラフを描く方法を紹介しました。

今回は、次のような条件での、『4群の血糖値のデータ』をアップロードして、箱髭図、または蜂群図（beesworm）を作成します。さらに、4群の平均値に差があるかどうか、つまりone-way ANOVAを行ってみたいと思います。また、4群の平均値に差があった場合に、各2群同士で平均値に有意な差があるかどうかの検定（Tukey）も行います。

1. 普通食、糖尿病薬投与無しのマウス（group A）の血糖値
2. 普通食、糖尿病薬投与ありのマウス（group B）の血糖値
3. 高脂肪食、糖尿病薬投与無しのマウス（group C）の血糖値
4. 高脂肪食、糖尿病薬投与ありのマウス（group D）の血糖値

Windows10 Pro (1803) (OSは何でもよいです)
Chrome
Googleアカウント

Google ColaboratoryでRプログラミングを行う方法については、Google Colaboratoryで3秒でR環境構築 2019年10月14日をご覧下さい。

tidyverse

今回、自前のデータとして、以下のエクセルファイルを用いてみます。

book3.xlsx （←クリックでダウンロードできます。）

Wild Typeマウスを、普通食vs高脂肪食、糖尿病薬（血糖降下薬）投与ありvs投与無しの4群に分けて飼育し、3ヶ月後の血糖値[mg/dL]を測定したデータです。（33行5列の仮のデータです。）

Rで統計する際の自前データの記載方法については、03.Excelでcsv形式でデータを保存をご覧頂けましたら幸いです。半角英数字とピリオドのみを用いて記載していただくのがお勧めです。なお、欠損データは半角英数字の大文字で『NA』と記載してください。

以下のテンプレートへ行き、Googleログインして、『ファイル』＞『ドライブにコピーを保存』を順にクリックすると、自分のipynbファイルとして、そのままRを利用することができます。

https://colab.research.google.com/drive/1Wbt_tHd9KFEEszzXchwuFrhyCNR_KgtC

今回は、191018_ggplot-boxplot-beesworm_001.ipynb という名前にしました。

https://qiita.com/wakuteka/items/95ac758070f6f4d89a96
@wakuteka
2016年12月26日に更新
R言語 標準データセットの私的まとめ

csvを取り込みたいとき

library(tidyverse)

setwd(エクセルファイルのある場所)
data <- read_csv("data.csv")
# 

エクセルファイルを取り込みたいとき

library(tidyverse)
library(readxl)

setwd(エクセルファイルのある場所)
data <- read_excel("data.xlsx", sheet = 1)

tribble()関数を用いる方法

library(tidyverse)

d <- tribble(
  ~id, ~highFatDiet, ~medicine, ~group,	~glucose,
  A1,	"nasi",	"nasi",	A,	76,
  A2,	"nasi",	"nasi",	A,	82,

)
d

箱ひげ図は、「箱」と「ひげ」を用いて、データの中央値、ばらつきを可視化するための図です。

『中央値』とは、例えば、11個の数値からなる集合があったとき、大きさが上から数えて6番目(下から数えても6番目の大きさ)の数値です。

同様に、『第1四分位数』は、下から数えて3番目の数の値です。

『第3四分位数』は、上から数えて3番目の数の値です。

https://udemy.benesse.co.jp/office-enhance/microsoft-office/box-and-whisker-plot.html
2018/06/25
箱ひげ図の概念から作り方まで、わかりやすく解説！
→箱ひげ図の説明と、エクセルでの描画方法がのっています。お勧めです。

https://mathwords.net/tyuouchi
中央値の意味と求め方、偶数の場合
最終更新日 2019/04/13

https://kazutan.github.io/fukuokaR11/intro_ggplot2.html#box_plot
ggplot2による可視化入門
fukuoka.R #11
kazutan
2018/9/15

beesworm蜂群図(ほうぐんず)とは、数値からなる集合があったとき、それを全て点などで表示したものです。

https://wikidocs.net/40935
(04) Exercise » geom_dotplot

# beesworm 蜂群図(ほうぐんず)
# https://wikidocs.net/40935
ggplot(airquality, aes(x = factor(Month), y = Wind)) +
  geom_dotplot(binaxis = "y", stackdir = "center")

https://ninosanpo.hatenablog.com/entry/2018/08/03/141738
2018-08-03
Rmarkdownで蜂群図(ほうぐんず)を書いてみよう！

https://qiita.com/uri/items/615653e83642d7e475de
@uri
2015年06月08日に更新
学会発表のためのggplot2の設定めも

https://gist.github.com/soh-i/6332270
ggplotの図を透明背景にする

https://note.mu/shigeruhanano/n/n1e1333f293b2
Rのggplot2でグラフを描く
Shigeru Hanano （花野 滋）
2018/10/02 15:58

目次：初めての医療統計RとEZR

前：06.Google Colaboratoryで時系列データで折れ線グラフを描画

次：