R中的双和计算

Question

我在将以下双和写入r中的函数时遇到问题：

其中I函数是指标函数，我们有n个观察值x_1 ... x_n。

我试过这样的事情

sum( (x <= t) / (sum((x<=x)))^2)

Answer 1

你可以这样做：

sumiI <- function(xi, x) sum(xi<=x)
sum(1/sapply(x[x<=t], FUN=sumI, x=x))