<div dir="ltr"><div>Hi Sir</div><div><br></div><div>Could you please help for the following error message.<br></div><div><br></div><div>Below is the minimum working example and error I got.</div><div><br></div><div>
<div>y <- matrix(1:10, nrow=5)</div><div>mean <- c(4, 5)</div><div>sigma <- matrix(c(0.01, 0.009, 0.009, 1.2), nrow =2, byrow = TRUE)</div><div>dmvnormC(y, mean, sigma, 0) # Function is below.<br></div><div><table class="gmail-GD15MCFCNTB gmail-ace_text-layer gmail-ace_line gmail-GD15MCFCAW" style="font-style:normal;font-variant-ligatures:normal;font-variant-caps:normal;font-variant-numeric:inherit;font-variant-east-asian:inherit;font-weight:400;font-stretch:inherit;font-size:12pt;line-height:1.2;font-family:"Lucida Console",monospace;padding-left:6px;padding-bottom:8px;border:medium none;outline:currentcolor none medium;white-space:pre-wrap;color:rgb(0,0,0);letter-spacing:normal;text-align:start;text-indent:0px;text-transform:none;word-spacing:0px;background-color:rgb(255,255,255);text-decoration-style:initial;text-decoration-color:initial;width:468.546px" cellspacing="0" cellpadding="0"><tbody><tr><td style="font-family:"Lucida Console",monospace;line-height:1.2;font-size:12pt;vertical-align:top" align="left"><pre tabindex="0" class="gmail-GD15MCFCEUB" id="gmail-rstudio_console_output" style="font-family:"Lucida Console",monospace;font-size:12pt;outline:currentcolor none medium;border:medium none;word-break:break-all;margin:0px;white-space:pre-wrap;line-height:1.2"><span><span class="gmail-GD15MCFCOTB gmail-ace_constant" style="color:rgb(197,6,11)">Error in dmvnormC(y, mean, sig, 0) :
Cannot convert object to a function: [type=NULL; target=CLOSXP, SPECIALSXP, or BUILTINSXP].
</span></span></pre></td></tr><tr><td style="font-family:"Lucida Console",monospace;line-height:1.2;font-size:12pt;vertical-align:top" align="left"><br></td></tr><tr><td style="font-family:"Lucida Console",monospace;line-height:1.2;font-size:12pt;vertical-align:top" align="left"><table style="width:468.788px" cellspacing="0" cellpadding="0"><tbody><tr><td rowspan="1" style="font-family:"Lucida Console",monospace;line-height:1.2;font-size:12pt;vertical-align:top" width="1" height="" align="left"></td></tr></tbody></table></td></tr></tbody></table></div><br></div><div>#include <RcppArmadillo.h><br>using namespace Rcpp;<br>using namespace RcppArmadillo;<br>using namespace arma;<br><br>//[[Rcpp::depends(RcppArmadillo)]]<br>//[[Rcpp::export]]<br><br>arma::vec dmvnormC(const arma::mat& x, const arma::vec& mean, const arma::mat& sigma, int give_log)<br>{<br> Environment pkg = Environment::namespace_env("mvtnorm");<br> Function dmvnorm = pkg["mvtnorm"];<br> int N = x.n_rows;<br> arma::rowvec f(N, fill::zeros);<br> for(int k=0; k<N; k++)<br> {<br> f(k) = Rcpp::as<double>(dmvnorm(x.row(k), mean, sigma, give_log));<br> }<br> return f;<br>}</div><div><br></div><div><br></div><div>Thank You<br></div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Thu, Dec 19, 2019 at 7:57 PM Shaami <<a href="mailto:nzshaam@gmail.com">nzshaam@gmail.com</a>> wrote:<br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div dir="ltr"><div>Dear Sir</div><div><br></div><div>It compiles successfully. But when I use it I get an error message.</div><div><br></div><div><br></div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Thu, Dec 19, 2019 at 7:41 PM Shaami <<a href="mailto:nzshaam@gmail.com" target="_blank">nzshaam@gmail.com</a>> wrote:<br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div dir="auto">Great. It worked.<div dir="auto"><br></div><div dir="auto">Thank you very much</div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Thu, 19 Dec 2019, 6:01 PM Qiang Kou, <<a href="mailto:qkou@qkou.info" target="_blank">qkou@qkou.info</a>> wrote:<br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div dir="ltr">Please try the line below:<div><br></div><div>f(k) = <b>Rcpp::as<double></b>(dmvnorm(x.row(k), mean, sigma, give_log));<br></div><div><br></div><div>Best,</div><div><br></div><div>KK</div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Wed, Dec 18, 2019 at 8:27 PM Shaami <<a href="mailto:nzshaam@gmail.com" rel="noreferrer" target="_blank">nzshaam@gmail.com</a>> wrote:<br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div dir="ltr"><div>Hi Dear</div><div><br></div><div>I want to use dmvnorm() function into an cpp function using RcppArmadillo I am trying to execute the following function. It gives me error <br></div><div><br></div><div>cannot convert 'SEXP' to 'double' in assignment</div><div><br></div><div>Could you please guide where I am wrong?</div><div><br></div><div>#include <RcppArmadillo.h><br>using namespace Rcpp;<br>using namespace RcppArmadillo;<br>using namespace arma;<br><br>//[[Rcpp::depends(RcppArmadillo)]]<br>//[[Rcpp::export]]<br><br>arma::rowvec dmvnormC(const arma::mat& x, const arma::vec& mean, const arma::mat& sigma, int give_log)<br>{<br> Environment pkg = Environment::namespace_env("mvtnorm");<br> Function dmvnorm = pkg["mvtnorm"];<br> int N = x.n_rows;<br> arma::rowvec f(N, fill::zeros);<br> for(int k=0; k<N; k++)<br> {<br> f(k) = dmvnorm(x.row(k), mean, sigma, give_log);<br> }<br> return f;<br>}</div><div><br></div><div>Thank you</div><div><br></div><div>Shaami<br></div></div>
_______________________________________________<br>
Rcpp-devel mailing list<br>
<a href="mailto:Rcpp-devel@lists.r-forge.r-project.org" rel="noreferrer" target="_blank">Rcpp-devel@lists.r-forge.r-project.org</a><br>
<a href="https://lists.r-forge.r-project.org/cgi-bin/mailman/listinfo/rcpp-devel" rel="noreferrer noreferrer" target="_blank">https://lists.r-forge.r-project.org/cgi-bin/mailman/listinfo/rcpp-devel</a></blockquote></div>
</blockquote></div>
</blockquote></div>
</blockquote></div>