asplit(行列)

asplit は行列を行ごと (MARGIN = 1) ないし列ごと (MARGIN = 2) のリストに変換する．

(M <- matrix(c(11, 21, 12, 22, 13, 23), 2, 3))
#>      [,1] [,2] [,3]
#> [1,]   11   12   13
#> [2,]   21   22   23

# 行ごとのリストに分割
(x <- asplit(M, 1))
#> [[1]]
#> [1] 11 12 13
#> 
#> [[2]]
#> [1] 21 22 23

# 行と列で分割しようとするとエラー
asplit(M, 1:2)
#> Error in array(newx[, i], d.call, dn.call): 'dims' cannot be of length 0

asplit(配列)

配列の場合は 3 以上の MARGIN を与えることができる．また， MARGIN = 1:2 といった具合に複数の次元で切ることもできる．

# 配列の容易
(A <- array(M, c(2, 2, 2)))
#> , , 1
#> 
#>      [,1] [,2]
#> [1,]   11   12
#> [2,]   21   22
#> 
#> , , 2
#> 
#>      [,1] [,2]
#> [1,]   13   11
#> [2,]   23   21

# 3次元目で分割
asplit(A, 3)
#> [[1]]
#>      [,1] [,2]
#> [1,]   11   12
#> [2,]   21   22
#> 
#> [[2]]
#>      [,1] [,2]
#> [1,]   13   11
#> [2,]   23   21

# 1:2 次元目で分割
(y <- asplit(A, 1:2))
#>      [,1]      [,2]     
#> [1,] Numeric,2 Numeric,2
#> [2,] Numeric,2 Numeric,2

猶， y の出力が一風変わっていることに注目されたい．実は，asplit はリストの配列を返り値に持つ． y の各要素は Numeric, 2 と表示されているので一見ベクトルであるが，実体は配列である．ある変数が配列であるかどうかは is.array が TRUE を返すことや dim が整数ベクトルを返すことで確認できる．

# y は配列なので，dimが数値を返す
dim(y)
#> [1] 2 2

# y の各要素も配列
lapply(y, dim)
#> [[1]]
#> [1] 2
#> 
#> [[2]]
#> [1] 2
#> 
#> [[3]]
#> [1] 2
#> 
#> [[4]]
#> [1] 2

先の行列の例においても，返り値は数値ベクトルのリストに見えるが，数値の一次元配列を内包したリストの一次元配列である．

# x は配列
dim(x)
#> [1] 2

# x の各要素も配列
lapply(x, dim)
#> [[1]]
#> [1] 3
#> 
#> [[2]]
#> [1] 3

array_tree(行列)

行列を行ごとに切ると，ベクトルを束ねたリストを返す．

M
#>      [,1] [,2] [,3]
#> [1,]   11   12   13
#> [2,]   21   22   23

str(array_tree(M, 1))
#> List of 2
#>  $ : num [1:3] 11 12 13
#>  $ : num [1:3] 21 22 23

# asplit は配列のリストを返すため，
# 要素の長さに (1d) と記載される
str(asplit(M, 1))
#> List of 2
#>  $ : num [1:3(1d)] 11 12 13
#>  $ : num [1:3(1d)] 21 22 23
#>  - attr(*, "dim")= int 2

2行3列の行列を行→列の順で切ると長さ3のリストを2つ束ねたリストを返す．

また逆順で列→行の順で切ると長さ2のリストを3つ束ねたリストを返す．

str(array_tree(M, 1:2))
#> List of 2
#>  $ :List of 3
#>   ..$ : num 11
#>   ..$ : num 12
#>   ..$ : num 13
#>  $ :List of 3
#>   ..$ : num 21
#>   ..$ : num 22
#>   ..$ : num 23

str(array_tree(M, 2:1))
#> List of 3
#>  $ :List of 2
#>   ..$ : num 11
#>   ..$ : num 21
#>  $ :List of 2
#>   ..$ : num 12
#>   ..$ : num 22
#>  $ :List of 2
#>   ..$ : num 13
#>   ..$ : num 23

# asplit では次元数と同じ長さの MARGIN に対してエラーを返す
asplit(M, 1:2)
#> Error in array(newx[, i], d.call, dn.call): 'dims' cannot be of length 0

margin 引数を省略すると， margin = 1:2 として扱う．

identical(array_tree(M), array_tree(M, 1:2))
#> [1] TRUE

array_branch(配列)

行列での挙動をそのまま配列に拡張しただけなので，特筆することはない，

A
#> , , 1
#> 
#>      [,1] [,2]
#> [1,]   11   12
#> [2,]   21   22
#> 
#> , , 2
#> 
#>      [,1] [,2]
#> [1,]   13   11
#> [2,]   23   21

# margin の長さが1の時
str(array_tree(A, 3))
#> List of 2
#>  $ : num [1:2, 1:2] 11 21 12 22
#>  $ : num [1:2, 1:2] 13 23 11 21

# margin の長さが2以上の時
str(array_tree(A, 1:2))
#> List of 2
#>  $ :List of 2
#>   ..$ : num [1:2] 11 13
#>   ..$ : num [1:2] 12 11
#>  $ :List of 2
#>   ..$ : num [1:2] 21 23
#>   ..$ : num [1:2] 22 21

# margin を省略した時
str(array_tree(A))
#> List of 2
#>  $ :List of 2
#>   ..$ :List of 2
#>   .. ..$ : num 11
#>   .. ..$ : num 13
#>   ..$ :List of 2
#>   .. ..$ : num 12
#>   .. ..$ : num 11
#>  $ :List of 2
#>   ..$ :List of 2
#>   .. ..$ : num 21
#>   .. ..$ : num 23
#>   ..$ :List of 2
#>   .. ..$ : num 22
#>   .. ..$ : num 21

array_tree(ベクトル)

ベクトルに対しては as.list に似た挙動を示す．

str(array_tree(1:3))
#> List of 3
#>  $ : int 1
#>  $ : int 2
#>  $ : int 3