我有一系列的1和0,我想计算更改次数。 e.g。
x <- rbinom(10, 1, 1/2)
> x
[1] 0 0 1 1 1 1 1 0 1 0
因此,我想计算(在R中)序列从1到0交替(或翻转)的次数。在上述顺序中,交替次数(手动计数)为4。
答案 0 :(得分:12)
您可以使用diff():
> x <- rbinom(10,1,1/2)
> x
[1] 0 0 0 1 1 1 1 0 1 0
> sum(diff(x)!=0)
[1] 4
答案 1 :(得分:6)
rle函数将计算向量中相同值的“运行”数。因此,此向量的长度(减1)为您提供了更改次数:
> x
[1] 0 0 0 1 1 1 1 0 1 0
> rle(x)
Run Length Encoding
lengths: int [1:5] 3 4 1 1 1
values : num [1:5] 0 1 0 1 0
> length(rle(x)$lengths)-1
[1] 4
可能比diff()方法更快或更慢,但如果你需要它也可以给你运行长度......
答案 2 :(得分:3)
它绝对不会在优雅方面超越差异,但另一种方式:
sum(x[-1] != head(x, n=-1))
在我的系统上,这似乎要快一点:
> x <- rbinom(1e6, 1, 0.5)
> system.time(replicate(100, sum(x[-1] != head(x, n=-1))))
user system elapsed
8.421 3.688 12.150
> system.time(replicate(100, sum(diff(x) != 0)))
user system elapsed
9.431 4.525 14.248
似乎应该有一个很好的分析解决方案来分配序列中不相等的相邻元素的数量。
答案 3 :(得分:0)
伪代码(序列是你的硬币翻转的数组):
variable count = 0
variable state = sequence[0]
iterate i from sequence[1] to sequence[n]
if (state not equal sequence[i])
state = 1 - state
count++
计数应该是你的结果