Pessoal estou criando uma função em R em que em determinada linha ocorre essa mensagem: $ operator is invalid for atomic vectors

O trecho do código é:

myintegrate = function(...) tryCatch(integrate(...), error=function(e) NA)
value_int = myintegrate(f=fdp,par=starts,lower=domain[1],upper=domain[2])$value
if(isTRUE(is.na(value_int))==TRUE) warning("Make sure that fdp is a probability density function. The integral in the domain specified is not convergent.")

[ ],
Pedro Rafael Diniz Marinho.