estou tentando integrar, sob toda a reta, a função:
integrand_1<-function(x){
exp(-x/phi-y*exp(-x)/(phi*2)-(x+sigma2/2)^2/(2*sigma2))}
Usando os seguintes dados:
df <- readRDS(url("https://rawcdn.githack.com/fsbmat/StackOverflow/db2ff43778b047eab9761cdedb16f70a08ec2f89/df.rds"))
phi <- 0.02216155
rho <- 0.6003019
sigma2 <- 0.02652626
y <- df$y
mu <- df$mu
Para isso, uso a função integrate, conforme abaixo:
integrate(Vectorize(integrand_1), -Inf, Inf)$value
No entanto, recebo a seguinte mensagem:
Error in integrate(Vectorize(integrand_1), -Inf, Inf) :
evaluation of function gave a result of wrong length
Tentei visualizar a função para ver se havia algo errado, usando esse código:
x <- seq(-10,10, length=300)
integrand_1(x)
is.numeric(integrand_1(x))
plot(x,integrand_1(x), xlim = c(2.3,2.6), type = "l")
O gráfico gerado foi:
Ou seja, não era para gerar o erro acima. Ou pelo menos não entendo o motivo do erro. Alguém, poderia me ajudar?