Есть ли способ, которым вы можете дифференцировать гауссову смесь в R? Чтобы получить гауссово распределение плотности смеси, я использовал следующий код:
dnorm_mix <- function(x, weights, means, sds) {
value <- 0
for (i in 1:length(weights)) {value <- value + weights[i]*dnorm(x, mean = means[i], sd = sds[i])}
return(value)
}
Может кто-нибудь помочь мне найти выражения для первого и второго производных этого? Я пытался использовать встроенную функцию 'производной' в R, но выдает ошибку 'dnorm_mix' отсутствует в таблице производных.
Спасибо!