【ggplot2】タグの記事一覧｜ドクターフント(Dr. Hund)

[R for beginners] Creating Error Bars for Bar Graphs and Line Graphs in ggplot2

brainblog — Fri, 24 Mar 2023 13:13:56 +0000

The ‘T’ symbol on bar graphs and line graphs represents error bars that extend above and below the data point, indicating standard error or standard deviation. Adding error bars in R’s ggplot2 is easy. Here’s a step-by-step guide:

Contents

Preparing the Data
Calculating the Mean, Standard Deviation, and Standard Error of the Data
- Use group_by and summarise_all functions in dplyr to Calculate
Drawing Error Bars in Bar Graphs
Drawing Error Bars on a Line Graph
Adding Error Bars to Grouped Data
- Adding Error Bars to Grouped Bar Graphs
- Drawing Error Bars on a Grouped Line Graph

Preparing the Data

We’ll be using tidyverse for this. If you have your own data, skip ahead.

# Set working directory
setwd("~/Rpractice/")

# Load tidyverse
library(tidyverse)

# Generate random data
dat <- list(
  X <- rnorm(50, 30, 10),
  Y <- rnorm(50, 50, 5),
  Z <- rnorm(50, 40, 15)
)
df <- data.frame(matrix(unlist(dat), nrow=50))
colnames(df) <- c("A","B","C")

# Transform the data from wide to long
df.long <- pivot_longer(df, cols = A:C, 
                        names_to = "Categories", 
                        values_to = "Values")

Now we have data that looks like this:

> head(df.long)
# A tibble: 6 x 2
  Categories Values
         
1 A            26.8
2 B            53.7
3 C            27.3
4 A            31.0
5 B            58.7
6 C            56.2
>

To add error bars to bar graphs and line graphs, you need the mean and standard error (or standard deviation) of the data. We use pivot_longer to transform the data from wide to long format.”

Calculating the Mean, Standard Deviation, and Standard Error of the Data

To add error bars, we need to calculate the mean, standard deviation, and standard error. If you already have this data, you can skip this section.

Calculating these values is easy with the dplyr package in tidyverse.

Use group_by and summarise_all functions in dplyr to Calculate

a <- group_by(df.long, Categories) %>% 
  summarise_all(list(mean = ~mean(.), 
                     sd = ~sd(.), 
                     se = ~sd(.)/sqrt(length(.))))

We can simplify this code further,

a <- group_by(df.long, Categories) %>% 
  summarise_all(list(mean = mean, 
                     sd = sd, 
                     se = ~sd/sqrt(length(.))))

Let’s take a look at the data for ‘a’

> a
# A tibble: 3 x 4
  Categories  mean    sd    se
          
1 A           30.1  9.93 1.40 
2 B           49.2  5.00 0.708
3 C           43.9 15.5  2.20

We have generated a distribution that looks like this. Since the data was generated randomly, the values may vary slightly if you follow the same steps.

Drawing Error Bars in Bar Graphs

We will use the calculated data to create a bar graph and specify the error bars using geom_errorbar(). First, let’s specify only ymin and ymax in geom_errorbar() and take a look at the graph.

ggplot(a, aes(x = Categories, y = mean, fill = Categories))+
  geom_bar(stat = "identity") +
  geom_errorbar(aes(ymin = mean - se, ymax = mean + se))

The width is too wide. Let’s adjust it using the width argument. We will narrow both the bar graph and error bars.

ggplot(a, aes(x = Categories, y = mean, fill = Categories))+
  geom_bar(stat = "identity", width = 0.6) +
  geom_errorbar(aes(ymin = mean - se, ymax = mean + se),
                width = .1)

It’s a good to adjust the width according to the size of the output image.

Drawing Error Bars on a Line Graph

The process for adding error bars to a line graph is the same as above. First, draw the line graph and then add geom_errorbar().

ggplot(a, aes(x = Categories, y = mean)) +
  geom_line(group = 1) +
  geom_point(size = 3) +
  geom_errorbar(aes(ymin = mean - se, ymax = mean + se),
                width = .1)

Adding Error Bars to Grouped Data

Let’s generate grouped data.

# Add ID column
data <- df.long %>% 
  tibble::rownames_to_column(var = "ID")

# Convert ID from character to numeric
data$ID <- as.numeric(data$ID)

# Assign 1 to ID <= 75 and 0 to ID >= 76 to create the "group" column
data <- mutate(data, group = ifelse(ID < 76, 1, 0))

# Calculate Mean, SD, and SE
b <- group_by(data, group, Categories) %>% 
  summarise_at(vars(Values), list(mean = ~mean(.), 
                     sd = ~sd(.), 
                     se = ~sd(.)/sqrt(length(.))))

# Convert group column values from numeric to character
b$group <- as.character(b$group)

# Check the data included in b
> head(b)
# A tibble: 6 x 5
# Groups:   group 
  group Categories  mean    sd    se
           
1 0     A           31.4 11.1  2.21 
2 0     B           51.0  4.31 0.862
3 0     C           33.9 13.5  2.70 
4 1     A           29.5 10.4  2.08 
5 1     B           49.5  4.02 0.804
6 1     C           36.5 15.8  3.16 
>

Now that we have grouped data, let’s draw a bar graph with error bars.

Adding Error Bars to Grouped Bar Graphs

By specifying position = "dodge", you can create a grouped bar chart like this.

ggplot(b, aes(x = Categories, y = mean, fill = group))+
  geom_bar(stat = "identity", width = 0.6, position = "dodge")

Let’s add error bars to this plot.

ggplot(b, aes(x = Categories, y = mean, fill = group))+
  geom_bar(stat = "identity", width = 0.6, position = "dodge")+
  geom_errorbar(aes(ymin = mean - se, ymax = mean + se), 
                    position = "dodge", width = .1)

It is shifted to the center. You can adjust the position with position = position_dodge().

ggplot(b, aes(x = Categories, y = mean, fill = group)) +
  geom_bar(stat = "identity", width = 0.6, position = "dodge") +
  geom_errorbar(aes(ymin = mean - se, ymax = mean + se),
                position = position_dodge(0.6), width = .1)

You’ve created a nice bar graph!

Drawing Error Bars on a Grouped Line Graph

The same can be done for a grouped line graph as well.

ggplot(b, aes(x = Categories, y = mean, group = group, color = group)) +
  geom_line() +
  geom_point(size = 3) +
  geom_errorbar(aes(ymin = mean - se, ymax = mean + se),
                width = .1, color = "black")

The error bars are on top and overlapping, making it difficult to see. Let’s move the error bars to the back first.

ggplot(b, aes(x = Categories, y = mean, group = group, color = group)) +
  geom_errorbar(aes(ymin = mean - se, ymax = mean + se),
                width = .1, color = "black") +
  geom_line() +
  geom_point(size = 3)

By shifting the position to the left or right using position_dodge(), overlapping of error bars can be avoided.

ggplot(b, aes(x = Categories, y = mean, group = group, color = group)) +
  geom_errorbar(aes(ymin = mean - se, ymax = mean + se),
                width = .1, color = "black", position = position_dodge(.2)) +
  geom_line(position = position_dodge(.2)) +
  geom_point(size = 3, position = position_dodge(.2))

Now you can add error bars to line graphs. I hope this was helpful.

R, ggplot2を使った論文に使えるグラフテクニック【まとめ】

brainblog — Fri, 24 Mar 2023 12:25:12 +0000

Rではキレイなビジュアルのグラフを無料で作ることができます。高いソフトウェアは不要で、論文に十分使える図ができます！特にggplot2はその使いやすさや見た目のキレイさからとても好評です。

様々なデータがあると思いますが、一番良く使うのは箱ひげ図や散布図ではないでしょうか。それらをよりわかりやすくビジュアライズしていきたい方向けの基礎的なテクニックです！

Contents

ボックスプロット(箱ひげ図)とドットプロットを描く
ggplot2でバイオリンプロットを描く
ggplot2を使って散布図を描く
- 散布図と近似直線を描く方法
複数データはfacet_gridとfacet_wrapを使う
geom_barで棒グラフを作図する
geom_errorbarで棒グラフと折れ線グラフにエラーバーをつける
geom_tileで相関行列のヒートマップをつくる
背景の一部(四角)に色をつけるgeom_rect
ROC曲線を描く

ボックスプロット(箱ひげ図)とドットプロットを描く

このグラフはドットプロットと呼ばれます。箱ひげ図の上にドットプロットを重ねています。

箱ひげ図(box plot)とドットプロット(dot plot)は実際のデータがどのような分布になっているのかみるのに有用です。これで外れ値などもよくわかりますし、誤魔化している感じがない正直なプロットです。

base機能とbeeswarmパッケージを用いた方法とggplot2を使う方法があります。ggplot2を使う場合にはggbeeswarmというパッケージを使うとよりきれいになります。

baseのboxplotとbeeswarmライブラリを用いて作図する

