我想用r实现梯度下降法，其中需要对损失函数求导。我将损失函数写为
loss<-expression(sum((y-(a[1]*1+a[2]*x+a[3]*x1+a[4]*x2+a[5]*x3)^2))

其中y是因变量,x,x1,x2,x3是自变量。

然后我对上式求导
for (j in 1:5) {
       d[j]<- D(loss,"a[j]")
        }

错误提示是微分表中没有sum函数，这是怎么回事？该如何解决？跪谢各位大侠！