写一段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)))