箱ひげ図は、baseと呼ばれるRに元々入っている機能を使ったものです。箱ひげ図はかんたんにかけますが、その後のドットプロットについてはbeeswarmというライブラリを使います。beeswarmは蜂群図とも呼ばれ、点を散らばらせることで視覚的に印象的にみせる効果があります。beeswarmライブラリで簡単に作図をすることができます。

x軸の並び替えもこちらに書いています。

https://brain-storm.space/beeswarm_r/698/

ggplot2を用いてボックスプロットとドットプロットを作図する

ggplot2を使うのに慣れた人向けの方法です。これは単純なドットプロットを描く方法で、いわゆるbeeswarm(蜂群図)ではありませんが、視覚的効果は十分にあります。

https://brain-storm.space/rggplot2_boxplot_dotplot/620/

ggbeeswarmを用いてggplot2のドットを散らばらせる

ggplot2でもbeeswarm(蜂群図)を作成することができます。これにはggplot2に加えて、ggbeeswarmというライブラリを使います。点の散らし方にもいくつかの方法があります。

https://brain-storm.space/ggbeeswarm/820/

ggplot2でバイオリンプロットを描く

バイオリンプロットも複数のデータを比較する時にとても便利です。とてもキレイな図になります。この中では、バイオリンプロットにドットプロットや箱ひげ図の重ね方を説明しています。バイオリンプロットはggplot2のgeom_violinで描くことができます。

https://brain-storm.space/violin-plot/871/

ggplot2を使って散布図を描く

相関を数字で出すよりも実際に散布図を出したほうがより直感的にもわかりやすくなります。これも正直な印象を与えることができます。　ggplot2を用いて、geom_pointで散布図を、geom_smoothまたはstat_smoothで近似曲線を表示できます。

散布図と近似直線を描く方法

https://brain-storm.space/rggplot_scatterplot/751/

複数データはfacet_gridとfacet_wrapを使う

複数データを一度に表示したいとき、一つ一つ出していたのでは大変です。そんなときはfacet_gridやfacet_wrapを使うと一度にできて、見栄えのあるグラフができます。

またタイトルラベルなどの変え方にも少しテクニックが必要なので、その説明も加えています。

https://brain-storm.space/r_facet/843/

geom_barで棒グラフを作図する

棒グラフは最も基本的なグラフの一つです。エクセルでも出来てしまうので、それほど困らないかもしれませんが、ggplot2でも簡単に出来てしまいます。次のリンクで紹介をしています。

https://brain-storm.space/bargraph/1025/

geom_errorbarで棒グラフと折れ線グラフにエラーバーをつける

棒グラフや折れ線グラフは平均値の大小や推移などを可視化するのによく使われます。データ解析ではそこにエラーバーをつけたりしますが、そのやり方も知らなくては作図ができません。

さらにRでの標準偏差や標準誤差の算出方法も説明しています。計算で標準誤差などを出せば簡単に作図ができます。

https://brain-storm.space/ggplot2_errorbar/895/

geom_tileで相関行列のヒートマップをつくる

相関行列をつくったとき、それをテーブルにしても十分ですが、それを視覚的にわかりやすくするとより印象的な図になります。一目でデータが理解出来るのは、大きなメリットになります。

ggplot2以外のいくつかのパッケージでも描くことはできますが、ggplot2を使った方法を説明しています。それを紹介したのは次のリンクになります。

https://brain-storm.space/corrheatmap/942/

ggplot2を使うときにデータを横長のデータから縦長のデータへ変換する必要があるときがあります。そんなときはpivot_longerを使うことで簡単に変換をすることができます。これは次の記事で詳細に説明をしています。

https://brain-storm.space/tidyr_gather_pivot_longer/784/

今回紹介した作図法だけでも、論文のグラフの多くは作図ができると思います。今後もさらなる方法を広めて行きたいと思っています。

背景の一部(四角)に色をつけるgeom_rect

背景に色をつけて、その部分を強調したり、領域を分けたいことがあります。そんな時はgeom_rectを使います。

https://brain-storm.space/geom_rect/1191/

ROC曲線を描く

ROC curve

ROC曲線を書くのは思ったよりも簡単です。Rの初心者でも短時間で書くことができます。AUC(Area under the curve)の比較もあっという間です。

https://brain-storm.space/roc-curve/1212/

今後もさらに追加していく予定です。お役に立ちましたら幸いです。

【R初心者向け】ggplot2で棒グラフと折れ線グラフにエラーバー(error bar)を書く

brainblog — Fri, 24 Mar 2023 12:20:55 +0000

棒グラフや折れ線グラフなどについている”T”が上下についているものがエラーバーです。エラーバーは標準誤差や標準偏差を示しています。

Rのggplot2を使えば簡単にできます。ステップバイステップで書いていきます。

Contents

データの準備する
データの平均、標準偏差、標準誤差を算出する
- dplyrのgroup_byとsummarise_allを使う
棒グラフでエラーバーを描く
折れ線グラフでエラーバーを描く
グループ化されたデータでのエラーバーを加える
- グループ化された棒グラフにエラーバーを追加する
- グループ化された折れ線グラフでのエラーバーを書く

データの準備する

tidyverseを利用します。自分のデータを使う場合には読み飛ばしてください。

# 作業ディレクトリのセット
setwd("~/Rpractice/")

# tidyverseの起動
library(tidyverse)

# ランダムデータの生成
dat <- list(
  X <- rnorm(50, 30, 10),
  Y <- rnorm(50, 50, 5),
  Z <- rnorm(50, 40, 15)
)
df <- data.frame(matrix(unlist(dat), nrow=50))
colnames(df) <- c("A","B","C")

#データの変形(横から縦へ)
df.long <- pivot_longer(df, cols = A:C, 
                        names_to = "Categories", 
                        values_to = "Values")

これでこのようなデータができました。

> head(df.long)
# A tibble: 6 x 2
  Categories Values
         
1 A            26.8
2 B            53.7
3 C            27.3
4 A            31.0
5 B            58.7
6 C            56.2
>

棒グラフと折れ線グラフにエラーバーを加えるには、平均値と標準誤差(または標準偏差)が必要です。pivot_longerで縦長データにしています。

https://brain-storm.space/tidyr_gather_pivot_longer/784/

データの平均、標準偏差、標準誤差を算出する

エラーバーのためには平均(mean)、標準偏差(sd)、標準誤差(se)を算出します。すでにこれらのデータがある場合にはここも読み飛ばしてください。

その算出方法も慣れれば簡単です！tidyverseに入っているdplyrを使っていきましょう。

dplyrのgroup_byとsummarise_allを使う

(summariseはアメリカ英語でsummarizeでもOKです)

a <- group_by(df.long, Categories) %>% 
  summarise_all(list(mean = ~mean(.), 
                     sd = ~sd(.), 
                     se = ~sd(.)/sqrt(length(.))))

これはもう少し簡素化できます。どちらを使っても大丈夫です。

a <- group_by(df.long, Categories) %>% 
  summarise_all(list(mean = mean, 
                     sd = sd, 
                     se = ~sd/sqrt(length(.))))

ではaのデータをみてみます。

> a
# A tibble: 3 x 4
  Categories  mean    sd    se
          
1 A           30.1  9.93 1.40 
2 B           49.2  5.00 0.708
3 C           43.9 15.5  2.20

このような分布になりました。データはランダムで生成しましたので、同じ手順をなぞってもこの数値は若干変化します。

棒グラフでエラーバーを描く

計算をしたデータを使って棒グラフをかきます。エラーバーはgeom_errorbarで指定します。まずはgeom_errorbarでyminとymaxのみを指定して、グラフをみてみましょう。

ggplot(a, aes(x = Categories, y = mean, fill = Categories))+
  geom_bar(stat = "identity") +
  geom_errorbar(aes(ymin = mean - se, ymax = mean + se))

横幅が広すぎますね。これをwidthで調整します。バーグラフ、エラーバーともに狭くしてみます。

ggplot(a, aes(x = Categories, y = mean, fill = Categories))+
  geom_bar(stat = "identity", width = 0.6) +
  geom_errorbar(aes(ymin = mean - se, ymax = mean + se),
                width = .1)

幅は出力する画像の大きさに合わせて適宜調整をするといいです。

折れ線グラフでエラーバーを描く

折れ線グラフも上と同じ要領です。折れ線グラフを描いて、geom_errorbarを付け足します。

ggplot(a, aes(x = Categories, y = mean)) +
  geom_line(group = 1) +
  geom_point(size = 3) +
  geom_errorbar(aes(ymin = mean - se, ymax = mean + se),
                width = .1)

グループ化されたデータでのエラーバーを加える

グループ化したデータを生成してみます。

# IDを追加する
data <- df.long %>% 
  tibble::rownames_to_column(var = "ID")

# IDを文字列から数値へ変換する
data$ID <- as.numeric(data$ID)

