具有矩阵乘法的Deriv函数

问题描述

如果定义的函数中有％*％，我们可以使用Deriv（）吗？

f <- function(alpha,beta,theta,rho,y,x,T,sig_u,sig_v,sig_uv){
  sigma = matrix(c(sig_u,sig_uv,sig_v),2,2)# Covariance matrix
  z = vec(cbind(y[-1],x[-1]))
  Z = kronecker(diag(2),cbind(1,x[-T]))
  B = as.vector(c(alpha,rho))
  res = z - Z %*% B
  cov_var = (kronecker(solve(sigma),diag(T-1)))
  r ~ T/2 * log(det(solve(sigma))) - 1/2 * t(res)  %*% cov_var %*% (res)
} 

Deriv(f(alpha=0,beta=0.5,theta=0,rho=0.999,y=1:10,x=1:10,T=10,sig_u =5,sig_v=1,sig_uv=0.2),c("alpha","beta","theta","rho"))

解决方法

暂无找到可以解决该程序问题的有效方法，小编努力寻找整理中！

如果你已经找到好的解决方法，欢迎将解决方案带上本链接一起发送给小编。

小编邮箱:dio#foxmail.com (将#修改为@）