[Code]
p <- c(1 / 2, 1 / 4, 1 / 4)
len <- c(1, 2, 2)
lengths_seen <- sort(unique(len))
prob_of_length <- vapply(
lengths_seen,
function(k) sum(p[len == k]),
numeric(1)
)
from_definition <- sum(lengths_seen * prob_of_length)
from_sum <- sum(p * len)
c(from_definition = from_definition, from_sum = from_sum)
#> from_definition from_sum
#> 1.5 1.5