# IDが75以下を1, 76以上を0と割り当てて、group列を作る
data <- mutate(data, group = ifelse(ID < 76, 1, 0))

# mean, SD, SEを算出
b <- group_by(data, group, Categories) %>% 
  summarise_at(vars(Values), list(mean = ~mean(.), 
                     sd = ~sd(.), 
                     se = ~sd(.)/sqrt(length(.))))

# group列の数値を文字列に変換する
b$group <- as.character(b$group)

# bに含まれているデータを確認する
> head(b)
# A tibble: 6 x 5
# Groups:   group 
  group Categories  mean    sd    se
           
1 0     A           31.4 11.1  2.21 
2 0     B           51.0  4.31 0.862
3 0     C           33.9 13.5  2.70 
4 1     A           29.5 10.4  2.08 
5 1     B           49.5  4.02 0.804
6 1     C           36.5 15.8  3.16 
>

これでグループ化されました。このデータでエラーバーの入った棒グラフを描いてみましょう。

グループ化された棒グラフにエラーバーを追加する

position = “dodge”をするとこのような棒グラフができます。

ggplot(b, aes(x = Categories, y = mean, fill = group))+
  geom_bar(stat = "identity", width = 0.6, position = "dodge")

これにエラーバーを加えていきます。

ggplot(b, aes(x = Categories, y = mean, fill = group))+
  geom_bar(stat = "identity", width = 0.6, position = "dodge")+
  geom_errorbar(aes(ymin = mean - se, ymax = mean + se), 
                    position = "dodge", width = .1)

真ん中に寄ってしまいました。この位置はposition = position_dodge()で調整できます。

ggplot(b, aes(x = Categories, y = mean, fill = group)) +
  geom_bar(stat = "identity", width = 0.6, position = "dodge") +
  geom_errorbar(aes(ymin = mean - se, ymax = mean + se),
                position = position_dodge(0.6), width = .1)

これでいい棒グラフができました！

グループ化された折れ線グラフでのエラーバーを書く

折れ線グラフも同じようにできます。

ggplot(b, aes(x = Categories, y = mean, group = group, color = group)) +
  geom_line() +
  geom_point(size = 3) +
  geom_errorbar(aes(ymin = mean - se, ymax = mean + se),
                width = .1, color = "black")

エラーバーが前面に出ていて重なっているので見にくくなっています。まずエラーバーを後ろに持っていきます。

ggplot(b, aes(x = Categories, y = mean, group = group, color = group)) +
  geom_errorbar(aes(ymin = mean - se, ymax = mean + se),
                width = .1, color = "black") +
  geom_line() +
  geom_point(size = 3)

エラーバーが重なりは左右にずらすことで解消ができます。さきほどのposition_dodgeを使います！

ggplot(b, aes(x = Categories, y = mean, group = group, color = group)) +
  geom_errorbar(aes(ymin = mean - se, ymax = mean + se),
                width = .1, color = "black", position = position_dodge(.2)) +
  geom_line(position = position_dodge(.2)) +
  geom_point(size = 3, position = position_dodge(.2))

これで折れ線グラフのエラーバーを加えることができました。お役に立ちましたら幸いです。

Violin Plot in R – How to Draw It –

brainblog — Fri, 24 Mar 2023 10:33:25 +0000

There are several ways to visualize data, and one of them is the violin plot! This method is useful for comparing multiple sets of data, and it has an appealing appearance.

First, you’ll need to prepare your data. While you’re at it, you can also check the data distribution using a dot plot.

Contents

Data Generation and Distribution Checking
Draw a Simple Violin Plot
Add Color to Violin Plot
- Add Color to the Borders
- Fill the interior of Violin Plot
Add Mean or Median to the Violin Plot
- Add Mean
- Add Median
Change the Degree of Smoothing
Overlaying a Box Plot

Data Generation and Distribution Checking

Now, let’s create some data. If you have your own data, feel free to use that instead.

# Set working directory
setwd("~/Rpractice/")
# load tidyverse and ggbeeswarm
library(tidyverse)
library(ggbeeswarm)

# generate sample data
dat <- list(
  X <- rnorm(100, 5, 10),
  Y <- rnorm(100, 20, 10),
  Z <- rnorm(100, 15, 15)
)

# format the data into a more usable shape using pivot_longer
df <- data.frame(matrix(unlist(dat), nrow=100))
colnames(df) <- c("A","B","C")
df.long <- pivot_longer(df, cols = A:C, names_to = "Categories", values_to = "Values")

# draw dotplot
ggplot(df.long, aes(x = Categories, y = Values))+
  geom_beeswarm(aes(color = Categories),
                size = 2,
                cex = 2,
                alpha = .8)+
  theme_classic()+
  theme(legend.position = "none")

The dot plot displays the data distribution and can be used to confirm that the data is appropriate for creating a violin plot.

Draw a Simple Violin Plot

To draw a violin plot using ggplot2, you can utilize the geom_violin() function. To create a clean and simple plot, set the background color to white using the theme_classic() function.

ggplot(df.long, aes(x = Categories, y = Values))+
  geom_violin()+
  theme_classic()

The ends of the violin plot may appear cut off. By overlaying a dot plot on top of the violin plot, you can address this issue. To do so, you can use either geom_dotplot() or geom_beeswarm(), which are both part of the ggplot2 package.

ggplot(df.long, aes(x = Categories, y = Values))+
  geom_violin()+
  geom_beeswarm(aes(color = Categories),
                size = 2,
                cex = 2,
                alpha = .8)+
  theme_classic()

You can see that the ends of this violin plot are cut off at the minimum and maximum values. If you don’t want to cut off the ends, you can use geom_violin(trim = FALSE) to specify this preference.

ggplot(df.long, aes(x = Categories, y = Values))+
  geom_violin(trim = FALSE)+
  geom_beeswarm(aes(color = Categories),
                size = 2,
                cex = 2,
                alpha = .8)+
  theme_classic()

The violin plot may stretch vertically up and down, even in areas where there are no data points.

Add Color to Violin Plot

Add Color to the Borders

To add color to the border of the violin plot, you can use aes(color = ) inside the geom_violin() function.

ggplot(df.long, aes(x = Categories, y = Values, color = Categories))+
  geom_violin()+
  theme_classic()

Fill the interior of Violin Plot

If you want to fill the interior of the violin plot with color, you can use aes(fill =) inside the geom_violin() function.

ggplot(df.long, aes(x = Categories, y = Values, fill = Categories))+
  geom_violin()+
  theme_classic()

There are several ways to change the color of the violin plot. One way is to use the scale_fill_brewer() function to specify the color scheme.

ggplot(df.long, aes(x = Categories, y = Values, fill = Categories))+
  geom_violin()+
  scale_fill_brewer(palette = "Set2")+
  theme_classic()

Add Mean or Median to the Violin Plot

To add the mean or median to the violin plot, you can use the stat_summary() function.

Add Mean

ggplot(df.long, aes(x = Categories, y = Values, color = Categories))+
  geom_violin()+
  stat_summary(fun = mean, geom = "point", 
               shape = 16, size = 2, color = "red")+
  theme_classic()

The shape parameter in stat_summary() is the same as pch in base R. Here’s a list of shapes that correspond to each numerical value:

Add Median

ggplot(df.long, aes(x = Categories, y = Values, color = Categories))+
  geom_violin()+
  stat_summary(fun = median, geom = "point", 
               shape = 3, size = 2, color = "red")+
  theme_classic()

Change the Degree of Smoothing

To change the degree of smoothing in the violin plot, you can use the adjust parameter inside the geom_violin() function. The default value for adjust is 1.

To decrease the degree of smoothing, you can set adjust to a smaller value. For example, to set adjust to 0.2, you can use the following code:

ggplot(df.long, aes(x = Categories, y = Values, fill = Categories))+
  geom_violin(adjust = .2)+
  theme_classic()

Conversely, if you want to increase the degree of smoothing in the violin plot, you can set the adjust parameter to a larger value.

ggplot(df.long, aes(x = Categories, y = Values, fill = Categories))+
  geom_violin(adjust = 2)+
  theme_classic()

If you increase the degree of smoothing too much, the violin plot can become overly smoothed and lose important details.

Overlaying a Box Plot

Overlaying a violin plot with a box plot is a common technique in data visualization, and it can be a powerful way to display data.

ggplot(df.long, aes(x = Categories, y = Values, fill = Categories))+
  geom_violin()+
  geom_boxplot(width = .1, fill = "white")+
  theme_classic()

To hide the outliers in the box plot when overlaying it with a violin plot, you can use the outlier.color parameter inside the geom_boxplot() function.

ggplot(df.long, aes(x = Categories, y = Values, fill = Categories))+
  geom_violin()+
  geom_boxplot(width = .1, fill = "white", outlier.color = NA)+
  theme_classic()

You can fill the box plot with black color and add a white circle at the median value.

