17

time2つの変数、およびを使用して、時系列の一次導関数(dpH / dtime)を計算しますpH

Rでこれを行うための関数はありますか、それともこれを行うために追加の関数を計算する必要がありますか?

4

3 に答える 3

19

pHtimeが単純なベクトルであると仮定すると、これを試してください。

library(pspline)
predict(sm.spline(time, pH), time, 1)
于 2012-12-29T15:01:51.247 に答える
12

Matt Lが提案したように、stats::deriv またはそれから始めたいと思うかもしれません。diff.ts私の教授がすべての学生に言ったことを覚えておいてください。数値微分は「誤差乗数」として知られています。

編集:明確にするために-彼が警告していたのは、データ内のノイズが導関数の推定値を大きく外す可能性があるということでした。積分はローパスフィルターであり、微分はハイパスフィルターであると言われています。したがって、重要なことは、導関数を計算する前にデータを平滑化することです。したがって、Gaborの優れた使用提案predict.spline。ただし、スプラインパラメータを変更すると、データがさまざまなレベルに滑らかになることに注意してください。常に結果を確認して、明らかなノイズは除去されているが、目的の機能は除去されていないことを確認してください。

于 2012-12-29T14:37:19.283 に答える
0

ここに「数値微分」へのリンクがあります。

http://en.wikipedia.org/wiki/Numerical_differentiation

テイラー級数展開に基づく方法を説明するリンクは次のとおりです。

http://ocw.usu.edu/civil_and_environmental_engineering/numerical_methods_in_civil_engineering/ODEsMatlab.pdf

于 2012-12-29T16:49:57.047 に答える