Please try the line below:
f(k) = *Rcpp::as<double>*(dmvnorm(x.row(k), mean, sigma, give_log));
Best,
KK
On Wed, Dec 18, 2019 at 8:27 PM Shaami <nzsh...@gmail.com> wrote:
> Hi Dear
>
> I want to use dmvnorm() function into an cpp function using RcppArmadillo
> I am trying to execute the following function. It gives me error
>
> cannot convert 'SEXP' to 'double' in assignment
>
> Could you please guide where I am wrong?
>
> #include <RcppArmadillo.h>
> using namespace Rcpp;
> using namespace RcppArmadillo;
> using namespace arma;
>
> //[[Rcpp::depends(RcppArmadillo)]]
> //[[Rcpp::export]]
>
> arma::rowvec dmvnormC(const arma::mat& x, const arma::vec& mean, const
> arma::mat& sigma, int give_log)
> {
> Environment pkg = Environment::namespace_env("mvtnorm");
> Function dmvnorm = pkg["mvtnorm"];
> int N = x.n_rows;
> arma::rowvec f(N, fill::zeros);
> for(int k=0; k<N; k++)
> {
> f(k) = dmvnorm(x.row(k), mean, sigma, give_log);
> }
> return f;
> }
>
> Thank you
>
> Shaami
> _______________________________________________
> Rcpp-devel mailing list
> Rcpp-devel@lists.r-forge.r-project.org
> https://lists.r-forge.r-project.org/cgi-bin/mailman/listinfo/rcpp-devel
_______________________________________________
Rcpp-devel mailing list
Rcpp-devel@lists.r-forge.r-project.org
https://lists.r-forge.r-project.org/cgi-bin/mailman/listinfo/rcpp-devel