ggplot(df.long, aes(x = Categories, y = Values))+
  geom_violin()+
  geom_boxplot(width = .1, fill = "black", outlier.color = NA) +
  stat_summary(fun = median, geom = "point", fill = "white", shape = 21, size = 3) +
  theme_classic()

With the information provided, I believe you can now create a violin plot. I hope this guidance proves helpful!

【R簡単テクニック】背景に四角の領域(背景)をいれる

brainblog — Mon, 06 Mar 2023 22:01:58 +0000

ggplotを使ったグラフの中に、四角の領域を入れる方法を紹介します。とても簡単です。ggplot2を使うためまずはライブラリのtidyverse(またはggplot2)を起動します。

library(tidyverse)

資格の領域を作るにはgeom_rectを使います。その領域をxmin, xmax, ymin, ymaxで指定します。その領域の色はデフォルトで灰色(grey20)になっています。

Contents

geom_rectの使い方
複数の領域を書きたい時
- この時はデータフレームで領域を指定すると楽になります
他のグラフの背景にする

geom_rectの使い方

まずは書いてみます

四角の領域を書くにはgeom_rectを使います。x軸の最小値(xmin)と最大値(xmax)、y軸の最小値yminと最大値ymaxをしています。すると次のようなグラフが描けます。

ggplot() +
　geom_rect(aes(xmin=1, xmax=3, ymin=2, ymax=5))

色を塗る

色をつけていきます。中の色はfillで指定します。今はレジェンドが邪魔なので消しておきます。

ggplot() +
  geom_rect(aes(xmin=1, xmax=3, ymin=2, ymax=5, fill='red')) +
  theme(legend.position = "none")

枠線の色を塗る

次にその領域の枠の色も指定します。color=’色’で枠線の色が入ります。

ggplot() +
geom_rect(aes(xmin=1, xmax=3, ymin=2, ymax=5, fill='red'), color='blue') +
  theme(legend.position = "none")

グラフ外まで領域が広がっている時

例えば、その図の外側までずっと続くような感じのグラフにしたい場合、Infを使うとできます。

ggplot()+
  geom_rect(aes(xmin=0, xmax=Inf, ymin=0, ymax=Inf, fill='red'), color='red') +
  scale_x_continuous(limits = c(-2, 6)) + 
  scale_y_continuous(limits = c(-2, 6)) +
  theme(legend.position = "none")

マイナス方向でずっと領域が続く時には-Infを使います。複数の領域を書く方法はまたあとで書きますが、-Infを使うと次のようになります。

透明度はalphaで1より小さく、そして小さい値ほど透明になります。

df <- data.frame(x1=c(-Inf, 0), x2=c(0, Inf), y1=c(-Inf, 0), y2=c(0, Inf), 
                group=c('a', 'b'))
ggplot() +
  geom_rect(data=df, aes(xmin=x1, xmax=x2, ymin=y1, ymax=y2, fill=group), 
           alpha=0.3) +
  scale_x_continuous(limits = c(-6, 6)) + 
  scale_y_continuous(limits = c(-6, 6)) +
  theme_minimal()

複数の領域を書きたい時

この時はデータフレームで領域を指定すると楽になります

複数の四角を作りたい時、データフレームにします。下の式ではx1が左辺、x2が右辺のx座標、y1が下辺、y2が上辺のy座標です。データフレームではさらにグループ分けもしています。

df=data.frame(x1=c(1,2, 3,5), x2=c(5, 8, 5,Inf), 
              y1=c(1, 2, 2, 4), y2=c(2, 3, 6,Inf), 
             group=c('a', 'a','b','b'))
ggplot() +
geom_rect(data=df, aes(xmin=x1, xmax=x2, ymin=y1, ymax=y2, fill=group), 
          color='black', alpha=0.3) +
scale_x_continuous(limits =c(0, 10)) +
scale_y_continuous(limits = c(0, 7))

scale_x_continuousとscale_y_continuousでx軸とy軸のどの数値まで入れるか指定しています。

他のグラフの背景にする

使い所としては他のグラフの背景にして目立たせるということでしょう。ggplot2は記述を＋で重ねて行くだけなので、とても簡単です。単純な直線のグラフの背景として書いてみます。geom_ablineが追加された斜めの線です。

rect <- data.frame(x1=c(-Inf, 0), x2=c(0, Inf), y1=c(-Inf, 0), y2=c(0, Inf), 
                group=c('a', 'b'))
ggplot()+
  geom_rect(data=rect, aes(xmin=x1, xmax=x2, ymin=y1, ymax=y2, fill=group), 
            alpha=0.3) +
  geom_abline(intercept=0, slope = 1) +
  scale_x_continuous(limits = c(-6, 6)) + 
  scale_y_continuous(limits = c(-6, 6)) +
  theme_minimal()

いかがだったでしょうか。わかってしまえば簡単です。ぜひ試してみてください。

【ggplot2】Rで棒グラフを作成する

brainblog — Thu, 01 Jul 2021 04:33:36 +0000

最もよく使われる棒グラフですが、その中でもいくつか使い分けが必要です。そのいくつかの種類の棒グラフの描き方を紹介します。

Contents

データの準備
棒グラフはgeom_bar()を使う
- 色を付ける
グループ化された棒グラフを作図する
エラーバーをつける
- 縦棒グラフ＋エラーバー
- 横棒グラフ＋エラーバー

データの準備

まずはデータの準備をします。デモ用のデータとして最もよく使われるiris(アヤメ)データセットを今回も使います。

# 作業ディレクトリ
setwd("~/Rpractice/")

# tidyverseの起動
library(tidyverse)

#irisデータをaに格納
a <- iris

#データの確認
> head(a)
  Sepal.Length Sepal.Width Petal.Length Petal.Width Species
1          5.1         3.5          1.4         0.2  setosa
2          4.9         3.0          1.4         0.2  setosa
3          4.7         3.2          1.3         0.2  setosa
4          4.6         3.1          1.5         0.2  setosa
5          5.0         3.6          1.4         0.2  setosa
6          5.4         3.9          1.7         0.4  setosa
>

それではこのirisデータを用いていくつかの棒グラフをつくっていきます。

棒グラフはgeom_bar()を使う

それではgeom_bar()を使っていきます。まずはアヤメデータの先頭にあるSepal.Lengthの値を使っていきます。

横軸にアヤメの種類(Species)、縦軸にSepal.Lengthをとります。値をy軸に取るときにはgeom_bar(stat = “identity”)とidentityを指定します。その他にはcountとbinを取ることが出来ます。

ggplot(a, aes(Species, Sepal.Length))+
  geom_bar(stat = "identity")

グラフを作成したら何が縦軸の値をとっているのか、しっかり確認することが必要です。元のデータと対比し、自分の出したい値が作図されているか確認しましょう。上の例では積み上げになっています。

色を付ける

塗りつぶしにはaesの中でfillを指定します。枠の色を指定したいときはcolourで指定します

ggplot(a, aes(Species, Sepal.Length, fill = Species))+
  geom_bar(stat = "identity")

個数を出したいときはgeom_bar(stat = “count”)と指定します。

ggplot(a, aes(x = Species, fill = Species))+
  geom_bar(stat = "count")

グループ化された棒グラフを作図する

まず最初に平均値を計算します。合わせて、SDやSEも計算をしておくとエラーバーを付け加えることができます。

# 縦長のグラフ
a.long <- pivot_longer(a, col = Sepal.Length:Petal.Width,
                         names_to = "iris.attr", 
                       values_to = "measures")

# グループごとの平均値などの計算
a.sum <- group_by(a.long, iris.attr, Species) %>% 
  summarize_all(list(mean = mean,
                     max = max,
                     sd = sd, 
                     se = ~sd/sqrt(length(.))))

積み上げ棒グラフ

geom_barでpositionを指定しなければ積み上げ棒グラフになります。

ggplot(a.sum, aes(x = iris.attr, y = mean,　fill = Species))+
  geom_bar(stat = "identity")

グループ別の棒グラフ

それぞれグループを横に並べた棒グラフはgeom_barでposition = “dodge”と指定します。

ggplot(a.sum, aes(x = iris.attr, y = mean,　fill = Species))+
  geom_bar(stat = "identity", position = "dodge")

幅を調整する

横幅についてはwidthで調整をします。width = 0.5にしてみます。

ggplot(a.long, aes(x = iris.attr, y =measures, fill = Species))+
  geom_bar(stat = "identity", position = "dodge" , width = 0.5)

100％積み上げ棒グラフ

実際の数値ではなく、割合で比較したい時に有用です。

geom_bar(position = “fill”)とすることで100％積み上げ棒グラフを作図できます。

ggplot(a.sum, aes(x = iris.attr, y = mean,　fill = Species))+
  geom_bar(stat = "identity", position = "fill")

横向き棒グラフ

横向きにするときにはcood_flip()を追加します。

