初心者用 R 言語講座演習ノート

 (x <- c(2,4,6,8,10))    # ベクトルの作成
 length(x)               # ベクトルの次元数、すなわちベクトル中の数値の個数
 x[3]                    # 3 番目の数値を抽出する
 x[c(2,4)]               # 2番目と 4 番目の数値を抽出
 x[2:4]                  # 2 番目から 4 番目までの数値を抽出

[1]  2  4  6  8 10
[1] 5
[1] 6
[1] 4 8
[1] 4 6 8

ベクトルを集めて行列を作成する

(x <- rbind(c(1,2,3), c(4,5,6)))
(x <- matrix(c(1, 4, 2, 5, 3, 6), nrow=2, ncol=3)) # nrow:行数， ncol:列数

     [,1] [,2] [,3]
[1,]    1    2    3
[2,]    4    5    6
     [,1] [,2] [,3]
[1,]    1    2    3
[2,]    4    5    6

x+1
x*2

y <- rbind(c(2, 4, 6), c(8, 10, 12))

x + y 
x * y

x[2,3]
x[2,]
x[,3]
t(x) # 転置

	 [,1] [,2] [,3]
    [1,]    2    3    4
    [2,]    5    6    7
	 [,1] [,2] [,3]
    [1,]    2    4    6
    [2,]    8   10   12
	 [,1] [,2] [,3]
    [1,]    3    6    9
    [2,]   12   15   18
	 [,1] [,2] [,3]
    [1,]    2    8   18
    [2,]   32   50   72
    [1] 6
    [1] 4 5 6
    [1] 3 6
	 [,1] [,2]
    [1,]    1    4
    [2,]    2    5
    [3,]    3    6

apply(x, 2, mean) # 行和，xの2次元目の各項の平均
apply(x, 1, mean) # 列和，xの1次元目の各項の平均

[1] 2.5 3.5 4.5
[1] 2 5

行列は 2 次元の数値の羅列ですが、 R では array(ベクトル、各次元の要素の個数)を使っ てさらに n 次元に拡張された配列を扱うことができます。

(x <- array(1:24, c(3,4,2))) # 3×4×2 の大きさの3次元配列を作成

    , , 1

	 [,1] [,2] [,3] [,4]
    [1,]    1    4    7   10
    [2,]    2    5    8   11
    [3,]    3    6    9   12

    , , 2

	 [,1] [,2] [,3] [,4]
    [1,]   13   16   19   22
    [2,]   14   17   20   23
    [3,]   15   18   21   24

課題5-1:

課題5-2: 上記結果に対して行、列それぞれの平均を求めましょう。

リストは様々な型のデータをひとまとめにして管理できるデータ型です。

(x <- list("Oh", "Nomura", "Kadota"))
(x <- list("Oh", "Nomura", "Kadota", c(868, 657, 567)))

x[[2]]
x[[4]]

    [[1]]
    [1] "Oh"

    [[2]]
    [1] "Nomura"

    [[3]]
    [1] "Kadota"
    [[1]]
    [1] "Oh"

    [[2]]
    [1] "Nomura"

    [[3]]
    [1] "Kadota"

    [[4]]
    [1] 868 657 567
    [1] "Nomura"
    [1] 868 657 567

各々の要素に名前を付けることも可能です。

(x <- list(First="Oh", Second="Nomura", Third="Kadota", Home_Runs=c(868, 657, 567)))

    $First
    [1] "Oh"

    $Second
    [1] "Nomura"

    $Third
    [1] "Kadota"

    $Home_Runs
    [1] 868 657 567

として、各々の要素に First, Second, Third, Home_Runs という名前が 付いていることを確認しましょう。

x[[4]]
x$Home_Runs

[1] 868 657 567
[1] 868 657 567

課題 6-1: 名前が Yamagata のベクトル(“Tsuruoka”, “Sakata”, “Yonezawa”)と名前が Niigata のベクトル(“Yuzawa”, “Niitsu”, “Shibata”)を要素にしたようなリストを作成し、そ の中味を表示し てみましょう。8

R にはデータフレームと呼ばれる型が存在します。

• リストの一種で、行列と同じ二次元配列です。

データフレームの各々の行は1つのサンプルに対する観測値、 各列はサン プルに対する1つの項目とみなすことができます。 これを使うと、表を簡単に表現するこ とができます。 以下の表は日本プロ野球の5人の元選手の名前、在籍チーム名、通算打数、 安打数、本塁打数を表したものです。

Home_Runs 2901 657 3085 504 2543 504 2786 868 2566 567

(x <- data.frame(
 row.names = c("Nomura", "Harimoto", "Kinugasa", "Oh", "Kadota"),
 Team = c("Hawks", "Flyers", "Carps", "Giants", "Hawks"),
 At_Bat = c(10472, 9666, 9404, 9250, 8868),
 Hits = c(2901, 3085, 2543, 2786, 2566),
 Home_Runs = c(657, 504, 504, 868, 567)
))

           Team At_Bat Hits Home_Runs
Nomura    Hawks  10472 2901       657
Harimoto Flyers   9666 3085       504
Kinugasa  Carps   9404 2543       504
Oh       Giants   9250 2786       868
Kadota    Hawks   8868 2566       567

x$Hits
x[ c(1,5), c(2,3,4) ] # 1、5 行目と 2、3、4 列目を抽出する

[1] 2901 3085 2543 2786 2566
       At_Bat Hits Home_Runs
Nomura  10472 2901       657
Kadota   8868 2566       567

attributes(x)

$names
[1] "Team"      "At_Bat"    "Hits"      "Home_Runs"

$class
[1] "data.frame"

$row.names
[1] "Nomura"   "Harimoto" "Kinugasa" "Oh"       "Kadota"

names(x)
row.names(x)
class(x)

[1] "Team"      "At_Bat"    "Hits"      "Home_Runs"
[1] "Nomura"   "Harimoto" "Kinugasa" "Oh"       "Kadota"
[1] "data.frame"