time
2つの変数、およびを使用して、時系列の一次導関数(dpH / dtime)を計算しますpH
。
Rでこれを行うための関数はありますか、それともこれを行うために追加の関数を計算する必要がありますか?
time
2つの変数、およびを使用して、時系列の一次導関数(dpH / dtime)を計算しますpH
。
Rでこれを行うための関数はありますか、それともこれを行うために追加の関数を計算する必要がありますか?
pH
とtime
が単純なベクトルであると仮定すると、これを試してください。
library(pspline)
predict(sm.spline(time, pH), time, 1)
Matt Lが提案したように、stats::deriv
またはそれから始めたいと思うかもしれません。diff.ts
私の教授がすべての学生に言ったことを覚えておいてください。数値微分は「誤差乗数」として知られています。
編集:明確にするために-彼が警告していたのは、データ内のノイズが導関数の推定値を大きく外す可能性があるということでした。積分はローパスフィルターであり、微分はハイパスフィルターであると言われています。したがって、重要なことは、導関数を計算する前にデータを平滑化することです。したがって、Gaborの優れた使用提案predict.spline
。ただし、スプラインパラメータを変更すると、データがさまざまなレベルに滑らかになることに注意してください。常に結果を確認して、明らかなノイズは除去されているが、目的の機能は除去されていないことを確認してください。
ここに「数値微分」へのリンクがあります。
http://en.wikipedia.org/wiki/Numerical_differentiation
テイラー級数展開に基づく方法を説明するリンクは次のとおりです。