ggplot(a.sum, aes(x = iris.attr, y =mean,　fill = Species))+
  geom_bar(stat = "identity", position = "dodge") +
  coord_flip()

凡例の順番を入れ替える

棒グラフの並びと凡例の並びが合わないときにはguides(fill = guide_legend(reverse = TRUE))を指定します。

ggplot(a.sum, aes(x = iris.attr, y =mean,　fill = Species))+
  geom_bar(stat = "identity", position = "dodge") +
  coord_flip()+
  guides(fill = guide_legend(reverse = TRUE))

色を変える

色を変える方法は色々とありますが、scale_fill_brewerで簡単にセット化されたものを使ってみます。

scale_fill_brewer(palette = “Set1”)と指定します。

ggplot(a.sum, aes(x = iris.attr, y = mean, fill = Species))+
  geom_bar(stat = "identity", position = "dodge" , width = 0.9)+
  scale_fill_brewer(palette = "Set1")

さらにscale_fill_brewer(palette = “Set2”)で色を変えてみます。

ggplot(a.sum, aes(x = iris.attr, y = mean, fill = Species))+
  geom_bar(stat = "identity", position = "dodge" , width = 0.9)+
  scale_fill_brewer(palette = "Set2")

エラーバーをつける

最後にエラーバーをつけます。

エラーバーをの付け方については、次のリンクで詳しく説明をしています。

https://brain-storm.space/ggplot2_errorbar/895/

縦棒グラフ＋エラーバー

ggplot(a.sum, aes(x = iris.attr, y = mean, fill = Species))+
  geom_bar(stat = "identity", position = "dodge" , width = 0.9)+
  scale_fill_brewer(palette = "Set1")+
    geom_errorbar(aes(ymin = mean - se, ymax = mean + se),
                  position = position_dodge(0.9), width = .2)

横棒グラフ＋エラーバー

ggplot(a.sum, aes(x = iris.attr, y = mean, fill = Species))+
  geom_bar(stat = "identity", position = "dodge" , width = 0.9)+
  scale_fill_brewer(palette = "Set1")+
  geom_errorbar(aes(ymin = mean - se, ymax = mean + se),
                position = position_dodge(0.9), width = .2)+
  guides(fill = guide_legend(reverse = TRUE))+
  coord_flip()

お役に立ちましたら幸いです。

【R】ggplot2で相関行列ヒートマップを描く

brainblog — Tue, 22 Jun 2021 08:09:26 +0000

相関行列の相関係数を表にすることも多いと思います。しかし視覚に訴求する強力なツールとしてRで図を作成することができます。

それの一つが相関行列です。

Rではいくつものパッケージがあるので、そちらを使うこともできますが、シンプルにggplot2を使った方法を紹介します。

Contents

データの準備をします
- mtcarsとは
- 相関行列の生成
geom_tileを使って相関行列のヒートマップを描く
- geoｍ_tile()を使う
- scale_fill_gradientで色の変更をする
ヒートマップの中に数値を表示する
軸ラベルの順番を変える
- y軸を反転させる
三角行列で表示をしたい時
- 上三角行列
- 下三角行列

データの準備をします

相関行列を使います。デモのためにこの記事ではmtcarsのデータを使ってみます。

mtcarsとは

mtcarsは1974年にMotor Trends US Magazineという雑誌から抽出された、各車種における燃費などの車性能が入ったデータです。

mpg:燃費、cyl:気筒数、disp:排気量などのデータが入っています。

> mtcars
                     mpg cyl  disp  hp drat    wt  qsec vs am gear carb
Mazda RX4           21.0   6 160.0 110 3.90 2.620 16.46  0  1    4    4
Mazda RX4 Wag       21.0   6 160.0 110 3.90 2.875 17.02  0  1    4    4
Datsun 710          22.8   4 108.0  93 3.85 2.320 18.61  1  1    4    1
Hornet 4 Drive      21.4   6 258.0 110 3.08 3.215 19.44  1  0    3    1
Hornet Sportabout   18.7   8 360.0 175 3.15 3.440 17.02  0  0    3    2
Valiant             18.1   6 225.0 105 2.76 3.460 20.22  1  0    3    1
Duster 360          14.3   8 360.0 245 3.21 3.570 15.84  0  0    3    4
Merc 240D           24.4   4 146.7  62 3.69 3.190 20.00  1  0    4    2
Merc 230            22.8   4 140.8  95 3.92 3.150 22.90  1  0    4    2
Merc 280            19.2   6 167.6 123 3.92 3.440 18.30  1  0    4    4
Merc 280C           17.8   6 167.6 123 3.92 3.440 18.90  1  0    4    4
Merc 450SE          16.4   8 275.8 180 3.07 4.070 17.40  0  0    3    3
Merc 450SL          17.3   8 275.8 180 3.07 3.730 17.60  0  0    3    3
Merc 450SLC         15.2   8 275.8 180 3.07 3.780 18.00  0  0    3    3
Cadillac Fleetwood  10.4   8 472.0 205 2.93 5.250 17.98  0  0    3    4
Lincoln Continental 10.4   8 460.0 215 3.00 5.424 17.82  0  0    3    4
Chrysler Imperial   14.7   8 440.0 230 3.23 5.345 17.42  0  0    3    4
Fiat 128            32.4   4  78.7  66 4.08 2.200 19.47  1  1    4    1
Honda Civic         30.4   4  75.7  52 4.93 1.615 18.52  1  1    4    2
Toyota Corolla      33.9   4  71.1  65 4.22 1.835 19.90  1  1    4    1
Toyota Corona       21.5   4 120.1  97 3.70 2.465 20.01  1  0    3    1
Dodge Challenger    15.5   8 318.0 150 2.76 3.520 16.87  0  0    3    2
AMC Javelin         15.2   8 304.0 150 3.15 3.435 17.30  0  0    3    2
Camaro Z28          13.3   8 350.0 245 3.73 3.840 15.41  0  0    3    4
Pontiac Firebird    19.2   8 400.0 175 3.08 3.845 17.05  0  0    3    2
Fiat X1-9           27.3   4  79.0  66 4.08 1.935 18.90  1  1    4    1
Porsche 914-2       26.0   4 120.3  91 4.43 2.140 16.70  0  1    5    2
Lotus Europa        30.4   4  95.1 113 3.77 1.513 16.90  1  1    5    2
Ford Pantera L      15.8   8 351.0 264 4.22 3.170 14.50  0  1    5    4
Ferrari Dino        19.7   6 145.0 175 3.62 2.770 15.50  0  1    5    6
Maserati Bora       15.0   8 301.0 335 3.54 3.570 14.60  0  1    5    8
Volvo 142E          21.4   4 121.0 109 4.11 2.780 18.60  1  1    4    2
>

相関行列の生成

相関行列を生成してみます。

mtcarsのデータからSpearmanの順位相関係数をみてみます。またggplot2で使いやすいようにデータを変形していきます。相関はcorで計算をしており、その中のmethodでspearmanと指定しています。

縦のデータへの変換はpivot_longerを用いています。pivot_longerは次の記事に詳しい説明があります。

https://brain-storm.space/tidyr_gather_pivot_longer/784/

# tidyverseの起動
library(tidyverse)

#mtcars
a <- mtcars

#相関行列(Spearman)
cormat <- cor(a, method = "spearman")

#データフレーム化
df <- data.frame(cormat) %>% 

#列名の抽出
tibble::rownames_to_column(var = "Y")

#縦のデータへ
df.long <- pivot_longer(df, 
                            cols = mpg:carb,
                            names_to = "X", 
                            values_to = "rho")

#データの確認
> head(df.long)
# A tibble: 6 x 3
  Y     X        rho
     
1 mpg   mpg    1    
2 mpg   cyl   -0.911
3 mpg   disp  -0.909
4 mpg   hp    -0.895
5 mpg   drat   0.651
6 mpg   wt    -0.886
>

以上のようなデータになりました。自分のデータで相関行列があればそれを縦のデータに変換して利用することができます。

geom_tileを使って相関行列のヒートマップを描く

geom_tile()を使うことで視覚化することができます。

次のようになります。

geoｍ_tile()を使う

ggplot(df.long, aes(x = X, y = Y, fill = rho)) + 
  geom_tile()

scale_fill_gradientで色の変更をする

色を変更してみます。よく見るのは値の高いところが赤、低いところが青ですよね？

scale_fill_gradient2を使ってみます。

ggplot(df.long, aes(x = X, y = Y, fill = rho)) + 
  geom_tile(color = "white") + 
  scale_fill_gradient2(low = "blue", high = "red", mid = "white", 
                       midpoint = 0, limit = c(-1, 1), 
                       name = "Rho", space = "Lab")

これで出来ました！

ところでここに数値を入れたい人もいるかもしれません。

ヒートマップの中に数値を表示する

数値を入れるにはgeom_text()で出来ます！

