R 에서 미분계수 구하는 법
개요
R 에서 미분계수를 구하기 위해선 numDeriv 패키지의 grad() 함수를 사용할 수 있다.
코드
예를 들어 $f(x) = x^2 + 4x + 1$ 과 $g(x) = e^{-x}$ 의 미분계수는 다음과 같이 구할 수 있다.
install.packages("numDeriv")
library(numDeriv)
f<-function(x) {x^2 + 4*x + 1}
g<-function(x) {exp(-x)}
grad(f,2)
grad(g,0)
실제로 계산해보면 $f ' (2) = 2 \cdot 2 + 4 = 8$ 이고 $g ' (0) = - e^{0} = -1$ 인 것을 확인할 수 있다.
참고로 스칼라 함수의 경우에도 x 옵션에 벡터를 넣어주면 그래디언트를 잘 계산해준다.