Rで散布図を描く（ggplotなし）

 以下のデータを視覚化する。これは、相関分析で使ったデータと同じ。今回のデータは計画行列ではない。架空のデータで、a列が握力、ｂ列が垂直飛びの結果だとする。

a	b
30.1	56.7
27.5	48.3
32.8	57.1
29.6	52.2
30.6	53.4
32.1	57.7
23.3	49.2
29.4	56.9
28.2	56.6
28.8	53.5
25.7	54.2
30.2	51
29.3	52.8
28	53.9
25.8	51.3
30.9	54.9
32.5	57.4
32.2	58.6
24.4	53.6
30.7	55.1

　このデータのa列をx軸とb列をy軸に割り当てて散布図を描いてみる。

　なお、学会発表や論文のために作図する場合は、より自由度が高く、美しい作図が可能なggplotパッケージを用いた方が良い。

　data1としてデータを読みこんでいる。

　コードは以下の通り。plotが関数の名前で () 内に必要な情報を入力していく。最初にx軸(data1$a)、そのあとにy軸(data1$b)を指定する。

plot(data1$a, data1$b)

　x軸とy軸のラベルを任意に変更する場合にはxlab=、ylab=で指定してやればよい。

plot(data1$a, data1$b, xlab = "a", ylab = "b")

　ここからは、ちょっと難しいので、初級っぽくないので流しても良い。

　複数の散布図を重ねることもできる。以下のように2つのグループの握力と垂直飛びのデータがあったとして、それを散布図で重ねて表示したいとする。

a	b	group
30.1	56.7	group_1
27.5	48.3	group_1
32.8	57.1	group_1
29.6	52.2	group_1
30.6	53.4	group_1
32.1	57.7	group_1
23.3	49.2	group_1
29.4	56.9	group_1
28.2	56.6	group_1
28.8	53.5	group_1
25.7	54.2	group_1
30.2	51	group_1
29.3	52.8	group_1
28	53.9	group_1
25.8	51.3	group_1
30.9	54.9	group_1
32.5	57.4	group_1
32.2	58.6	group_1
24.4	53.6	group_1
30.7	55.1	group_1
37.2	62.2	group_2
31	49.7	group_2
35.8	58.4	group_2
30	54.3	group_2
32.2	53	group_2
32.3	59.7	group_2
23.9	52.5	group_2
33.2	58.5	group_2
31.4	60.1	group_2
29.3	53.3	group_2
28	53.7	group_2
34.1	50.8	group_2
29.6	56.4	group_2
30.6	56.4	group_2
28.5	54.3	group_2
34.6	57.3	group_2
35.9	56.6	group_2
33.8	60.4	group_2
26	53.6	group_2
31.3	55.6	group_2

　コードは以下の通り。上のデータをdata2として読み込んでいる。ドットに色を付けて区別させたいので、色がわかりやすいドットの形に変更する。pch=で変更する。19番で塗りつぶしのドットを出力できる。また、col=で色分けを指定する。factor()として色分けしたい列の名前を入れる。

plot(data2$a, data2$b, xlab = "a", ylab = "b", pch = 19, col = factor(data2$group))

ただし、これだとどの色がどのグループかわからない。そこで、plot関数を実行した後に、legends関数を実行する。コードと図は以下の通り。 () 内は最初に凡例を表示する場所、凡例として表示するデータ、凡例に表示するドットの形、凡例の色、を指定していく。表示場所は左上にしたいので"topleft"とする。凡例を表示するためにはlegend=levels(factor())と入力し、factorの () 内に列名data2$groupを指定する。pch=は19番としてドットの形を指定する。col=factor()で色を指定する列data2$groupを指定する（ただし、factorをlevelsで囲む）。

plot(data2$a, data2$b, xlab = "a", ylab = "b", pch = 19, col = factor(data2$group))
legend("topleft", legend = levels(factor(data2$group)), pch = 19, col = factor(levels(factor(data2$group))))