df.long %>% 
  ggplot(aes(x = X, y = Y, fill = rho)) + 
  geom_tile(color = "white") + 
  scale_fill_gradient2(low = "blue", high = "red", mid = "white", 
                       midpoint = 0, limit = c(-1, 1), name = "Rho", space = "Lab")+
  geom_text(aes(label = sprintf("%0.2f", rho)),
            color = "black", size = 2.5)

軸ラベルの順番を変える

軸ラベルがアルファベット順になっていることに気がついたでしょうか。

これが不便に感じるときがあります。そういうときはfactor(data, levels = )で軸ラベルの順番を指定すれば解決です。

lbs <- df[,1]
df.long$X <- factor(df.long$X, levels = lbs)
df.long$Y <- factor(df.long$Y, levels = lbs)

df.long %>% 
  ggplot(aes(x = X, y = Y, fill = rho)) + 
  geom_tile(color = "white") + 
  scale_fill_gradient2(low = "blue", high = "red", mid = "white", 
                       midpoint = 0, limit = c(-1, 1), name = "Rho", space = "Lab")+
  geom_text(aes(label = sprintf("%0.2f", rho)),
            color = "black", size = 2.5)

順番はdata <- factor(data, levels = c(“mpg”, “cyl”, “disp”, ….))というように一つ一つ指定することもできます。

y軸を反転させる

上のグラフは左下からラベルが並んでいます。y軸は上からはじまって欲しい場合があります。

そういうときはy軸の並びを反転させます。

scale_y_discrete(limits = rev(levels( )))で反転できます。

df.long %>% 
  ggplot(aes(x = X, y = Y, fill = rho)) + 
  geom_tile(color = "white") + 
  scale_fill_gradient2(low = "blue", high = "red", mid = "white", 
                       midpoint = 0, limit = c(-1, 1), name = "Rho", space = "Lab")+
  geom_text(aes(label = sprintf("%0.2f", rho)),
            color = "black", size = 2.5)+
  scale_y_discrete(limits = rev(levels(df.long$Y)))

三角行列で表示をしたい時

上の図では対角線で対称になっていますが、余計な情報だと思うかもしれません。それであれば三角行列にしてみます。

ただし、ggplotで三角行列にするというより、データ自体を三角行列になるようにしていきます。

上三角行列

上三角行列にするには下の部分をNAにして、その部分を空欄にします。

# 上三角行列にするためのfunction
# 下にはNAを代入する
get_upper_tri <- function(cormat){
  cormat[lower.tri(cormat)]<- NA
  return(cormat)
}

# 実際のデータに適応
upper_tri <- get_upper_tri(cormat)

# データの変形
df_upper_tri <- data.frame(upper_tri) %>% 
  tibble::rownames_to_column(var = "Y")

df_upper_tri.long <- pivot_longer(df_upper_tri, 
                        cols = mpg:carb,
                        names_to = "X", 
                        values_to = "rho")

# ラベルの並び替え
lbs2 <- df_upper_tri[,1]
df_upper_tri.long$X <- factor(df_upper_tri.long$X, levels = lbs2)
df_upper_tri.long$Y <- factor(df_upper_tri.long$Y, levels = lbs2)

# 作図
df_upper_tri.long %>% 
  ggplot(aes(x = X, y = Y, fill = rho)) + 
  geom_tile(color = "white") + 
  scale_fill_gradient2(low = "blue", high = "red", mid = "white", 
                       midpoint = 0, limit = c(-1, 1), 
                       name = "Rho", space = "Lab",
                       na.value = "white")+
  geom_text(aes(label = ifelse(is.na(rho) | rho %in% "",NA, sprintf("%0.2f", rho))),
            color = "black", size = 4.5)+
  scale_y_discrete(limits = rev(levels(df_upper_tri.long$Y)))+
  theme_classic()

下三角行列

下三角も同様に出来ます。

# 下三角行列にするためのfunction
# 上にはNAを代入する
get_lower_tri <- function(cormat){
  cormat[upper.tri(cormat)]<- NA
  return(cormat)
}

# 実際のデータに適応
lower_tri <- get_lower_tri(cormat)

# データの変形
df_lower_tri <- data.frame(lower_tri) %>% 
  tibble::rownames_to_column(var = "Y")

df_lower_tri.long <- pivot_longer(df_lower_tri, 
                                  cols = mpg:carb,
                                  names_to = "X", 
                                  values_to = "rho")

# ラベルの並び替え
lbs3 <- df_lower_tri[,1]
df_lower_tri.long$X <- factor(df_lower_tri.long$X, levels = lbs3)
df_lower_tri.long$Y <- factor(df_lower_tri.long$Y, levels = lbs3)

# 作図
df_lower_tri.long %>% 
  ggplot(aes(x = X, y = Y, fill = rho)) + 
  geom_tile(color = "white") + 
  scale_fill_gradient2(low = "blue", high = "red", mid = "white", 
                       midpoint = 0, limit = c(-1, 1), 
                       name = "Rho", space = "Lab",
                       na.value = "white")+
  geom_text(aes(label = ifelse(is.na(rho) | rho %in% "",NA, sprintf("%0.2f", rho))),
            color = "black", size = 2.5)+
  scale_y_discrete(limits = rev(levels(df_lower_tri.long$Y)))+
  theme_classic()

ifelseの使い方が初心者には難しいかもしれませんが、慣れれば特定の閾値以下は表示しない、などのテクニックなども使えます。

便利なグラフなので、覚えて損はありません。お役に立ちましたら幸いです。

【R初心者向け】ggplot2でバイオリンプロット(violin plot)を描く

brainblog — Mon, 31 May 2021 06:50:12 +0000

データを可視化するときにいくつかの方法があります。その一つがバイオリンプロットです。

バイオリンプロットは複数のデータを比較する時に便利な方法です。バイオリンプロットを書いているだけで、「この人出来る」というような印象があるかもしれません。是非チャンスがあれば使ってみてください。

それでは書き方にいきます。ついでにドットプロットでデータの分布を確認してみます。

Contents

データ生成と分布の確認
バイオリンプロットを描く
バイオリンに色をつける
- 枠に色を付ける
- バイオリンの中に色を塗る
平均値や中央値を入れる
- 平均値をマークする
- 中央値をマークする
平滑化の度合いを変更する
ボックスプロットと重ね合わせる

データ生成と分布の確認

データを作ってみます。自分のデータがある場合にはそちらを使ってみてください。

#作業ディレクトリのセットとtidyverse、ggbeeswarmの起動
#ggbeeswarmはドットプロットを描くためです
setwd("~/Rpractice/")
library(tidyverse)
library(ggbeeswarm)

#データを生成します
dat <- list(
  X <- rnorm(100, 5, 10),
  Y <- rnorm(100, 20, 10),
  Z <- rnorm(100, 15, 15)
)

#データを整えます
#単に使いやすいデータの形にするだけです
df <- data.frame(matrix(unlist(dat), nrow=100))
colnames(df) <- c("A","B","C")
df.long <- pivot_longer(df, cols = A:C, names_to = "Categories", values_to = "Values")

#ドットプロットです
ggplot(df.long, aes(x = Categories, y = Values))+
  geom_beeswarm(aes(color = Categories),
                size = 2,
                cex = 2,
                alpha = .8)+
  theme_classic()+
  theme(legend.position = "none")

ドットプロット、ggbeeswarmについてはこちらからどうぞ。またデータをpivot_longerで縦長に変換しています。

https://brain-storm.space/ggbeeswarm/820/ https://brain-storm.space/tidyr_gather_pivot_longer/784/

それでは、バイオリンプロットを描いていきましょう！

バイオリンプロットを描く

ggplot2でバイオリンプロットはgeom_violin()を使います！

theme_classic()で背景を白くします。

ggplot(df.long, aes(x = Categories, y = Values))+
  geom_violin()+
  theme_classic()

これがバイオリンプロットの基本形になります。バイオリンプロットの両端が切れていますね。この上にドットプロットを重ねてみると、その答えがわかります。

ggplot(df.long, aes(x = Categories, y = Values))+
  geom_violin()+
  geom_beeswarm(aes(color = Categories),
                size = 2,
                cex = 2,
                alpha = .8)+
  theme_classic()

このバイオリンプロットの両端は最大値と最小値でカットされていることがわかります。これをカットしないのはgeom_violin(trim = FALSE)で指定します。

ドットプロットも重ねてみます。

ggplot(df.long, aes(x = Categories, y = Values))+
  geom_violin(trim = FALSE)+
  geom_beeswarm(aes(color = Categories),
                size = 2,
                cex = 2,
                alpha = .8)+
  theme_classic()

このように点のないところまで上下に伸びますので要注意です。

バイオリンに色をつける

バイオリンに色を付ける場合、枠に色をつけるか、中を塗りつぶすかで指定の仕方を変えます。

枠に色を付ける

枠に色をつけるのはcolorで指定します。

ggplot(df.long, aes(x = Categories, y = Values, color = Categories))+
  geom_violin()+
  theme_classic()

