Rにより連続したデータ数をカウントしたいです。

Question

ID毎の形質の”0”と”1”の連続した長さとその回数をカウントしたいです。
データ：

ID	形質1	形質2
a	0	1
a	0	1
a	1	1
a	0	0
a	1	0
a	1	0
b	0	1
b	1	1
b	1	1
b	0	1
b	1	0
b	1	1

欲しい結果は
|ID|値|連続長さ|形質1についての回数|形質2についての回数|
|:--|:--:|--:|--:|
|a|0|1|1|0
|a|0|2|1|0
|a|0|3|0|1
|a|0|4|0|0
|b|1|1|0|1
|b|1|2|2|0
|b|1|3|0|0
|b|1|4|0|1

表の見方は（2行目）IDがaのとき”0”が2連続の長さで、形質1では1回、形質2では0回になります。
このようにカウントする方法を知りたいです。
よろしくお願いいたします。

Answer 1

「連続した長さとその回数」を得るにはRle関数が使えます。

> Rle(c(1,1,0,0,0,1,1,1,0,0,0,0))
numeric-Rle of length 12 with 4 runs
  Lengths: 2 3 3 4
  Values : 1 0 1 0

データフレームのコラムごと／値ごとにRleを呼び出すようなループを組めば良いと思います。

Comments

[orgil]

ありがとうございました。
大変助かりました。
b<- cbind(ID=c(rep("a",6), rep("b", 6)), 形質1=c(0,0,1,0,1,1,0,1,1,0,1,1),形質2=c(1,1,1,0,0,0,1,1,1,1,0,1))
> b
ID 形質1 形質2
[1,] "a" "0" "1"
[2,] "a" "0" "1"
[3,] "a" "1" "1"
[4,] "a" "0" "0"
[5,] "a" "1" "0"
[6,] "a" "1" "0"
[7,] "b" "0" "1"
[8,] "b" "1" "1"
[9,] "b" "1" "1"
[10,] "b" "0" "1"
[11,] "b" "1" "0"
[12,] "b" "1" "1"
> tapply(b[,2], b[,1], rle)
$`a`
Run Length Encoding
lengths: int [1:4] 2 1 1 2
values : chr [1:4] "0" "1" "0" "1"

$b
Run Length Encoding
lengths: int [1:4] 1 2 1 2
values : chr [1:4] "0" "1" "0" "1"

さらに、質問したいですがよろしくお願いいたします。

今回はID毎にグループ化してカウントすることができました。
さらに、IDと形質2毎にグループ化して形質1をカウントしたいです。
例：IDがa、形質２が１のときの形質２をカウントしたい、結果は長さ：2、1で、値：０、１になります。
よろしくお願いいたします。

[KojiDoi]

質問はつぎ足さないで新しい記事に書いてください。

[orgil]

了解いたしました。
ありがとうございました。