1

背景： 我正在使用这个Ryacas包，试图找出一个MLE估计量的大样本方差的符号表达式。

为此，我需要Hessian似然函数的（逆）矩阵。我没有mathematica（而且在线版本对我来说似乎太麻烦了），因此我正在尝试使用作为计算机代数系统Ryacas接口的包。YACAS

问题： 但是，我似乎无法弄清楚如何计算Hessian。在此处使用指南：https ://cran.r-project.org/web/packages/Ryacas/vignettes/high-level.html给了我一个错误。这是一个最小的可重现示例（取自该指南）

L <- yac_symbol("x^2 * (y/4) - a*(3*x + 3*y/2 - 45)")
Hessian(L)
Error in Hessian(L) : could not find function "Hessian"

当我以另一种方式尝试时，即使用（新？）界面

L <- yac_symbol("x^2 * (y/4) - a*(3*x + 3*y/2 - 45)")
y_fn(L, "HessianMatrix")

我似乎也无法得到可用的答案，但只是HessianMatrix((x^2*y)/4-a*(3*x+(3*y)/2-45))

有人知道如何解决这个问题吗？将不胜感激！

谢谢

4

3 回答 3

2 
library(Ryacas)

yac_str("HessianMatrix(x^2 * (y/4) - a*(3*x + 3*y/2 - 45), {x,y})")
# "{{y/2,x/2},{x/2,0}}"

yac_str("PrettyForm(HessianMatrix(x^2 * (y/4) - a*(3*x + 3*y/2 - 45), {x,y}))") %>% cat
# /              \
# | / y \ / x \  |
# | | - | | - |  |
# | \ 2 / \ 2 /  |
# |              |
# | / x \ ( 0 )  |
# | | - |        |
# | \ 2 /        |
# \              /

参考：yacas.readthedocs

于 2019-10-28T08:05:23.790 回答
1

以下是您可以使用 Ryacas v1.1.1 执行的操作：

> library(Ryacas)
> packageVersion("Ryacas")
[1] ‘1.1.1’
> x <- ysym("x")
> y <- ysym("y")
> a <- ysym("a")
> L <- x^2 * (y/4) - a*(3*x + 3*y/2 - 45)
> H <- Hessian(L, c("x", "y", "a"))
> H
{{   y/2,    x/2,     -3},
 {   x/2,      0, (-3)/2},
 {    -3, (-3)/2,      0}} 
> as_r(H)
expression(rbind(c(y/2, x/2, -3), c(x/2, 0, -3/2), c(-3, -3/2, 
    0)))
> eval(as_r(H), list(x = 2, y = 2, a = 2))
     [,1] [,2] [,3]
[1,]    1  1.0 -3.0
[2,]    1  0.0 -1.5
[3,]   -3 -1.5  0.0
于 2019-11-07T22:59:11.187 回答
0

错误Error in Hessian(L) : could not find function "Hessian"是由于Hessian函数不是Ryacas包的一部分，而是numDerivR 中包的一部分。

Hessian下面还有R 中函数的文档：

https://www.rdocumentation.org/packages/numDeriv/versions/2016.8-1.1/topics/hessian

我希望这能解决你的问题。所以请先安装numDeriv软件包，然后再使用该Hessian功能。

我也尝试使用上面的例子，结果如下：

library(Ryacas)
library(numDeriv)
L <- yac_symbol("x^2 * (y/4) - a*(3*x + 3*y/2 - 45)")
Hessian(L)

输出：

{}

输出是空集，但我想你可以参考Hessian函数的文档，它肯定会有所帮助。

于 2019-10-27T18:16:02.700 回答