重ね合わせの時にバイオリンだけ色をつけたい、などでcolorをgeom_violinのaesで指定もできます。

ggplot(df.long, aes(x = Categories, y = Values))+
  geom_violin(aes(color = Categories))+
  theme_classic()

出てくる図は上と同じです。慣れた人には当たり前なんですが、初心者のうちは意外に迷うものです。

バイオリンの中に色を塗る

　バイオリンの中を塗りつぶす時には、fillを使います。

ggplot(df.long, aes(x = Categories, y = Values, fill = Categories))+
  geom_violin()+
  theme_classic()

色を変えるときにはいくつか方法があります。今回はscale_fill_brewerで指定をしてみます。

ggplot(df.long, aes(x = Categories, y = Values, fill = Categories))+
  geom_violin()+
  scale_fill_brewer(palette = "Set2")+
  theme_classic()

平均値や中央値を入れる

バイオリンプロットに平均値や中央値を入れるのは、stat_summaryを使います。

平均値をマークする

平均値をマークしてみましょう。

ggplot(df.long, aes(x = Categories, y = Values, color = Categories))+
  geom_violin()+
  stat_summary(fun = mean, geom = "point", 
               shape =16, size = 2, color = "red")+
  theme_classic()

平均値を赤い点にしました。stat_summary()の中のshapeとcolorで形と色を指定をしています。

stat_summaryのshapeはpchと同じです。

数値に対応する図形はこのようになっています！

中央値をマークする

今度は中央値です。

ボックスプロットを重ねれば必ずしも必要ないかもしれませんが、先程と同じようにstat_summaryを使うことができます。shapeで十字マーク(shape = 3)にしています。

ggplot(df.long, aes(x = Categories, y = Values, color = Categories))+
  geom_violin()+
  stat_summary(fun = median, geom = "point", 
               shape = 3, size = 2, color = "red")+
  theme_classic()

平滑化の度合いを変更する

平滑化の度合いを変更したいときにはgeom_violinの中でadjustを指定します。

デフォルトはadjust = 1です。

まずはadjustの値を小さくしてみます。adjustを0.2にしてみましょう！

ggplot(df.long, aes(x = Categories, y = Values, fill = Categories))+
  geom_violin(adjust = .2)+
  theme_classic()

細かい並々がたくさんできました。反対にadjustを大きくしてみます。

ggplot(df.long, aes(x = Categories, y = Values, fill = Categories))+
  geom_violin(adjust = 2)+
  theme_classic()

今度は赤の群にあったくびれがなくなってしまいました。

ボックスプロットと重ね合わせる

ボックスプロットとの重ね合わせはよく使われます。

ボックスプロットを重ねることでその威力を発揮します。

ggplot(df.long, aes(x = Categories, y = Values, fill = Categories))+
  geom_violin()+
  geom_boxplot(width = .1, fill = "white")+
  theme_classic()

外れ値を表示しないようにするにはoutlier.color = NAを指定します。

ggplot(df.long, aes(x = Categories, y = Values, fill = Categories))+
  geom_violin()+
  geom_boxplot(width = .1, fill = "white", outlier.color = NA)+
  theme_classic()

ボックスプロットを黒く塗って、中央値を白丸にする方法です。白黒印刷物に投稿する場合は役に立ちます。

ggplot(df.long, aes(x = Categories, y = Values))+
  geom_violin()+
  geom_boxplot(width = .1, fill = "black", outlier.color = NA) +
  stat_summary(fun = median, geom = "point", fill = "white", shape = 21, size = 3) +
  theme_classic()

これでバイオリンプロットも描けると思います。お役に立ちましたら幸いです。

【Rテクニック】facet_gripとfacet_wrapの使い方とタイトルラベルを変えたい！

brainblog — Sat, 29 May 2021 03:29:46 +0000

facet_gripとfacet_wrapは複数グラフを一度に作成できるのでとても便利です！

またビジュアルとしてとても有用で、魅力的なグラフになります。

タイトルラベルを変えるのは少しテクニックがありますので、それを含めて紹介します。

それではいってみましょう！

Contents

facet_gridとfacet_wrapの使い方
ファセットラベルの変更の仕方
- labellerの設定をする
- タイトル、サブタイトルをつける！

facet_gridとfacet_wrapの使い方

facet_gridとfacet_wrapはどちらもデータをグループごとに並べるものになります。いつも便利なiris dataを用いていきます。

facet_gridを使ってみる！

まずはデータの準備をします。

tidyverseを起動しています。

#作業ディレクトリのセット
setwd("~/Rpractice/")

#libraryの起動
library(tidyverse)

#irisデータ
a <- iris

> head(a)
  Sepal.Length Sepal.Width Petal.Length Petal.Width Species
1          5.1         3.5          1.4         0.2  setosa
2          4.9         3.0          1.4         0.2  setosa
3          4.7         3.2          1.3         0.2  setosa
4          4.6         3.1          1.5         0.2  setosa
5          5.0         3.6          1.4         0.2  setosa
6          5.4         3.9          1.7         0.4  setosa
>

まずはfacet_gridの使い方を見てみましょう！

まずは水平方向です。

facet_gridのカッコの中を(. ~ Species)で指定しています。

ggplot(a, aes(x = Sepal.Length, 
              y = Sepal.Width, 
              color = Species)) +
  geom_point(size = 3, alpha = .5)+
  geom_smooth(method = lm, se = TRUE)+
  facet_grid(. ~ Species)

次は垂直方向です！

これは反対にfacet_grid(Species ~ .)で指定しています。

ggplot(a, aes(x = Sepal.Length, y = Sepal.Width, color = Species)) +
  geom_point(size = 3, alpha = .5)+
  geom_smooth(method = lm, se = TRUE)+
  facet_grid(Species ~ .)

このような感じになりました。

今回は１行または１列で表示しましたが、それぞれ垂直方向、水平方向を指定して複数の行列で表示することもできます。

facet_wrapを使ってみる！

ggplot(a, aes(x = Sepal.Length, y = Sepal.Width, color = Species)) +
  geom_point(size = 3, alpha = .5)+
  geom_smooth(method = lm, se = TRUE)+
  facet_wrap(~ Species)  #ドットは不要です

２行にしてみましょう。

facet_wrapのカッコの中のncolまたはnrowで指定できます。

ggplot(a, aes(x = Sepal.Length, y = Sepal.Width, color = Species)) +
  geom_point(size = 3, alpha = .5)+
  geom_smooth(method = lm, se = TRUE)+
  facet_wrap(~ Species, nrow = 2)

スケールを別々にする方法

スケールが合わさっていると極端に偏ったグラフが出来てしまうことがあります。

そのときはscales = “free”を設定します。

ただし、逆に見にくくなることもあるので注意をしてください。

ファセットラベルの変更の仕方

これを変更する一つの方法はデータを変更してしまうことです。

もう一つの方法として、labellerを指定します！

早速やり方を見ていきましょう。

labellerの設定をする

facet_gridまたはfacet_wrapでlabeller = labellar()で指定をすることができます。

下記の方法は少しテクニックが必要ですが、真似をするだけで変更できます。

Att.labs <- c("Att:Setosa", 
              "Att:Versicolor", 
              "Att:Virginica")
names(Att.labs) <- c("setosa",
                      "versicolor",
                      "virginica")

ggplot(a, aes(x = Sepal.Length, y = Sepal.Width, color = Species)) +
  geom_point(size = 3, alpha = .5)+
  geom_smooth(method = lm, se = FALSE)+
  facet_wrap(~ Species,  scales = "free", 
             labeller = labeller(Species = Att.labs))

またタイトルやサブタイトルもつけることができます

タイトル、サブタイトルをつける！

タイトルのラベルはlabs(title = )でしていできます。

また同様にlabsでサブタイトルやx軸、y軸のラベルを指定することができます。

ggplot(a, aes(x = Sepal.Length, y = Sepal.Width, color = Species)) +
  geom_point(size = 3, alpha = .5)+
  geom_smooth(method = lm, se = FALSE)+
  facet_wrap(~ Species,  scales = "free", 
             labeller = labeller(Species = Att.labs)
  )+
  labs(title = "Iris data analysis", subtitle = "Sepal measurement",
       x = "Length",
       y = "Width")

フォントサイズを整えます。またtheme()でその大きさをそれぞれ指定できます。

レジェンドもあんまり役に立っていないので、なくしちゃいましょう。その時はlegend.position = “none”とすることで表示をなくします。

