エントリーの編集
エントリーの編集は全ユーザーに共通の機能です。
必ずガイドラインを一読の上ご利用ください。
記事へのコメント0件
- 注目コメント
- 新着コメント
このエントリーにコメントしてみましょう。
注目コメント算出アルゴリズムの一部にLINEヤフー株式会社の「建設的コメント順位付けモデルAPI」を使用しています
- バナー広告なし
- ミュート機能あり
- ダークモード搭載
関連記事
Rによる数値微分 - 餡子付゛録゛
Rはderiv()関数で記号/微分ができるものの、和や積、つまりsum()やprod()があると微分できません。以下... Rはderiv()関数で記号/微分ができるものの、和や積、つまりsum()やprod()があると微分できません。以下のようなコードを実行すると、『関数 'sum' は導関数の表中にありません』のようにエラーメッセージが出ます。 # ベクトルxを代入 x <- (1:10*23)%%7 # 変数bを2とする b <- 2 # sum()を含む関数を設定 e <- expression(exp(b*sum(x))) # 微分してみるとエラーが出る deriv(e, "b") 解決策は調べた限りは三つあります。まずは微分してしまいましょう。 # 一階微分 sum(x)*exp(b*sum(x)) # 二階微分 (sum(x)^2)*exp(b*sum(x)) 1.520372e+30と5.017226e+31が得られましたが、複雑な式になると便利ではありません。こういうときは数値微分を使います。