对于两个给定的向量，不一定是相同的长度，$x=(x_1,\dots,x_N)$和$s=(s_1,\dots,s_M)$我想在 R 中尽可能有效地计算以下术语。对于$j=1, \dots, N$

下面是一些非常低效的测试代码。但是，我看不出如何使用矢量化样式或其他聪明的技巧来加快速度。

x <- rnorm(1000)
s <- rnorm(10,2,0.2)

l <- rep(0,length(s))
output <- rep(0,length(x))
for(j in 1:length(x)){
  xj <- x[j]
  for(i in 1:length(s)){
    s.i.dropped <- s[-i]
    l[i] <- prod((xj-s.i.dropped)/(s[i]-s.i.dropped))
  } 
  output[j] <- sum(l)
}