ggplot(a, aes(x = Sepal.Length, y = Sepal.Width, 
              color = Species)) +
  geom_point(size = 3, alpha = .5)+
  geom_smooth(method = lm, se = FALSE)+
  facet_wrap(. ~ Species,  scales = "free", 
             labeller = labeller(Species = Att.labs))+
  labs(title = "Iris data analysis", 
       subtitle = "Sepal measurement",
       x = "Length", y = "Width")+
  theme(axis.text.x = element_text(vjust = 1, size = 16, 
                                   hjust = 1, face = "bold"),
        axis.text.y = element_text(hjust = 1, size = 16),
        axis.title.x = element_text(size = 20),
        axis.title.y = element_text(size = 20),
        legend.title = element_text(size=16),
        legend.text = element_text(size=16),
        strip.text.x = element_text(size = 20),
        plot.title = element_text(size=22),
        plot.subtitle = element_text(size=20),
        legend.position = "none")

スケールメモリを統一したほうが実際は見やすいことが多いです。

最後にスケールメモリを統一してみましょう。

theme_classic()で背景を白くしてみます。

ggplot(a, aes(x = Sepal.Length, y = Sepal.Width, 
              color = Species)) +
  geom_point(size = 3, alpha = .5)+
  geom_smooth(method = lm, se = FALSE)+
  facet_wrap(. ~ Species,
             labeller = labeller(Species = Att.labs))+
  labs(title = "Iris data analysis", 
       subtitle = "Sepal measurement",
       x = "Length", y = "Width")+
  theme_classic() +
  theme(axis.text.x = element_text(vjust = 1, size = 16, 
                                   hjust = 1, face = "bold"),
        axis.text.y = element_text(hjust = 1, size = 16),
        axis.title.x = element_text(size = 20),
        axis.title.y = element_text(size = 20),
        legend.title = element_text(size=16),
        legend.text = element_text(size=16),
        strip.text.x = element_text(size = 20),
        plot.title = element_text(size=22),
        plot.subtitle = element_text(size=20),
        legend.position = "none")

それぞれを比較するにはとても見やすくなったと思います。

facet_gridとfacet_wrapはとても便利ですよね！

お役に立ちましたら幸いです。

【R初心者向け】ggplot2でもbeeswarmで作図したい！ggbeeswarmを使った方法

brainblog — Mon, 24 May 2021 23:00:00 +0000

ggplot2でドットプロットを紹介しましたが、ggplot2でbeeswarmに出来ないでしょうか、という質問をいただきました。

今回はggplot2でbeeswarmを描く方法を紹介します。ggbeeswarmというパッケージを使えばできます。

過去にggplot2を使ったドットプロットの書き方と、ggplot2を使わずbeeswarmパッケージを使った方法を紹介しています。それについてはこちらを参考にしてください。

https://brain-storm.space/rggplot2_boxplot_dotplot/620/ https://brain-storm.space/beeswarm_r/698/

それでは早速ggplot2でbeeswarmを書いていきましょう。

Contents

データ準備とパッケージのインストール
ggbeeswarmのgeom_beeswarmを使う
methodによる違い

データ準備とパッケージのインストール

まずはtidyverseとggbeeswarmをいれていきます。ggplot2はtidyverseの中に入っています。

まずは作業ディレクトリの設定

setwd("~/Rpractice/")

ggbeeswarmをCRANからインストールをします。

tidyverseやgglot2はすでにインストールされている前提ですすめていきますね。

install.packages("ggbeeswarm")

それではパッケージを起動します

library(tidyverse)
library(ggbeeswarm)

アヤメデータを使います

アヤメはアヤメ科の植物ですが、英語でiris(アイリス)といいます

ギリシャ語の虹(イリス)に由来するそうです

#irisデータをaに格納
a <- iris
#ボックスプロットを描画
ggplot(a, aes(x = Species, y = Sepal.Length, fill = Species))+
  geom_boxplot(fill = "white")+
  theme_classic()

ggbeeswarmのgeom_beeswarmを使う

ggbeeswarmを使っていないドットプロット

以前の記事でggplotを用いて書いた図です。ggplotでのボックスプロットはこちらを参考にしてください。

https://brain-storm.space/rggplot2_boxplot_dotplot/620/

ggplot(a, aes(x = Species, y = Sepal.Length, fill = Species))+
  geom_boxplot(fill = "white")+
  geom_dotplot(binaxis = "y", binwidth = 0.1, stackdir = "center")+
  theme_classic()+
  theme(axis.text.x = element_text(angle = 50, vjust = 1, size = 16, hjust = 1, face = "bold"),
        axis.text.y = element_text(hjust = 1, size = 16),
        axis.title.x = element_text(size = 20),
        axis.title.y = element_text(size = 20),
        legend.title = element_text(size =16),
        legend.text = element_text(size =16))

それではggbeeswarmを使っていきたいと思います。

geom_beeswarmを使う

geom_beeswarmを使ってみましょう。

ggplot(a, aes(x=Species, y=Sepal.Length))+
  geom_boxplot(fill="white")+
  geom_beeswarm()+
  theme_classic()

いまいちですね。色をつけてみます。

ggplot(a, aes(x=Species, y=Sepal.Length))+
  geom_boxplot(fill="white")+
  geom_beeswarm(aes(color = Species))+
  theme_classic()

点が小さいし、点どうしが重なっています！これでは使い物になりませんので、大きくして、もう少し点をちらしてみます。

ggplot(a, aes(x = Species, y = Sepal.Length))+
  geom_boxplot(fill = "white")+
  geom_beeswarm(aes(color = Species), size = 3, cex = 3)+
  theme_classic()

まあまあいいですが、このグラフではbeeswarmの良さが見えてこないかもしれません。下でもう少し別のセットを使って示してみます。

geom_quasirandomを使う

今度はggbeeswarmパッケージに入っているgeom_quasirandomを使ってみます

ggplot(a, aes(x = Species, y = Sepal.Length))+
  geom_boxplot(fill = "white")+
  geom_quasirandom()+
  theme_classic()

結構点がバラバラに表示をされました

色と点の大きさを変えてみましょう

ggplot(a, aes(x = Species, y = Sepal.Length))+
  geom_boxplot(fill = "white")+
  geom_quasirandom(aes(color = Species), size = 3)+
  theme_classic()

methodによる違い

点の数は多いほうが、ggbeeswarmの特徴がよくわかります。rnormを使って、さらに点の数が多い分布を作ってみます。

dat <- list(
  X <- rnorm(150, 10, 10),
  Y <- rnorm(150, 25, 10),
  Z <- rnorm(150, 20, 15)
)
df <- data.frame(matrix(unlist(dat), nrow=150))
colnames(df) <- c("A","B","C")
df.long <- pivot_longer(df, cols = A:C, names_to = "Categories", values_to = "Values")

ではいってみましょう。pivot_longerの使い方はこちらから。

https://brain-storm.space/tidyr_gather_pivot_longer/784/

geom_beeswarm

ggplot(df.long, aes(x=Categories, y = Values))+
  geom_boxplot(fill="white")+
  geom_beeswarm(aes(color = Categories),
                size = 3,
                cex = 2,
                alpha =.5)+
  theme_classic()

重なりを表示するためにalphaを指定して半透明にしています。

beeswarmの良さがみえてきました。十分発表や論文に使えるレベルだと思います。

quasirandom

geom_quasirandomでいくつかのmethodを指定することができます。

まずはデフォルトです。

ggplot(df.long, aes(x=Categories, y = Values))+
  geom_boxplot(fill="white")+
  geom_quasirandom(aes(color = Categories),
                   size = 3, 
                   alpha =.5)+
  theme_classic()

次は他のいくつかの方法を示していきます。

pseudorandom

ggplot(df.long, aes(x=Categories, y = Values))+
  geom_boxplot(fill="white")+
  geom_quasirandom(aes(color = Categories),
                   method ='pseudorandom', 
                   size = 3, 
                   alpha =.5)+
  theme_classic()

smiley

わかりやすい名前ですよね。

ggplot(df.long, aes(x=Categories, y = Values))+
  geom_boxplot(fill="white")+
  geom_quasirandom(aes(color = Categories),
                   method ='smiley', 
                   size = 3, 
                   alpha =.5)+
  theme_classic()

口角が上がっているようにみえます。

frowney

smileyとは反対に口角が下がっているような図になります。

ggplot(df.long, aes(x=Categories, y = Values))+
  geom_boxplot(fill="white")+
  geom_quasirandom(aes(color = Categories),
                   method ='frowney', 
                   size = 3, 
                   alpha =.5)+
  theme_classic()

おもしろい名前ですよね。

tukey

turkeyという方法です。

ggplot(df.long, aes(x=Categories, y = Values))+
  geom_boxplot(fill="white")+
  geom_quasirandom(aes(color = Categories),
                   method ='tukey', 
                   size = 3, 
                   alpha =.5)+
  theme_classic()

tukeyDense

最後にtukeyDenseを紹介します。

ggplot(df.long, aes(x=Categories, y = Values))+
  geom_boxplot(fill="white")+
  geom_quasirandom(aes(color = Categories),
                   method ='tukeyDense', 
                   size = 3, 
                   alpha =.5)+
  theme_classic()

ggplotをベースにしたggbeeswarmはキレイですね。お役に立ちましたら幸いです。