我想写一个函数来计算方差 - 协方差矩阵或二阶导数矩阵相对于下面函数的v
,我怎么能在R中做到这一点?我可以得到对角线值,但我不知道如何计算R中的协方差值。
x = c(1, 3, 5, 7, 10)
y = c(1.4, 3.5, 7.1, 8, 9)
v = list(a = .1, b= .5)
f = function(x, y, v) exp(v$a *x + v$b * y)
z = matrix(c(x, y), ncol=2)
var = colSums(z^2 * f(x, y, v))
var
# [1] 31363.68 29849.59