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))