R中的多變量導數

Question

有人可以用R中的簡單數學運算來幫助我嗎？

我要做的就是：

我有一個功能

f <-(a1*x1+a2*x2+…+an*xn~x1&…&xn)

和該函數的派生

df<-D(f(x1…xn)~…)

我需要填補“ ...”缺口以獲得多變量導數的結果嗎？ 問題是我無法處理向量c（x1，x2，x3 ... xn）的值

結果，我需要在每個步驟上寫一個“ for”循環，以獲取派生結果

Answer 1

你有固定的n嗎？ 如果是這樣，則可以使用parse來轉換表達式中的字符串：

x = c('x1','x2','x3')
a = c('a1','a2','a3')
expr=parse(text=paste(paste(x, a, sep="*"), collapse="+"))

#> sapply(x, function(u) D(expr,u))
#$x1
#a1

#$x2
#a2

#$x3
#a3