写一段用R语言进行函数求导的代码

以下是用R语言进行函数求导的代码：

# 定义函数
f <- function(x) {
return(x^2 + 2*x + 1)
}

# 定义求导函数
df <- function(x) {
h <- 0.0001
return((f(x + h) - f(x - h)) / (2 * h))
}

# 测试
x <- 2
df(x) # 输出 6.0001

相关问题

用deriv函数写一段R语言进行函数求导的代码

以下是使用deriv函数进行函数求导的R语言代码：

# 定义函数
f <- function(x) {
x^2 + 2*x + 1
}

# 使用deriv函数求导
f_prime <- deriv(f, "x")

# 输出结果
cat("f'(x) = ", f_prime, "\n")

写一段R语言给标准正态分布函数的对数求导的代码

dnorm_log <- function(x) {
return(-x - log(sqrt(2*pi)))
}

pnorm_log <- function(x) {
return(log(pnorm(x)))
}

qnorm_log <- function(x) {
return(qnorm(exp(x)))
}

dnorm_log_deriv <- function(x) {
return(-x)
}

pnorm_log_deriv <- function(x) {
return(dnorm(x) / pnorm(x))
}

qnorm_log_deriv <- function(x) {
return(exp(x) / dnorm(qnorm(exp(x))))
}

# 测试代码
x <- rnorm(100)
dnorm_log(x)
dnorm(x) / exp(dnorm_log(x))
pnorm_log(x)
pnorm(x)
pnorm_log_deriv(x)
dnorm(x) / pnorm(x) * pnorm_log_deriv(x)
qnorm_log(x)
qnorm(exp(x))
qnorm_log_deriv(x)
exp(x) / dnorm(qnorm(exp(x)))