我的左侧公式如下,我将其展开到右侧,+
表示如果b
_ i大于x
,则表示b_i - x
a = c(1, 2)
b = c(10, 30)
x = c(5, 15, 50, 100)
> colSums(a * pmax(outer(b, x, "-"), 0)^2)
[1] 1275 450 0 0
,否则为零。我可以写下左侧公式的代码如下,但是如何得到R中二次右侧公式的系数?
\ sum_ {i = 1} ^ {2} a_i(b_i - x)^ 2 _ {+} = \ sum_ {i = 1} ^ {2} a_i(b_i ^ 2 I(b_i> = x) + x ^ 2 I(b_i> = x) - 2 b_i x I(b_i> = x))
{{1}}