程序员经验分享-hwhale

如何计算 R 中样条函数的导数？

2023-11-24

R 可以使用 splines 库中的 splinefun() 生成样条函数。但是，我需要评估该函数的一阶和二阶导数。有没有办法做到这一点？

例如

library(splines)
x <- 1:10
y <- sin(pi/x) #just an example
f_of_x <- splinefun(x,y)

如何计算 x 向量的 f'(x)？


这很容易做到，因为计算函数导数的能力是内置于函数中的！

 f_of_x(x, deriv = 1)

谢谢R核心！

本文内容由网友自发贡献，版权归原作者所有，本站不承担相应法律责任。如您发现有涉嫌抄袭侵权的内容，请联系:hwhale#tublm.com(使用前将#替换为@)

r

numericalmethods

Spline

如何计算 R 中样条函数的导数？ 的相关文章

  • 跳过 data.table 中的 NA

    我想用data table但想跳过计算j部分如果by对应于缺失 NA 这是一个示例 data table library data table DT lt data table y 10 g c 1 1 1 2 2 2 2 2 NA NA

  • data.table `:=` 具有动态输入（现有列）和输出（新列名称）的赋值表达式

    Note 我在这个问题中遇到的确切问题不适用于最新版本的数据表 如果您想做标题中描述的事情 请查看软件包常见问题解答中的相应问题 1 6 好的 但我事先不知道表达方式 我如何以编程方式传递它们 https rawgit com wiki R

  • 在编译 html 笔记本中找不到对象

    使用 Rstudio 编译 html 笔记本 功能 在某些脚本上我收到对象未找到错误 当我逐行或按块运行脚本时 它运行良好并在控制台中生成输出 当我在同一块上运行 编译 html 笔记本 命令时 我得到的只是对象未找到错误 它似乎与哪个数据

  • 在 ggplot2 中对多个变量进行排序

    我正在尝试将变量分组并按降序排序 mydf region airport value MIA FLL 0 244587909 MIA PBI 0 824144687 MIA MIA 0 484907626 NYC EWR 0 7310755

  • 有没有比“[”更快的方法来对稀疏矩阵进行子集化？

    我是 seqMeta 包的维护者 正在寻找如何加速将大矩阵多次分割成小块的瓶颈的想法 背景 seqMeta 包用于分析遗传数据 所以你有一组受试者 n subject 和一些遗传标记 n snps 这导致 n subject x n snp

  • PowerBI：如何保存R脚本的结果？

    是否可以在 Power BI Desktop 中实现以下场景 将数据从 Excel 文件加载到多个表 使用 R 脚本从多个数据源进行计算 将计算结果存储到 Power BI pbix 中的新表 这个想法是使用 Power BI Deskto

  • R中引用和表达式的区别

    gettypes function x paramx substitute x print typeof paramx print class paramx print mode paramx gettypes expression sin

  • 检查相等性

    我想检查数据集的相等性 数据集看起来像这样 Equips lt c 1 1 1 2 2 2 3 3 3 3 3 3 3 4 4 4 4 4 4 5 5 5 5 5 5 5 6 7 8 Notifs lt c 10 10 20 55 63 6

  • R：如何检查向量元素是否相同

    gt dput ranks structure c NA NA 2L 2L NA NA NA 2L 1L Names c LANG1 LANG2 LANG3 LANG4 LANG5 LANG6 LANG7 LANG8 LANG9 gt ra

  • 在 R 中，如何将整个命令行放入二元运算符的 sys.call() 中？

    我为 R 编写了一个二元运算符函数 即名称类似于 X 这样就不必输入 X a b 我可以使用更方便的语法a X b 我的目标是有一个包装器 lt 它会执行一些操作 例如记录对该环境中的对象所做的操作并检查 受保护 属性 该属性会在覆盖该对象

  • r 选择数据框中某个值以下和之后的值

    我有一个问题如何从表中选择某些值 我有一个包含时间和值的表 我想在特定时间之后获取下面的行 示例 Data Frame Time Value 02 51 0 08033405 05 30 0 43456738 09 45 0 3605207

  • 使用 ff 包导入文本文件

    我有一个包含 450 万行和 90 列的文本文件要导入到 R 中 read table我得到了cannot allocate vector of size 错误消息 因此我尝试使用 ff 包导入 然后再对数据进行子集化以提取我感兴趣的观察结

  • 如何在 R 中建立具有某些条件的排列[重复]

    这个问题在这里已经有答案了 我是 R 新手 有点困惑 假设我有一个向量 c 1 2 3 4 5 6 我想生成具有四个元素的排列 每个排列应涉及 1 和 5 谢谢 您可以使用permutations https www rdocumentat

  • 根据“df1”的一个变量（“df1$var1”）在“df1”中创建一个变量，并根据“df1$var1”创建一个可更改的“df2”变量

    我有数据框df1总结了一段时间内鱼的深度 df1 Site告诉你鱼所在的地点 df1 Ind告诉你个人和df1 Depth告诉您鱼在特定位置的深度df1 Datetime 另一方面 我有df2总结了一段时间内 每三小时 从表面到 39 米深

  • R 将向量重塑为多列

    假设我在 R 中有一个向量 如下所示 d lt seq 1 100 我想将这个向量重塑为 10x10 矩阵 这样我就可以得到以下数据 1 2 3 10 1 2 3 10 11 12 13 20 21 22 23 30 91 92 93 10

  • R文本文件和文本挖掘...如何加载数据

    我正在使用 R 包tm我想做一些文本挖掘 这是一个文档 被视为一个词袋 我不明白有关如何加载文本文件并创建必要的对象以开始使用诸如 之类的功能的文档 stemDocument x language map IETF Language x 所

  • 使用函数保存 ggplot

    我想创建一个函数来保存绘图 来自ggplot 这是一个数据框 creating data frame music lt c Blues Hip hop Jazz Metal Rock number lt c 8 7 4 6 11 df mu

  • 如何在 ggplot2 图中有两个源标题？ [复制]

    这个问题在这里已经有答案了 我正在尝试在 ggplot2 图中添加第二个标题 与这位经济学家制作此图的方式类似 这是我制作的一个基本图 我知道如何在右下角添加一个标题 但如何在左下角添加另一个标题 ggplot mtcars aes mpg

  • 如何在R中对多个时间序列应用dtw算法？

    Problem 我有不同车辆速度的时间序列 我的最终目标是根据不同车辆在时间上的速度相似性对它们进行聚类 因此 我基本上需要生成一个距离矩阵 其中每个单元格包含一对车辆速度时间序列之间的距离 我想使用动态时间扭曲 dtw 作为距离度量 所以

  • data.table 和 pmin 带有 na.rm=TRUE 参数

    我正在尝试使用 pmin 函数和 data table 计算行间的最小值 类似于帖子data table 中的逐行操作和更新 https stackoverflow com questions 16860200 row by row ope

随机推荐

热门标签