我有以下公式:
其中T
是转置的,我想计算上述函数相对于x的一阶导数。我怎么能在R?
x = c(1, 2)
y = cbind(y1, y2)
y1 = runif(n = 100, -2, 2)
y2 = rbinom(n = 100, 1, 0.5)
最后我需要四个值,这是x的导数。
答案 0 :(得分:0)
有一个包numDeriv来计算数值导数。所以你会这样做:
library("numDeriv")
lapply(x, function(z) grad(f, z))
计算它
什么是T
?如果T
是常数,则可以通过分析计算导数,结果为
直接计算。