<div dir="ltr">Hi, I just tried to calculate the density of a normal variable using dnorm. I found that if the input is of type NumericVector, then it was fine. But if the input is of type double, it yielded error in compilation. For example:<div><br></div><div> // [[Rcpp::export]]</div><div>NumericVector test5(NumericVector x)</div><div>{</div><div>  NumericVector y = dnorm(x,5.0,2.0,true);</div><div>  return(y);</div><div>}</div><div><br></div><div>and</div><div><br></div><div><div>// [[Rcpp::export]]</div><div>double test6(double x)</div><div>{</div><div>  double y = dnorm(x,5.0,2.0,true);</div><div>  return(y);</div><div>}</div></div><div><br></div><div>May I ask why and how can I calculate this density (and call other common functions like dbeta, dt, beta, gamma, lbeta, lgamma, etc...) on double input?</div><div><br></div><div>Many thanks,</div><div>Eridk</div></div>