使用马赛克包中的 D 函数计算 R 中的导数

Question

我安装了 R 中可用的马赛克包，通过键入 install.packages('mosaic') 来计算导数。 我使用 makeFun 创建了一个函数，然后我尝试按以下方式计算导数

y1 <- makeFun(a +b *x  ~x, a=2, b=2)
dy1.dx <- D(a + b* x ~ x, a=2, b=2)
dy1.dx

但控制台显示

" dy1.dx <- D(a+b*x~x, a=2, b=2)
Error in D(a + b * x ~ x, a = 2, b = 2) : unused arguments (a = 2, b = 2)"

我该如何纠正？

Answer 1

使用基本 R 的stats包。 无需安装其他软件包，如mosaic 。

D将给出表达式的导数，因此使用expression()函数，我们创建一个表达式并将其传递给D函数。

然后eval将评估表达式，而substitute将替换表达式中a和b值。

得到表达式关于 x 的导数：

stats::D(expression(a + b * x), "x")
# b

在用导数中的值替换后计算表达式。 b被替换为值 2。

eval(substitute( stats::D(expression(a + b * x), "x"), list(a=2, b = 2) ))
# [1] 2

另一个例子：

stats::D(expression(a + a*b * x), "x")
# a * b
eval(substitute( D(expression(a + a*b * x), "x"), list(a=3, b = 2) ))
# 3 * 2
eval(eval(substitute( D(expression(a + a*b * x), "x"), list(a=3, b = 2) )))
# 6

Answer 2

您的示例中使用的D来自library(mosaicCalc)而不是来自 base stat。 安装并调用库。 您的功能正常工作。

require(mosaicCalc)
#> Loading required package: mosaicCalc
#> Warning: package 'mosaicCalc' was built under R version 3.6.3
#> Loading required package: mosaicCore
#> Warning: package 'mosaicCore' was built under R version 3.6.3
#> Registered S3 method overwritten by 'mosaic':
#>   method                           from   
#>   fortify.SpatialPolygonsDataFrame ggplot2
#> 
#> Attaching package: 'mosaicCalc'
#> The following object is masked from 'package:stats':
#> 
#>     D
dy1.dx <- D(a + b * x ~ x, a = 2, b = 2)
dy1.dx
#> function (x, a = 2, b = 2) 
#> b